Jmat.Real.erf

Average Error: 13.4 → 0.2

Time: 16.6s

Precision: binary64

Cost: 133700

Math

\[1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|} \]

↓

\[\begin{array}{l} t_0 := \mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\\ t_1 := \frac{1}{1 + 0.3275911 \cdot \left|x\right|}\\ \mathbf{if}\;\left(t_1 \cdot \left(0.254829592 + t_1 \cdot \left(-0.284496736 + t_1 \cdot \left(1.421413741 + t_1 \cdot \left(-1.453152027 + t_1 \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-x \cdot x} \leq 0.995:\\ \;\;\;\;\mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{t_0}, -0.254829592 + \frac{0.284496736 + \frac{-1.421413741 + \left(\frac{-1.061405429}{{t_0}^{2}} + \frac{1.453152027}{t_0}\right)}{t_0}}{t_0}, 1\right)\\ \mathbf{else}:\\ \;\;\;\;10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \sqrt{x \cdot \left(x \cdot 1.2732557730789702\right)}\right)\right)\\ \end{array} \]

FPCore

(FPCore (x)
 :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)))))))

↓

(FPCore (x)
 :precision binary64
 (let* ((t_0 (fma 0.3275911 (fabs x) 1.0))
        (t_1 (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x))))))
   (if (<=
        (*
         (*
          t_1
          (+
           0.254829592
           (*
            t_1
            (+
             -0.284496736
             (*
              t_1
              (+ 1.421413741 (* t_1 (+ -1.453152027 (* t_1 1.061405429)))))))))
         (exp (- (* x x))))
        0.995)
     (fma
      (/ (pow (exp x) (- x)) t_0)
      (+
       -0.254829592
       (/
        (+
         0.284496736
         (/
          (+
           -1.421413741
           (+ (/ -1.061405429 (pow t_0 2.0)) (/ 1.453152027 t_0)))
          t_0))
        t_0))
      1.0)
     (+
      1e-9
      (+
       (* -0.00011824294398844343 (pow x 2.0))
       (+
        (* -0.37545125292247583 (pow x 3.0))
        (sqrt (* x (* x 1.2732557730789702)))))))))

C

double code(double x) {
	return 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))));
}

↓

double code(double x) {
	double t_0 = fma(0.3275911, fabs(x), 1.0);
	double t_1 = 1.0 / (1.0 + (0.3275911 * fabs(x)));
	double tmp;
	if (((t_1 * (0.254829592 + (t_1 * (-0.284496736 + (t_1 * (1.421413741 + (t_1 * (-1.453152027 + (t_1 * 1.061405429))))))))) * exp(-(x * x))) <= 0.995) {
		tmp = fma((pow(exp(x), -x) / t_0), (-0.254829592 + ((0.284496736 + ((-1.421413741 + ((-1.061405429 / pow(t_0, 2.0)) + (1.453152027 / t_0))) / t_0)) / t_0)), 1.0);
	} else {
		tmp = 1e-9 + ((-0.00011824294398844343 * pow(x, 2.0)) + ((-0.37545125292247583 * pow(x, 3.0)) + sqrt((x * (x * 1.2732557730789702)))));
	}
	return tmp;
}

Julia

function code(x)
	return Float64(1.0 - Float64(Float64(Float64(1.0 / Float64(1.0 + Float64(0.3275911 * abs(x)))) * Float64(0.254829592 + Float64(Float64(1.0 / Float64(1.0 + Float64(0.3275911 * abs(x)))) * Float64(-0.284496736 + Float64(Float64(1.0 / Float64(1.0 + Float64(0.3275911 * abs(x)))) * Float64(1.421413741 + Float64(Float64(1.0 / Float64(1.0 + Float64(0.3275911 * abs(x)))) * Float64(-1.453152027 + Float64(Float64(1.0 / Float64(1.0 + Float64(0.3275911 * abs(x)))) * 1.061405429))))))))) * exp(Float64(-Float64(abs(x) * abs(x))))))
end

↓

function code(x)
	t_0 = fma(0.3275911, abs(x), 1.0)
	t_1 = Float64(1.0 / Float64(1.0 + Float64(0.3275911 * abs(x))))
	tmp = 0.0
	if (Float64(Float64(t_1 * Float64(0.254829592 + Float64(t_1 * Float64(-0.284496736 + Float64(t_1 * Float64(1.421413741 + Float64(t_1 * Float64(-1.453152027 + Float64(t_1 * 1.061405429))))))))) * exp(Float64(-Float64(x * x)))) <= 0.995)
		tmp = fma(Float64((exp(x) ^ Float64(-x)) / t_0), Float64(-0.254829592 + Float64(Float64(0.284496736 + Float64(Float64(-1.421413741 + Float64(Float64(-1.061405429 / (t_0 ^ 2.0)) + Float64(1.453152027 / t_0))) / t_0)) / t_0)), 1.0);
	else
		tmp = Float64(1e-9 + Float64(Float64(-0.00011824294398844343 * (x ^ 2.0)) + Float64(Float64(-0.37545125292247583 * (x ^ 3.0)) + sqrt(Float64(x * Float64(x * 1.2732557730789702))))));
	end
	return tmp
end

Wolfram

code[x_] := N[(1.0 - N[(N[(N[(1.0 / N[(1.0 + N[(0.3275911 * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(0.254829592 + N[(N[(1.0 / N[(1.0 + N[(0.3275911 * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(-0.284496736 + N[(N[(1.0 / N[(1.0 + N[(0.3275911 * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.421413741 + N[(N[(1.0 / N[(1.0 + N[(0.3275911 * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(-1.453152027 + N[(N[(1.0 / N[(1.0 + N[(0.3275911 * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 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]

↓

code[x_] := Block[{t$95$0 = N[(0.3275911 * N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(1.0 / N[(1.0 + N[(0.3275911 * N[Abs[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(t$95$1 * N[(0.254829592 + N[(t$95$1 * N[(-0.284496736 + N[(t$95$1 * N[(1.421413741 + N[(t$95$1 * N[(-1.453152027 + N[(t$95$1 * 1.061405429), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Exp[(-N[(x * x), $MachinePrecision])], $MachinePrecision]), $MachinePrecision], 0.995], N[(N[(N[Power[N[Exp[x], $MachinePrecision], (-x)], $MachinePrecision] / t$95$0), $MachinePrecision] * N[(-0.254829592 + N[(N[(0.284496736 + N[(N[(-1.421413741 + N[(N[(-1.061405429 / N[Power[t$95$0, 2.0], $MachinePrecision]), $MachinePrecision] + N[(1.453152027 / t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision], N[(1e-9 + N[(N[(-0.00011824294398844343 * N[Power[x, 2.0], $MachinePrecision]), $MachinePrecision] + N[(N[(-0.37545125292247583 * N[Power[x, 3.0], $MachinePrecision]), $MachinePrecision] + N[Sqrt[N[(x * N[(x * 1.2732557730789702), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]

Derivation

1. Split input into 2 regimes

2. if (*.f64 (*.f64 (/.f64 1 (+.f64 1 (*.f64 3275911/10000000 (fabs.f64 x)))) (+.f64 31853699/125000000 (*.f64 (/.f64 1 (+.f64 1 (*.f64 3275911/10000000 (fabs.f64 x)))) (+.f64 -8890523/31250000 (*.f64 (/.f64 1 (+.f64 1 (*.f64 3275911/10000000 (fabs.f64 x)))) (+.f64 1421413741/1000000000 (*.f64 (/.f64 1 (+.f64 1 (*.f64 3275911/10000000 (fabs.f64 x)))) (+.f64 -1453152027/1000000000 (*.f64 (/.f64 1 (+.f64 1 (*.f64 3275911/10000000 (fabs.f64 x)))) 1061405429/1000000000))))))))) (exp.f64 (neg.f64 (*.f64 (fabs.f64 x) (fabs.f64 x))))) < 0.994999999999999996

   1. Initial program 0.0

      \[1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|} \]

   2. Simplified 0.0

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right)} \]
      Proof

      [Start] 0.0	\[ 1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|} \]
      cancel-sign-sub-inv [=>] 0.0	\[ \color{blue}{1 + \left(-\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|}} \]
      +-commutative [=>] 0.0	\[ \color{blue}{\left(-\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|} + 1} \]

   3. Taylor expanded in x around 0 0.0

      \[\leadsto \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{\color{blue}{\left(1.421413741 + 1.061405429 \cdot \frac{1}{{\left(0.3275911 \cdot \left|x\right| + 1\right)}^{2}}\right) - 1.453152027 \cdot \frac{1}{0.3275911 \cdot \left|x\right| + 1}}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]

   4. Simplified 0.0

      \[\leadsto \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{\color{blue}{1.421413741 + \left(\frac{1.061405429}{{\left(\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\right)}^{2}} + \frac{-1.453152027}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]
      Proof

      [Start] 0.0	\[ \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{\left(1.421413741 + 1.061405429 \cdot \frac{1}{{\left(0.3275911 \cdot \left|x\right| + 1\right)}^{2}}\right) - 1.453152027 \cdot \frac{1}{0.3275911 \cdot \left|x\right| + 1}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]
      associate--l+ [=>] 0.0	\[ \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{\color{blue}{1.421413741 + \left(1.061405429 \cdot \frac{1}{{\left(0.3275911 \cdot \left|x\right| + 1\right)}^{2}} - 1.453152027 \cdot \frac{1}{0.3275911 \cdot \left|x\right| + 1}\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]
      +-commutative [=>] 0.0	\[ \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{1.421413741 + \left(1.061405429 \cdot \frac{1}{{\color{blue}{\left(1 + 0.3275911 \cdot \left|x\right|\right)}}^{2}} - 1.453152027 \cdot \frac{1}{0.3275911 \cdot \left|x\right| + 1}\right)}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]
      +-commutative [=>] 0.0	\[ \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{1.421413741 + \left(1.061405429 \cdot \frac{1}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{2}} - 1.453152027 \cdot \frac{1}{\color{blue}{1 + 0.3275911 \cdot \left|x\right|}}\right)}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]
      cancel-sign-sub-inv [=>] 0.0	\[ \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{1.421413741 + \color{blue}{\left(1.061405429 \cdot \frac{1}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{2}} + \left(-1.453152027\right) \cdot \frac{1}{1 + 0.3275911 \cdot \left|x\right|}\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]
      associate-*r/ [=>] 0.0	\[ \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{1.421413741 + \left(\color{blue}{\frac{1.061405429 \cdot 1}{{\left(1 + 0.3275911 \cdot \left|x\right|\right)}^{2}}} + \left(-1.453152027\right) \cdot \frac{1}{1 + 0.3275911 \cdot \left|x\right|}\right)}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]
      +-commutative [<=] 0.0	\[ \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{1.421413741 + \left(\frac{1.061405429 \cdot 1}{{\color{blue}{\left(0.3275911 \cdot \left|x\right| + 1\right)}}^{2}} + \left(-1.453152027\right) \cdot \frac{1}{1 + 0.3275911 \cdot \left|x\right|}\right)}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]
      metadata-eval [=>] 0.0	\[ \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{1.421413741 + \left(\frac{\color{blue}{1.061405429}}{{\left(0.3275911 \cdot \left|x\right| + 1\right)}^{2}} + \left(-1.453152027\right) \cdot \frac{1}{1 + 0.3275911 \cdot \left|x\right|}\right)}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]
      fma-def [=>] 0.0	\[ \mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 - \frac{1.421413741 + \left(\frac{1.061405429}{{\color{blue}{\left(\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\right)}}^{2}} + \left(-1.453152027\right) \cdot \frac{1}{1 + 0.3275911 \cdot \left|x\right|}\right)}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right) \]

   if 0.994999999999999996 < (*.f64 (*.f64 (/.f64 1 (+.f64 1 (*.f64 3275911/10000000 (fabs.f64 x)))) (+.f64 31853699/125000000 (*.f64 (/.f64 1 (+.f64 1 (*.f64 3275911/10000000 (fabs.f64 x)))) (+.f64 -8890523/31250000 (*.f64 (/.f64 1 (+.f64 1 (*.f64 3275911/10000000 (fabs.f64 x)))) (+.f64 1421413741/1000000000 (*.f64 (/.f64 1 (+.f64 1 (*.f64 3275911/10000000 (fabs.f64 x)))) (+.f64 -1453152027/1000000000 (*.f64 (/.f64 1 (+.f64 1 (*.f64 3275911/10000000 (fabs.f64 x)))) 1061405429/1000000000))))))))) (exp.f64 (neg.f64 (*.f64 (fabs.f64 x) (fabs.f64 x)))))

   1. Initial program 26.9

      \[1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|} \]

   2. Simplified 26.9

      \[\leadsto \color{blue}{1 - \frac{\frac{0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{{\left(e^{x}\right)}^{x}}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}} \]
      Proof

      [Start] 26.9	\[ 1 - \left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|} \]
      associate-*l/ [=>] 26.9	\[ 1 - \color{blue}{\frac{1 \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)}{1 + 0.3275911 \cdot \left|x\right|}} \cdot e^{-\left|x\right| \cdot \left|x\right|} \]
      associate-*l/ [=>] 26.9	\[ 1 - \color{blue}{\frac{\left(1 \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-\left|x\right| \cdot \left|x\right|}}{1 + 0.3275911 \cdot \left|x\right|}} \]

   3. Applied egg-rr 27.6

      \[\leadsto \color{blue}{1 \cdot \left(1 - \frac{0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right) \cdot {\left(e^{x}\right)}^{x}}\right)} \]

   4. Simplified 27.6

      \[\leadsto \color{blue}{1 - \frac{0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{e^{x \cdot x} \cdot \mathsf{fma}\left(0.3275911, x, 1\right)}} \]
      Proof

      [Start] 27.6	\[ 1 \cdot \left(1 - \frac{0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right) \cdot {\left(e^{x}\right)}^{x}}\right) \]
      *-lft-identity [=>] 27.6	\[ \color{blue}{1 - \frac{0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right) \cdot {\left(e^{x}\right)}^{x}}} \]
      *-commutative [=>] 27.6	\[ 1 - \frac{0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\color{blue}{{\left(e^{x}\right)}^{x} \cdot \mathsf{fma}\left(0.3275911, x, 1\right)}} \]
      exp-prod [<=] 27.6	\[ 1 - \frac{0.254829592 + \frac{-0.284496736 + \frac{1.421413741 + \frac{-1.453152027 + \frac{1.061405429}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\mathsf{fma}\left(0.3275911, x, 1\right)}}{\color{blue}{e^{x \cdot x}} \cdot \mathsf{fma}\left(0.3275911, x, 1\right)} \]

   5. Taylor expanded in x around 0 1.6

      \[\leadsto \color{blue}{10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + 1.128386358070218 \cdot x\right)\right)} \]

   6. Applied egg-rr 0.3

      \[\leadsto 10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \color{blue}{\sqrt{1.2732557730789702 \cdot \left(x \cdot x\right)}}\right)\right) \]

   7. Simplified 0.3

      \[\leadsto 10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \color{blue}{\sqrt{x \cdot \left(x \cdot 1.2732557730789702\right)}}\right)\right) \]
      Proof

      [Start] 0.3	\[ 10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \sqrt{1.2732557730789702 \cdot \left(x \cdot x\right)}\right)\right) \]
      *-commutative [=>] 0.3	\[ 10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \sqrt{\color{blue}{\left(x \cdot x\right) \cdot 1.2732557730789702}}\right)\right) \]
      associate-*l* [=>] 0.3	\[ 10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \sqrt{\color{blue}{x \cdot \left(x \cdot 1.2732557730789702\right)}}\right)\right) \]

3. Recombined 2 regimes into one program.

4. Final simplification 0.2

   \[\leadsto \begin{array}{l} \mathbf{if}\;\left(\frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(0.254829592 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-0.284496736 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(1.421413741 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot \left(-1.453152027 + \frac{1}{1 + 0.3275911 \cdot \left|x\right|} \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-x \cdot x} \leq 0.995:\\ \;\;\;\;\mathsf{fma}\left(\frac{{\left(e^{x}\right)}^{\left(-x\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, -0.254829592 + \frac{0.284496736 + \frac{-1.421413741 + \left(\frac{-1.061405429}{{\left(\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)\right)}^{2}} + \frac{1.453152027}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}\right)}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}}{\mathsf{fma}\left(0.3275911, \left|x\right|, 1\right)}, 1\right)\\ \mathbf{else}:\\ \;\;\;\;10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \sqrt{x \cdot \left(x \cdot 1.2732557730789702\right)}\right)\right)\\ \end{array} \]

Alternatives

Alternative 1
Error	0.2
Cost	96004

\[\begin{array}{l} t_0 := 1 + 0.3275911 \cdot \left|x\right|\\ t_1 := \frac{1}{t_0}\\ \mathbf{if}\;\left(t_1 \cdot \left(0.254829592 + t_1 \cdot \left(-0.284496736 + t_1 \cdot \left(1.421413741 + t_1 \cdot \left(-1.453152027 + t_1 \cdot 1.061405429\right)\right)\right)\right)\right) \cdot e^{-x \cdot x} \leq 0.995:\\ \;\;\;\;\log \left(e^{1 + \frac{\frac{-0.254829592 + \frac{0.284496736 + \frac{-1.421413741 + \frac{1.453152027 + \frac{-1.061405429}{t_0}}{t_0}}{t_0}}{t_0}}{e^{x \cdot x}}}{t_0}}\right)\\ \mathbf{else}:\\ \;\;\;\;10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \sqrt{x \cdot \left(x \cdot 1.2732557730789702\right)}\right)\right)\\ \end{array} \]

Alternative 2
Error	0.1
Cost	61316

\[\begin{array}{l} t_0 := 1 + 0.3275911 \cdot \left|x\right|\\ \mathbf{if}\;\left|x\right| \leq 2 \cdot 10^{-6}:\\ \;\;\;\;10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \sqrt{x \cdot \left(x \cdot 1.2732557730789702\right)}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;1 + \frac{\frac{-0.254829592 + \frac{0.284496736 + \frac{-1.421413741 + \frac{\frac{2.111650813574209 + {t_0}^{-2} \cdot -1.126581484710674}{\frac{1.061405429}{t_0} + 1.453152027}}{t_0}}{t_0}}{t_0}}{e^{x \cdot x}}}{t_0}\\ \end{array} \]

Alternative 3
Error	0.2
Cost	41412

\[\begin{array}{l} t_0 := 1 + 0.3275911 \cdot \left|x\right|\\ \mathbf{if}\;x \leq -7.4 \cdot 10^{-6}:\\ \;\;\;\;1 + \frac{\frac{-0.254829592 + \frac{0.284496736 + \frac{-1.421413741 + \frac{1.453152027 + \frac{-1.061405429}{t_0}}{t_0}}{t_0}}{t_0}}{e^{x \cdot x}}}{t_0}\\ \mathbf{elif}\;x \leq 1.05:\\ \;\;\;\;10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \sqrt{x \cdot \left(x \cdot 1.2732557730789702\right)}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;1\\ \end{array} \]

Alternative 4
Error	0.4
Cost	26820

\[\begin{array}{l} \mathbf{if}\;\left|x\right| \leq 0.01:\\ \;\;\;\;10^{-9} + \left(-0.00011824294398844343 \cdot {x}^{2} + \left(-0.37545125292247583 \cdot {x}^{3} + \sqrt{x \cdot \left(x \cdot 1.2732557730789702\right)}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;1\\ \end{array} \]

Alternative 5
Error	0.4
Cost	20232

\[\begin{array}{l} \mathbf{if}\;x \leq -1.1:\\ \;\;\;\;1\\ \mathbf{elif}\;x \leq 1.05:\\ \;\;\;\;10^{-9} + \left|\mathsf{fma}\left(-0.37545125292247583, {x}^{3}, x \cdot \left(1.128386358070218 + x \cdot -0.00011824294398844343\right)\right)\right|\\ \mathbf{else}:\\ \;\;\;\;1\\ \end{array} \]

Alternative 6
Error	0.5
Cost	19848

\[\begin{array}{l} \mathbf{if}\;x \leq -0.88:\\ \;\;\;\;1\\ \mathbf{elif}\;x \leq 0.88:\\ \;\;\;\;10^{-9} + \sqrt{{\left(\sqrt[3]{x \cdot 1.128386358070218}\right)}^{6}}\\ \mathbf{else}:\\ \;\;\;\;1\\ \end{array} \]

Alternative 7
Error	0.5
Cost	7112

\[\begin{array}{l} \mathbf{if}\;x \leq -0.88:\\ \;\;\;\;1\\ \mathbf{elif}\;x \leq 0.88:\\ \;\;\;\;10^{-9} + \sqrt{x \cdot \left(x \cdot 1.2732557730789702\right)}\\ \mathbf{else}:\\ \;\;\;\;1\\ \end{array} \]

Alternative 8
Error	1.1
Cost	840

\[\begin{array}{l} \mathbf{if}\;x \leq -8.8 \cdot 10^{-10}:\\ \;\;\;\;1\\ \mathbf{elif}\;x \leq 0.88:\\ \;\;\;\;10^{-9} + x \cdot \left(1.128386358070218 + x \cdot -0.00011824294398844343\right)\\ \mathbf{else}:\\ \;\;\;\;1\\ \end{array} \]

Alternative 9
Error	1.1
Cost	584

\[\begin{array}{l} \mathbf{if}\;x \leq -8.8 \cdot 10^{-10}:\\ \;\;\;\;1\\ \mathbf{elif}\;x \leq 0.88:\\ \;\;\;\;10^{-9} + x \cdot 1.128386358070218\\ \mathbf{else}:\\ \;\;\;\;1\\ \end{array} \]

Alternative 10
Error	1.6
Cost	328

\[\begin{array}{l} \mathbf{if}\;x \leq -2.8 \cdot 10^{-5}:\\ \;\;\;\;1\\ \mathbf{elif}\;x \leq 2.8 \cdot 10^{-5}:\\ \;\;\;\;10^{-9}\\ \mathbf{else}:\\ \;\;\;\;1\\ \end{array} \]

Alternative 11
Error	29.8
Cost	64

\[10^{-9} \]

Reproduce

herbie shell --seed 2023038 
(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)))))))