
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* (fabs x) (fabs x)) (fabs x)))
(t_1 (* (* t_0 (fabs x)) (fabs x))))
(fabs
(*
(/ 1.0 (sqrt (PI)))
(+
(+ (+ (* 2.0 (fabs x)) (* (/ 2.0 3.0) t_0)) (* (/ 1.0 5.0) t_1))
(* (/ 1.0 21.0) (* (* t_1 (fabs x)) (fabs x))))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\\
t_1 := \left(t\_0 \cdot \left|x\right|\right) \cdot \left|x\right|\\
\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot t\_0\right) + \frac{1}{5} \cdot t\_1\right) + \frac{1}{21} \cdot \left(\left(t\_1 \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 5 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* (fabs x) (fabs x)) (fabs x)))
(t_1 (* (* t_0 (fabs x)) (fabs x))))
(fabs
(*
(/ 1.0 (sqrt (PI)))
(+
(+ (+ (* 2.0 (fabs x)) (* (/ 2.0 3.0) t_0)) (* (/ 1.0 5.0) t_1))
(* (/ 1.0 21.0) (* (* t_1 (fabs x)) (fabs x))))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\\
t_1 := \left(t\_0 \cdot \left|x\right|\right) \cdot \left|x\right|\\
\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot t\_0\right) + \frac{1}{5} \cdot t\_1\right) + \frac{1}{21} \cdot \left(\left(t\_1 \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|
\end{array}
\end{array}
x_m = (fabs.f64 x)
(FPCore (x_m)
:precision binary64
(fabs
(*
(/ 1.0 (sqrt (PI)))
(+
(+
(* (fabs x_m) (fma (* x_m x_m) 0.6666666666666666 2.0))
(* (/ 1.0 5.0) (fabs (* (* (* (* x_m x_m) x_m) x_m) x_m))))
(* (* 0.047619047619047616 (* (* (pow x_m 4.0) x_m) x_m)) x_m)))))\begin{array}{l}
x_m = \left|x\right|
\\
\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left|x\_m\right| \cdot \mathsf{fma}\left(x\_m \cdot x\_m, 0.6666666666666666, 2\right) + \frac{1}{5} \cdot \left|\left(\left(\left(x\_m \cdot x\_m\right) \cdot x\_m\right) \cdot x\_m\right) \cdot x\_m\right|\right) + \left(0.047619047619047616 \cdot \left(\left({x\_m}^{4} \cdot x\_m\right) \cdot x\_m\right)\right) \cdot x\_m\right)\right|
\end{array}
Initial program 99.9%
lift-+.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-fabs.f64N/A
distribute-rgt-outN/A
lower-*.f64N/A
*-commutativeN/A
lower-fma.f6499.9
Applied rewrites99.9%
Applied rewrites72.8%
Final simplification72.8%
x_m = (fabs.f64 x)
(FPCore (x_m)
:precision binary64
(let* ((t_0 (* (* x_m x_m) x_m)))
(fabs
(*
(/ 1.0 (sqrt (PI)))
(+
(+
(* (fma (* x_m x_m) 0.6666666666666666 2.0) x_m)
(* 0.2 (fabs (* (* t_0 x_m) x_m))))
(* (/ 1.0 21.0) (* (* t_0 t_0) (fabs x_m))))))))\begin{array}{l}
x_m = \left|x\right|
\\
\begin{array}{l}
t_0 := \left(x\_m \cdot x\_m\right) \cdot x\_m\\
\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\mathsf{fma}\left(x\_m \cdot x\_m, 0.6666666666666666, 2\right) \cdot x\_m + 0.2 \cdot \left|\left(t\_0 \cdot x\_m\right) \cdot x\_m\right|\right) + \frac{1}{21} \cdot \left(\left(t\_0 \cdot t\_0\right) \cdot \left|x\_m\right|\right)\right)\right|
\end{array}
\end{array}
Initial program 99.9%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
swap-sqrN/A
lift-fabs.f64N/A
lift-fabs.f64N/A
sqr-absN/A
unswap-sqrN/A
Applied rewrites99.9%
lift-/.f64N/A
metadata-eval99.9
Applied rewrites99.9%
lift-+.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
sqr-abs-revN/A
pow2N/A
associate-*r*N/A
distribute-rgt-inN/A
*-commutativeN/A
pow2N/A
Applied rewrites80.4%
Final simplification80.4%
x_m = (fabs.f64 x)
(FPCore (x_m)
:precision binary64
(fabs
(*
(/ 1.0 (sqrt (PI)))
(+
(+
(* (fabs x_m) (fma (* x_m x_m) 0.6666666666666666 2.0))
(* 0.2 (fabs (* (* (* (* x_m x_m) x_m) x_m) x_m))))
(*
0.047619047619047616
(* (fabs (* (* (* (* x_m x_m) (* x_m x_m)) x_m) x_m)) (fabs x_m)))))))\begin{array}{l}
x_m = \left|x\right|
\\
\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left|x\_m\right| \cdot \mathsf{fma}\left(x\_m \cdot x\_m, 0.6666666666666666, 2\right) + 0.2 \cdot \left|\left(\left(\left(x\_m \cdot x\_m\right) \cdot x\_m\right) \cdot x\_m\right) \cdot x\_m\right|\right) + 0.047619047619047616 \cdot \left(\left|\left(\left(\left(x\_m \cdot x\_m\right) \cdot \left(x\_m \cdot x\_m\right)\right) \cdot x\_m\right) \cdot x\_m\right| \cdot \left|x\_m\right|\right)\right)\right|
\end{array}
Initial program 99.9%
lift-+.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-fabs.f64N/A
distribute-rgt-outN/A
lower-*.f64N/A
*-commutativeN/A
lower-fma.f6499.9
Applied rewrites99.9%
lift-/.f64N/A
metadata-eval99.9
Applied rewrites99.9%
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
sqr-abs-revN/A
pow2N/A
associate-*l*N/A
sqr-abs-revN/A
pow2N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6499.9
Applied rewrites99.9%
lift-/.f64N/A
metadata-eval99.9
Applied rewrites99.9%
Final simplification99.9%
x_m = (fabs.f64 x)
(FPCore (x_m)
:precision binary64
(let* ((t_0 (fabs (* (* (* (* x_m x_m) x_m) x_m) x_m))))
(fabs
(*
(/ 1.0 (sqrt (PI)))
(+
(+ (* x_m 2.0) (* 0.2 t_0))
(* 0.047619047619047616 (* (* t_0 (fabs x_m)) (fabs x_m))))))))\begin{array}{l}
x_m = \left|x\right|
\\
\begin{array}{l}
t_0 := \left|\left(\left(\left(x\_m \cdot x\_m\right) \cdot x\_m\right) \cdot x\_m\right) \cdot x\_m\right|\\
\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(x\_m \cdot 2 + 0.2 \cdot t\_0\right) + 0.047619047619047616 \cdot \left(\left(t\_0 \cdot \left|x\_m\right|\right) \cdot \left|x\_m\right|\right)\right)\right|
\end{array}
\end{array}
Initial program 99.9%
lift-+.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-fabs.f64N/A
distribute-rgt-outN/A
lower-*.f64N/A
*-commutativeN/A
lower-fma.f6499.9
Applied rewrites99.9%
lift-/.f64N/A
metadata-eval99.9
Applied rewrites99.9%
Taylor expanded in x around 0
lift-*.f64N/A
+-commutativeN/A
lift-*.f64N/A
pow2N/A
*-commutativeN/A
distribute-rgt-inN/A
associate-*r*N/A
metadata-evalN/A
pow2N/A
sqr-abs-revN/A
*-commutativeN/A
lower-*.f64N/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrt98.8
Applied rewrites98.8%
lift-/.f64N/A
metadata-eval98.8
Applied rewrites98.8%
Final simplification98.8%
x_m = (fabs.f64 x) (FPCore (x_m) :precision binary64 (fabs (* (/ 1.0 (sqrt (PI))) (* (* x_m (* 0.6666666666666666 x_m)) x_m))))
\begin{array}{l}
x_m = \left|x\right|
\\
\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(x\_m \cdot \left(0.6666666666666666 \cdot x\_m\right)\right) \cdot x\_m\right)\right|
\end{array}
Initial program 99.9%
lift-+.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-fabs.f64N/A
distribute-rgt-outN/A
lower-*.f64N/A
*-commutativeN/A
lower-fma.f6499.9
Applied rewrites99.9%
Taylor expanded in x around inf
Applied rewrites23.9%
lift-*.f64N/A
lift-pow.f64N/A
*-commutativeN/A
pow3N/A
pow2N/A
rem-square-sqrtN/A
sqrt-unprodN/A
rem-sqrt-square-revN/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrt23.9
Applied rewrites23.9%
metadata-evalN/A
lift-*.f64N/A
lift-*.f64N/A
sqr-abs-revN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
lower-*.f64N/A
metadata-evalN/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrt23.9
Applied rewrites23.9%
herbie shell --seed 2025050
(FPCore (x)
:name "Jmat.Real.erfi, branch x less than or equal to 0.5"
:precision binary64
:pre (<= x 0.5)
(fabs (* (/ 1.0 (sqrt (PI))) (+ (+ (+ (* 2.0 (fabs x)) (* (/ 2.0 3.0) (* (* (fabs x) (fabs x)) (fabs x)))) (* (/ 1.0 5.0) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))))))