


Bits error versus z
Results
Initial program 1.8
Simplified0.6
Simplified0.6
rmApplied exp-sum0.6
Applied sub-neg0.6
Applied unpow-prod-up0.6
Applied times-frac0.6
Simplified0.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2019008 +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: 3.0m)Debug log
1 calls. Slowest were:
| 121.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 | (* (pow (- (- 1 z) (- -6 0.5)) (- (- 1 z) (- 1 0.5))) (/ 1 (exp (- (- 1 z) (- -6 0.5))))) |
4 calls. Slowest were:
| 29.0ms | (* (pow (- (- 1 z) (- -6 0.5)) (- (- 1 z) (- 1 0.5))) (/ 1 (exp (- (- 1 z) (- -6 0.5))))) |
| 9.0ms | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
| 8.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
4 calls. Slowest were:
| 84.0ms | (* (pow (- (- 1 z) (- -6 0.5)) (- (- 1 z) (- 1 0.5))) (/ 1 (exp (- (- 1 z) (- -6 0.5))))) |
| 31.0ms | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
| 20.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 0.0ms | (sqrt (* 2 PI)) |
65 calls. Slowest were:
| 1.2s | (pow (+ (sqrt (- 1 z)) (sqrt (- -6 0.5))) (- (- 1 z) (- 1 0.5))) |
| 936.0ms | (pow (- (- 1 z) (- -6 0.5)) (fma (sqrt (- 1 z)) (sqrt (- 1 z)) (- (* (sqrt (- 1 0.5)) (sqrt (- 1 0.5)))))) |
| 726.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.1b
Found 4 expressions with local error:
| 1.0b | (sqrt (* 2 PI)) |
| 1.0b | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 1.0b | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 1.0b | (/ (pow (+ (- 1 z) (+ 6 0.5)) (- (- 1 z) (- 1 0.5))) (exp (+ (- 1 z) (+ 6 0.5)))) |
4 calls. Slowest were:
| 36.0ms | (/ (pow (+ (- 1 z) (+ 6 0.5)) (- (- 1 z) (- 1 0.5))) (exp (+ (- 1 z) (+ 6 0.5)))) |
| 8.0ms | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 5.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
4 calls. Slowest were:
| 104.0ms | (/ (pow (+ (- 1 z) (+ 6 0.5)) (- (- 1 z) (- 1 0.5))) (exp (+ (- 1 z) (+ 6 0.5)))) |
| 28.0ms | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 25.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 0.0ms | (sqrt (* 2 PI)) |
225 calls. Slowest were:
| 1.3s | (/ (pow (sqrt (+ (- 1 z) (+ 6 0.5))) (- (- 1 z) (- 1 0.5))) (exp (- 1 z))) |
| 1.3s | (/ (pow 1 (- (- 1 z) (- 1 0.5))) 1) |
| 1.2s | (pow (* (cbrt (+ (- 1 z) (+ 6 0.5))) (cbrt (+ (- 1 z) (+ 6 0.5)))) (- (- 1 z) (- 1 0.5))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 1.0b | (sqrt (* 2 PI)) |
| 1.0b | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 1.0b | (/ (pow (+ 0.5 (- 7 z)) (- 1 z)) (/ E (exp z))) |
| 1.0b | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
4 calls. Slowest were:
| 12.0ms | (/ (pow (+ 0.5 (- 7 z)) (- 1 z)) (/ E (exp z))) |
| 10.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 6.0ms | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
4 calls. Slowest were:
| 95.0ms | (/ (pow (+ 0.5 (- 7 z)) (- 1 z)) (/ E (exp z))) |
| 36.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 26.0ms | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 0.0ms | (sqrt (* 2 PI)) |
108 calls. Slowest were:
| 573.0ms | (- (* (log (+ 0.5 (- 7 z))) (- 1 z)) (log (/ E (exp z)))) |
| 537.0ms | (- (* (log (+ 0.5 (- 7 z))) (- 1 z)) (log (/ E (exp z)))) |
| 475.0ms | (* (+ (+ (/ 1.5056327351493116e-07 (- 8 z)) (/ 9.984369578019572e-06 (- (- 1 z) -6))) (+ (+ (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) (+ (/ 676.5203681218851 (- 1 z)) 0.9999999999998099)) (+ (/ -1259.1392167224028 (- 2 z)) (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z)))))) (* (* (/ (pow (+ 0.5 (- 7 z)) (- 1 z)) (/ E (exp z))) (/ (pow (+ (- 7 z) 0.5) (+ -1 0.5)) (exp (+ 6 0.5)))) (* (sqrt (* 2 PI)) (/ PI (sin (* PI z)))))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 1.0b | (sqrt (* 2 PI)) |
| 1.0b | (* (cbrt (+ (- 7 z) 0.5)) (cbrt (+ (- 7 z) 0.5))) |
| 1.0b | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 1.0b | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
4 calls. Slowest were:
| 17.0ms | (* (cbrt (+ (- 7 z) 0.5)) (cbrt (+ (- 7 z) 0.5))) |
| 16.0ms | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 10.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
4 calls. Slowest were:
| 47.0ms | (* (cbrt (+ (- 7 z) 0.5)) (cbrt (+ (- 7 z) 0.5))) |
| 34.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 33.0ms | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 0.0ms | (sqrt (* 2 PI)) |
59 calls. Slowest were:
| 640.0ms | (* (cbrt (- (* (- 7 z) (- 7 z)) (* 0.5 0.5))) (cbrt (- (* (- 7 z) (- 7 z)) (* 0.5 0.5)))) |
| 567.0ms | (cbrt (+ (* (- 7 z) (- 7 z)) (- (* 0.5 0.5) (* (- 7 z) 0.5)))) |
| 545.0ms | (* (cbrt (+ (* (- 7 z) (- 7 z)) (- (* 0.5 0.5) (* (- 7 z) 0.5)))) (cbrt (+ (* (- 7 z) (- 7 z)) (- (* 0.5 0.5) (* (- 7 z) 0.5))))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 0.0b
0% (0.5b remaining)
Error of 0.6b against oracle of 0.1b and baseline of 0.6b