Average Error: 1.5 → 0.6
Time: 1.7m
Precision: 64
Internal Precision: 128
\[\left(\frac{1}{\sqrt{\pi}} \cdot e^{\left|x\right| \cdot \left|x\right|}\right) \cdot \left(\left(\left(\frac{1}{\left|x\right|} + \frac{1}{2} \cdot \left(\left(\frac{1}{\left|x\right|} \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right)\right) + \frac{3}{4} \cdot \left(\left(\left(\left(\frac{1}{\left|x\right|} \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right)\right) + \frac{15}{8} \cdot \left(\left(\left(\left(\left(\left(\frac{1}{\left|x\right|} \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right)\right)\]
\[(\left(\frac{e^{\left|x\right| \cdot \left|x\right|}}{\sqrt{\pi}}\right) \cdot \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{4} \cdot \frac{3}{4}\right) \cdot \left(\frac{1}{\left|x\right|}\right) + \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{3}\right) \cdot \frac{1}{2} + \left(\frac{1}{\left|x\right|}\right))_*\right))_*\right) + \left(\left(e^{\left|x\right| \cdot \left|x\right|} \cdot \left(\frac{15}{8} \cdot \frac{1}{\sqrt{\pi}}\right)\right) \cdot \left({\left(\left|x\right|\right)}^{-4} \cdot {\left(\left|x\right|\right)}^{-3}\right)\right))_*\]

Error

Bits error versus x

Derivation

  1. Initial program 1.5

    \[\left(\frac{1}{\sqrt{\pi}} \cdot e^{\left|x\right| \cdot \left|x\right|}\right) \cdot \left(\left(\left(\frac{1}{\left|x\right|} + \frac{1}{2} \cdot \left(\left(\frac{1}{\left|x\right|} \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right)\right) + \frac{3}{4} \cdot \left(\left(\left(\left(\frac{1}{\left|x\right|} \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right)\right) + \frac{15}{8} \cdot \left(\left(\left(\left(\left(\left(\frac{1}{\left|x\right|} \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right) \cdot \frac{1}{\left|x\right|}\right)\right)\]
  2. Initial simplification1.4

    \[\leadsto (\left(\frac{e^{\left|x\right| \cdot \left|x\right|}}{\sqrt{\pi}}\right) \cdot \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{4} \cdot \frac{3}{4}\right) \cdot \left(\frac{1}{\left|x\right|}\right) + \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{3}\right) \cdot \frac{1}{2} + \left(\frac{1}{\left|x\right|}\right))_*\right))_*\right) + \left(\left({\left(\frac{1}{\left|x\right|}\right)}^{3} \cdot {\left(\frac{1}{\left|x\right|}\right)}^{4}\right) \cdot \left(\left(\frac{15}{8} \cdot \frac{1}{\sqrt{\pi}}\right) \cdot e^{\left|x\right| \cdot \left|x\right|}\right)\right))_*\]
  3. Using strategy rm
  4. Applied inv-pow1.4

    \[\leadsto (\left(\frac{e^{\left|x\right| \cdot \left|x\right|}}{\sqrt{\pi}}\right) \cdot \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{4} \cdot \frac{3}{4}\right) \cdot \left(\frac{1}{\left|x\right|}\right) + \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{3}\right) \cdot \frac{1}{2} + \left(\frac{1}{\left|x\right|}\right))_*\right))_*\right) + \left(\left({\left(\frac{1}{\left|x\right|}\right)}^{3} \cdot {\color{blue}{\left({\left(\left|x\right|\right)}^{-1}\right)}}^{4}\right) \cdot \left(\left(\frac{15}{8} \cdot \frac{1}{\sqrt{\pi}}\right) \cdot e^{\left|x\right| \cdot \left|x\right|}\right)\right))_*\]
  5. Applied pow-pow0.9

    \[\leadsto (\left(\frac{e^{\left|x\right| \cdot \left|x\right|}}{\sqrt{\pi}}\right) \cdot \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{4} \cdot \frac{3}{4}\right) \cdot \left(\frac{1}{\left|x\right|}\right) + \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{3}\right) \cdot \frac{1}{2} + \left(\frac{1}{\left|x\right|}\right))_*\right))_*\right) + \left(\left({\left(\frac{1}{\left|x\right|}\right)}^{3} \cdot \color{blue}{{\left(\left|x\right|\right)}^{\left(-1 \cdot 4\right)}}\right) \cdot \left(\left(\frac{15}{8} \cdot \frac{1}{\sqrt{\pi}}\right) \cdot e^{\left|x\right| \cdot \left|x\right|}\right)\right))_*\]
  6. Using strategy rm
  7. Applied inv-pow0.9

    \[\leadsto (\left(\frac{e^{\left|x\right| \cdot \left|x\right|}}{\sqrt{\pi}}\right) \cdot \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{4} \cdot \frac{3}{4}\right) \cdot \left(\frac{1}{\left|x\right|}\right) + \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{3}\right) \cdot \frac{1}{2} + \left(\frac{1}{\left|x\right|}\right))_*\right))_*\right) + \left(\left({\color{blue}{\left({\left(\left|x\right|\right)}^{-1}\right)}}^{3} \cdot {\left(\left|x\right|\right)}^{\left(-1 \cdot 4\right)}\right) \cdot \left(\left(\frac{15}{8} \cdot \frac{1}{\sqrt{\pi}}\right) \cdot e^{\left|x\right| \cdot \left|x\right|}\right)\right))_*\]
  8. Applied pow-pow0.6

    \[\leadsto (\left(\frac{e^{\left|x\right| \cdot \left|x\right|}}{\sqrt{\pi}}\right) \cdot \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{4} \cdot \frac{3}{4}\right) \cdot \left(\frac{1}{\left|x\right|}\right) + \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{3}\right) \cdot \frac{1}{2} + \left(\frac{1}{\left|x\right|}\right))_*\right))_*\right) + \left(\left(\color{blue}{{\left(\left|x\right|\right)}^{\left(-1 \cdot 3\right)}} \cdot {\left(\left|x\right|\right)}^{\left(-1 \cdot 4\right)}\right) \cdot \left(\left(\frac{15}{8} \cdot \frac{1}{\sqrt{\pi}}\right) \cdot e^{\left|x\right| \cdot \left|x\right|}\right)\right))_*\]
  9. Final simplification0.6

    \[\leadsto (\left(\frac{e^{\left|x\right| \cdot \left|x\right|}}{\sqrt{\pi}}\right) \cdot \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{4} \cdot \frac{3}{4}\right) \cdot \left(\frac{1}{\left|x\right|}\right) + \left((\left({\left(\frac{1}{\left|x\right|}\right)}^{3}\right) \cdot \frac{1}{2} + \left(\frac{1}{\left|x\right|}\right))_*\right))_*\right) + \left(\left(e^{\left|x\right| \cdot \left|x\right|} \cdot \left(\frac{15}{8} \cdot \frac{1}{\sqrt{\pi}}\right)\right) \cdot \left({\left(\left|x\right|\right)}^{-4} \cdot {\left(\left|x\right|\right)}^{-3}\right)\right))_*\]

Reproduce

herbie shell --seed 2018360 +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)))))))

Details

Time bar (total: 16.2s)Debug log

start3.4s

Algorithm
intervals

setup365.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 1.2b

localize79.0ms

Local error

Found 4 expressions with local error:

1.0b
(sqrt PI)
1.0b
(sqrt PI)
0.9b
(pow (/ 1 (fabs x)) 4)
0.9b
(pow (/ 1 (fabs x)) 4)

rewrite7.0ms

Algorithm
rewrite-expression-head
Counts
4 → 58
Calls

4 calls. Slowest were:

3.0ms
(pow (/ 1 (fabs x)) 4)
3.0ms
(pow (/ 1 (fabs x)) 4)
0.0ms
(sqrt PI)

series48.0ms

Counts
4 → 8
Calls

4 calls. Slowest were:

30.0ms
(pow (/ 1 (fabs x)) 4)
18.0ms
(pow (/ 1 (fabs x)) 4)
0.0ms
(sqrt PI)
0.0ms
(sqrt PI)

simplify1.1s

Counts
26 → 66
Calls

26 calls. Slowest were:

518.0ms
(fma (/ (exp (* (fabs x) (fabs x))) (sqrt PI)) (fma (* (pow (/ 1 (fabs x)) 4) 3/4) (/ 1 (fabs x)) (fma (pow (/ 1 (fabs x)) 3) 1/2 (/ 1 (fabs x)))) (* (* (pow (/ 1 (fabs x)) 3) (pow (/ 1 (fabs x)) 4)) (* (* 15/8 (/ 1 (sqrt PI))) (exp (* (fabs x) (fabs x))))))
475.0ms
(fma (/ (exp (* (fabs x) (fabs x))) (sqrt PI)) (fma (* (pow (/ 1 (fabs x)) 4) 3/4) (/ 1 (fabs x)) (fma (pow (/ 1 (fabs x)) 3) 1/2 (/ 1 (fabs x)))) (* (* (pow (/ 1 (fabs x)) 3) (pow (/ 1 (fabs x)) 4)) (* (* 15/8 (/ 1 (sqrt PI))) (exp (* (fabs x) (fabs x))))))
19.0ms
(* (log (/ 1 (fabs x))) 4)

prune1.9s

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 0.5b

localize27.0ms

Local error

Found 4 expressions with local error:

1.0b
(sqrt PI)
1.0b
(sqrt PI)
0.9b
(pow (/ 1 (fabs x)) 4)
0.7b
(pow (/ 1 (fabs x)) 3)

rewrite4.0ms

Algorithm
rewrite-expression-head
Counts
4 → 66
Calls

4 calls. Slowest were:

1.0ms
(pow (/ 1 (fabs x)) 4)
1.0ms
(pow (/ 1 (fabs x)) 3)
0.0ms
(sqrt PI)

series52.0ms

Counts
4 → 8
Calls

4 calls. Slowest were:

29.0ms
(pow (/ 1 (fabs x)) 3)
22.0ms
(pow (/ 1 (fabs x)) 4)
0.0ms
(sqrt PI)
0.0ms
(sqrt PI)

simplify258.0ms

Counts
32 → 74
Calls

32 calls. Slowest were:

66.0ms
(pow (* (cbrt (/ 1 (fabs x))) (cbrt (/ 1 (fabs x)))) 3)
49.0ms
(pow (* (cbrt (/ 1 (fabs x))) (cbrt (/ 1 (fabs x)))) 3)
16.0ms
(* (/ 1 (fabs x)) (/ 1 (fabs x)))

prune2.5s

Pruning

12 alts after pruning (12 fresh and 0 done)

Merged error: 0.1b

localize44.0ms

Local error

Found 4 expressions with local error:

1.0b
(sqrt PI)
1.0b
(sqrt PI)
0.9b
(pow (/ 1 (fabs x)) 4)
0.7b
(pow (/ 1 (fabs x)) 3)

rewrite7.0ms

Algorithm
rewrite-expression-head
Counts
4 → 66
Calls

4 calls. Slowest were:

3.0ms
(pow (/ 1 (fabs x)) 3)
3.0ms
(pow (/ 1 (fabs x)) 4)
0.0ms
(sqrt PI)

series54.0ms

Counts
4 → 8
Calls

4 calls. Slowest were:

33.0ms
(pow (/ 1 (fabs x)) 4)
21.0ms
(pow (/ 1 (fabs x)) 3)
0.0ms
(sqrt PI)
0.0ms
(sqrt PI)

simplify220.0ms

Counts
32 → 74
Calls

32 calls. Slowest were:

58.0ms
(pow (* (cbrt (/ 1 (fabs x))) (cbrt (/ 1 (fabs x)))) 3)
36.0ms
(pow (* (cbrt (/ 1 (fabs x))) (cbrt (/ 1 (fabs x)))) 3)
19.0ms
(* (log (/ 1 (fabs x))) 3)

prune2.6s

Pruning

14 alts after pruning (13 fresh and 1 done)

Merged error: 0.1b

localize21.0ms

Local error

Found 4 expressions with local error:

1.0b
(sqrt PI)
1.0b
(sqrt PI)
0.9b
(pow (/ 1 (fabs x)) 4)
0.7b
(pow (/ 1 (fabs x)) 3)

rewrite7.0ms

Algorithm
rewrite-expression-head
Counts
4 → 66
Calls

4 calls. Slowest were:

3.0ms
(pow (/ 1 (fabs x)) 4)
3.0ms
(pow (/ 1 (fabs x)) 3)
0.0ms
(sqrt PI)

series45.0ms

Counts
4 → 8
Calls

4 calls. Slowest were:

26.0ms
(pow (/ 1 (fabs x)) 4)
19.0ms
(pow (/ 1 (fabs x)) 3)
0.0ms
(sqrt PI)
0.0ms
(sqrt PI)

simplify235.0ms

Counts
32 → 74
Calls

32 calls. Slowest were:

38.0ms
(pow (* (cbrt (/ 1 (fabs x))) (cbrt (/ 1 (fabs x)))) 3)
37.0ms
(pow (* (cbrt (/ 1 (fabs x))) (cbrt (/ 1 (fabs x)))) 3)
19.0ms
(* (/ 1 (fabs x)) (/ 1 (fabs x)))

prune2.4s

Pruning

13 alts after pruning (12 fresh and 1 done)

Merged error: 0.1b

regimes749.0ms

Accuracy

0% (0.4b remaining)

Error of 0.6b against oracle of 0.2b and baseline of 0.6b

bsearch9.0ms