


Bits error versus z
Initial program 1.8
Simplified0.7
rmApplied add-exp-log1.0
rmApplied frac-add1.2
Simplified1.2
Simplified1.2
rmApplied sqrt-prod0.9
rmApplied add-sqr-sqrt0.9
Applied exp-prod0.9
Simplified0.9
Final simplification0.9
herbie shell --seed 2019004 +o rules:numerics
(FPCore (z)
:name "Jmat.Real.gamma, branch z less than 0.5"
(* (/ PI (sin (* PI z))) (* (* (* (sqrt (* PI 2)) (pow (+ (+ (- (- 1 z) 1) 7) 0.5) (+ (- (- 1 z) 1) 0.5))) (exp (- (+ (+ (- (- 1 z) 1) 7) 0.5)))) (+ (+ (+ (+ (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (+ (- (- 1 z) 1) 1))) (/ -1259.1392167224028 (+ (- (- 1 z) 1) 2))) (/ 771.3234287776531 (+ (- (- 1 z) 1) 3))) (/ -176.6150291621406 (+ (- (- 1 z) 1) 4))) (/ 12.507343278686905 (+ (- (- 1 z) 1) 5))) (/ -0.13857109526572012 (+ (- (- 1 z) 1) 6))) (/ 9.984369578019572e-06 (+ (- (- 1 z) 1) 7))) (/ 1.5056327351493116e-07 (+ (- (- 1 z) 1) 8))))))
Time bar (total: 1.5m)Debug log
1 calls. Slowest were:
| 111.0ms | (* (/ PI (sin (* PI z))) (* (* (* (sqrt (* PI 2)) (pow (+ (+ (- (- 1 z) 1) 7) 0.5) (+ (- (- 1 z) 1) 0.5))) (exp (- (+ (+ (- (- 1 z) 1) 7) 0.5)))) (+ (+ (+ (+ (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (+ (- (- 1 z) 1) 1))) (/ -1259.1392167224028 (+ (- (- 1 z) 1) 2))) (/ 771.3234287776531 (+ (- (- 1 z) 1) 3))) (/ -176.6150291621406 (+ (- (- 1 z) 1) 4))) (/ 12.507343278686905 (+ (- (- 1 z) 1) 5))) (/ -0.13857109526572012 (+ (- (- 1 z) 1) 6))) (/ 9.984369578019572e-06 (+ (- (- 1 z) 1) 7))) (/ 1.5056327351493116e-07 (+ (- (- 1 z) 1) 8))))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 0.6b
Found 4 expressions with local error:
| 2.6b | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 1.0b | (sqrt (* 2 PI)) |
| 1.0b | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
| 1.0b | (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) |
4 calls. Slowest were:
| 14.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 13.0ms | (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) |
| 10.0ms | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
4 calls. Slowest were:
| 40.0ms | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
| 32.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 31.0ms | (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) |
| 0.0ms | (sqrt (* 2 PI)) |
35 calls. Slowest were:
| 1.6s | (+ (* -0.13857109526572012 (- 6 (+ z 1))) (* (+ (- 1 z) 5) 12.507343278686905)) |
| 675.0ms | (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) (+ (- 1 z) 1)) |
| 647.0ms | (+ (* (+ (pow 0.9999999999998099 3) (pow (/ 676.5203681218851 (- 1 z)) 3)) (+ (- 1 z) 1)) (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) -1259.1392167224028)) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0.2b
Found 4 expressions with local error:
| 2.6b | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 1.5b | (exp (log (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))))) |
| 1.0b | (sqrt (* 2 PI)) |
| 1.0b | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
4 calls. Slowest were:
| 14.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 12.0ms | (exp (log (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))))) |
| 10.0ms | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
4 calls. Slowest were:
| 32.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 32.0ms | (exp (log (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))))) |
| 24.0ms | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
| 0.0ms | (sqrt (* 2 PI)) |
39 calls. Slowest were:
| 767.0ms | (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) (+ (- 1 z) 1)) |
| 706.0ms | (exp (log (+ (* (+ (pow 0.9999999999998099 3) (pow (/ 676.5203681218851 (- 1 z)) 3)) (+ (- 1 z) 1)) (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) -1259.1392167224028)))) |
| 692.0ms | (+ (* (+ (pow 0.9999999999998099 3) (pow (/ 676.5203681218851 (- 1 z)) 3)) (+ (- 1 z) 1)) (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) -1259.1392167224028)) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 0.2b
Found 4 expressions with local error:
| 2.6b | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 1.5b | (exp (log (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))))) |
| 1.0b | (sqrt (* 2 PI)) |
| 1.0b | (/ (fma 771.3234287776531 (- 4 z) (* (- 3 z) -176.6150291621406)) (* (- 4 z) (- 3 z))) |
4 calls. Slowest were:
| 10.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 6.0ms | (exp (log (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))))) |
| 6.0ms | (/ (fma 771.3234287776531 (- 4 z) (* (- 3 z) -176.6150291621406)) (* (- 4 z) (- 3 z))) |
4 calls. Slowest were:
| 55.0ms | (/ (fma 771.3234287776531 (- 4 z) (* (- 3 z) -176.6150291621406)) (* (- 4 z) (- 3 z))) |
| 27.0ms | (exp (log (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))))) |
| 19.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 0.0ms | (sqrt (* 2 PI)) |
53 calls. Slowest were:
| 788.0ms | (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) (+ (- 1 z) 1)) |
| 746.0ms | (+ (* (+ (pow 0.9999999999998099 3) (pow (/ 676.5203681218851 (- 1 z)) 3)) (+ (- 1 z) 1)) (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) -1259.1392167224028)) |
| 731.0ms | (exp (log (+ (* (+ (pow 0.9999999999998099 3) (pow (/ 676.5203681218851 (- 1 z)) 3)) (+ (- 1 z) 1)) (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) -1259.1392167224028)))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 2.6b | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 1.5b | (exp (log (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))))) |
| 1.0b | (sqrt PI) |
| 1.0b | (/ (fma 771.3234287776531 (- 4 z) (* (- 3 z) -176.6150291621406)) (* (- 4 z) (- 3 z))) |
4 calls. Slowest were:
| 7.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 6.0ms | (exp (log (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))))) |
| 6.0ms | (/ (fma 771.3234287776531 (- 4 z) (* (- 3 z) -176.6150291621406)) (* (- 4 z) (- 3 z))) |
4 calls. Slowest were:
| 34.0ms | (/ (fma 771.3234287776531 (- 4 z) (* (- 3 z) -176.6150291621406)) (* (- 4 z) (- 3 z))) |
| 23.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 21.0ms | (exp (log (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))))) |
| 0.0ms | (sqrt PI) |
52 calls. Slowest were:
| 769.0ms | (exp (log (+ (* (+ (pow 0.9999999999998099 3) (pow (/ 676.5203681218851 (- 1 z)) 3)) (+ (- 1 z) 1)) (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) -1259.1392167224028)))) |
| 750.0ms | (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) (+ (- 1 z) 1)) |
| 599.0ms | (+ (* (+ (pow 0.9999999999998099 3) (pow (/ 676.5203681218851 (- 1 z)) 3)) (+ (- 1 z) 1)) (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) -1259.1392167224028)) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 0.1b
0% (0.7b remaining)
Error of 0.9b against oracle of 0.1b and baseline of 0.9b