
(FPCore (x)
:precision binary64
(let* ((t_0 (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x))))))
(-
1.0
(*
(*
t_0
(+
0.254829592
(*
t_0
(+
-0.284496736
(*
t_0
(+ 1.421413741 (* t_0 (+ -1.453152027 (* t_0 1.061405429)))))))))
(exp (- (* (fabs x) (fabs x))))))))
double code(double x) {
double t_0 = 1.0 / (1.0 + (0.3275911 * fabs(x)));
return 1.0 - ((t_0 * (0.254829592 + (t_0 * (-0.284496736 + (t_0 * (1.421413741 + (t_0 * (-1.453152027 + (t_0 * 1.061405429))))))))) * exp(-(fabs(x) * fabs(x))));
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
t_0 = 1.0d0 / (1.0d0 + (0.3275911d0 * abs(x)))
code = 1.0d0 - ((t_0 * (0.254829592d0 + (t_0 * ((-0.284496736d0) + (t_0 * (1.421413741d0 + (t_0 * ((-1.453152027d0) + (t_0 * 1.061405429d0))))))))) * exp(-(abs(x) * abs(x))))
end function
public static double code(double x) {
double t_0 = 1.0 / (1.0 + (0.3275911 * Math.abs(x)));
return 1.0 - ((t_0 * (0.254829592 + (t_0 * (-0.284496736 + (t_0 * (1.421413741 + (t_0 * (-1.453152027 + (t_0 * 1.061405429))))))))) * Math.exp(-(Math.abs(x) * Math.abs(x))));
}
def code(x): t_0 = 1.0 / (1.0 + (0.3275911 * math.fabs(x))) return 1.0 - ((t_0 * (0.254829592 + (t_0 * (-0.284496736 + (t_0 * (1.421413741 + (t_0 * (-1.453152027 + (t_0 * 1.061405429))))))))) * math.exp(-(math.fabs(x) * math.fabs(x))))
function code(x) t_0 = Float64(1.0 / Float64(1.0 + Float64(0.3275911 * abs(x)))) return Float64(1.0 - Float64(Float64(t_0 * Float64(0.254829592 + Float64(t_0 * Float64(-0.284496736 + Float64(t_0 * Float64(1.421413741 + Float64(t_0 * Float64(-1.453152027 + Float64(t_0 * 1.061405429))))))))) * exp(Float64(-Float64(abs(x) * abs(x)))))) end
function tmp = code(x) t_0 = 1.0 / (1.0 + (0.3275911 * abs(x))); tmp = 1.0 - ((t_0 * (0.254829592 + (t_0 * (-0.284496736 + (t_0 * (1.421413741 + (t_0 * (-1.453152027 + (t_0 * 1.061405429))))))))) * exp(-(abs(x) * abs(x)))); end
code[x_] := Block[{t$95$0 = N[(1.0 / N[(1.0 + N[(0.3275911 * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(1.0 - N[(N[(t$95$0 * N[(0.254829592 + N[(t$95$0 * N[(-0.284496736 + N[(t$95$0 * N[(1.421413741 + N[(t$95$0 * N[(-1.453152027 + N[(t$95$0 * 1.061405429), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Exp[(-N[(N[Abs[x], $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision])], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{1}{1 + 0.3275911 \cdot \left|x\right|}\\
1 - \left(t\_0 \cdot \left(0.254829592 + t\_0 \cdot \left(-0.284496736 + t\_0 \cdot \left(1.421413741 + t\_0 \cdot \left(-1.453152027 + t\_0 \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 10 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x)
:precision binary64
(let* ((t_0 (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x))))))
(-
1.0
(*
(*
t_0
(+
0.254829592
(*
t_0
(+
-0.284496736
(*
t_0
(+ 1.421413741 (* t_0 (+ -1.453152027 (* t_0 1.061405429)))))))))
(exp (- (* (fabs x) (fabs x))))))))
double code(double x) {
double t_0 = 1.0 / (1.0 + (0.3275911 * fabs(x)));
return 1.0 - ((t_0 * (0.254829592 + (t_0 * (-0.284496736 + (t_0 * (1.421413741 + (t_0 * (-1.453152027 + (t_0 * 1.061405429))))))))) * exp(-(fabs(x) * fabs(x))));
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
t_0 = 1.0d0 / (1.0d0 + (0.3275911d0 * abs(x)))
code = 1.0d0 - ((t_0 * (0.254829592d0 + (t_0 * ((-0.284496736d0) + (t_0 * (1.421413741d0 + (t_0 * ((-1.453152027d0) + (t_0 * 1.061405429d0))))))))) * exp(-(abs(x) * abs(x))))
end function
public static double code(double x) {
double t_0 = 1.0 / (1.0 + (0.3275911 * Math.abs(x)));
return 1.0 - ((t_0 * (0.254829592 + (t_0 * (-0.284496736 + (t_0 * (1.421413741 + (t_0 * (-1.453152027 + (t_0 * 1.061405429))))))))) * Math.exp(-(Math.abs(x) * Math.abs(x))));
}
def code(x): t_0 = 1.0 / (1.0 + (0.3275911 * math.fabs(x))) return 1.0 - ((t_0 * (0.254829592 + (t_0 * (-0.284496736 + (t_0 * (1.421413741 + (t_0 * (-1.453152027 + (t_0 * 1.061405429))))))))) * math.exp(-(math.fabs(x) * math.fabs(x))))
function code(x) t_0 = Float64(1.0 / Float64(1.0 + Float64(0.3275911 * abs(x)))) return Float64(1.0 - Float64(Float64(t_0 * Float64(0.254829592 + Float64(t_0 * Float64(-0.284496736 + Float64(t_0 * Float64(1.421413741 + Float64(t_0 * Float64(-1.453152027 + Float64(t_0 * 1.061405429))))))))) * exp(Float64(-Float64(abs(x) * abs(x)))))) end
function tmp = code(x) t_0 = 1.0 / (1.0 + (0.3275911 * abs(x))); tmp = 1.0 - ((t_0 * (0.254829592 + (t_0 * (-0.284496736 + (t_0 * (1.421413741 + (t_0 * (-1.453152027 + (t_0 * 1.061405429))))))))) * exp(-(abs(x) * abs(x)))); end
code[x_] := Block[{t$95$0 = N[(1.0 / N[(1.0 + N[(0.3275911 * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(1.0 - N[(N[(t$95$0 * N[(0.254829592 + N[(t$95$0 * N[(-0.284496736 + N[(t$95$0 * N[(1.421413741 + N[(t$95$0 * N[(-1.453152027 + N[(t$95$0 * 1.061405429), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Exp[(-N[(N[Abs[x], $MachinePrecision] * N[Abs[x], $MachinePrecision]), $MachinePrecision])], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{1}{1 + 0.3275911 \cdot \left|x\right|}\\
1 - \left(t\_0 \cdot \left(0.254829592 + t\_0 \cdot \left(-0.284496736 + t\_0 \cdot \left(1.421413741 + t\_0 \cdot \left(-1.453152027 + t\_0 \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|}
\end{array}
\end{array}
(FPCore (x)
:precision binary64
(let* ((t_0 (fma 0.3275911 (fabs x) 1.0))
(t_1
(+
0.254829592
(/
(+
-0.284496736
(/
(+ 1.421413741 (/ (+ -1.453152027 (/ 1.061405429 t_0)) t_0))
t_0))
t_0)))
(t_2 (/ t_1 (* t_0 (exp (* x x)))))
(t_3 (pow t_2 2.0)))
(/
(+ (/ 1.0 (+ 1.0 t_3)) (/ (pow t_2 4.0) (- -1.0 t_3)))
(fma t_1 (/ (exp (- (* x x))) t_0) 1.0))))
double code(double x) {
double t_0 = fma(0.3275911, fabs(x), 1.0);
double t_1 = 0.254829592 + ((-0.284496736 + ((1.421413741 + ((-1.453152027 + (1.061405429 / t_0)) / t_0)) / t_0)) / t_0);
double t_2 = t_1 / (t_0 * exp((x * x)));
double t_3 = pow(t_2, 2.0);
return ((1.0 / (1.0 + t_3)) + (pow(t_2, 4.0) / (-1.0 - t_3))) / fma(t_1, (exp(-(x * x)) / t_0), 1.0);
}
function code(x) t_0 = fma(0.3275911, abs(x), 1.0) t_1 = Float64(0.254829592 + Float64(Float64(-0.284496736 + Float64(Float64(1.421413741 + Float64(Float64(-1.453152027 + Float64(1.061405429 / t_0)) / t_0)) / t_0)) / t_0)) t_2 = Float64(t_1 / Float64(t_0 * exp(Float64(x * x)))) t_3 = t_2 ^ 2.0 return Float64(Float64(Float64(1.0 / Float64(1.0 + t_3)) + Float64((t_2 ^ 4.0) / Float64(-1.0 - t_3))) / fma(t_1, Float64(exp(Float64(-Float64(x * x))) / t_0), 1.0)) end
code[x_] := Block[{t$95$0 = N[(0.3275911 * N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(0.254829592 + N[(N[(-0.284496736 + N[(N[(1.421413741 + N[(N[(-1.453152027 + N[(1.061405429 / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 / N[(t$95$0 * N[Exp[N[(x * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[Power[t$95$2, 2.0], $MachinePrecision]}, N[(N[(N[(1.0 / N[(1.0 + t$95$3), $MachinePrecision]), $MachinePrecision] + N[(N[Power[t$95$2, 4.0], $MachinePrecision] / N[(-1.0 - t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$1 * N[(N[Exp[(-N[(x * x), $MachinePrecision])], $MachinePrecision] / t$95$0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\\
t_1 := 0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{t\_0}}{t\_0}}{t\_0}}{t\_0}\\
t_2 := \frac{t\_1}{t\_0 \cdot e^{x \cdot x}}\\
t_3 := {t\_2}^{2}\\
\frac{\frac{1}{1 + t\_3} + \frac{{t\_2}^{4}}{-1 - t\_3}}{\mathsf{fma}\left(t\_1, \frac{e^{-x \cdot x}}{t\_0}, 1\right)}
\end{array}
\end{array}
Initial program 78.1%
Applied rewrites78.1%
lift--.f64N/A
flip--N/A
Applied rewrites78.1%
lift--.f64N/A
flip--N/A
Applied rewrites86.1%
Final simplification86.1%
(FPCore (x)
:precision binary64
(let* ((t_0 (fma 0.3275911 (fabs x) 1.0))
(t_1 (* t_0 (exp (* x x))))
(t_2 (/ 1.061405429 t_0))
(t_3
(+
0.254829592
(/
(+
-0.284496736
(/ (+ 1.421413741 (/ (+ -1.453152027 t_2) t_0)) t_0))
t_0))))
(/
(+ 1.0 (/ -1.0 (/ (pow t_1 3.0) (pow t_3 3.0))))
(fma
(/
(+
0.254829592
(/
(+
-0.284496736
(/
(+
1.421413741
(/
(/
(+ -3.0685496600615605 (/ 1.1957597040827899 (* t_0 (* t_0 t_0))))
(fma t_2 (- t_2 -1.453152027) 2.111650813574209))
t_0))
t_0))
t_0))
t_1)
(fma t_3 (/ (exp (- (* x x))) t_0) 1.0)
1.0))))
double code(double x) {
double t_0 = fma(0.3275911, fabs(x), 1.0);
double t_1 = t_0 * exp((x * x));
double t_2 = 1.061405429 / t_0;
double t_3 = 0.254829592 + ((-0.284496736 + ((1.421413741 + ((-1.453152027 + t_2) / t_0)) / t_0)) / t_0);
return (1.0 + (-1.0 / (pow(t_1, 3.0) / pow(t_3, 3.0)))) / fma(((0.254829592 + ((-0.284496736 + ((1.421413741 + (((-3.0685496600615605 + (1.1957597040827899 / (t_0 * (t_0 * t_0)))) / fma(t_2, (t_2 - -1.453152027), 2.111650813574209)) / t_0)) / t_0)) / t_0)) / t_1), fma(t_3, (exp(-(x * x)) / t_0), 1.0), 1.0);
}
function code(x) t_0 = fma(0.3275911, abs(x), 1.0) t_1 = Float64(t_0 * exp(Float64(x * x))) t_2 = Float64(1.061405429 / t_0) t_3 = Float64(0.254829592 + Float64(Float64(-0.284496736 + Float64(Float64(1.421413741 + Float64(Float64(-1.453152027 + t_2) / t_0)) / t_0)) / t_0)) return Float64(Float64(1.0 + Float64(-1.0 / Float64((t_1 ^ 3.0) / (t_3 ^ 3.0)))) / fma(Float64(Float64(0.254829592 + Float64(Float64(-0.284496736 + Float64(Float64(1.421413741 + Float64(Float64(Float64(-3.0685496600615605 + Float64(1.1957597040827899 / Float64(t_0 * Float64(t_0 * t_0)))) / fma(t_2, Float64(t_2 - -1.453152027), 2.111650813574209)) / t_0)) / t_0)) / t_0)) / t_1), fma(t_3, Float64(exp(Float64(-Float64(x * x))) / t_0), 1.0), 1.0)) end
code[x_] := Block[{t$95$0 = N[(0.3275911 * N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 * N[Exp[N[(x * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(1.061405429 / t$95$0), $MachinePrecision]}, Block[{t$95$3 = N[(0.254829592 + N[(N[(-0.284496736 + N[(N[(1.421413741 + N[(N[(-1.453152027 + t$95$2), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision]}, N[(N[(1.0 + N[(-1.0 / N[(N[Power[t$95$1, 3.0], $MachinePrecision] / N[Power[t$95$3, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(0.254829592 + N[(N[(-0.284496736 + N[(N[(1.421413741 + N[(N[(N[(-3.0685496600615605 + N[(1.1957597040827899 / N[(t$95$0 * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$2 * N[(t$95$2 - -1.453152027), $MachinePrecision] + 2.111650813574209), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision] * N[(t$95$3 * N[(N[Exp[(-N[(x * x), $MachinePrecision])], $MachinePrecision] / t$95$0), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\\
t_1 := t\_0 \cdot e^{x \cdot x}\\
t_2 := \frac{1.061405429}{t\_0}\\
t_3 := 0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + t\_2}{t\_0}}{t\_0}}{t\_0}\\
\frac{1 + \frac{-1}{\frac{{t\_1}^{3}}{{t\_3}^{3}}}}{\mathsf{fma}\left(\frac{0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{\frac{-3.0685496600615605 + \frac{1.1957597040827899}{t\_0 \cdot \left(t\_0 \cdot t\_0\right)}}{\mathsf{fma}\left(t\_2, t\_2 - -1.453152027, 2.111650813574209\right)}}{t\_0}}{t\_0}}{t\_0}}{t\_1}, \mathsf{fma}\left(t\_3, \frac{e^{-x \cdot x}}{t\_0}, 1\right), 1\right)}
\end{array}
\end{array}
Initial program 78.1%
Applied rewrites78.1%
lift--.f64N/A
flip3--N/A
Applied rewrites78.1%
lift-pow.f64N/A
lift-/.f64N/A
cube-divN/A
clear-numN/A
lower-/.f64N/A
Applied rewrites79.3%
lift-+.f64N/A
flip3-+N/A
lower-/.f64N/A
lower-+.f64N/A
metadata-evalN/A
lift-/.f64N/A
cube-divN/A
cube-unmultN/A
lift-*.f64N/A
lift-*.f64N/A
lower-/.f64N/A
metadata-evalN/A
Applied rewrites79.3%
Final simplification79.3%
(FPCore (x)
:precision binary64
(let* ((t_0 (fma 0.3275911 (fabs x) 1.0))
(t_1
(+
0.254829592
(/
(+
-0.284496736
(/
(+ 1.421413741 (/ (+ -1.453152027 (/ 1.061405429 t_0)) t_0))
t_0))
t_0)))
(t_2 (* t_0 (exp (* x x)))))
(/
(+ 1.0 (/ -1.0 (/ (pow t_2 3.0) (pow t_1 3.0))))
(fma (/ t_1 t_2) (fma t_1 (/ (exp (- (* x x))) t_0) 1.0) 1.0))))
double code(double x) {
double t_0 = fma(0.3275911, fabs(x), 1.0);
double t_1 = 0.254829592 + ((-0.284496736 + ((1.421413741 + ((-1.453152027 + (1.061405429 / t_0)) / t_0)) / t_0)) / t_0);
double t_2 = t_0 * exp((x * x));
return (1.0 + (-1.0 / (pow(t_2, 3.0) / pow(t_1, 3.0)))) / fma((t_1 / t_2), fma(t_1, (exp(-(x * x)) / t_0), 1.0), 1.0);
}
function code(x) t_0 = fma(0.3275911, abs(x), 1.0) t_1 = Float64(0.254829592 + Float64(Float64(-0.284496736 + Float64(Float64(1.421413741 + Float64(Float64(-1.453152027 + Float64(1.061405429 / t_0)) / t_0)) / t_0)) / t_0)) t_2 = Float64(t_0 * exp(Float64(x * x))) return Float64(Float64(1.0 + Float64(-1.0 / Float64((t_2 ^ 3.0) / (t_1 ^ 3.0)))) / fma(Float64(t_1 / t_2), fma(t_1, Float64(exp(Float64(-Float64(x * x))) / t_0), 1.0), 1.0)) end
code[x_] := Block[{t$95$0 = N[(0.3275911 * N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(0.254829592 + N[(N[(-0.284496736 + N[(N[(1.421413741 + N[(N[(-1.453152027 + N[(1.061405429 / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$0 * N[Exp[N[(x * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(N[(1.0 + N[(-1.0 / N[(N[Power[t$95$2, 3.0], $MachinePrecision] / N[Power[t$95$1, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$1 / t$95$2), $MachinePrecision] * N[(t$95$1 * N[(N[Exp[(-N[(x * x), $MachinePrecision])], $MachinePrecision] / t$95$0), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\\
t_1 := 0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{t\_0}}{t\_0}}{t\_0}}{t\_0}\\
t_2 := t\_0 \cdot e^{x \cdot x}\\
\frac{1 + \frac{-1}{\frac{{t\_2}^{3}}{{t\_1}^{3}}}}{\mathsf{fma}\left(\frac{t\_1}{t\_2}, \mathsf{fma}\left(t\_1, \frac{e^{-x \cdot x}}{t\_0}, 1\right), 1\right)}
\end{array}
\end{array}
Initial program 78.1%
Applied rewrites78.1%
lift--.f64N/A
flip3--N/A
Applied rewrites78.1%
lift-pow.f64N/A
lift-/.f64N/A
cube-divN/A
clear-numN/A
lower-/.f64N/A
Applied rewrites79.3%
Final simplification79.3%
(FPCore (x)
:precision binary64
(let* ((t_0 (fma 0.3275911 (fabs x) 1.0))
(t_1
(+
0.254829592
(/
(+
-0.284496736
(/
(+ 1.421413741 (/ (+ -1.453152027 (/ 1.061405429 t_0)) t_0))
t_0))
t_0))))
(/
(+ 1.0 (/ -1.0 (/ (* t_0 (* t_0 t_0)) (pow t_1 3.0))))
(fma
(/ t_1 (* t_0 (exp (* x x))))
(fma t_1 (/ (exp (- (* x x))) t_0) 1.0)
1.0))))
double code(double x) {
double t_0 = fma(0.3275911, fabs(x), 1.0);
double t_1 = 0.254829592 + ((-0.284496736 + ((1.421413741 + ((-1.453152027 + (1.061405429 / t_0)) / t_0)) / t_0)) / t_0);
return (1.0 + (-1.0 / ((t_0 * (t_0 * t_0)) / pow(t_1, 3.0)))) / fma((t_1 / (t_0 * exp((x * x)))), fma(t_1, (exp(-(x * x)) / t_0), 1.0), 1.0);
}
function code(x) t_0 = fma(0.3275911, abs(x), 1.0) t_1 = Float64(0.254829592 + Float64(Float64(-0.284496736 + Float64(Float64(1.421413741 + Float64(Float64(-1.453152027 + Float64(1.061405429 / t_0)) / t_0)) / t_0)) / t_0)) return Float64(Float64(1.0 + Float64(-1.0 / Float64(Float64(t_0 * Float64(t_0 * t_0)) / (t_1 ^ 3.0)))) / fma(Float64(t_1 / Float64(t_0 * exp(Float64(x * x)))), fma(t_1, Float64(exp(Float64(-Float64(x * x))) / t_0), 1.0), 1.0)) end
code[x_] := Block[{t$95$0 = N[(0.3275911 * N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(0.254829592 + N[(N[(-0.284496736 + N[(N[(1.421413741 + N[(N[(-1.453152027 + N[(1.061405429 / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision]}, N[(N[(1.0 + N[(-1.0 / N[(N[(t$95$0 * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision] / N[Power[t$95$1, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$1 / N[(t$95$0 * N[Exp[N[(x * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$1 * N[(N[Exp[(-N[(x * x), $MachinePrecision])], $MachinePrecision] / t$95$0), $MachinePrecision] + 1.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\\
t_1 := 0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{t\_0}}{t\_0}}{t\_0}}{t\_0}\\
\frac{1 + \frac{-1}{\frac{t\_0 \cdot \left(t\_0 \cdot t\_0\right)}{{t\_1}^{3}}}}{\mathsf{fma}\left(\frac{t\_1}{t\_0 \cdot e^{x \cdot x}}, \mathsf{fma}\left(t\_1, \frac{e^{-x \cdot x}}{t\_0}, 1\right), 1\right)}
\end{array}
\end{array}
Initial program 78.1%
Applied rewrites78.1%
lift--.f64N/A
flip3--N/A
Applied rewrites78.1%
lift-pow.f64N/A
lift-/.f64N/A
cube-divN/A
clear-numN/A
lower-/.f64N/A
Applied rewrites79.3%
Taylor expanded in x around 0
cube-multN/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-fabs.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-fabs.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-fabs.f6478.8
Applied rewrites78.8%
Final simplification78.8%
(FPCore (x)
:precision binary64
(let* ((t_0 (fma 0.3275911 (fabs x) 1.0))
(t_1
(/
(+
0.254829592
(/
(+
(/ (+ -1.453152027 (/ 1.061405429 t_0)) (* t_0 t_0))
(+ -0.284496736 (/ 1.421413741 t_0)))
t_0))
(* t_0 (exp (* x x))))))
(/ (- 1.0 (pow t_1 2.0)) (+ 1.0 t_1))))
double code(double x) {
double t_0 = fma(0.3275911, fabs(x), 1.0);
double t_1 = (0.254829592 + ((((-1.453152027 + (1.061405429 / t_0)) / (t_0 * t_0)) + (-0.284496736 + (1.421413741 / t_0))) / t_0)) / (t_0 * exp((x * x)));
return (1.0 - pow(t_1, 2.0)) / (1.0 + t_1);
}
function code(x) t_0 = fma(0.3275911, abs(x), 1.0) t_1 = Float64(Float64(0.254829592 + Float64(Float64(Float64(Float64(-1.453152027 + Float64(1.061405429 / t_0)) / Float64(t_0 * t_0)) + Float64(-0.284496736 + Float64(1.421413741 / t_0))) / t_0)) / Float64(t_0 * exp(Float64(x * x)))) return Float64(Float64(1.0 - (t_1 ^ 2.0)) / Float64(1.0 + t_1)) end
code[x_] := Block[{t$95$0 = N[(0.3275911 * N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(N[(0.254829592 + N[(N[(N[(N[(-1.453152027 + N[(1.061405429 / t$95$0), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision] + N[(-0.284496736 + N[(1.421413741 / t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * N[Exp[N[(x * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(N[(1.0 - N[Power[t$95$1, 2.0], $MachinePrecision]), $MachinePrecision] / N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\\
t_1 := \frac{0.254829592 + \frac{\frac{-1.453152027 + \frac{1.061405429}{t\_0}}{t\_0 \cdot t\_0} + \left(-0.284496736 + \frac{1.421413741}{t\_0}\right)}{t\_0}}{t\_0 \cdot e^{x \cdot x}}\\
\frac{1 - {t\_1}^{2}}{1 + t\_1}
\end{array}
\end{array}
Initial program 78.1%
lift-*.f64N/A
lift-+.f64N/A
distribute-lft-inN/A
+-commutativeN/A
*-commutativeN/A
lower-+.f64N/A
Applied rewrites78.1%
Applied rewrites78.2%
Final simplification78.2%
(FPCore (x)
:precision binary64
(let* ((t_0 (fma 0.3275911 (fabs x) 1.0)) (t_1 (* t_0 t_0)))
(-
1.0
(/
(+
0.254829592
(/
(+
(+ (/ 1.421413741 t_0) (/ 1.061405429 (* t_0 t_1)))
(+ -0.284496736 (/ -1.453152027 t_1)))
t_0))
(* t_0 (exp (* x x)))))))
double code(double x) {
double t_0 = fma(0.3275911, fabs(x), 1.0);
double t_1 = t_0 * t_0;
return 1.0 - ((0.254829592 + ((((1.421413741 / t_0) + (1.061405429 / (t_0 * t_1))) + (-0.284496736 + (-1.453152027 / t_1))) / t_0)) / (t_0 * exp((x * x))));
}
function code(x) t_0 = fma(0.3275911, abs(x), 1.0) t_1 = Float64(t_0 * t_0) return Float64(1.0 - Float64(Float64(0.254829592 + Float64(Float64(Float64(Float64(1.421413741 / t_0) + Float64(1.061405429 / Float64(t_0 * t_1))) + Float64(-0.284496736 + Float64(-1.453152027 / t_1))) / t_0)) / Float64(t_0 * exp(Float64(x * x))))) end
code[x_] := Block[{t$95$0 = N[(0.3275911 * N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 * t$95$0), $MachinePrecision]}, N[(1.0 - N[(N[(0.254829592 + N[(N[(N[(N[(1.421413741 / t$95$0), $MachinePrecision] + N[(1.061405429 / N[(t$95$0 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-0.284496736 + N[(-1.453152027 / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * N[Exp[N[(x * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\\
t_1 := t\_0 \cdot t\_0\\
1 - \frac{0.254829592 + \frac{\left(\frac{1.421413741}{t\_0} + \frac{1.061405429}{t\_0 \cdot t\_1}\right) + \left(-0.284496736 + \frac{-1.453152027}{t\_1}\right)}{t\_0}}{t\_0 \cdot e^{x \cdot x}}
\end{array}
\end{array}
Initial program 78.1%
Applied rewrites78.1%
Taylor expanded in x around 0
lower-/.f64N/A
Applied rewrites78.1%
Taylor expanded in x around 0
lower-/.f64N/A
Applied rewrites78.1%
Final simplification78.1%
(FPCore (x)
:precision binary64
(let* ((t_0 (fma 0.3275911 (fabs x) 1.0)))
(-
1.0
(/
(+
0.254829592
(/
(+
(/ (+ -1.453152027 (/ 1.061405429 t_0)) (* t_0 t_0))
(+ -0.284496736 (/ 1.421413741 t_0)))
t_0))
(* t_0 (exp (* x x)))))))
double code(double x) {
double t_0 = fma(0.3275911, fabs(x), 1.0);
return 1.0 - ((0.254829592 + ((((-1.453152027 + (1.061405429 / t_0)) / (t_0 * t_0)) + (-0.284496736 + (1.421413741 / t_0))) / t_0)) / (t_0 * exp((x * x))));
}
function code(x) t_0 = fma(0.3275911, abs(x), 1.0) return Float64(1.0 - Float64(Float64(0.254829592 + Float64(Float64(Float64(Float64(-1.453152027 + Float64(1.061405429 / t_0)) / Float64(t_0 * t_0)) + Float64(-0.284496736 + Float64(1.421413741 / t_0))) / t_0)) / Float64(t_0 * exp(Float64(x * x))))) end
code[x_] := Block[{t$95$0 = N[(0.3275911 * N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, N[(1.0 - N[(N[(0.254829592 + N[(N[(N[(N[(-1.453152027 + N[(1.061405429 / t$95$0), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision] + N[(-0.284496736 + N[(1.421413741 / t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * N[Exp[N[(x * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\\
1 - \frac{0.254829592 + \frac{\frac{-1.453152027 + \frac{1.061405429}{t\_0}}{t\_0 \cdot t\_0} + \left(-0.284496736 + \frac{1.421413741}{t\_0}\right)}{t\_0}}{t\_0 \cdot e^{x \cdot x}}
\end{array}
\end{array}
Initial program 78.1%
lift-*.f64N/A
lift-+.f64N/A
distribute-lft-inN/A
+-commutativeN/A
*-commutativeN/A
lower-+.f64N/A
Applied rewrites78.1%
Applied rewrites78.1%
Final simplification78.1%
(FPCore (x)
:precision binary64
(let* ((t_0 (fma 0.3275911 (fabs x) 1.0)))
(-
1.0
(/
(+
0.254829592
(/
(-
(+
(/ (+ -1.453152027 (/ 1.061405429 t_0)) (* t_0 t_0))
(/ 1.421413741 t_0))
0.284496736)
t_0))
(* t_0 (exp (* x x)))))))
double code(double x) {
double t_0 = fma(0.3275911, fabs(x), 1.0);
return 1.0 - ((0.254829592 + (((((-1.453152027 + (1.061405429 / t_0)) / (t_0 * t_0)) + (1.421413741 / t_0)) - 0.284496736) / t_0)) / (t_0 * exp((x * x))));
}
function code(x) t_0 = fma(0.3275911, abs(x), 1.0) return Float64(1.0 - Float64(Float64(0.254829592 + Float64(Float64(Float64(Float64(Float64(-1.453152027 + Float64(1.061405429 / t_0)) / Float64(t_0 * t_0)) + Float64(1.421413741 / t_0)) - 0.284496736) / t_0)) / Float64(t_0 * exp(Float64(x * x))))) end
code[x_] := Block[{t$95$0 = N[(0.3275911 * N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, N[(1.0 - N[(N[(0.254829592 + N[(N[(N[(N[(N[(-1.453152027 + N[(1.061405429 / t$95$0), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision] + N[(1.421413741 / t$95$0), $MachinePrecision]), $MachinePrecision] - 0.284496736), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * N[Exp[N[(x * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\\
1 - \frac{0.254829592 + \frac{\left(\frac{-1.453152027 + \frac{1.061405429}{t\_0}}{t\_0 \cdot t\_0} + \frac{1.421413741}{t\_0}\right) - 0.284496736}{t\_0}}{t\_0 \cdot e^{x \cdot x}}
\end{array}
\end{array}
Initial program 78.1%
Applied rewrites78.1%
Taylor expanded in x around 0
lower-/.f64N/A
Applied rewrites78.1%
Final simplification78.1%
(FPCore (x)
:precision binary64
(let* ((t_0 (fma 0.3275911 (fabs x) 1.0)))
(-
1.0
(/
(+
0.254829592
(/
(+
-0.284496736
(/ (+ 1.421413741 (/ (+ -1.453152027 (/ 1.061405429 t_0)) t_0)) t_0))
t_0))
(* t_0 (exp (* x x)))))))
double code(double x) {
double t_0 = fma(0.3275911, fabs(x), 1.0);
return 1.0 - ((0.254829592 + ((-0.284496736 + ((1.421413741 + ((-1.453152027 + (1.061405429 / t_0)) / t_0)) / t_0)) / t_0)) / (t_0 * exp((x * x))));
}
function code(x) t_0 = fma(0.3275911, abs(x), 1.0) return Float64(1.0 - Float64(Float64(0.254829592 + Float64(Float64(-0.284496736 + Float64(Float64(1.421413741 + Float64(Float64(-1.453152027 + Float64(1.061405429 / t_0)) / t_0)) / t_0)) / t_0)) / Float64(t_0 * exp(Float64(x * x))))) end
code[x_] := Block[{t$95$0 = N[(0.3275911 * N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, N[(1.0 - N[(N[(0.254829592 + N[(N[(-0.284496736 + N[(N[(1.421413741 + N[(N[(-1.453152027 + N[(1.061405429 / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * N[Exp[N[(x * x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\\
1 - \frac{0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{t\_0}}{t\_0}}{t\_0}}{t\_0}}{t\_0 \cdot e^{x \cdot x}}
\end{array}
\end{array}
Initial program 78.1%
Applied rewrites78.1%
(FPCore (x) :precision binary64 1.0)
double code(double x) {
return 1.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0
end function
public static double code(double x) {
return 1.0;
}
def code(x): return 1.0
function code(x) return 1.0 end
function tmp = code(x) tmp = 1.0; end
code[x_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 78.1%
Applied rewrites78.1%
Taylor expanded in x around inf
Applied rewrites53.5%
herbie shell --seed 2024233
(FPCore (x)
:name "Jmat.Real.erf"
:precision binary64
(- 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))