


Bits error versus x
Initial program 1.5
Simplified1.1
rmApplied inv-pow1.1
Applied pow-pow0.6
rmApplied pow10.6
Applied pow-flip0.6
Applied inv-pow0.6
Applied pow-prod-up0.6
Simplified0.6
Simplified1.0
Taylor expanded around 0 0.5
Simplified0.5
Final simplification0.5
herbie shell --seed 2019010 +o rules:numerics
(FPCore (x)
:name "Jmat.Real.erfi, branch x greater than or equal to 5"
(* (* (/ 1 (sqrt PI)) (exp (* (fabs x) (fabs x)))) (+ (+ (+ (/ 1 (fabs x)) (* (/ 1 2) (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))))) (* (/ 3 4) (* (* (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))))) (* (/ 15 8) (* (* (* (* (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x)))))))
Time bar (total: 1.2m)Debug log
| 643.0ms | (* (* (/ 1 (sqrt PI)) (exp (* (fabs x) (fabs x)))) (+ (+ (+ (/ 1 (fabs x)) (* (/ 1 2) (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))))) (* (/ 3 4) (* (* (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))))) (* (/ 15 8) (* (* (* (* (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x)))))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 0.9b
Found 4 expressions with local error:
| 1.0b | (sqrt PI) |
| 1.0b | (sqrt PI) |
| 0.9b | (pow (/ 1 (fabs x)) 4) |
| 0.5b | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 17× | pow1 |
| 14× | add-exp-log |
| 7× | inv-pow |
| 7× | pow-flip |
| 7× | add-cube-cbrt |
| 7× | *-un-lft-identity |
| 7× | add-sqr-sqrt |
| 6× | add-cbrt-cube |
| 5× | pow-prod-up |
| 5× | rec-exp |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | prod-exp |
| 4× | associate-*l* |
| 4× | associate-*r* |
| 4× | unpow-prod-down |
| 4× | expm1-log1p-u |
| 3× | div-inv |
| 3× | pow-prod-down |
| 3× | pow-pow |
| 2× | pow-exp |
| 2× | pow1/2 |
| 1× | cbrt-unprod |
| 1× | *-commutative |
| 1× | un-div-inv |
| 1× | associate-*r/ |
| 1× | associate-*l/ |
| 1× | pow-plus |
| 1× | pow-to-exp |
| 1× | frac-times |
| 1× | pow2 |
| 9.0ms | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 3.0ms | (pow (/ 1 (fabs x)) 4) |
| 0.0ms | (sqrt PI) |
| 0.0ms | (sqrt PI) |
| 32.0ms | (pow (/ 1 (fabs x)) 4) |
| 31.0ms | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 0.0ms | (sqrt PI) |
| 0.0ms | (sqrt PI) |
| 217.0ms | (fma (/ (exp (* (fabs x) (fabs x))) (sqrt PI)) (fma (* (* 3/4 (/ 1 (fabs x))) (* (/ 1 (fabs x)) (/ 1 (fabs x)))) (* (/ 1 (fabs x)) (/ 1 (fabs x))) (fma (* 1/2 (/ 1 (fabs x))) (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x)))) (* (* (/ 1 (pow (fabs x) 3)) (pow (/ 1 (fabs x)) 4)) (* (exp (* (fabs x) (fabs x))) (* (/ 1 (sqrt PI)) 15/8)))) |
| 195.0ms | (* (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x)))) |
| 136.0ms | (fma (/ (exp (* (fabs x) (fabs x))) (sqrt PI)) (fma (* (* 3/4 (/ 1 (fabs x))) (* (/ 1 (fabs x)) (/ 1 (fabs x)))) (* (/ 1 (fabs x)) (/ 1 (fabs x))) (fma (* 1/2 (/ 1 (fabs x))) (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x)))) (* (* (/ 1 (pow (fabs x) 3)) (pow (/ 1 (fabs x)) 4)) (* (exp (* (fabs x) (fabs x))) (* (/ 1 (sqrt PI)) 15/8)))) |
| 26.0ms | (+ (log (/ 1 (fabs x))) (log (/ 1 (fabs x)))) |
| 20.0ms | (+ (- (log (fabs x))) (log (/ 1 (fabs x)))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 1.0b | (sqrt PI) |
| 1.0b | (sqrt PI) |
| 0.5b | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 0.5b | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 26× | pow1 |
| 20× | add-exp-log |
| 12× | inv-pow |
| 12× | pow-flip |
| 10× | pow-prod-up |
| 8× | add-cube-cbrt |
| 8× | prod-exp |
| 8× | associate-*l* |
| 8× | associate-*r* |
| 8× | add-cbrt-cube |
| 8× | *-un-lft-identity |
| 8× | rec-exp |
| 8× | add-sqr-sqrt |
| 6× | pow-prod-down |
| 4× | add-log-exp |
| 4× | div-inv |
| 4× | log1p-expm1-u |
| 4× | expm1-log1p-u |
| 2× | cbrt-unprod |
| 2× | *-commutative |
| 2× | un-div-inv |
| 2× | associate-*r/ |
| 2× | associate-*l/ |
| 2× | pow1/2 |
| 2× | pow-plus |
| 2× | frac-times |
| 2× | pow2 |
| 5.0ms | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 4.0ms | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 0.0ms | (sqrt PI) |
| 0.0ms | (sqrt PI) |
| 24.0ms | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 17.0ms | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 0.0ms | (sqrt PI) |
| 0.0ms | (sqrt PI) |
| 526.0ms | (* (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x)))) |
| 171.0ms | (* (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x)))) |
| 40.0ms | (+ (- (log (fabs x))) (log (/ 1 (fabs x)))) |
| 31.0ms | (+ (log (/ 1 (fabs x))) (log (/ 1 (fabs x)))) |
| 26.0ms | (+ (- (log (fabs x))) (log (/ 1 (fabs x)))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 1.0b | (sqrt PI) |
| 1.0b | (sqrt PI) |
| 0.5b | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 0.5b | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 26× | pow1 |
| 20× | add-exp-log |
| 12× | inv-pow |
| 12× | pow-flip |
| 10× | pow-prod-up |
| 8× | add-cube-cbrt |
| 8× | prod-exp |
| 8× | associate-*l* |
| 8× | associate-*r* |
| 8× | add-cbrt-cube |
| 8× | *-un-lft-identity |
| 8× | rec-exp |
| 8× | add-sqr-sqrt |
| 6× | pow-prod-down |
| 4× | add-log-exp |
| 4× | div-inv |
| 4× | log1p-expm1-u |
| 4× | expm1-log1p-u |
| 2× | cbrt-unprod |
| 2× | *-commutative |
| 2× | un-div-inv |
| 2× | associate-*r/ |
| 2× | associate-*l/ |
| 2× | pow1/2 |
| 2× | pow-plus |
| 2× | frac-times |
| 2× | pow2 |
| 9.0ms | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 9.0ms | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 0.0ms | (sqrt PI) |
| 0.0ms | (sqrt PI) |
| 24.0ms | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 23.0ms | (* (/ 1 (fabs x)) (/ 1 (fabs x))) |
| 0.0ms | (sqrt PI) |
| 0.0ms | (sqrt PI) |
| 225.0ms | (* (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x)))) |
| 202.0ms | (fma (/ (exp (* (fabs x) (fabs x))) (sqrt PI)) (fma (* (* 3/4 (/ 1 (fabs x))) (* (/ 1 (fabs x)) (/ 1 (fabs x)))) (pow (fabs x) -2) (fma (* 1/2 (/ 1 (fabs x))) (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x)))) (* (* (/ 1 (pow (fabs x) 3)) (pow (fabs x) (* -1 4))) (* (exp (* (fabs x) (fabs x))) (* (/ 1 (sqrt PI)) 15/8)))) |
| 183.0ms | (* (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x)))) |
| 130.0ms | (fma (/ (exp (* (fabs x) (fabs x))) (sqrt PI)) (fma (* (* 3/4 (/ 1 (fabs x))) (* (/ 1 (fabs x)) (/ 1 (fabs x)))) (pow (fabs x) -2) (fma (* 1/2 (/ 1 (fabs x))) (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x)))) (* (* (/ 1 (pow (fabs x) 3)) (pow (fabs x) (* -1 4))) (* (exp (* (fabs x) (fabs x))) (* (/ 1 (sqrt PI)) 15/8)))) |
| 32.0ms | (+ (log (/ 1 (fabs x))) (log (/ 1 (fabs x)))) |
13 alts after pruning (12 fresh and 1 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 1.0b | (sqrt PI) |
| 1.0b | (sqrt PI) |
| 0.4b | (* (/ (/ 1 (fabs x)) (* (fabs x) (fabs x))) (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4))) |
| 0.4b | (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4)) |
| 83× | add-exp-log |
| 46× | prod-exp |
| 27× | add-cbrt-cube |
| 16× | div-exp |
| 12× | cbrt-unprod |
| 12× | associate-*l* |
| 9× | add-cube-cbrt |
| 9× | *-un-lft-identity |
| 9× | add-sqr-sqrt |
| 8× | rec-exp |
| 7× | associate-*r* |
| 6× | pow-exp |
| 6× | pow1 |
| 6× | pow-to-exp |
| 4× | add-log-exp |
| 4× | log1p-expm1-u |
| 4× | associate-*l/ |
| 4× | times-frac |
| 4× | cbrt-undiv |
| 4× | expm1-log1p-u |
| 3× | div-inv |
| 3× | unpow-prod-down |
| 2× | *-commutative |
| 2× | pow1/2 |
| 1× | associate-*r/ |
| 1× | pow-prod-down |
| 1× | frac-times |
| 29.0ms | (* (/ (/ 1 (fabs x)) (* (fabs x) (fabs x))) (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4))) |
| 9.0ms | (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4)) |
| 0.0ms | (sqrt PI) |
| 0.0ms | (sqrt PI) |
| 202.0ms | (* (/ (/ 1 (fabs x)) (* (fabs x) (fabs x))) (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4))) |
| 141.0ms | (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4)) |
| 0.0ms | (sqrt PI) |
| 0.0ms | (sqrt PI) |
| 597.0ms | (+ (- (- (log (fabs x))) (log (* (fabs x) (fabs x)))) (log (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4)))) |
| 515.0ms | (+ (- (- (log (fabs x))) (+ (log (fabs x)) (log (fabs x)))) (log (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4)))) |
| 322.0ms | (+ (log (/ (/ 1 (fabs x)) (* (fabs x) (fabs x)))) (log (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4)))) |
| 311.0ms | (* (/ (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (* (* (* (fabs x) (fabs x)) (fabs x)) (* (* (fabs x) (fabs x)) (fabs x)))) (* (* (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4)) (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4))) (* (/ 15/8 (sqrt PI)) (pow (fabs x) -4)))) |
| 276.0ms | (* (/ (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (* (* (* (fabs x) (fabs x)) (fabs x)) (* (* (fabs x) (fabs x)) (fabs x)))) (* (* (* (/ 15/8 (sqrt PI)) (/ 15/8 (sqrt PI))) (/ 15/8 (sqrt PI))) (* (* (pow (fabs x) -4) (pow (fabs x) -4)) (pow (fabs x) -4)))) |
10 alts after pruning (9 fresh and 1 done)
Merged error: 0b
0% (0.4b remaining)
Error of 0.5b against oracle of 0.1b and baseline of 0.5b