


Bits error versus z
Initial program 1.8
Initial simplification0.9
Simplified0.9
Taylor expanded around 0 1.2
Simplified1.2
Final simplification1.2
herbie shell --seed 2018360 +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
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 | (* (* (sqrt (* 2 PI)) (pow (+ (- 8 (+ z 1)) 0.5) (- (- 1 z) (- 1 0.5)))) (/ 1 (exp (+ (- 8 (+ z 1)) 0.5)))) |
| 0.9b | (+ (+ (/ -0.13857109526572012 (- 7 (+ z 1))) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (+ (- 1 z) 7)) (/ 9.984369578019572e-06 (- 8 (+ z 1))))) |
4 calls. Slowest were:
| 71.0ms | (+ (+ (/ -0.13857109526572012 (- 7 (+ z 1))) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (+ (- 1 z) 7)) (/ 9.984369578019572e-06 (- 8 (+ z 1))))) |
| 66.0ms | (* (* (sqrt (* 2 PI)) (pow (+ (- 8 (+ z 1)) 0.5) (- (- 1 z) (- 1 0.5)))) (/ 1 (exp (+ (- 8 (+ z 1)) 0.5)))) |
| 14.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
4 calls. Slowest were:
| 187.0ms | (* (* (sqrt (* 2 PI)) (pow (+ (- 8 (+ z 1)) 0.5) (- (- 1 z) (- 1 0.5)))) (/ 1 (exp (+ (- 8 (+ z 1)) 0.5)))) |
| 55.0ms | (+ (+ (/ -0.13857109526572012 (- 7 (+ z 1))) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (+ (- 1 z) 7)) (/ 9.984369578019572e-06 (- 8 (+ z 1))))) |
| 22.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 0.0ms | (sqrt (* 2 PI)) |
71 calls. Slowest were:
| 1.0s | (+ (+ (/ -0.13857109526572012 (- 7 (+ z 1))) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (+ (- 1 z) 7)) (/ 9.984369578019572e-06 (- 8 (+ z 1))))) |
| 1.0s | (+ (+ (/ -0.13857109526572012 (- 7 (+ z 1))) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (+ (- 1 z) 7)) (/ 9.984369578019572e-06 (- 8 (+ z 1))))) |
| 1.0s | (+ (+ (/ -0.13857109526572012 (- 7 (+ z 1))) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (+ (- 1 z) 7)) (/ 9.984369578019572e-06 (- 8 (+ z 1))))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 1.0b | (sqrt (* 2 PI)) |
| 0.9b | (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) |
| 0.9b | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 0.9b | (* (* (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) (pow (+ (- 7 z) 0.5) (- 0 (- z 0.5)))) (fma (+ (+ (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) (/ -176.6150291621406 (- 4 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (+ (/ -1259.1392167224028 (- 2 z)) (/ 771.3234287776531 (- 3 z))))) (/ PI (sin (* z PI))) (* (+ (/ 9.984369578019572e-06 (- 7 z)) (/ 1.5056327351493116e-07 (- 8 z))) (/ PI (sin (* z PI)))))) |
4 calls. Slowest were:
| 64.0ms | (* (* (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) (pow (+ (- 7 z) 0.5) (- 0 (- z 0.5)))) (fma (+ (+ (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) (/ -176.6150291621406 (- 4 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (+ (/ -1259.1392167224028 (- 2 z)) (/ 771.3234287776531 (- 3 z))))) (/ PI (sin (* z PI))) (* (+ (/ 9.984369578019572e-06 (- 7 z)) (/ 1.5056327351493116e-07 (- 8 z))) (/ PI (sin (* z PI)))))) |
| 15.0ms | (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) |
| 5.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
4 calls. Slowest were:
| 387.0ms | (* (* (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) (pow (+ (- 7 z) 0.5) (- 0 (- z 0.5)))) (fma (+ (+ (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) (/ -176.6150291621406 (- 4 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (+ (/ -1259.1392167224028 (- 2 z)) (/ 771.3234287776531 (- 3 z))))) (/ PI (sin (* z PI))) (* (+ (/ 9.984369578019572e-06 (- 7 z)) (/ 1.5056327351493116e-07 (- 8 z))) (/ PI (sin (* z PI)))))) |
| 84.0ms | (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) |
| 33.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 0.0ms | (sqrt (* 2 PI)) |
84 calls. Slowest were:
| 561.0ms | (+ (+ (log (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5)))) (* (log (+ (- 7 z) 0.5)) (- 0 (- z 0.5)))) (log (fma (+ (+ (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) (/ -176.6150291621406 (- 4 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (+ (/ -1259.1392167224028 (- 2 z)) (/ 771.3234287776531 (- 3 z))))) (/ PI (sin (* z PI))) (* (+ (/ 9.984369578019572e-06 (- 7 z)) (/ 1.5056327351493116e-07 (- 8 z))) (/ PI (sin (* z PI))))))) |
| 543.0ms | (* (* (/ (* (* (sqrt (* 2 PI)) (sqrt (* 2 PI))) (sqrt (* 2 PI))) (* (* (exp (+ (- 7 z) 0.5)) (exp (+ (- 7 z) 0.5))) (exp (+ (- 7 z) 0.5)))) (* (* (pow (+ (- 7 z) 0.5) (- 0 (- z 0.5))) (pow (+ (- 7 z) 0.5) (- 0 (- z 0.5)))) (pow (+ (- 7 z) 0.5) (- 0 (- z 0.5))))) (* (* (fma (+ (+ (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) (/ -176.6150291621406 (- 4 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (+ (/ -1259.1392167224028 (- 2 z)) (/ 771.3234287776531 (- 3 z))))) (/ PI (sin (* z PI))) (* (+ (/ 9.984369578019572e-06 (- 7 z)) (/ 1.5056327351493116e-07 (- 8 z))) (/ PI (sin (* z PI))))) (fma (+ (+ (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) (/ -176.6150291621406 (- 4 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (+ (/ -1259.1392167224028 (- 2 z)) (/ 771.3234287776531 (- 3 z))))) (/ PI (sin (* z PI))) (* (+ (/ 9.984369578019572e-06 (- 7 z)) (/ 1.5056327351493116e-07 (- 8 z))) (/ PI (sin (* z PI)))))) (fma (+ (+ (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) (/ -176.6150291621406 (- 4 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (+ (/ -1259.1392167224028 (- 2 z)) (/ 771.3234287776531 (- 3 z))))) (/ PI (sin (* z PI))) (* (+ (/ 9.984369578019572e-06 (- 7 z)) (/ 1.5056327351493116e-07 (- 8 z))) (/ PI (sin (* z PI))))))) |
| 459.0ms | (+ (* 1/2 (* (/ (* (sqrt 2) (pow z 2)) (exp 7.5)) (sqrt PI))) (+ (* (/ (sqrt 2) (exp 7.5)) (sqrt PI)) (* (/ (* (sqrt 2) z) (exp 7.5)) (sqrt PI)))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 1.6b | (* (* 682.7204804053209 (/ (sqrt 2) (exp 7.5))) (* (pow (pow 7.5 1.0) 0.5) (sqrt PI))) |
| 1.0b | (/ (sqrt 2) (exp 7.5)) |
| 1.0b | (/ (sqrt 2) (exp 7.5)) |
| 1.0b | (* (pow (/ 1 (pow 7.5 3.0)) 0.5) (sqrt PI)) |
4 calls. Slowest were:
| 17.0ms | (* (* 682.7204804053209 (/ (sqrt 2) (exp 7.5))) (* (pow (pow 7.5 1.0) 0.5) (sqrt PI))) |
| 10.0ms | (* (pow (/ 1 (pow 7.5 3.0)) 0.5) (sqrt PI)) |
| 5.0ms | (/ (sqrt 2) (exp 7.5)) |
4 calls. Slowest were:
| 0.0ms | (* (* 682.7204804053209 (/ (sqrt 2) (exp 7.5))) (* (pow (pow 7.5 1.0) 0.5) (sqrt PI))) |
| 0.0ms | (* (pow (/ 1 (pow 7.5 3.0)) 0.5) (sqrt PI)) |
| 0.0ms | (/ (sqrt 2) (exp 7.5)) |
| 0.0ms | (/ (sqrt 2) (exp 7.5)) |
79 calls. Slowest were:
| 426.0ms | (+ (log (* 682.7204804053209 (/ (sqrt 2) (exp 7.5)))) (+ (log (pow (pow 7.5 1.0) 0.5)) (log (sqrt PI)))) |
| 410.0ms | (+ (log (* 682.7204804053209 (/ (sqrt 2) (exp 7.5)))) (log (* (pow (pow 7.5 1.0) 0.5) (sqrt PI)))) |
| 407.0ms | (+ (log (* 682.7204804053209 (/ (sqrt 2) (exp 7.5)))) (+ (* (log (pow 7.5 1.0)) 0.5) (log (sqrt PI)))) |
10 alts after pruning (9 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 1.0b | (sqrt (* 2 PI)) |
| 1.0b | (sqrt (* 2 PI)) |
| 0.9b | (/ (sqrt (* 2 PI)) (exp (+ 7 (- 0.5 z)))) |
| 0.9b | (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) |
4 calls. Slowest were:
| 22.0ms | (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) |
| 11.0ms | (/ (sqrt (* 2 PI)) (exp (+ 7 (- 0.5 z)))) |
| 1.0ms | (sqrt (* 2 PI)) |
4 calls. Slowest were:
| 93.0ms | (/ (sqrt (* 2 PI)) (exp (+ 7 (- 0.5 z)))) |
| 77.0ms | (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) |
| 0.0ms | (sqrt (* 2 PI)) |
| 0.0ms | (sqrt (* 2 PI)) |
94 calls. Slowest were:
| 668.0ms | (+ (* (* (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) (pow (+ (- 7 z) 0.5) (- 0 (- z 0.5)))) (* (+ (+ (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) (/ -176.6150291621406 (- 4 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (+ (/ -1259.1392167224028 (- 2 z)) (/ 771.3234287776531 (- 3 z))))) (/ PI (sin (* z PI))))) (* (* (/ PI (sin (* z PI))) (/ (sqrt (* 2 PI)) (exp (+ 7 (- 0.5 z))))) (* (pow (+ 7 (- 0.5 z)) (- 0.5 z)) (+ (/ 9.984369578019572e-06 (- 7 z)) (/ 1.5056327351493116e-07 (- 8 z)))))) |
| 484.0ms | (+ (* (* (/ (sqrt (* 2 PI)) (exp (+ (- 7 z) 0.5))) (pow (+ (- 7 z) 0.5) (- 0 (- z 0.5)))) (* (+ (+ (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) (/ -176.6150291621406 (- 4 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (+ (/ -1259.1392167224028 (- 2 z)) (/ 771.3234287776531 (- 3 z))))) (/ PI (sin (* z PI))))) (* (* (/ PI (sin (* z PI))) (/ (sqrt (* 2 PI)) (exp (+ 7 (- 0.5 z))))) (* (pow (+ 7 (- 0.5 z)) (- 0.5 z)) (+ (/ 9.984369578019572e-06 (- 7 z)) (/ 1.5056327351493116e-07 (- 8 z)))))) |
| 366.0ms | (+ (* 1/2 (* (/ (* (sqrt 2) (pow z 2)) (exp 7.5)) (sqrt PI))) (+ (* (/ (sqrt 2) (exp 7.5)) (sqrt PI)) (* (/ (* (sqrt 2) z) (exp 7.5)) (sqrt PI)))) |
10 alts after pruning (8 fresh and 2 done)
Merged error: 0.0b
0% (1.1b remaining)
Error of 1.2b against oracle of 0.1b and baseline of 1.2b