


Bits error versus c0



Bits error versus w



Bits error versus h



Bits error versus D



Bits error versus d



Bits error versus M
Results
Initial program 57.8
Initial simplification55.1
Taylor expanded around inf 33.6
Final simplification33.6
herbie shell --seed 2018362 +o rules:numerics
(FPCore (c0 w h D d M)
:name "Henrywood and Agarwal, Equation (13)"
(* (/ c0 (* 2 w)) (+ (/ (* c0 (* d d)) (* (* w h) (* D D))) (sqrt (- (* (/ (* c0 (* d d)) (* (* w h) (* D D))) (/ (* c0 (* d d)) (* (* w h) (* D D)))) (* M M))))))
Time bar (total: 41.5s)Debug log
2 alts after pruning (2 fresh and 0 done)
Merged error: 53.5b
Found 4 expressions with local error:
| 46.9b | (fma (/ c0 (* w 2)) (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* M (- M)))) (* (/ c0 (* w 2)) (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)))) |
| 13.7b | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 13.7b | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 13.7b | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
4 calls. Slowest were:
| 13.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 8.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 8.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
4 calls. Slowest were:
| 434.0ms | (fma (/ c0 (* w 2)) (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* M (- M)))) (* (/ c0 (* w 2)) (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)))) |
| 63.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 60.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 52.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
72 calls. Slowest were:
| 479.0ms | (/ (* c0 (pow d 2)) (* w (* (pow D 2) h))) |
| 446.0ms | (/ (* c0 (pow d 2)) (* w (* (pow D 2) h))) |
| 429.0ms | (* (* (* (* (/ d D) (/ d D)) (* (/ d D) (/ d D))) (* (/ d D) (/ d D))) (* (* (/ (/ c0 h) w) (/ (/ c0 h) w)) (/ (/ c0 h) w))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 29.2b
Found 0 expressions with local error:
0 calls. Slowest were:
0 calls. Slowest were:
0 calls. Slowest were:
10 alts after pruning (9 fresh and 1 done)
Merged error: 29.2b
Found 4 expressions with local error:
| 46.5b | (fma (/ c0 (* w 2)) (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M)))) (* (/ c0 (* w 2)) (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)))) |
| 13.7b | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 13.7b | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 13.2b | (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M)))) |
4 calls. Slowest were:
| 18.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 17.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 2.0ms | (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M)))) |
4 calls. Slowest were:
| 439.0ms | (fma (/ c0 (* w 2)) (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M)))) (* (/ c0 (* w 2)) (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)))) |
| 248.0ms | (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M)))) |
| 46.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 46.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
56 calls. Slowest were:
| 396.0ms | (+ (+ (log (/ d D)) (log (/ d D))) (log (/ (/ c0 h) w))) |
| 394.0ms | (/ (* c0 (pow d 2)) (* w (* (pow D 2) h))) |
| 367.0ms | (* (* d (/ d D)) (/ (/ c0 h) w)) |
10 alts after pruning (8 fresh and 2 done)
Merged error: 29.1b
Found 4 expressions with local error:
| 29.0b | (+ (* (/ c0 (* w 2)) (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M))))) (* (/ c0 (* w 2)) (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)))) |
| 13.7b | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 13.7b | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 13.2b | (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M)))) |
4 calls. Slowest were:
| 12.0ms | (+ (* (/ c0 (* w 2)) (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M))))) (* (/ c0 (* w 2)) (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)))) |
| 9.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 8.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
4 calls. Slowest were:
| 509.0ms | (+ (* (/ c0 (* w 2)) (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M))))) (* (/ c0 (* w 2)) (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)))) |
| 176.0ms | (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M)))) |
| 69.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
| 59.0ms | (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) |
90 calls. Slowest were:
| 589.0ms | (+ (* (* c0 (sqrt (fma (* (* (/ d D) (/ d D)) (/ (/ c0 h) w)) (* (/ d D) (* (/ d D) (/ (/ c0 h) w))) (* M (- M))))) (* D D)) (* (* w 2) (* (/ c0 (* w 2)) (* (* d d) (/ (/ c0 h) w))))) |
| 424.0ms | (/ (* c0 (pow d 2)) (* w (* (pow D 2) h))) |
| 401.0ms | (+ (+ (log (/ d D)) (log (/ d D))) (log (/ (/ c0 h) w))) |
14 alts after pruning (12 fresh and 2 done)
Merged error: 28.2b
0% (6.6b remaining)
Error of 33.6b against oracle of 27.0b and baseline of 33.6b