
(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))))))))
double code(double x) {
double t_0 = (fabs(x) * fabs(x)) * fabs(x);
double t_1 = (t_0 * fabs(x)) * fabs(x);
return fabs(((1.0 / sqrt(((double) M_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))))));
}
public static double code(double x) {
double t_0 = (Math.abs(x) * Math.abs(x)) * Math.abs(x);
double t_1 = (t_0 * Math.abs(x)) * Math.abs(x);
return Math.abs(((1.0 / Math.sqrt(Math.PI)) * ((((2.0 * Math.abs(x)) + ((2.0 / 3.0) * t_0)) + ((1.0 / 5.0) * t_1)) + ((1.0 / 21.0) * ((t_1 * Math.abs(x)) * Math.abs(x))))));
}
def code(x): t_0 = (math.fabs(x) * math.fabs(x)) * math.fabs(x) t_1 = (t_0 * math.fabs(x)) * math.fabs(x) return math.fabs(((1.0 / math.sqrt(math.pi)) * ((((2.0 * math.fabs(x)) + ((2.0 / 3.0) * t_0)) + ((1.0 / 5.0) * t_1)) + ((1.0 / 21.0) * ((t_1 * math.fabs(x)) * math.fabs(x))))))
function code(x) t_0 = Float64(Float64(abs(x) * abs(x)) * abs(x)) t_1 = Float64(Float64(t_0 * abs(x)) * abs(x)) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(Float64(Float64(Float64(2.0 * abs(x)) + Float64(Float64(2.0 / 3.0) * t_0)) + Float64(Float64(1.0 / 5.0) * t_1)) + Float64(Float64(1.0 / 21.0) * Float64(Float64(t_1 * abs(x)) * abs(x)))))) end
function tmp = code(x) t_0 = (abs(x) * abs(x)) * abs(x); t_1 = (t_0 * abs(x)) * abs(x); tmp = abs(((1.0 / sqrt(pi)) * ((((2.0 * abs(x)) + ((2.0 / 3.0) * t_0)) + ((1.0 / 5.0) * t_1)) + ((1.0 / 21.0) * ((t_1 * abs(x)) * abs(x)))))); end
code[x_] := Block[{t$95$0 = N[(N[(N[Abs[x], $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(t$95$0 * N[Abs[x], $MachinePrecision]), $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision]}, N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(2.0 * N[Abs[x], $MachinePrecision]), $MachinePrecision] + N[(N[(2.0 / 3.0), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision] + N[(N[(1.0 / 5.0), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] + N[(N[(1.0 / 21.0), $MachinePrecision] * N[(N[(t$95$1 * N[Abs[x], $MachinePrecision]), $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\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{\pi}} \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}
Herbie found 12 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))))))))
double code(double x) {
double t_0 = (fabs(x) * fabs(x)) * fabs(x);
double t_1 = (t_0 * fabs(x)) * fabs(x);
return fabs(((1.0 / sqrt(((double) M_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))))));
}
public static double code(double x) {
double t_0 = (Math.abs(x) * Math.abs(x)) * Math.abs(x);
double t_1 = (t_0 * Math.abs(x)) * Math.abs(x);
return Math.abs(((1.0 / Math.sqrt(Math.PI)) * ((((2.0 * Math.abs(x)) + ((2.0 / 3.0) * t_0)) + ((1.0 / 5.0) * t_1)) + ((1.0 / 21.0) * ((t_1 * Math.abs(x)) * Math.abs(x))))));
}
def code(x): t_0 = (math.fabs(x) * math.fabs(x)) * math.fabs(x) t_1 = (t_0 * math.fabs(x)) * math.fabs(x) return math.fabs(((1.0 / math.sqrt(math.pi)) * ((((2.0 * math.fabs(x)) + ((2.0 / 3.0) * t_0)) + ((1.0 / 5.0) * t_1)) + ((1.0 / 21.0) * ((t_1 * math.fabs(x)) * math.fabs(x))))))
function code(x) t_0 = Float64(Float64(abs(x) * abs(x)) * abs(x)) t_1 = Float64(Float64(t_0 * abs(x)) * abs(x)) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(Float64(Float64(Float64(2.0 * abs(x)) + Float64(Float64(2.0 / 3.0) * t_0)) + Float64(Float64(1.0 / 5.0) * t_1)) + Float64(Float64(1.0 / 21.0) * Float64(Float64(t_1 * abs(x)) * abs(x)))))) end
function tmp = code(x) t_0 = (abs(x) * abs(x)) * abs(x); t_1 = (t_0 * abs(x)) * abs(x); tmp = abs(((1.0 / sqrt(pi)) * ((((2.0 * abs(x)) + ((2.0 / 3.0) * t_0)) + ((1.0 / 5.0) * t_1)) + ((1.0 / 21.0) * ((t_1 * abs(x)) * abs(x)))))); end
code[x_] := Block[{t$95$0 = N[(N[(N[Abs[x], $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(t$95$0 * N[Abs[x], $MachinePrecision]), $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision]}, N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(2.0 * N[Abs[x], $MachinePrecision]), $MachinePrecision] + N[(N[(2.0 / 3.0), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision] + N[(N[(1.0 / 5.0), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] + N[(N[(1.0 / 21.0), $MachinePrecision] * N[(N[(t$95$1 * N[Abs[x], $MachinePrecision]), $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\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{\pi}} \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}
(FPCore (x)
:precision binary64
(fabs
(*
(*
(+
(/
(+
(fma
(pow x 7.0)
0.047619047619047616
(fma (pow x 5.0) 0.2 (* 0.6666666666666666 (pow x 3.0))))
x)
x)
1.0)
(pow (sqrt PI) -1.0))
x)))
double code(double x) {
return fabs((((((fma(pow(x, 7.0), 0.047619047619047616, fma(pow(x, 5.0), 0.2, (0.6666666666666666 * pow(x, 3.0)))) + x) / x) + 1.0) * pow(sqrt(((double) M_PI)), -1.0)) * x));
}
function code(x) return abs(Float64(Float64(Float64(Float64(Float64(fma((x ^ 7.0), 0.047619047619047616, fma((x ^ 5.0), 0.2, Float64(0.6666666666666666 * (x ^ 3.0)))) + x) / x) + 1.0) * (sqrt(pi) ^ -1.0)) * x)) end
code[x_] := N[Abs[N[(N[(N[(N[(N[(N[(N[Power[x, 7.0], $MachinePrecision] * 0.047619047619047616 + N[(N[Power[x, 5.0], $MachinePrecision] * 0.2 + N[(0.6666666666666666 * N[Power[x, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] / x), $MachinePrecision] + 1.0), $MachinePrecision] * N[Power[N[Sqrt[Pi], $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\left|\left(\left(\frac{\mathsf{fma}\left({x}^{7}, 0.047619047619047616, \mathsf{fma}\left({x}^{5}, 0.2, 0.6666666666666666 \cdot {x}^{3}\right)\right) + x}{x} + 1\right) \cdot {\left(\sqrt{\pi}\right)}^{-1}\right) \cdot x\right|
\end{array}
Initial program 99.9%
lift-*.f64N/A
count-2-revN/A
lift-fabs.f64N/A
rem-sqrt-square-revN/A
sqrt-prodN/A
lower-fma.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f6435.1
Applied rewrites35.1%
Taylor expanded in x around inf
Applied rewrites99.9%
(FPCore (x)
:precision binary64
(fabs
(*
(/ 1.0 (sqrt PI))
(+
(+
(+ (* 2.0 (fabs x)) (* (/ 2.0 3.0) (* (* x x) (fabs x))))
(* (/ 1.0 5.0) (fabs (* (* (* (* x x) x) x) x))))
(* (pow (fabs x) 7.0) 0.047619047619047616)))))
double code(double x) {
return fabs(((1.0 / sqrt(((double) M_PI))) * ((((2.0 * fabs(x)) + ((2.0 / 3.0) * ((x * x) * fabs(x)))) + ((1.0 / 5.0) * fabs(((((x * x) * x) * x) * x)))) + (pow(fabs(x), 7.0) * 0.047619047619047616))));
}
public static double code(double x) {
return Math.abs(((1.0 / Math.sqrt(Math.PI)) * ((((2.0 * Math.abs(x)) + ((2.0 / 3.0) * ((x * x) * Math.abs(x)))) + ((1.0 / 5.0) * Math.abs(((((x * x) * x) * x) * x)))) + (Math.pow(Math.abs(x), 7.0) * 0.047619047619047616))));
}
def code(x): return math.fabs(((1.0 / math.sqrt(math.pi)) * ((((2.0 * math.fabs(x)) + ((2.0 / 3.0) * ((x * x) * math.fabs(x)))) + ((1.0 / 5.0) * math.fabs(((((x * x) * x) * x) * x)))) + (math.pow(math.fabs(x), 7.0) * 0.047619047619047616))))
function code(x) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(Float64(Float64(Float64(2.0 * abs(x)) + Float64(Float64(2.0 / 3.0) * Float64(Float64(x * x) * abs(x)))) + Float64(Float64(1.0 / 5.0) * abs(Float64(Float64(Float64(Float64(x * x) * x) * x) * x)))) + Float64((abs(x) ^ 7.0) * 0.047619047619047616)))) end
function tmp = code(x) tmp = abs(((1.0 / sqrt(pi)) * ((((2.0 * abs(x)) + ((2.0 / 3.0) * ((x * x) * abs(x)))) + ((1.0 / 5.0) * abs(((((x * x) * x) * x) * x)))) + ((abs(x) ^ 7.0) * 0.047619047619047616)))); end
code[x_] := N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(2.0 * N[Abs[x], $MachinePrecision]), $MachinePrecision] + N[(N[(2.0 / 3.0), $MachinePrecision] * N[(N[(x * x), $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(1.0 / 5.0), $MachinePrecision] * N[Abs[N[(N[(N[(N[(x * x), $MachinePrecision] * x), $MachinePrecision] * x), $MachinePrecision] * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[Power[N[Abs[x], $MachinePrecision], 7.0], $MachinePrecision] * 0.047619047619047616), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(x \cdot x\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left|\left(\left(\left(x \cdot x\right) \cdot x\right) \cdot x\right) \cdot x\right|\right) + {\left(\left|x\right|\right)}^{7} \cdot 0.047619047619047616\right)\right|
\end{array}
Initial program 99.9%
Taylor expanded in x around 0
metadata-evalN/A
lift-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-fabs.f64N/A
lower-pow.f6499.9
lift-/.f64N/A
metadata-eval99.9
Applied rewrites99.9%
Final simplification99.9%
(FPCore (x)
:precision binary64
(fabs
(*
(/ 1.0 (sqrt PI))
(+
(+
(* (fma (* x x) 0.6666666666666666 2.0) x)
(* (/ 1.0 5.0) (fabs (* (* (* (* x x) x) x) x))))
(* (pow (fabs x) 7.0) 0.047619047619047616)))))
double code(double x) {
return fabs(((1.0 / sqrt(((double) M_PI))) * (((fma((x * x), 0.6666666666666666, 2.0) * x) + ((1.0 / 5.0) * fabs(((((x * x) * x) * x) * x)))) + (pow(fabs(x), 7.0) * 0.047619047619047616))));
}
function code(x) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(Float64(Float64(fma(Float64(x * x), 0.6666666666666666, 2.0) * x) + Float64(Float64(1.0 / 5.0) * abs(Float64(Float64(Float64(Float64(x * x) * x) * x) * x)))) + Float64((abs(x) ^ 7.0) * 0.047619047619047616)))) end
code[x_] := N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(N[(x * x), $MachinePrecision] * 0.6666666666666666 + 2.0), $MachinePrecision] * x), $MachinePrecision] + N[(N[(1.0 / 5.0), $MachinePrecision] * N[Abs[N[(N[(N[(N[(x * x), $MachinePrecision] * x), $MachinePrecision] * x), $MachinePrecision] * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[Power[N[Abs[x], $MachinePrecision], 7.0], $MachinePrecision] * 0.047619047619047616), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x + \frac{1}{5} \cdot \left|\left(\left(\left(x \cdot x\right) \cdot x\right) \cdot x\right) \cdot x\right|\right) + {\left(\left|x\right|\right)}^{7} \cdot 0.047619047619047616\right)\right|
\end{array}
Initial program 99.9%
Taylor expanded in x around 0
metadata-evalN/A
lift-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-fabs.f64N/A
lower-pow.f6499.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
metadata-evalN/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
Applied rewrites76.6%
Final simplification76.6%
(FPCore (x)
:precision binary64
(fabs
(*
(*
(fma
(fma (fma x (* 0.047619047619047616 x) 0.2) (* x x) 0.6666666666666666)
(* x x)
2.0)
(pow (sqrt PI) -1.0))
x)))
double code(double x) {
return fabs(((fma(fma(fma(x, (0.047619047619047616 * x), 0.2), (x * x), 0.6666666666666666), (x * x), 2.0) * pow(sqrt(((double) M_PI)), -1.0)) * x));
}
function code(x) return abs(Float64(Float64(fma(fma(fma(x, Float64(0.047619047619047616 * x), 0.2), Float64(x * x), 0.6666666666666666), Float64(x * x), 2.0) * (sqrt(pi) ^ -1.0)) * x)) end
code[x_] := N[Abs[N[(N[(N[(N[(N[(x * N[(0.047619047619047616 * x), $MachinePrecision] + 0.2), $MachinePrecision] * N[(x * x), $MachinePrecision] + 0.6666666666666666), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * N[Power[N[Sqrt[Pi], $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\left|\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(x, 0.047619047619047616 \cdot x, 0.2\right), x \cdot x, 0.6666666666666666\right), x \cdot x, 2\right) \cdot {\left(\sqrt{\pi}\right)}^{-1}\right) \cdot x\right|
\end{array}
Initial program 99.9%
lift-*.f64N/A
count-2-revN/A
lift-fabs.f64N/A
rem-sqrt-square-revN/A
sqrt-prodN/A
lower-fma.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f6435.1
Applied rewrites35.1%
Taylor expanded in x around inf
Applied rewrites99.9%
Taylor expanded in x around 0
metadata-evalN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
metadata-evalN/A
metadata-evalN/A
lift-fma.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
rem-square-sqrtN/A
sqrt-unprodN/A
rem-sqrt-square-revN/A
lower-fma.f64N/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
lower-*.f64N/A
metadata-evalN/A
metadata-eval99.9
Applied rewrites99.9%
(FPCore (x)
:precision binary64
(fabs
(*
(*
(pow PI -0.5)
(fma
(fma (fma (* x x) 0.047619047619047616 0.2) (* x x) 0.6666666666666666)
(* x x)
2.0))
x)))
double code(double x) {
return fabs(((pow(((double) M_PI), -0.5) * fma(fma(fma((x * x), 0.047619047619047616, 0.2), (x * x), 0.6666666666666666), (x * x), 2.0)) * x));
}
function code(x) return abs(Float64(Float64((pi ^ -0.5) * fma(fma(fma(Float64(x * x), 0.047619047619047616, 0.2), Float64(x * x), 0.6666666666666666), Float64(x * x), 2.0)) * x)) end
code[x_] := N[Abs[N[(N[(N[Power[Pi, -0.5], $MachinePrecision] * N[(N[(N[(N[(x * x), $MachinePrecision] * 0.047619047619047616 + 0.2), $MachinePrecision] * N[(x * x), $MachinePrecision] + 0.6666666666666666), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\left|\left({\pi}^{-0.5} \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(x \cdot x, 0.047619047619047616, 0.2\right), x \cdot x, 0.6666666666666666\right), x \cdot x, 2\right)\right) \cdot x\right|
\end{array}
Initial program 99.9%
lift-*.f64N/A
count-2-revN/A
lift-fabs.f64N/A
rem-sqrt-square-revN/A
sqrt-prodN/A
lower-fma.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f6435.1
Applied rewrites35.1%
Taylor expanded in x around inf
Applied rewrites99.9%
Taylor expanded in x around 0
metadata-evalN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Applied rewrites99.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) x)))
(fabs
(*
(/ 1.0 (sqrt PI))
(+
(* (fma (fma (* 0.2 x) x 0.6666666666666666) (* x x) 2.0) x)
(* (/ 1.0 21.0) (* (* t_0 t_0) (fabs x))))))))
double code(double x) {
double t_0 = (x * x) * x;
return fabs(((1.0 / sqrt(((double) M_PI))) * ((fma(fma((0.2 * x), x, 0.6666666666666666), (x * x), 2.0) * x) + ((1.0 / 21.0) * ((t_0 * t_0) * fabs(x))))));
}
function code(x) t_0 = Float64(Float64(x * x) * x) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(Float64(fma(fma(Float64(0.2 * x), x, 0.6666666666666666), Float64(x * x), 2.0) * x) + Float64(Float64(1.0 / 21.0) * Float64(Float64(t_0 * t_0) * abs(x)))))) end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * x), $MachinePrecision]}, N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(N[(0.2 * x), $MachinePrecision] * x + 0.6666666666666666), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * x), $MachinePrecision] + N[(N[(1.0 / 21.0), $MachinePrecision] * N[(N[(t$95$0 * t$95$0), $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot x\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(0.2 \cdot x, x, 0.6666666666666666\right), x \cdot x, 2\right) \cdot x + \frac{1}{21} \cdot \left(\left(t\_0 \cdot t\_0\right) \cdot \left|x\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
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-/.f64N/A
metadata-evalN/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
Applied rewrites76.6%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
pow3N/A
pow2N/A
associate-*l*N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6476.6
Applied rewrites76.6%
Taylor expanded in x around 0
metadata-evalN/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites72.2%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) x)))
(fabs
(*
(/ 1.0 (sqrt PI))
(+
(* (fma (* x x) 0.6666666666666666 2.0) x)
(* (/ 1.0 21.0) (* (* t_0 t_0) (fabs x))))))))
double code(double x) {
double t_0 = (x * x) * x;
return fabs(((1.0 / sqrt(((double) M_PI))) * ((fma((x * x), 0.6666666666666666, 2.0) * x) + ((1.0 / 21.0) * ((t_0 * t_0) * fabs(x))))));
}
function code(x) t_0 = Float64(Float64(x * x) * x) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(Float64(fma(Float64(x * x), 0.6666666666666666, 2.0) * x) + Float64(Float64(1.0 / 21.0) * Float64(Float64(t_0 * t_0) * abs(x)))))) end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * x), $MachinePrecision]}, N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(x * x), $MachinePrecision] * 0.6666666666666666 + 2.0), $MachinePrecision] * x), $MachinePrecision] + N[(N[(1.0 / 21.0), $MachinePrecision] * N[(N[(t$95$0 * t$95$0), $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot x\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x + \frac{1}{21} \cdot \left(\left(t\_0 \cdot t\_0\right) \cdot \left|x\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
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-/.f64N/A
metadata-evalN/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
Applied rewrites76.6%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
pow3N/A
pow2N/A
associate-*l*N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6476.6
Applied rewrites76.6%
Taylor expanded in x around 0
*-commutativeN/A
+-commutativeN/A
*-commutativeN/A
pow2N/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f6476.2
Applied rewrites76.2%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (* x x) x)))
(fabs
(*
(/ 1.0 (sqrt PI))
(+ (* 2.0 x) (* (/ 1.0 21.0) (* (* t_0 t_0) (fabs x))))))))
double code(double x) {
double t_0 = (x * x) * x;
return fabs(((1.0 / sqrt(((double) M_PI))) * ((2.0 * x) + ((1.0 / 21.0) * ((t_0 * t_0) * fabs(x))))));
}
public static double code(double x) {
double t_0 = (x * x) * x;
return Math.abs(((1.0 / Math.sqrt(Math.PI)) * ((2.0 * x) + ((1.0 / 21.0) * ((t_0 * t_0) * Math.abs(x))))));
}
def code(x): t_0 = (x * x) * x return math.fabs(((1.0 / math.sqrt(math.pi)) * ((2.0 * x) + ((1.0 / 21.0) * ((t_0 * t_0) * math.fabs(x))))))
function code(x) t_0 = Float64(Float64(x * x) * x) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(Float64(2.0 * x) + Float64(Float64(1.0 / 21.0) * Float64(Float64(t_0 * t_0) * abs(x)))))) end
function tmp = code(x) t_0 = (x * x) * x; tmp = abs(((1.0 / sqrt(pi)) * ((2.0 * x) + ((1.0 / 21.0) * ((t_0 * t_0) * abs(x)))))); end
code[x_] := Block[{t$95$0 = N[(N[(x * x), $MachinePrecision] * x), $MachinePrecision]}, N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * x), $MachinePrecision] + N[(N[(1.0 / 21.0), $MachinePrecision] * N[(N[(t$95$0 * t$95$0), $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot x\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(2 \cdot x + \frac{1}{21} \cdot \left(\left(t\_0 \cdot t\_0\right) \cdot \left|x\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
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-/.f64N/A
metadata-evalN/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
Applied rewrites76.6%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
pow3N/A
pow2N/A
associate-*l*N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6476.6
Applied rewrites76.6%
Taylor expanded in x around 0
lower-*.f6499.2
Applied rewrites99.2%
(FPCore (x) :precision binary64 (fabs (* (/ 1.0 (sqrt PI)) (* (fma (fma (* 0.2 x) x 0.6666666666666666) (* x x) 2.0) x))))
double code(double x) {
return fabs(((1.0 / sqrt(((double) M_PI))) * (fma(fma((0.2 * x), x, 0.6666666666666666), (x * x), 2.0) * x)));
}
function code(x) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(fma(fma(Float64(0.2 * x), x, 0.6666666666666666), Float64(x * x), 2.0) * x))) end
code[x_] := N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(0.2 * x), $MachinePrecision] * x + 0.6666666666666666), $MachinePrecision] * N[(x * x), $MachinePrecision] + 2.0), $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(0.2 \cdot x, x, 0.6666666666666666\right), x \cdot x, 2\right) \cdot x\right)\right|
\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
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-/.f64N/A
metadata-evalN/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
Applied rewrites76.6%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
pow3N/A
pow2N/A
associate-*l*N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6476.6
Applied rewrites76.6%
Taylor expanded in x around 0
metadata-evalN/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites93.6%
(FPCore (x) :precision binary64 (fabs (* (/ 1.0 (sqrt PI)) (* (fma (* x x) 0.6666666666666666 2.0) x))))
double code(double x) {
return fabs(((1.0 / sqrt(((double) M_PI))) * (fma((x * x), 0.6666666666666666, 2.0) * x)));
}
function code(x) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(fma(Float64(x * x), 0.6666666666666666, 2.0) * x))) end
code[x_] := N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(x * x), $MachinePrecision] * 0.6666666666666666 + 2.0), $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x\right)\right|
\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
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-/.f64N/A
metadata-evalN/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
Applied rewrites76.6%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
pow3N/A
pow2N/A
associate-*l*N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6476.6
Applied rewrites76.6%
Taylor expanded in x around 0
*-commutativeN/A
+-commutativeN/A
*-commutativeN/A
pow2N/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f6489.4
Applied rewrites89.4%
(FPCore (x) :precision binary64 (fabs (* (/ 1.0 (sqrt PI)) (* 2.0 x))))
double code(double x) {
return fabs(((1.0 / sqrt(((double) M_PI))) * (2.0 * x)));
}
public static double code(double x) {
return Math.abs(((1.0 / Math.sqrt(Math.PI)) * (2.0 * x)));
}
def code(x): return math.fabs(((1.0 / math.sqrt(math.pi)) * (2.0 * x)))
function code(x) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(2.0 * x))) end
function tmp = code(x) tmp = abs(((1.0 / sqrt(pi)) * (2.0 * x))); end
code[x_] := N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(2.0 * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(2 \cdot x\right)\right|
\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
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-/.f64N/A
metadata-evalN/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
Applied rewrites76.6%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
rem-sqrt-square-revN/A
sqrt-unprodN/A
rem-square-sqrtN/A
pow3N/A
pow2N/A
associate-*l*N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6476.6
Applied rewrites76.6%
Taylor expanded in x around 0
lower-*.f6467.5
Applied rewrites67.5%
(FPCore (x) :precision binary64 (fabs (/ x (sqrt PI))))
double code(double x) {
return fabs((x / sqrt(((double) M_PI))));
}
public static double code(double x) {
return Math.abs((x / Math.sqrt(Math.PI)));
}
def code(x): return math.fabs((x / math.sqrt(math.pi)))
function code(x) return abs(Float64(x / sqrt(pi))) end
function tmp = code(x) tmp = abs((x / sqrt(pi))); end
code[x_] := N[Abs[N[(x / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\left|\frac{x}{\sqrt{\pi}}\right|
\end{array}
Initial program 99.9%
lift-*.f64N/A
count-2-revN/A
lift-fabs.f64N/A
rem-sqrt-square-revN/A
sqrt-prodN/A
lower-fma.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f6435.1
Applied rewrites35.1%
Taylor expanded in x around inf
Applied rewrites14.2%
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
lower-/.f64N/A
Applied rewrites14.2%
Applied rewrites14.2%
herbie shell --seed 2025084
(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)))))))