


Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 15.1
Initial simplification15.1
rmApplied cos-sum0.3
rmApplied *-un-lft-identity0.3
Applied times-frac0.3
Simplified0.3
rmApplied expm1-log1p-u0.3
rmApplied div-inv0.4
Applied associate-*r*0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2018360 +o rules:numerics
(FPCore (r a b)
:name "r*sin(b)/cos(a+b), B"
(* r (/ (sin b) (cos (+ a b)))))
Time bar (total: 17.1s)Debug log
2 alts after pruning (2 fresh and 0 done)
Merged error: 15.6b
Found 3 expressions with local error:
| 2.3b | (cos (+ b a)) |
| 0.2b | (/ (* r (sin b)) (cos (+ b a))) |
| 0.1b | (* r (sin b)) |
3 calls. Slowest were:
| 8.0ms | (/ (* r (sin b)) (cos (+ b a))) |
| 5.0ms | (cos (+ b a)) |
| 2.0ms | (* r (sin b)) |
3 calls. Slowest were:
| 437.0ms | (/ (* r (sin b)) (cos (+ b a))) |
| 51.0ms | (* r (sin b)) |
| 35.0ms | (cos (+ b a)) |
24 calls. Slowest were:
| 580.0ms | (/ (* (* (* r (sin b)) (* r (sin b))) (* r (sin b))) (* (* (cos (+ b a)) (cos (+ b a))) (cos (+ b a)))) |
| 283.0ms | (+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b))) |
| 131.0ms | (- 1 (+ (* 1/2 (pow b 2)) (* a b))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 0.2b
Found 4 expressions with local error:
| 0.2b | (/ (* r (sin b)) (- (* (cos b) (cos a)) (* (sin b) (sin a)))) |
| 0.2b | (* (sin b) (sin a)) |
| 0.1b | (- (* (cos b) (cos a)) (* (sin b) (sin a))) |
| 0.1b | (* (cos b) (cos a)) |
4 calls. Slowest were:
| 10.0ms | (/ (* r (sin b)) (- (* (cos b) (cos a)) (* (sin b) (sin a)))) |
| 7.0ms | (- (* (cos b) (cos a)) (* (sin b) (sin a))) |
| 3.0ms | (* (sin b) (sin a)) |
4 calls. Slowest were:
| 261.0ms | (/ (* r (sin b)) (- (* (cos b) (cos a)) (* (sin b) (sin a)))) |
| 111.0ms | (- (* (cos b) (cos a)) (* (sin b) (sin a))) |
| 75.0ms | (* (sin b) (sin a)) |
| 36.0ms | (* (cos b) (cos a)) |
51 calls. Slowest were:
| 406.0ms | (- (* (+ (cos (+ b a)) (cos (- b a))) 2) (* 2 (- (cos (- b a)) (cos (+ b a))))) |
| 321.0ms | (- (+ (cos (+ b a)) (cos (- b a))) (- (cos (- b a)) (cos (+ b a)))) |
| 256.0ms | (+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.2b | (/ (sin b) (- (* (cos b) (cos a)) (* (sin b) (sin a)))) |
| 0.2b | (* (sin b) (sin a)) |
| 0.2b | (* r (/ (sin b) (- (* (cos b) (cos a)) (* (sin b) (sin a))))) |
| 0.1b | (- (* (cos b) (cos a)) (* (sin b) (sin a))) |
4 calls. Slowest were:
| 18.0ms | (/ (sin b) (- (* (cos b) (cos a)) (* (sin b) (sin a)))) |
| 11.0ms | (* r (/ (sin b) (- (* (cos b) (cos a)) (* (sin b) (sin a))))) |
| 7.0ms | (- (* (cos b) (cos a)) (* (sin b) (sin a))) |
4 calls. Slowest were:
| 237.0ms | (* r (/ (sin b) (- (* (cos b) (cos a)) (* (sin b) (sin a))))) |
| 166.0ms | (/ (sin b) (- (* (cos b) (cos a)) (* (sin b) (sin a)))) |
| 91.0ms | (- (* (cos b) (cos a)) (* (sin b) (sin a))) |
| 66.0ms | (* (sin b) (sin a)) |
73 calls. Slowest were:
| 426.0ms | (- (* (+ (cos (+ b a)) (cos (- b a))) 2) (* 2 (- (cos (- b a)) (cos (+ b a))))) |
| 377.0ms | (- (+ (cos (+ b a)) (cos (- b a))) (- (cos (- b a)) (cos (+ b a)))) |
| 249.0ms | (+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.2b | (/ (sin b) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a)))))) |
| 0.2b | (* (sin b) (sin a)) |
| 0.2b | (* r (/ (sin b) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))))) |
| 0.1b | (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))) |
4 calls. Slowest were:
| 17.0ms | (/ (sin b) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a)))))) |
| 14.0ms | (* r (/ (sin b) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))))) |
| 7.0ms | (* (sin b) (sin a)) |
4 calls. Slowest were:
| 292.0ms | (* r (/ (sin b) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))))) |
| 151.0ms | (/ (sin b) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a)))))) |
| 149.0ms | (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))) |
| 59.0ms | (* (sin b) (sin a)) |
72 calls. Slowest were:
| 263.0ms | (+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b))) |
| 179.0ms | (+ (* (* (cos b) (cos a)) (* (cos b) (cos a))) (+ (* (expm1 (log1p (* (sin b) (sin a)))) (expm1 (log1p (* (sin b) (sin a))))) (* (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))))) |
| 167.0ms | (/ (* (* (sin b) (sin b)) (sin b)) (* (* (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a)))))) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 0.0b
0% (0.3b remaining)
Error of 0.4b against oracle of 0.1b and baseline of 0.4b