


Bits error versus d



Bits error versus h



Bits error versus l



Bits error versus M



Bits error versus D
Results
if l < -2.1332837374323508e-27Initial program 24.8
rmApplied associate-*r/26.2
rmApplied add-cube-cbrt26.2
Applied associate-*r*26.2
Taylor expanded around -inf 22.4
Simplified19.1
if -2.1332837374323508e-27 < l < -9.2206008893997e-311Initial program 27.3
rmApplied associate-*r/22.0
Taylor expanded around -inf 16.6
Simplified12.6
if -9.2206008893997e-311 < l < 6.506805187540822e-39Initial program 28.1
rmApplied associate-*r/22.6
rmApplied div-inv22.6
Applied unpow-prod-down12.2
Simplified12.2
if 6.506805187540822e-39 < l Initial program 24.0
rmApplied div-inv24.0
Applied unpow-prod-down17.3
Simplified17.3
Final simplification16.4
herbie shell --seed 2019004
(FPCore (d h l M D)
:name "Henrywood and Agarwal, Equation (12)"
(* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)))))
Time bar (total: 1.1m)Debug log
1 calls. Slowest were:
| 148.0ms | (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 26.1b
Found 4 expressions with local error:
| 12.9b | (pow (/ d l) (/ 1 2)) |
| 12.8b | (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)) |
| 12.4b | (pow (/ d h) (/ 1 2)) |
| 7.5b | (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)))) |
4 calls. Slowest were:
| 16.0ms | (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)))) |
| 13.0ms | (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)) |
| 3.0ms | (pow (/ d l) (/ 1 2)) |
4 calls. Slowest were:
| 169.0ms | (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)))) |
| 63.0ms | (pow (/ d l) (/ 1 2)) |
| 61.0ms | (pow (/ d h) (/ 1 2)) |
| 53.0ms | (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)) |
71 calls. Slowest were:
| 718.0ms | (+ (+ (log (/ 1 2)) (log (pow (/ (* M D) (* 2 d)) 2))) (log (/ h l))) |
| 628.0ms | (* (* 1 (pow (/ (* M D) (* 2 d)) 2)) h) |
| 367.0ms | (* (* (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2))) (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2))) (* (* (/ h l) (/ h l)) (/ h l))) |
18 alts after pruning (18 fresh and 0 done)
Merged error: 12.1b
Found 4 expressions with local error:
| 12.9b | (pow (/ d l) (/ 1 2)) |
| 12.4b | (pow (/ d h) (/ 1 2)) |
| 7.5b | (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l))) |
| 7.1b | (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l) |
4 calls. Slowest were:
| 16.0ms | (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l))) |
| 12.0ms | (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l) |
| 1.0ms | (pow (/ d l) (/ 1 2)) |
4 calls. Slowest were:
| 257.0ms | (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l))) |
| 86.0ms | (pow (/ d l) (/ 1 2)) |
| 67.0ms | (pow (/ d h) (/ 1 2)) |
| 53.0ms | (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l) |
57 calls. Slowest were:
| 751.0ms | (* (* 1 (pow (/ (* M D) (* 2 d)) 2)) h) |
| 591.0ms | (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) |
| 195.0ms | (+ (+ (* (log (/ d h)) (/ 1 2)) (* (log (/ d l)) (/ 1 2))) (log (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l)))) |
17 alts after pruning (17 fresh and 0 done)
Merged error: 9.6b
Found 4 expressions with local error:
| 12.9b | (pow (/ d l) (/ 1 2)) |
| 12.4b | (pow (/ d h) (/ 1 2)) |
| 7.1b | (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l) |
| 7.1b | (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l) |
4 calls. Slowest were:
| 15.0ms | (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l) |
| 12.0ms | (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l) |
| 2.0ms | (pow (/ d l) (/ 1 2)) |
4 calls. Slowest were:
| 101.0ms | (pow (/ d l) (/ 1 2)) |
| 91.0ms | (pow (/ d h) (/ 1 2)) |
| 69.0ms | (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l) |
| 61.0ms | (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l) |
36 calls. Slowest were:
| 984.0ms | (* (* 1 (pow (/ (* M D) (* 2 d)) 2)) h) |
| 886.0ms | (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) |
| 773.0ms | (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) |
17 alts after pruning (17 fresh and 0 done)
Merged error: 9.6b
Found 4 expressions with local error:
| 12.9b | (pow (/ d l) (/ 1 2)) |
| 12.4b | (pow (/ d h) (/ 1 2)) |
| 12.4b | (pow (/ d h) (/ 1 2)) |
| 7.5b | (* (* (* (sqrt (pow (/ d h) (/ 1 2))) (sqrt (pow (/ d h) (/ 1 2)))) (pow (/ d l) (/ 1 2))) (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l))) |
4 calls. Slowest were:
| 53.0ms | (* (* (* (sqrt (pow (/ d h) (/ 1 2))) (sqrt (pow (/ d h) (/ 1 2)))) (pow (/ d l) (/ 1 2))) (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l))) |
| 2.0ms | (pow (/ d l) (/ 1 2)) |
| 1.0ms | (pow (/ d h) (/ 1 2)) |
4 calls. Slowest were:
| 243.0ms | (* (* (* (sqrt (pow (/ d h) (/ 1 2))) (sqrt (pow (/ d h) (/ 1 2)))) (pow (/ d l) (/ 1 2))) (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l))) |
| 89.0ms | (pow (/ d h) (/ 1 2)) |
| 74.0ms | (pow (/ d h) (/ 1 2)) |
| 72.0ms | (pow (/ d l) (/ 1 2)) |
63 calls. Slowest were:
| 298.0ms | (* (* (* (* (* (sqrt (pow (/ d h) (/ 1 2))) (sqrt (pow (/ d h) (/ 1 2)))) (* (sqrt (pow (/ d h) (/ 1 2))) (sqrt (pow (/ d h) (/ 1 2))))) (* (sqrt (pow (/ d h) (/ 1 2))) (sqrt (pow (/ d h) (/ 1 2))))) (* (* (pow (/ d l) (/ 1 2)) (pow (/ d l) (/ 1 2))) (pow (/ d l) (/ 1 2)))) (* (* (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l)) (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l))) (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l)))) |
| 298.0ms | (* (* (* (* (* (sqrt (pow (/ d h) (/ 1 2))) (sqrt (pow (/ d h) (/ 1 2)))) (sqrt (pow (/ d h) (/ 1 2)))) (* (* (sqrt (pow (/ d h) (/ 1 2))) (sqrt (pow (/ d h) (/ 1 2)))) (sqrt (pow (/ d h) (/ 1 2))))) (* (* (pow (/ d l) (/ 1 2)) (pow (/ d l) (/ 1 2))) (pow (/ d l) (/ 1 2)))) (* (* (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l)) (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l))) (- 1 (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l)))) |
| 218.0ms | (* (- (/ (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) h) l)) (* (* (sqrt (pow (/ d h) (/ 1 2))) (sqrt (pow (/ d h) (/ 1 2)))) (pow (/ d l) (/ 1 2)))) |
18 alts after pruning (18 fresh and 0 done)
Merged error: 9.6b
53% (7.3b remaining)
Error of 16.4b against oracle of 9.1b and baseline of 24.7b