


Bits error versus a1



Bits error versus a2



Bits error versus th
Initial program 0.5
Simplified0.5
rmApplied associate-*l/0.5
rmApplied add-sqr-sqrt0.5
Applied associate-*r*0.5
Simplified0.5
rmApplied add-cube-cbrt0.5
Applied associate-/r*0.5
Final simplification0.5
herbie shell --seed 2019004 +o rules:numerics
(FPCore (a1 a2 th)
:name "Migdal et al, Equation (64)"
(+ (* (/ (cos th) (sqrt 2)) (* a1 a1)) (* (/ (cos th) (sqrt 2)) (* a2 a2))))
Time bar (total: 44.3s)Debug log
1 calls. Slowest were:
| 99.0ms | (+ (* (/ (cos th) (sqrt 2)) (* a1 a1)) (* (/ (cos th) (sqrt 2)) (* a2 a2))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 0.5b
Found 2 expressions with local error:
| 0.7b | (/ (cos th) (sqrt 2)) |
| 0.4b | (* (/ (cos th) (sqrt 2)) (fma a1 a1 (* a2 a2))) |
2 calls. Slowest were:
| 13.0ms | (* (/ (cos th) (sqrt 2)) (fma a1 a1 (* a2 a2))) |
| 7.0ms | (/ (cos th) (sqrt 2)) |
2 calls. Slowest were:
| 180.0ms | (* (/ (cos th) (sqrt 2)) (fma a1 a1 (* a2 a2))) |
| 49.0ms | (/ (cos th) (sqrt 2)) |
54 calls. Slowest were:
| 291.0ms | (/ (* (* (cos th) (cos th)) (cos th)) (* (* (sqrt 2) (sqrt 2)) (sqrt 2))) |
| 124.0ms | (- (+ (/ 1 (sqrt 2)) (* 1/24 (/ (pow th 4) (sqrt 2)))) (* 1/2 (/ (pow th 2) (sqrt 2)))) |
| 117.0ms | (* (/ (* (* (cos th) (cos th)) (cos th)) (* (* (sqrt 2) (sqrt 2)) (sqrt 2))) (* (* (fma a1 a1 (* a2 a2)) (fma a1 a1 (* a2 a2))) (fma a1 a1 (* a2 a2)))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 0.0b
Found 2 expressions with local error:
| 0.4b | (/ (* (cos th) (fma a1 a1 (* a2 a2))) (sqrt 2)) |
| 0.2b | (* (cos th) (fma a1 a1 (* a2 a2))) |
2 calls. Slowest were:
| 9.0ms | (* (cos th) (fma a1 a1 (* a2 a2))) |
| 8.0ms | (/ (* (cos th) (fma a1 a1 (* a2 a2))) (sqrt 2)) |
2 calls. Slowest were:
| 147.0ms | (/ (* (cos th) (fma a1 a1 (* a2 a2))) (sqrt 2)) |
| 124.0ms | (* (cos th) (fma a1 a1 (* a2 a2))) |
30 calls. Slowest were:
| 498.0ms | (* (* (* (cos th) (cos th)) (cos th)) (* (* (fma a1 a1 (* a2 a2)) (fma a1 a1 (* a2 a2))) (fma a1 a1 (* a2 a2)))) |
| 103.0ms | (/ (* (* (* (cos th) (cos th)) (cos th)) (* (* (fma a1 a1 (* a2 a2)) (fma a1 a1 (* a2 a2))) (fma a1 a1 (* a2 a2)))) (* (* (sqrt 2) (sqrt 2)) (sqrt 2))) |
| 98.0ms | (/ (* (* (* (cos th) (fma a1 a1 (* a2 a2))) (* (cos th) (fma a1 a1 (* a2 a2)))) (* (cos th) (fma a1 a1 (* a2 a2)))) (* (* (sqrt 2) (sqrt 2)) (sqrt 2))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 5.0b | (sqrt (fma a1 a1 (* a2 a2))) |
| 0.4b | (/ (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) (sqrt 2)) |
| 0.1b | (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) |
| 0.1b | (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) |
4 calls. Slowest were:
| 13.0ms | (/ (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) (sqrt 2)) |
| 11.0ms | (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) |
| 5.0ms | (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) |
4 calls. Slowest were:
| 221.0ms | (/ (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) (sqrt 2)) |
| 154.0ms | (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) |
| 107.0ms | (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) |
| 35.0ms | (sqrt (fma a1 a1 (* a2 a2))) |
53 calls. Slowest were:
| 1.1s | (* (* (* (* (cos th) (cos th)) (cos th)) (* (* (sqrt (fma a1 a1 (* a2 a2))) (sqrt (fma a1 a1 (* a2 a2)))) (sqrt (fma a1 a1 (* a2 a2))))) (* (* (hypot a1 a2) (hypot a1 a2)) (hypot a1 a2))) |
| 278.0ms | (* (* (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (* (cos th) (sqrt (fma a1 a1 (* a2 a2))))) (* (cos th) (sqrt (fma a1 a1 (* a2 a2))))) (* (* (hypot a1 a2) (hypot a1 a2)) (hypot a1 a2))) |
| 221.0ms | (- (log (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2))) (log (sqrt 2))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 0b
Found 4 expressions with local error:
| 5.0b | (sqrt (fma a1 a1 (* a2 a2))) |
| 0.3b | (/ (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) |
| 0.2b | (/ (/ (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (cbrt (sqrt 2))) |
| 0.1b | (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) |
4 calls. Slowest were:
| 29.0ms | (/ (/ (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (cbrt (sqrt 2))) |
| 16.0ms | (/ (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) |
| 9.0ms | (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) |
4 calls. Slowest were:
| 868.0ms | (/ (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) |
| 262.0ms | (/ (/ (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) (* (cbrt (sqrt 2)) (cbrt (sqrt 2)))) (cbrt (sqrt 2))) |
| 101.0ms | (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) |
| 28.0ms | (sqrt (fma a1 a1 (* a2 a2))) |
128 calls. Slowest were:
| 714.0ms | (/ (* (* (* (* (cos th) (cos th)) (cos th)) (* (* (sqrt (fma a1 a1 (* a2 a2))) (sqrt (fma a1 a1 (* a2 a2)))) (sqrt (fma a1 a1 (* a2 a2))))) (* (* (hypot a1 a2) (hypot a1 a2)) (hypot a1 a2))) (* (sqrt 2) (sqrt 2))) |
| 663.0ms | (- (- (log (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2))) (log (* (cbrt (sqrt 2)) (cbrt (sqrt 2))))) (log (cbrt (sqrt 2)))) |
| 565.0ms | (/ (/ (* (* (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2)) (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2))) (* (* (cos th) (sqrt (fma a1 a1 (* a2 a2)))) (hypot a1 a2))) (* (sqrt 2) (sqrt 2))) (sqrt 2)) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 0b
0% (0.4b remaining)
Error of 0.5b against oracle of 0.0b and baseline of 0.5b