\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)\frac{1}{\sqrt{\sqrt{\pi}}} \cdot \left(\left(1 + \left(\log \left(e^{\frac{0.5 + \frac{0.75}{x \cdot x}}{x \cdot x}}\right) + \frac{1.875}{{x}^{6}}\right)\right) \cdot \frac{\frac{{\left(e^{x}\right)}^{x}}{\sqrt{\sqrt{\pi}}}}{\left|x\right|}\right)(FPCore (x)
:precision binary64
(*
(* (/ 1.0 (sqrt PI)) (exp (* (fabs x) (fabs x))))
(+
(+
(+
(/ 1.0 (fabs x))
(*
(/ 1.0 2.0)
(* (* (/ 1.0 (fabs x)) (/ 1.0 (fabs x))) (/ 1.0 (fabs x)))))
(*
(/ 3.0 4.0)
(*
(*
(* (* (/ 1.0 (fabs x)) (/ 1.0 (fabs x))) (/ 1.0 (fabs x)))
(/ 1.0 (fabs x)))
(/ 1.0 (fabs x)))))
(*
(/ 15.0 8.0)
(*
(*
(*
(*
(* (* (/ 1.0 (fabs x)) (/ 1.0 (fabs x))) (/ 1.0 (fabs x)))
(/ 1.0 (fabs x)))
(/ 1.0 (fabs x)))
(/ 1.0 (fabs x)))
(/ 1.0 (fabs x)))))))(FPCore (x)
:precision binary64
(*
(/ 1.0 (sqrt (sqrt PI)))
(*
(+
1.0
(+ (log (exp (/ (+ 0.5 (/ 0.75 (* x x))) (* x x)))) (/ 1.875 (pow x 6.0))))
(/ (/ (pow (exp x) x) (sqrt (sqrt PI))) (fabs x)))))double code(double x) {
return ((1.0 / sqrt((double) M_PI)) * exp(fabs(x) * fabs(x))) * ((((1.0 / fabs(x)) + ((1.0 / 2.0) * (((1.0 / fabs(x)) * (1.0 / fabs(x))) * (1.0 / fabs(x))))) + ((3.0 / 4.0) * (((((1.0 / fabs(x)) * (1.0 / fabs(x))) * (1.0 / fabs(x))) * (1.0 / fabs(x))) * (1.0 / fabs(x))))) + ((15.0 / 8.0) * (((((((1.0 / fabs(x)) * (1.0 / fabs(x))) * (1.0 / fabs(x))) * (1.0 / fabs(x))) * (1.0 / fabs(x))) * (1.0 / fabs(x))) * (1.0 / fabs(x)))));
}
double code(double x) {
return (1.0 / sqrt(sqrt((double) M_PI))) * ((1.0 + (log(exp((0.5 + (0.75 / (x * x))) / (x * x))) + (1.875 / pow(x, 6.0)))) * ((pow(exp(x), x) / sqrt(sqrt((double) M_PI))) / fabs(x)));
}














Bits error versus x
Results
| Alternative 1 | |
|---|---|
| Error | 1.2 |
| Cost | 53184 |
| Alternative 2 | |
|---|---|
| Error | 1.2 |
| Cost | 53056 |
| Alternative 3 | |
|---|---|
| Error | 1.2 |
| Cost | 53056 |
| Alternative 4 | |
|---|---|
| Error | 1.3 |
| Cost | 33600 |
| Alternative 5 | |
|---|---|
| Error | 2.7 |
| Cost | 27264 |
| Alternative 6 | |
|---|---|
| Error | 2.7 |
| Cost | 27264 |
| Alternative 7 | |
|---|---|
| Error | 43.8 |
| Cost | 26880 |
| Alternative 8 | |
|---|---|
| Error | 47.5 |
| Cost | 32833 |
| Alternative 9 | |
|---|---|
| Error | 55.5 |
| Cost | 20992 |
| Alternative 10 | |
|---|---|
| Error | 56.2 |
| Cost | 20736 |
| Alternative 11 | |
|---|---|
| Error | 56.6 |
| Cost | 64 |


Initial program 2.8
Simplified2.7
rmApplied *-un-lft-identity_binary64_21242.7
Applied add-sqr-sqrt_binary64_21462.6
Applied *-un-lft-identity_binary64_21242.6
Applied times-frac_binary64_21302.6
Applied times-frac_binary64_21302.6
Applied associate-*l*_binary64_20652.6
Simplified1.2
rmApplied add-log-exp_binary64_21631.2
Simplified1.2
Final simplification1.2
herbie shell --seed 2021065
(FPCore (x)
:name "Jmat.Real.erfi, branch x greater than or equal to 5"
:precision binary64
:pre (>= x 0.5)
(* (* (/ 1.0 (sqrt PI)) (exp (* (fabs x) (fabs x)))) (+ (+ (+ (/ 1.0 (fabs x)) (* (/ 1.0 2.0) (* (* (/ 1.0 (fabs x)) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))))) (* (/ 3.0 4.0) (* (* (* (* (/ 1.0 (fabs x)) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))))) (* (/ 15.0 8.0) (* (* (* (* (* (* (/ 1.0 (fabs x)) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x)))))))