


Bits error versus z
Results
Initial program 1.8
Simplified0.7
Simplified0.7
rmApplied exp-sum0.7
Applied *-un-lft-identity0.7
Applied times-frac0.6
Simplified0.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2019022 +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
| 589.0ms | 110× | body | 10240 | exit |
| 148.0ms | 42× | body | 1280 | nan |
| 137.0ms | 259× | body | 80 | valid |
| 106.0ms | 41× | body | 640 | nan |
| 35.0ms | 21× | body | 320 | nan |
| 11.0ms | 8× | body | 160 | nan |
| 11.0ms | 17× | body | 80 | nan |
| 131.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.7b
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 | (+ (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 1.0b | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
| 25× | *-un-lft-identity |
| 20× | add-log-exp |
| 18× | frac-add |
| 13× | fma-def |
| 11× | distribute-lft-out |
| 10× | flip-+ |
| 10× | flip3-+ |
| 10× | sum-log |
| 7× | add-cube-cbrt |
| 7× | add-sqr-sqrt |
| 4× | log1p-expm1-u |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 4× | expm1-log1p-u |
| 3× | +-commutative |
| 2× | associate-/r/ |
| 2× | associate-+l+ |
| 1× | flip-- |
| 1× | div-inv |
| 1× | sqrt-prod |
| 1× | pow1/2 |
| 1× | flip3-- |
| 1× | associate-+r+ |
| 41.0ms | (+ (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 13.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 9.0ms | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
| 1.0ms | (sqrt (* 2 PI)) |
| 69.0ms | (+ (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 35.0ms | (+ (/ 771.3234287776531 (- 1 (+ z -2))) (/ -176.6150291621406 (- 5 (+ z 1)))) |
| 33.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 0.0ms | (sqrt (* 2 PI)) |
| 1.1s | (* (* (+ (- 1 z) 5) (- 6 (+ z 1))) (- (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 804.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)) |
| 753.0ms | (* (- (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) (* (- 9 (+ z 1)) (- (- 1 z) -6))) |
| 719.0ms | (* (+ (* 0.9999999999998099 0.9999999999998099) (- (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (* 0.9999999999998099 (/ 676.5203681218851 (- 1 z))))) (+ (- 1 z) 1)) |
| 705.0ms | (* (exp (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1))))) (exp (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6))))) |
10 alts after pruning (10 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.0b | (sqrt (* 2 PI)) |
| 1.0b | (+ (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 1.0b | (+ (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 37× | *-un-lft-identity |
| 32× | frac-add |
| 30× | add-log-exp |
| 19× | distribute-lft-out |
| 17× | sum-log |
| 16× | flip-+ |
| 16× | flip3-+ |
| 11× | fma-def |
| 7× | add-cube-cbrt |
| 7× | add-sqr-sqrt |
| 4× | log1p-expm1-u |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 4× | expm1-log1p-u |
| 3× | associate-+l+ |
| 3× | +-commutative |
| 2× | associate-+r+ |
| 1× | sqrt-prod |
| 1× | pow1/2 |
| 63.0ms | (+ (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 54.0ms | (+ (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 6.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 1.0ms | (sqrt (* 2 PI)) |
| 52.0ms | (+ (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 50.0ms | (+ (+ (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 29.0ms | (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1 z))) (/ -1259.1392167224028 (+ (- 1 z) 1))) |
| 0.0ms | (sqrt (* 2 PI)) |
| 796.0ms | (* (* (exp (/ -0.13857109526572012 (+ (- 1 z) 5))) (exp (/ 12.507343278686905 (- 6 (+ z 1))))) (* (exp (/ 1.5056327351493116e-07 (- 9 (+ z 1)))) (exp (/ 9.984369578019572e-06 (- (- 1 z) -6))))) |
| 780.0ms | (* (* (exp (/ -0.13857109526572012 (+ (- 1 z) 5))) (exp (/ 12.507343278686905 (- 6 (+ z 1))))) (* (exp (/ 1.5056327351493116e-07 (- 9 (+ z 1)))) (exp (/ 9.984369578019572e-06 (- (- 1 z) -6))))) |
| 765.0ms | (* (* (+ (- 1 z) 5) (- 6 (+ z 1))) (- (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 747.0ms | (* (* (+ (- 1 z) 5) (- 6 (+ z 1))) (- (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))) |
| 743.0ms | (+ (* (+ (pow (/ -0.13857109526572012 (+ (- 1 z) 5)) 3) (pow (/ 12.507343278686905 (- 6 (+ z 1))) 3)) (+ (* (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 1.5056327351493116e-07 (- 9 (+ z 1)))) (- (* (/ 9.984369578019572e-06 (- (- 1 z) -6)) (/ 9.984369578019572e-06 (- (- 1 z) -6))) (* (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6)))))) (* (+ (* (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ -0.13857109526572012 (+ (- 1 z) 5))) (- (* (/ 12.507343278686905 (- 6 (+ z 1))) (/ 12.507343278686905 (- 6 (+ z 1)))) (* (/ -0.13857109526572012 (+ (- 1 z) 5)) (/ 12.507343278686905 (- 6 (+ z 1)))))) (+ (pow (/ 1.5056327351493116e-07 (- 9 (+ z 1))) 3) (pow (/ 9.984369578019572e-06 (- (- 1 z) -6)) 3)))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0.2b
Found 4 expressions with local error:
| 1.0b | (sqrt (* 2 PI)) |
| 1.0b | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 1.0b | (/ (pow (+ (- 1 z) (+ 6 0.5)) (- (- 1 z) (- 1 0.5))) (exp (+ (- 1 z) (+ 6 0.5)))) |
| 1.0b | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 84× | times-frac |
| 82× | *-un-lft-identity |
| 70× | unpow-prod-up |
| 68× | add-cube-cbrt |
| 68× | add-sqr-sqrt |
| 45× | prod-diff |
| 22× | exp-sum |
| 21× | associate-/l* |
| 20× | unpow-prod-down |
| 12× | fma-def |
| 10× | associate--l+ |
| 10× | sub-neg |
| 8× | add-log-exp |
| 7× | distribute-lft-out |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 5× | associate--r- |
| 5× | associate-/r/ |
| 4× | log1p-expm1-u |
| 4× | associate-/r* |
| 4× | pow1 |
| 4× | expm1-log1p-u |
| 3× | div-inv |
| 3× | div-exp |
| 2× | flip-- |
| 2× | flip-+ |
| 2× | frac-add |
| 2× | flip3-- |
| 2× | flip3-+ |
| 2× | sum-log |
| 2× | +-commutative |
| 1× | associate-/l/ |
| 1× | pow-sub |
| 1× | sqrt-prod |
| 1× | associate-+l- |
| 1× | pow-exp |
| 1× | exp-diff |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | cbrt-undiv |
| 1× | pow-to-exp |
| 41.0ms | (/ (pow (+ (- 1 z) (+ 6 0.5)) (- (- 1 z) (- 1 0.5))) (exp (+ (- 1 z) (+ 6 0.5)))) |
| 10.0ms | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 5.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 1.0ms | (sqrt (* 2 PI)) |
| 76.0ms | (/ (pow (+ (- 1 z) (+ 6 0.5)) (- (- 1 z) (- 1 0.5))) (exp (+ (- 1 z) (+ 6 0.5)))) |
| 38.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 28.0ms | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 0.0ms | (sqrt (* 2 PI)) |
| 1.2s | (/ (pow 1 (- (- 1 z) (- 1 0.5))) 1) |
| 1.2s | (/ (pow 1 (- (- 1 z) (- 1 0.5))) 1) |
| 1.1s | (/ (sqrt (pow (+ (- 1 z) (+ 6 0.5)) (- (- 1 z) (- 1 0.5)))) (exp (+ 6 0.5))) |
| 1.1s | (/ (pow (sqrt (+ (- 1 z) (+ 6 0.5))) (- (- 1 z) (- 1 0.5))) (exp (- 1 z))) |
| 966.0ms | (/ (* (cbrt (pow (+ (- 1 z) (+ 6 0.5)) (- (- 1 z) (- 1 0.5)))) (cbrt (pow (+ (- 1 z) (+ 6 0.5)) (- (- 1 z) (- 1 0.5))))) 1) |
9 alts after pruning (9 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 | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 1.0b | (+ (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) (+ (/ 676.5203681218851 (- 1 z)) 0.9999999999998099)) |
| 18× | *-un-lft-identity |
| 16× | fma-def |
| 13× | add-log-exp |
| 10× | frac-add |
| 8× | flip-+ |
| 8× | flip3-+ |
| 7× | add-cube-cbrt |
| 7× | add-sqr-sqrt |
| 6× | distribute-lft-out |
| 5× | sum-log |
| 4× | log1p-expm1-u |
| 4× | associate-/r/ |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 4× | expm1-log1p-u |
| 3× | +-commutative |
| 2× | flip-- |
| 2× | div-inv |
| 2× | flip3-- |
| 1× | sqrt-prod |
| 1× | pow1/2 |
| 1× | associate-+r+ |
| 1× | associate-+l+ |
| 63.0ms | (+ (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) (+ (/ 676.5203681218851 (- 1 z)) 0.9999999999998099)) |
| 9.0ms | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 9.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 1.0ms | (sqrt (* 2 PI)) |
| 36.0ms | (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) |
| 35.0ms | (+ (/ -0.13857109526572012 (- 6 z)) (/ 12.507343278686905 (- 5 z))) |
| 34.0ms | (+ (+ (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) (+ (/ 676.5203681218851 (- 1 z)) 0.9999999999998099)) |
| 0.0ms | (sqrt (* 2 PI)) |
| 453.0ms | (- (+ (* 4116.070392404628 (/ 1 (pow z 3))) (+ (* 1607.5101696843972 (/ 1 (pow z 2))) (* 594.7083996155126 (/ 1 z))))) |
| 407.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)))))) (* (* (exp (+ -1 z)) (/ (pow (+ (- 7 z) 0.5) (- 0.5 z)) (exp (+ 6 0.5)))) (* (sqrt (* 2 PI)) (/ PI (sin (* PI z)))))) |
| 349.0ms | (- (+ (* 4116.070392404628 (/ 1 (pow z 3))) (+ (* 1607.5101696843972 (/ 1 (pow z 2))) (* 594.7083996155126 (/ 1 z))))) |
| 335.0ms | (* (- (/ -176.6150291621406 (- 4 z)) (/ 771.3234287776531 (- 3 z))) (+ (* (/ 676.5203681218851 (- 1 z)) (/ 676.5203681218851 (- 1 z))) (- (* 0.9999999999998099 0.9999999999998099) (* (/ 676.5203681218851 (- 1 z)) 0.9999999999998099)))) |
| 293.0ms | (- (+ (* 307.69502253760675 (/ 1 (pow z 3))) (+ (* 61.7052898218402 (/ 1 (pow z 2))) (* 12.368772183421186 (/ 1 z))))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 0.0b
0% (0.6b remaining)
Error of 0.6b against oracle of 0.1b and baseline of 0.6b
| 22.8s | 3832× | body | 10240 | exit |
| 6.6s | 1504× | body | 1280 | nan |
| 6.0s | 8099× | body | 80 | valid |
| 3.4s | 1206× | body | 640 | nan |
| 1.3s | 643× | body | 320 | nan |
| 379.0ms | 319× | body | 160 | nan |
| 231.0ms | 301× | body | 80 | nan |
| 1.0ms | 1× | body | 160 | valid |