
(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}
Sampling outcomes in binary64 precision:
Herbie found 3 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}
x_m = (fabs.f64 x)
(FPCore (x_m)
:precision binary64
(fabs
(*
(/ 1.0 (sqrt PI))
(+
(+
(+ (* 2.0 (fabs x_m)) (* x_m (* (* x_m x_m) 0.6666666666666666)))
(* (* x_m (* x_m (* x_m (* x_m x_m)))) 0.2))
(*
(/ 1.0 21.0)
(*
(fabs x_m)
(*
(fabs x_m)
(*
(fabs x_m)
(* (fabs x_m) (* (* x_m x_m) (* (* x_m x_m) (/ 1.0 x_m))))))))))))x_m = fabs(x);
double code(double x_m) {
return fabs(((1.0 / sqrt(((double) M_PI))) * ((((2.0 * fabs(x_m)) + (x_m * ((x_m * x_m) * 0.6666666666666666))) + ((x_m * (x_m * (x_m * (x_m * x_m)))) * 0.2)) + ((1.0 / 21.0) * (fabs(x_m) * (fabs(x_m) * (fabs(x_m) * (fabs(x_m) * ((x_m * x_m) * ((x_m * x_m) * (1.0 / x_m)))))))))));
}
x_m = Math.abs(x);
public static double code(double x_m) {
return Math.abs(((1.0 / Math.sqrt(Math.PI)) * ((((2.0 * Math.abs(x_m)) + (x_m * ((x_m * x_m) * 0.6666666666666666))) + ((x_m * (x_m * (x_m * (x_m * x_m)))) * 0.2)) + ((1.0 / 21.0) * (Math.abs(x_m) * (Math.abs(x_m) * (Math.abs(x_m) * (Math.abs(x_m) * ((x_m * x_m) * ((x_m * x_m) * (1.0 / x_m)))))))))));
}
x_m = math.fabs(x) def code(x_m): return math.fabs(((1.0 / math.sqrt(math.pi)) * ((((2.0 * math.fabs(x_m)) + (x_m * ((x_m * x_m) * 0.6666666666666666))) + ((x_m * (x_m * (x_m * (x_m * x_m)))) * 0.2)) + ((1.0 / 21.0) * (math.fabs(x_m) * (math.fabs(x_m) * (math.fabs(x_m) * (math.fabs(x_m) * ((x_m * x_m) * ((x_m * x_m) * (1.0 / x_m)))))))))))
x_m = abs(x) function code(x_m) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(Float64(Float64(Float64(2.0 * abs(x_m)) + Float64(x_m * Float64(Float64(x_m * x_m) * 0.6666666666666666))) + Float64(Float64(x_m * Float64(x_m * Float64(x_m * Float64(x_m * x_m)))) * 0.2)) + Float64(Float64(1.0 / 21.0) * Float64(abs(x_m) * Float64(abs(x_m) * Float64(abs(x_m) * Float64(abs(x_m) * Float64(Float64(x_m * x_m) * Float64(Float64(x_m * x_m) * Float64(1.0 / x_m))))))))))) end
x_m = abs(x); function tmp = code(x_m) tmp = abs(((1.0 / sqrt(pi)) * ((((2.0 * abs(x_m)) + (x_m * ((x_m * x_m) * 0.6666666666666666))) + ((x_m * (x_m * (x_m * (x_m * x_m)))) * 0.2)) + ((1.0 / 21.0) * (abs(x_m) * (abs(x_m) * (abs(x_m) * (abs(x_m) * ((x_m * x_m) * ((x_m * x_m) * (1.0 / x_m))))))))))); end
x_m = N[Abs[x], $MachinePrecision] code[x$95$m_] := N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(2.0 * N[Abs[x$95$m], $MachinePrecision]), $MachinePrecision] + N[(x$95$m * N[(N[(x$95$m * x$95$m), $MachinePrecision] * 0.6666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(x$95$m * N[(x$95$m * N[(x$95$m * N[(x$95$m * x$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.2), $MachinePrecision]), $MachinePrecision] + N[(N[(1.0 / 21.0), $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(N[(x$95$m * x$95$m), $MachinePrecision] * N[(N[(x$95$m * x$95$m), $MachinePrecision] * N[(1.0 / x$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
x_m = \left|x\right|
\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\_m\right| + x\_m \cdot \left(\left(x\_m \cdot x\_m\right) \cdot 0.6666666666666666\right)\right) + \left(x\_m \cdot \left(x\_m \cdot \left(x\_m \cdot \left(x\_m \cdot x\_m\right)\right)\right)\right) \cdot 0.2\right) + \frac{1}{21} \cdot \left(\left|x\_m\right| \cdot \left(\left|x\_m\right| \cdot \left(\left|x\_m\right| \cdot \left(\left|x\_m\right| \cdot \left(\left(x\_m \cdot x\_m\right) \cdot \left(\left(x\_m \cdot x\_m\right) \cdot \frac{1}{x\_m}\right)\right)\right)\right)\right)\right)\right)\right|
\end{array}
Initial program 99.9%
lift-/.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.9
Applied egg-rr73.2%
Applied egg-rr31.8%
lift-/.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
sqr-powN/A
unpow-prod-downN/A
lift-fabs.f64N/A
lift-fabs.f64N/A
sqr-absN/A
unpow-prod-downN/A
sqr-powN/A
pow3N/A
lift-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6431.8
lift-/.f64N/A
metadata-eval31.8
Applied egg-rr31.8%
Applied egg-rr99.4%
Final simplification99.4%
x_m = (fabs.f64 x)
(FPCore (x_m)
:precision binary64
(let* ((t_0 (* x_m (* x_m x_m))))
(fabs
(*
(/ 1.0 (sqrt PI))
(+
(+
(* (* x_m (* x_m t_0)) 0.2)
(+ (* 2.0 (fabs x_m)) (* 0.6666666666666666 t_0)))
(*
(/ 1.0 21.0)
(*
(fabs x_m)
(*
(fabs x_m)
(* (fabs x_m) (* (fabs x_m) (* x_m (/ x_m (/ 1.0 x_m)))))))))))))x_m = fabs(x);
double code(double x_m) {
double t_0 = x_m * (x_m * x_m);
return fabs(((1.0 / sqrt(((double) M_PI))) * ((((x_m * (x_m * t_0)) * 0.2) + ((2.0 * fabs(x_m)) + (0.6666666666666666 * t_0))) + ((1.0 / 21.0) * (fabs(x_m) * (fabs(x_m) * (fabs(x_m) * (fabs(x_m) * (x_m * (x_m / (1.0 / x_m)))))))))));
}
x_m = Math.abs(x);
public static double code(double x_m) {
double t_0 = x_m * (x_m * x_m);
return Math.abs(((1.0 / Math.sqrt(Math.PI)) * ((((x_m * (x_m * t_0)) * 0.2) + ((2.0 * Math.abs(x_m)) + (0.6666666666666666 * t_0))) + ((1.0 / 21.0) * (Math.abs(x_m) * (Math.abs(x_m) * (Math.abs(x_m) * (Math.abs(x_m) * (x_m * (x_m / (1.0 / x_m)))))))))));
}
x_m = math.fabs(x) def code(x_m): t_0 = x_m * (x_m * x_m) return math.fabs(((1.0 / math.sqrt(math.pi)) * ((((x_m * (x_m * t_0)) * 0.2) + ((2.0 * math.fabs(x_m)) + (0.6666666666666666 * t_0))) + ((1.0 / 21.0) * (math.fabs(x_m) * (math.fabs(x_m) * (math.fabs(x_m) * (math.fabs(x_m) * (x_m * (x_m / (1.0 / x_m)))))))))))
x_m = abs(x) function code(x_m) t_0 = Float64(x_m * Float64(x_m * x_m)) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(Float64(Float64(Float64(x_m * Float64(x_m * t_0)) * 0.2) + Float64(Float64(2.0 * abs(x_m)) + Float64(0.6666666666666666 * t_0))) + Float64(Float64(1.0 / 21.0) * Float64(abs(x_m) * Float64(abs(x_m) * Float64(abs(x_m) * Float64(abs(x_m) * Float64(x_m * Float64(x_m / Float64(1.0 / x_m))))))))))) end
x_m = abs(x); function tmp = code(x_m) t_0 = x_m * (x_m * x_m); tmp = abs(((1.0 / sqrt(pi)) * ((((x_m * (x_m * t_0)) * 0.2) + ((2.0 * abs(x_m)) + (0.6666666666666666 * t_0))) + ((1.0 / 21.0) * (abs(x_m) * (abs(x_m) * (abs(x_m) * (abs(x_m) * (x_m * (x_m / (1.0 / x_m))))))))))); end
x_m = N[Abs[x], $MachinePrecision]
code[x$95$m_] := Block[{t$95$0 = N[(x$95$m * N[(x$95$m * x$95$m), $MachinePrecision]), $MachinePrecision]}, N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(x$95$m * N[(x$95$m * t$95$0), $MachinePrecision]), $MachinePrecision] * 0.2), $MachinePrecision] + N[(N[(2.0 * N[Abs[x$95$m], $MachinePrecision]), $MachinePrecision] + N[(0.6666666666666666 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(1.0 / 21.0), $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(x$95$m * N[(x$95$m / N[(1.0 / x$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
x_m = \left|x\right|
\\
\begin{array}{l}
t_0 := x\_m \cdot \left(x\_m \cdot x\_m\right)\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(x\_m \cdot \left(x\_m \cdot t\_0\right)\right) \cdot 0.2 + \left(2 \cdot \left|x\_m\right| + 0.6666666666666666 \cdot t\_0\right)\right) + \frac{1}{21} \cdot \left(\left|x\_m\right| \cdot \left(\left|x\_m\right| \cdot \left(\left|x\_m\right| \cdot \left(\left|x\_m\right| \cdot \left(x\_m \cdot \frac{x\_m}{\frac{1}{x\_m}}\right)\right)\right)\right)\right)\right)\right|
\end{array}
\end{array}
Initial program 99.9%
lift-/.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.9
Applied egg-rr73.2%
lift-/.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
sqr-powN/A
unpow-prod-downN/A
lift-fabs.f64N/A
lift-fabs.f64N/A
sqr-absN/A
unpow-prod-downN/A
sqr-powN/A
cube-unmultN/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f6473.1
lift-/.f64N/A
metadata-eval73.1
Applied egg-rr73.1%
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
sqr-powN/A
unpow-prod-downN/A
lift-fabs.f64N/A
lift-fabs.f64N/A
sqr-absN/A
unpow-prod-downN/A
sqr-powN/A
metadata-evalN/A
pow-divN/A
pow2N/A
inv-powN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6499.4
Applied egg-rr99.4%
Final simplification99.4%
x_m = (fabs.f64 x)
(FPCore (x_m)
:precision binary64
(let* ((t_0 (* x_m (* x_m x_m))))
(fabs
(*
(/ 1.0 (sqrt PI))
(+
(+
(* (* x_m (* x_m t_0)) 0.2)
(+ (* 2.0 (fabs x_m)) (* 0.6666666666666666 t_0)))
(*
(/ 1.0 21.0)
(* (fabs x_m) (* (fabs x_m) (* (fabs x_m) (* (fabs x_m) t_0))))))))))x_m = fabs(x);
double code(double x_m) {
double t_0 = x_m * (x_m * x_m);
return fabs(((1.0 / sqrt(((double) M_PI))) * ((((x_m * (x_m * t_0)) * 0.2) + ((2.0 * fabs(x_m)) + (0.6666666666666666 * t_0))) + ((1.0 / 21.0) * (fabs(x_m) * (fabs(x_m) * (fabs(x_m) * (fabs(x_m) * t_0))))))));
}
x_m = Math.abs(x);
public static double code(double x_m) {
double t_0 = x_m * (x_m * x_m);
return Math.abs(((1.0 / Math.sqrt(Math.PI)) * ((((x_m * (x_m * t_0)) * 0.2) + ((2.0 * Math.abs(x_m)) + (0.6666666666666666 * t_0))) + ((1.0 / 21.0) * (Math.abs(x_m) * (Math.abs(x_m) * (Math.abs(x_m) * (Math.abs(x_m) * t_0))))))));
}
x_m = math.fabs(x) def code(x_m): t_0 = x_m * (x_m * x_m) return math.fabs(((1.0 / math.sqrt(math.pi)) * ((((x_m * (x_m * t_0)) * 0.2) + ((2.0 * math.fabs(x_m)) + (0.6666666666666666 * t_0))) + ((1.0 / 21.0) * (math.fabs(x_m) * (math.fabs(x_m) * (math.fabs(x_m) * (math.fabs(x_m) * t_0))))))))
x_m = abs(x) function code(x_m) t_0 = Float64(x_m * Float64(x_m * x_m)) return abs(Float64(Float64(1.0 / sqrt(pi)) * Float64(Float64(Float64(Float64(x_m * Float64(x_m * t_0)) * 0.2) + Float64(Float64(2.0 * abs(x_m)) + Float64(0.6666666666666666 * t_0))) + Float64(Float64(1.0 / 21.0) * Float64(abs(x_m) * Float64(abs(x_m) * Float64(abs(x_m) * Float64(abs(x_m) * t_0)))))))) end
x_m = abs(x); function tmp = code(x_m) t_0 = x_m * (x_m * x_m); tmp = abs(((1.0 / sqrt(pi)) * ((((x_m * (x_m * t_0)) * 0.2) + ((2.0 * abs(x_m)) + (0.6666666666666666 * t_0))) + ((1.0 / 21.0) * (abs(x_m) * (abs(x_m) * (abs(x_m) * (abs(x_m) * t_0)))))))); end
x_m = N[Abs[x], $MachinePrecision]
code[x$95$m_] := Block[{t$95$0 = N[(x$95$m * N[(x$95$m * x$95$m), $MachinePrecision]), $MachinePrecision]}, N[Abs[N[(N[(1.0 / N[Sqrt[Pi], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(x$95$m * N[(x$95$m * t$95$0), $MachinePrecision]), $MachinePrecision] * 0.2), $MachinePrecision] + N[(N[(2.0 * N[Abs[x$95$m], $MachinePrecision]), $MachinePrecision] + N[(0.6666666666666666 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(1.0 / 21.0), $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * N[(N[Abs[x$95$m], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
x_m = \left|x\right|
\\
\begin{array}{l}
t_0 := x\_m \cdot \left(x\_m \cdot x\_m\right)\\
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(x\_m \cdot \left(x\_m \cdot t\_0\right)\right) \cdot 0.2 + \left(2 \cdot \left|x\_m\right| + 0.6666666666666666 \cdot t\_0\right)\right) + \frac{1}{21} \cdot \left(\left|x\_m\right| \cdot \left(\left|x\_m\right| \cdot \left(\left|x\_m\right| \cdot \left(\left|x\_m\right| \cdot t\_0\right)\right)\right)\right)\right)\right|
\end{array}
\end{array}
Initial program 99.9%
lift-/.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
lift-fabs.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.9
Applied egg-rr73.2%
lift-/.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
sqr-powN/A
unpow-prod-downN/A
lift-fabs.f64N/A
lift-fabs.f64N/A
sqr-absN/A
unpow-prod-downN/A
sqr-powN/A
cube-unmultN/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f6473.1
lift-/.f64N/A
metadata-eval73.1
Applied egg-rr73.1%
lift-fabs.f64N/A
lift-fabs.f64N/A
lift-fabs.f64N/A
pow3N/A
sqr-powN/A
unpow-prod-downN/A
lift-fabs.f64N/A
lift-fabs.f64N/A
sqr-absN/A
unpow-prod-downN/A
sqr-powN/A
pow3N/A
lift-*.f64N/A
lower-*.f6499.4
Applied egg-rr99.4%
Final simplification99.4%
herbie shell --seed 2024219
(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)))))))