Jmat.Real.erfi, branch x less than or equal to 0.5

Percentage Accurate: 99.8% → 99.8%
Time: 9.7s
Alternatives: 9
Speedup: 0.5×

Specification

?
\[x \leq 0.5\]
\[\begin{array}{l} \\ \begin{array}{l} t_0 := \left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\\ t_1 := \left(t\_0 \cdot \left|x\right|\right) \cdot \left|x\right|\\ \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot t\_0\right) + \frac{1}{5} \cdot t\_1\right) + \frac{1}{21} \cdot \left(\left(t\_1 \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \end{array} \end{array} \]
(FPCore (x)
 :precision binary64
 (let* ((t_0 (* (* (fabs x) (fabs x)) (fabs x)))
        (t_1 (* (* t_0 (fabs x)) (fabs x))))
   (fabs
    (*
     (/ 1.0 (sqrt (PI)))
     (+
      (+ (+ (* 2.0 (fabs x)) (* (/ 2.0 3.0) t_0)) (* (/ 1.0 5.0) t_1))
      (* (/ 1.0 21.0) (* (* t_1 (fabs x)) (fabs x))))))))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\\
t_1 := \left(t\_0 \cdot \left|x\right|\right) \cdot \left|x\right|\\
\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot t\_0\right) + \frac{1}{5} \cdot t\_1\right) + \frac{1}{21} \cdot \left(\left(t\_1 \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|
\end{array}
\end{array}

Sampling outcomes in binary64 precision:

Local Percentage Accuracy vs ?

The average percentage accuracy by input value. Horizontal axis shows value of an input variable; the variable is choosen in the title. Vertical axis is accuracy; higher is better. Red represent the original program, while blue represents Herbie's suggestion. These can be toggled with buttons below the plot. The line is an average while dots represent individual samples.

Accuracy vs Speed?

Herbie found 9 alternatives:

AlternativeAccuracySpeedup
The accuracy (vertical axis) and speed (horizontal axis) of each alternatives. Up and to the right is better. The red square shows the initial program, and each blue circle shows an alternative.The line shows the best available speed-accuracy tradeoffs.

Initial Program: 99.8% accurate, 1.0× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\\ t_1 := \left(t\_0 \cdot \left|x\right|\right) \cdot \left|x\right|\\ \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot t\_0\right) + \frac{1}{5} \cdot t\_1\right) + \frac{1}{21} \cdot \left(\left(t\_1 \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \end{array} \end{array} \]
(FPCore (x)
 :precision binary64
 (let* ((t_0 (* (* (fabs x) (fabs x)) (fabs x)))
        (t_1 (* (* t_0 (fabs x)) (fabs x))))
   (fabs
    (*
     (/ 1.0 (sqrt (PI)))
     (+
      (+ (+ (* 2.0 (fabs x)) (* (/ 2.0 3.0) t_0)) (* (/ 1.0 5.0) t_1))
      (* (/ 1.0 21.0) (* (* t_1 (fabs x)) (fabs x))))))))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\\
t_1 := \left(t\_0 \cdot \left|x\right|\right) \cdot \left|x\right|\\
\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot t\_0\right) + \frac{1}{5} \cdot t\_1\right) + \frac{1}{21} \cdot \left(\left(t\_1 \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|
\end{array}
\end{array}

Alternative 1: 99.8% accurate, 0.3× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := {\mathsf{PI}\left(\right)}^{-0.5}\\ \left|\mathsf{fma}\left(\mathsf{fma}\left(0.2, {x}^{5}, \mathsf{fma}\left(0.6666666666666666, {x}^{3}, 2 \cdot x\right)\right), t\_0, \left(t\_0 \cdot 0.047619047619047616\right) \cdot {x}^{7}\right)\right| \end{array} \end{array} \]
(FPCore (x)
 :precision binary64
 (let* ((t_0 (pow (PI) -0.5)))
   (fabs
    (fma
     (fma 0.2 (pow x 5.0) (fma 0.6666666666666666 (pow x 3.0) (* 2.0 x)))
     t_0
     (* (* t_0 0.047619047619047616) (pow x 7.0))))))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := {\mathsf{PI}\left(\right)}^{-0.5}\\
\left|\mathsf{fma}\left(\mathsf{fma}\left(0.2, {x}^{5}, \mathsf{fma}\left(0.6666666666666666, {x}^{3}, 2 \cdot x\right)\right), t\_0, \left(t\_0 \cdot 0.047619047619047616\right) \cdot {x}^{7}\right)\right|
\end{array}
\end{array}
Derivation
  1. Initial program 99.9%

    \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  2. Add Preprocessing
  3. Applied rewrites99.9%

    \[\leadsto \left|\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(0.2, {x}^{5}, \mathsf{fma}\left(0.6666666666666666, {x}^{3}, 2 \cdot x\right)\right), {\mathsf{PI}\left(\right)}^{-0.5}, \left({\mathsf{PI}\left(\right)}^{-0.5} \cdot 0.047619047619047616\right) \cdot {x}^{7}\right)}\right| \]
  4. Add Preprocessing

Alternative 2: 99.8% accurate, 0.4× speedup?

\[\begin{array}{l} \\ \left|{\left(\sqrt{\mathsf{PI}\left(\right)}\right)}^{-1} \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) + {5}^{-1} \cdot \left|\left(\left(\left(x \cdot x\right) \cdot x\right) \cdot x\right) \cdot x\right|\right) + {21}^{-1} \cdot \left(\left({x}^{5} \cdot x\right) \cdot \left|x\right|\right)\right)\right| \end{array} \]
(FPCore (x)
 :precision binary64
 (fabs
  (*
   (pow (sqrt (PI)) -1.0)
   (+
    (+
     (+ (* 2.0 (fabs x)) (* (/ 2.0 3.0) (* (* x x) (fabs x))))
     (* (pow 5.0 -1.0) (fabs (* (* (* (* x x) x) x) x))))
    (* (pow 21.0 -1.0) (* (* (pow x 5.0) x) (fabs x)))))))
\begin{array}{l}

\\
\left|{\left(\sqrt{\mathsf{PI}\left(\right)}\right)}^{-1} \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) + {5}^{-1} \cdot \left|\left(\left(\left(x \cdot x\right) \cdot x\right) \cdot x\right) \cdot x\right|\right) + {21}^{-1} \cdot \left(\left({x}^{5} \cdot x\right) \cdot \left|x\right|\right)\right)\right|
\end{array}
Derivation
  1. Initial program 99.9%

    \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  2. Add Preprocessing
  3. Step-by-step derivation
    1. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\color{blue}{\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    2. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\color{blue}{\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    3. associate-*l*N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\color{blue}{\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \left|x\right|\right)\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    4. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\color{blue}{\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)} \cdot \left(\left|x\right| \cdot \left|x\right|\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    5. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\color{blue}{\left(\left|x\right| \cdot \left|x\right|\right)} \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \left|x\right|\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    6. pow3N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\color{blue}{{\left(\left|x\right|\right)}^{3}} \cdot \left(\left|x\right| \cdot \left|x\right|\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    7. pow2N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left({\left(\left|x\right|\right)}^{3} \cdot \color{blue}{{\left(\left|x\right|\right)}^{2}}\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    8. pow-prod-upN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\color{blue}{{\left(\left|x\right|\right)}^{\left(3 + 2\right)}} \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    9. metadata-evalN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left({\left(\left|x\right|\right)}^{\color{blue}{5}} \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    10. lower-pow.f6499.9

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\color{blue}{{\left(\left|x\right|\right)}^{5}} \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    11. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left({\color{blue}{\left(\left|x\right|\right)}}^{5} \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    12. rem-sqrt-square-revN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left({\color{blue}{\left(\sqrt{x \cdot x}\right)}}^{5} \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    13. sqrt-prodN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left({\color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)}}^{5} \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    14. rem-square-sqrt78.1

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left({\color{blue}{x}}^{5} \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    15. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left({x}^{5} \cdot \color{blue}{\left|x\right|}\right) \cdot \left|x\right|\right)\right)\right| \]
    16. rem-sqrt-square-revN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left({x}^{5} \cdot \color{blue}{\sqrt{x \cdot x}}\right) \cdot \left|x\right|\right)\right)\right| \]
    17. sqrt-prodN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left({x}^{5} \cdot \color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)}\right) \cdot \left|x\right|\right)\right)\right| \]
    18. rem-square-sqrt99.9

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left({x}^{5} \cdot \color{blue}{x}\right) \cdot \left|x\right|\right)\right)\right| \]
  4. Applied rewrites99.9%

    \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\color{blue}{\left({x}^{5} \cdot x\right)} \cdot \left|x\right|\right)\right)\right| \]
  5. Final simplification99.9%

    \[\leadsto \left|{\left(\sqrt{\mathsf{PI}\left(\right)}\right)}^{-1} \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) + {5}^{-1} \cdot \left|\left(\left(\left(x \cdot x\right) \cdot x\right) \cdot x\right) \cdot x\right|\right) + {21}^{-1} \cdot \left(\left({x}^{5} \cdot x\right) \cdot \left|x\right|\right)\right)\right| \]
  6. Add Preprocessing

Alternative 3: 75.9% accurate, 0.4× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\ \left|{\left(\sqrt{\mathsf{PI}\left(\right)}\right)}^{-1} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot 0.6666666666666666\right)\right) + {5}^{-1} \cdot \left(t\_0 \cdot \left|x\right|\right)\right) + {21}^{-1} \cdot \left(\left|\left(t\_0 \cdot x\right) \cdot x\right| \cdot \left|x\right|\right)\right)\right| \end{array} \end{array} \]
(FPCore (x)
 :precision binary64
 (let* ((t_0 (* (* x x) (* x x))))
   (fabs
    (*
     (pow (sqrt (PI)) -1.0)
     (+
      (+
       (+ (* 2.0 (fabs x)) (* (* x x) (* x 0.6666666666666666)))
       (* (pow 5.0 -1.0) (* t_0 (fabs x))))
      (* (pow 21.0 -1.0) (* (fabs (* (* t_0 x) x)) (fabs x))))))))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \left(x \cdot x\right) \cdot \left(x \cdot x\right)\\
\left|{\left(\sqrt{\mathsf{PI}\left(\right)}\right)}^{-1} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot 0.6666666666666666\right)\right) + {5}^{-1} \cdot \left(t\_0 \cdot \left|x\right|\right)\right) + {21}^{-1} \cdot \left(\left|\left(t\_0 \cdot x\right) \cdot x\right| \cdot \left|x\right|\right)\right)\right|
\end{array}
\end{array}
Derivation
  1. Initial program 99.9%

    \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  2. Add Preprocessing
  3. Step-by-step derivation
    1. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    2. *-commutativeN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \frac{2}{3}}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    3. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)} \cdot \frac{2}{3}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    4. associate-*l*N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    5. lower-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    6. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(\color{blue}{\left|x\right|} \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    7. rem-sqrt-square-revN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(\color{blue}{\sqrt{x \cdot x}} \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    8. sqrt-prodN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(\color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)} \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    9. rem-square-sqrtN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(\color{blue}{x} \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    10. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot \color{blue}{\left|x\right|}\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    11. rem-sqrt-square-revN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot \color{blue}{\sqrt{x \cdot x}}\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    12. sqrt-prodN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot \color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)}\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    13. rem-square-sqrtN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot \color{blue}{x}\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    14. lower-*.f6499.9

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \color{blue}{\left(\left|x\right| \cdot \frac{2}{3}\right)}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    15. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(\color{blue}{\left|x\right|} \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    16. rem-sqrt-square-revN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(\color{blue}{\sqrt{x \cdot x}} \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    17. sqrt-prodN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(\color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)} \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    18. rem-square-sqrt81.2

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(\color{blue}{x} \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    19. lift-/.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \color{blue}{\frac{2}{3}}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    20. metadata-eval81.2

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \color{blue}{0.6666666666666666}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  4. Applied rewrites81.2%

    \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\left(x \cdot x\right) \cdot \left(x \cdot 0.6666666666666666\right)}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  5. Step-by-step derivation
    1. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\color{blue}{\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    2. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\color{blue}{\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    3. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\color{blue}{\left(\left|x\right| \cdot \left|x\right|\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    4. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\color{blue}{\left|x\right|} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    5. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \color{blue}{\left|x\right|}\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    6. sqr-abs-revN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\color{blue}{\left(x \cdot x\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    7. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\color{blue}{\left(x \cdot x\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    8. associate-*l*N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\color{blue}{\left(\left(x \cdot x\right) \cdot \left(\left|x\right| \cdot \left|x\right|\right)\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    9. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(\color{blue}{\left|x\right|} \cdot \left|x\right|\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    10. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(\left|x\right| \cdot \color{blue}{\left|x\right|}\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    11. sqr-abs-revN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \color{blue}{\left(x \cdot x\right)}\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    12. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \color{blue}{\left(x \cdot x\right)}\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    13. lower-*.f6481.2

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot 0.6666666666666666\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\color{blue}{\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  6. Applied rewrites81.2%

    \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot 0.6666666666666666\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\color{blue}{\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  7. Step-by-step derivation
    1. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\color{blue}{\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)} \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    2. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\color{blue}{\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    3. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\color{blue}{\left(\left|x\right| \cdot \left|x\right|\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    4. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\color{blue}{\left|x\right|} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    5. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \color{blue}{\left|x\right|}\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    6. sqr-abs-revN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\color{blue}{\left(x \cdot x\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    7. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\color{blue}{\left(x \cdot x\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    8. associate-*l*N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\color{blue}{\left(\left(x \cdot x\right) \cdot \left(\left|x\right| \cdot \left|x\right|\right)\right)} \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    9. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\color{blue}{\left|x\right|} \cdot \left|x\right|\right)\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    10. lift-fabs.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left|x\right| \cdot \color{blue}{\left|x\right|}\right)\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    11. sqr-abs-revN/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(x \cdot x\right) \cdot \color{blue}{\left(x \cdot x\right)}\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    12. lift-*.f64N/A

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(x \cdot x\right) \cdot \color{blue}{\left(x \cdot x\right)}\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    13. lift-*.f6481.2

      \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot 0.6666666666666666\right)\right) + \frac{1}{5} \cdot \left(\color{blue}{\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)} \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  8. Applied rewrites81.2%

    \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot 0.6666666666666666\right)\right) + \frac{1}{5} \cdot \left(\color{blue}{\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)} \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  9. Final simplification81.2%

    \[\leadsto \left|{\left(\sqrt{\mathsf{PI}\left(\right)}\right)}^{-1} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot 0.6666666666666666\right)\right) + {5}^{-1} \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left|x\right|\right)\right) + {21}^{-1} \cdot \left(\left|\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot x\right) \cdot x\right| \cdot \left|x\right|\right)\right)\right| \]
  10. Add Preprocessing

Alternative 4: 99.8% accurate, 0.5× speedup?

\[\begin{array}{l} \\ \left|{\left(\sqrt{\mathsf{PI}\left(\right)}\right)}^{-1} \cdot \mathsf{fma}\left({x}^{5}, 0.2, \mathsf{fma}\left(x, \mathsf{fma}\left(0.6666666666666666, x \cdot x, 2\right), {x}^{7} \cdot 0.047619047619047616\right)\right)\right| \end{array} \]
(FPCore (x)
 :precision binary64
 (fabs
  (*
   (pow (sqrt (PI)) -1.0)
   (fma
    (pow x 5.0)
    0.2
    (fma
     x
     (fma 0.6666666666666666 (* x x) 2.0)
     (* (pow x 7.0) 0.047619047619047616))))))
\begin{array}{l}

\\
\left|{\left(\sqrt{\mathsf{PI}\left(\right)}\right)}^{-1} \cdot \mathsf{fma}\left({x}^{5}, 0.2, \mathsf{fma}\left(x, \mathsf{fma}\left(0.6666666666666666, x \cdot x, 2\right), {x}^{7} \cdot 0.047619047619047616\right)\right)\right|
\end{array}
Derivation
  1. Initial program 99.9%

    \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  2. Add Preprocessing
  3. Applied rewrites99.9%

    \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \color{blue}{\mathsf{fma}\left({x}^{5}, 0.2, \mathsf{fma}\left(x, \mathsf{fma}\left(0.6666666666666666, x \cdot x, 2\right), {x}^{7} \cdot 0.047619047619047616\right)\right)}\right| \]
  4. Final simplification99.9%

    \[\leadsto \left|{\left(\sqrt{\mathsf{PI}\left(\right)}\right)}^{-1} \cdot \mathsf{fma}\left({x}^{5}, 0.2, \mathsf{fma}\left(x, \mathsf{fma}\left(0.6666666666666666, x \cdot x, 2\right), {x}^{7} \cdot 0.047619047619047616\right)\right)\right| \]
  5. Add Preprocessing

Alternative 5: 88.9% accurate, 1.3× speedup?

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq 2.2:\\ \;\;\;\;\left|\left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right|\\ \mathbf{else}:\\ \;\;\;\;{x}^{7} \cdot \frac{0.047619047619047616}{\sqrt{\mathsf{PI}\left(\right)}}\\ \end{array} \end{array} \]
(FPCore (x)
 :precision binary64
 (if (<= x 2.2)
   (fabs (* (* (fma (* x x) 0.6666666666666666 2.0) x) (sqrt (pow (PI) -1.0))))
   (* (pow x 7.0) (/ 0.047619047619047616 (sqrt (PI))))))
\begin{array}{l}

\\
\begin{array}{l}
\mathbf{if}\;x \leq 2.2:\\
\;\;\;\;\left|\left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right|\\

\mathbf{else}:\\
\;\;\;\;{x}^{7} \cdot \frac{0.047619047619047616}{\sqrt{\mathsf{PI}\left(\right)}}\\


\end{array}
\end{array}
Derivation
  1. Split input into 2 regimes
  2. if x < 2.2000000000000002

    1. Initial program 99.9%

      \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    2. Add Preprocessing
    3. Applied rewrites99.9%

      \[\leadsto \left|\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(0.2, {x}^{5}, \mathsf{fma}\left(0.6666666666666666, {x}^{3}, 2 \cdot x\right)\right), {\mathsf{PI}\left(\right)}^{-0.5}, \left({\mathsf{PI}\left(\right)}^{-0.5} \cdot 0.047619047619047616\right) \cdot {x}^{7}\right)}\right| \]
    4. Taylor expanded in x around 0

      \[\leadsto \left|\color{blue}{x \cdot \left(\frac{2}{3} \cdot \left({x}^{2} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right) + 2 \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right)}\right| \]
    5. Applied rewrites91.7%

      \[\leadsto \left|\color{blue}{\left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]

    if 2.2000000000000002 < x

    1. Initial program 99.9%

      \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    2. Add Preprocessing
    3. Applied rewrites99.9%

      \[\leadsto \left|\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(0.2, {x}^{5}, \mathsf{fma}\left(0.6666666666666666, {x}^{3}, 2 \cdot x\right)\right), {\mathsf{PI}\left(\right)}^{-0.5}, \left({\mathsf{PI}\left(\right)}^{-0.5} \cdot 0.047619047619047616\right) \cdot {x}^{7}\right)}\right| \]
    4. Taylor expanded in x around 0

      \[\leadsto \left|\color{blue}{2 \cdot \left(x \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right)}\right| \]
    5. Step-by-step derivation
      1. associate-*r*N/A

        \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      2. lower-*.f64N/A

        \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      3. lower-*.f64N/A

        \[\leadsto \left|\color{blue}{\left(2 \cdot x\right)} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right| \]
      4. rem-square-sqrtN/A

        \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
      5. lower-sqrt.f64N/A

        \[\leadsto \left|\left(2 \cdot x\right) \cdot \color{blue}{\sqrt{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
      6. rem-square-sqrtN/A

        \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      7. lower-/.f64N/A

        \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      8. lower-PI.f6474.2

        \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\color{blue}{\mathsf{PI}\left(\right)}}}\right| \]
    6. Applied rewrites74.2%

      \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
    7. Taylor expanded in x around inf

      \[\leadsto \left|\color{blue}{\frac{1}{21} \cdot \left({x}^{7} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right)}\right| \]
    8. Step-by-step derivation
      1. associate-*r*N/A

        \[\leadsto \left|\color{blue}{\left(\frac{1}{21} \cdot {x}^{7}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      2. lower-*.f64N/A

        \[\leadsto \left|\color{blue}{\left(\frac{1}{21} \cdot {x}^{7}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      3. *-commutativeN/A

        \[\leadsto \left|\color{blue}{\left({x}^{7} \cdot \frac{1}{21}\right)} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right| \]
      4. lower-*.f64N/A

        \[\leadsto \left|\color{blue}{\left({x}^{7} \cdot \frac{1}{21}\right)} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right| \]
      5. lower-pow.f64N/A

        \[\leadsto \left|\left(\color{blue}{{x}^{7}} \cdot \frac{1}{21}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right| \]
      6. rem-square-sqrtN/A

        \[\leadsto \left|\left({x}^{7} \cdot \frac{1}{21}\right) \cdot \sqrt{\color{blue}{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
      7. lower-sqrt.f64N/A

        \[\leadsto \left|\left({x}^{7} \cdot \frac{1}{21}\right) \cdot \color{blue}{\sqrt{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
      8. rem-square-sqrtN/A

        \[\leadsto \left|\left({x}^{7} \cdot \frac{1}{21}\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      9. lower-/.f64N/A

        \[\leadsto \left|\left({x}^{7} \cdot \frac{1}{21}\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      10. lower-PI.f6430.7

        \[\leadsto \left|\left({x}^{7} \cdot 0.047619047619047616\right) \cdot \sqrt{\frac{1}{\color{blue}{\mathsf{PI}\left(\right)}}}\right| \]
    9. Applied rewrites30.7%

      \[\leadsto \left|\color{blue}{\left({x}^{7} \cdot 0.047619047619047616\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
    10. Step-by-step derivation
      1. lift-fabs.f64N/A

        \[\leadsto \color{blue}{\left|\left({x}^{7} \cdot \frac{1}{21}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right|} \]
      2. rem-sqrt-square-revN/A

        \[\leadsto \color{blue}{\sqrt{\left(\left({x}^{7} \cdot \frac{1}{21}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right) \cdot \left(\left({x}^{7} \cdot \frac{1}{21}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right)}} \]
      3. sqrt-prodN/A

        \[\leadsto \color{blue}{\sqrt{\left({x}^{7} \cdot \frac{1}{21}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}} \cdot \sqrt{\left({x}^{7} \cdot \frac{1}{21}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}} \]
      4. rem-square-sqrt4.0

        \[\leadsto \color{blue}{\left({x}^{7} \cdot 0.047619047619047616\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}} \]
    11. Applied rewrites4.0%

      \[\leadsto \color{blue}{{x}^{7} \cdot \frac{0.047619047619047616}{\sqrt{\mathsf{PI}\left(\right)}}} \]
  3. Recombined 2 regimes into one program.
  4. Final simplification91.7%

    \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq 2.2:\\ \;\;\;\;\left|\left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right|\\ \mathbf{else}:\\ \;\;\;\;{x}^{7} \cdot \frac{0.047619047619047616}{\sqrt{\mathsf{PI}\left(\right)}}\\ \end{array} \]
  5. Add Preprocessing

Alternative 6: 88.9% accurate, 1.4× speedup?

\[\begin{array}{l} \\ \left|\left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right| \end{array} \]
(FPCore (x)
 :precision binary64
 (fabs (* (* (fma (* x x) 0.6666666666666666 2.0) x) (sqrt (pow (PI) -1.0)))))
\begin{array}{l}

\\
\left|\left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right|
\end{array}
Derivation
  1. Initial program 99.9%

    \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
  2. Add Preprocessing
  3. Applied rewrites99.9%

    \[\leadsto \left|\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(0.2, {x}^{5}, \mathsf{fma}\left(0.6666666666666666, {x}^{3}, 2 \cdot x\right)\right), {\mathsf{PI}\left(\right)}^{-0.5}, \left({\mathsf{PI}\left(\right)}^{-0.5} \cdot 0.047619047619047616\right) \cdot {x}^{7}\right)}\right| \]
  4. Taylor expanded in x around 0

    \[\leadsto \left|\color{blue}{x \cdot \left(\frac{2}{3} \cdot \left({x}^{2} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right) + 2 \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right)}\right| \]
  5. Applied rewrites91.7%

    \[\leadsto \left|\color{blue}{\left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
  6. Final simplification91.7%

    \[\leadsto \left|\left(\mathsf{fma}\left(x \cdot x, 0.6666666666666666, 2\right) \cdot x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right| \]
  7. Add Preprocessing

Alternative 7: 66.5% accurate, 1.4× speedup?

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq 1.75:\\ \;\;\;\;\left|\left(x + x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right|\\ \mathbf{else}:\\ \;\;\;\;\left|x \cdot \left(\left(x \cdot x\right) \cdot \frac{0.6666666666666666}{\sqrt{\mathsf{PI}\left(\right)}}\right)\right|\\ \end{array} \end{array} \]
(FPCore (x)
 :precision binary64
 (if (<= x 1.75)
   (fabs (* (+ x x) (sqrt (pow (PI) -1.0))))
   (fabs (* x (* (* x x) (/ 0.6666666666666666 (sqrt (PI))))))))
\begin{array}{l}

\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.75:\\
\;\;\;\;\left|\left(x + x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right|\\

\mathbf{else}:\\
\;\;\;\;\left|x \cdot \left(\left(x \cdot x\right) \cdot \frac{0.6666666666666666}{\sqrt{\mathsf{PI}\left(\right)}}\right)\right|\\


\end{array}
\end{array}
Derivation
  1. Split input into 2 regimes
  2. if x < 1.75

    1. Initial program 99.9%

      \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
    2. Add Preprocessing
    3. Applied rewrites99.9%

      \[\leadsto \left|\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(0.2, {x}^{5}, \mathsf{fma}\left(0.6666666666666666, {x}^{3}, 2 \cdot x\right)\right), {\mathsf{PI}\left(\right)}^{-0.5}, \left({\mathsf{PI}\left(\right)}^{-0.5} \cdot 0.047619047619047616\right) \cdot {x}^{7}\right)}\right| \]
    4. Taylor expanded in x around 0

      \[\leadsto \left|\color{blue}{2 \cdot \left(x \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right)}\right| \]
    5. Step-by-step derivation
      1. associate-*r*N/A

        \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      2. lower-*.f64N/A

        \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      3. lower-*.f64N/A

        \[\leadsto \left|\color{blue}{\left(2 \cdot x\right)} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right| \]
      4. rem-square-sqrtN/A

        \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
      5. lower-sqrt.f64N/A

        \[\leadsto \left|\left(2 \cdot x\right) \cdot \color{blue}{\sqrt{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
      6. rem-square-sqrtN/A

        \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      7. lower-/.f64N/A

        \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      8. lower-PI.f6474.2

        \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\color{blue}{\mathsf{PI}\left(\right)}}}\right| \]
    6. Applied rewrites74.2%

      \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
    7. Step-by-step derivation
      1. Applied rewrites74.2%

        \[\leadsto \left|\left(x + x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]

      if 1.75 < x

      1. Initial program 99.9%

        \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
      2. Add Preprocessing
      3. Step-by-step derivation
        1. lift-*.f64N/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        2. *-commutativeN/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \frac{2}{3}}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        3. lift-*.f64N/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)} \cdot \frac{2}{3}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        4. associate-*l*N/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        5. lower-*.f64N/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        6. lift-fabs.f64N/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(\color{blue}{\left|x\right|} \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        7. rem-sqrt-square-revN/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(\color{blue}{\sqrt{x \cdot x}} \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        8. sqrt-prodN/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(\color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)} \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        9. rem-square-sqrtN/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(\color{blue}{x} \cdot \left|x\right|\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        10. lift-fabs.f64N/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot \color{blue}{\left|x\right|}\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        11. rem-sqrt-square-revN/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot \color{blue}{\sqrt{x \cdot x}}\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        12. sqrt-prodN/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot \color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)}\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        13. rem-square-sqrtN/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot \color{blue}{x}\right) \cdot \left(\left|x\right| \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        14. lower-*.f6499.9

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \color{blue}{\left(\left|x\right| \cdot \frac{2}{3}\right)}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        15. lift-fabs.f64N/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(\color{blue}{\left|x\right|} \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        16. rem-sqrt-square-revN/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(\color{blue}{\sqrt{x \cdot x}} \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        17. sqrt-prodN/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(\color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)} \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        18. rem-square-sqrt81.2

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(\color{blue}{x} \cdot \frac{2}{3}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        19. lift-/.f64N/A

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \color{blue}{\frac{2}{3}}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        20. metadata-eval81.2

          \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \left(x \cdot x\right) \cdot \left(x \cdot \color{blue}{0.6666666666666666}\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
      4. Applied rewrites81.2%

        \[\leadsto \left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \color{blue}{\left(x \cdot x\right) \cdot \left(x \cdot 0.6666666666666666\right)}\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
      5. Taylor expanded in x around inf

        \[\leadsto \left|\color{blue}{\frac{2}{3} \cdot \left({x}^{3} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right)}\right| \]
      6. Step-by-step derivation
        1. associate-*r*N/A

          \[\leadsto \left|\color{blue}{\left(\frac{2}{3} \cdot {x}^{3}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
        2. lower-*.f64N/A

          \[\leadsto \left|\color{blue}{\left(\frac{2}{3} \cdot {x}^{3}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
        3. *-commutativeN/A

          \[\leadsto \left|\color{blue}{\left({x}^{3} \cdot \frac{2}{3}\right)} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right| \]
        4. lower-*.f64N/A

          \[\leadsto \left|\color{blue}{\left({x}^{3} \cdot \frac{2}{3}\right)} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right| \]
        5. lower-pow.f64N/A

          \[\leadsto \left|\left(\color{blue}{{x}^{3}} \cdot \frac{2}{3}\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right| \]
        6. rem-square-sqrtN/A

          \[\leadsto \left|\left({x}^{3} \cdot \frac{2}{3}\right) \cdot \sqrt{\color{blue}{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
        7. lower-sqrt.f64N/A

          \[\leadsto \left|\left({x}^{3} \cdot \frac{2}{3}\right) \cdot \color{blue}{\sqrt{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
        8. rem-square-sqrtN/A

          \[\leadsto \left|\left({x}^{3} \cdot \frac{2}{3}\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
        9. lower-/.f64N/A

          \[\leadsto \left|\left({x}^{3} \cdot \frac{2}{3}\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
        10. lower-PI.f6423.5

          \[\leadsto \left|\left({x}^{3} \cdot 0.6666666666666666\right) \cdot \sqrt{\frac{1}{\color{blue}{\mathsf{PI}\left(\right)}}}\right| \]
      7. Applied rewrites23.5%

        \[\leadsto \left|\color{blue}{\left({x}^{3} \cdot 0.6666666666666666\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      8. Step-by-step derivation
        1. Applied rewrites23.5%

          \[\leadsto \left|x \cdot \color{blue}{\left(\left(x \cdot x\right) \cdot \frac{0.6666666666666666}{\sqrt{\mathsf{PI}\left(\right)}}\right)}\right| \]
      9. Recombined 2 regimes into one program.
      10. Final simplification74.2%

        \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq 1.75:\\ \;\;\;\;\left|\left(x + x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right|\\ \mathbf{else}:\\ \;\;\;\;\left|x \cdot \left(\left(x \cdot x\right) \cdot \frac{0.6666666666666666}{\sqrt{\mathsf{PI}\left(\right)}}\right)\right|\\ \end{array} \]
      11. Add Preprocessing

      Alternative 8: 66.5% accurate, 1.5× speedup?

      \[\begin{array}{l} \\ \left|\left(x + x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right| \end{array} \]
      (FPCore (x) :precision binary64 (fabs (* (+ x x) (sqrt (pow (PI) -1.0)))))
      \begin{array}{l}
      
      \\
      \left|\left(x + x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right|
      \end{array}
      
      Derivation
      1. Initial program 99.9%

        \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
      2. Add Preprocessing
      3. Applied rewrites99.9%

        \[\leadsto \left|\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(0.2, {x}^{5}, \mathsf{fma}\left(0.6666666666666666, {x}^{3}, 2 \cdot x\right)\right), {\mathsf{PI}\left(\right)}^{-0.5}, \left({\mathsf{PI}\left(\right)}^{-0.5} \cdot 0.047619047619047616\right) \cdot {x}^{7}\right)}\right| \]
      4. Taylor expanded in x around 0

        \[\leadsto \left|\color{blue}{2 \cdot \left(x \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right)}\right| \]
      5. Step-by-step derivation
        1. associate-*r*N/A

          \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
        2. lower-*.f64N/A

          \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
        3. lower-*.f64N/A

          \[\leadsto \left|\color{blue}{\left(2 \cdot x\right)} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right| \]
        4. rem-square-sqrtN/A

          \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
        5. lower-sqrt.f64N/A

          \[\leadsto \left|\left(2 \cdot x\right) \cdot \color{blue}{\sqrt{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
        6. rem-square-sqrtN/A

          \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
        7. lower-/.f64N/A

          \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
        8. lower-PI.f6474.2

          \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\color{blue}{\mathsf{PI}\left(\right)}}}\right| \]
      6. Applied rewrites74.2%

        \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
      7. Step-by-step derivation
        1. Applied rewrites74.2%

          \[\leadsto \left|\left(x + x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
        2. Final simplification74.2%

          \[\leadsto \left|\left(x + x\right) \cdot \sqrt{{\mathsf{PI}\left(\right)}^{-1}}\right| \]
        3. Add Preprocessing

        Alternative 9: 4.1% accurate, 183.0× speedup?

        \[\begin{array}{l} \\ 0 \end{array} \]
        (FPCore (x) :precision binary64 0.0)
        double code(double x) {
        	return 0.0;
        }
        
        module fmin_fmax_functions
            implicit none
            private
            public fmax
            public fmin
        
            interface fmax
                module procedure fmax88
                module procedure fmax44
                module procedure fmax84
                module procedure fmax48
            end interface
            interface fmin
                module procedure fmin88
                module procedure fmin44
                module procedure fmin84
                module procedure fmin48
            end interface
        contains
            real(8) function fmax88(x, y) result (res)
                real(8), intent (in) :: x
                real(8), intent (in) :: y
                res = merge(y, merge(x, max(x, y), y /= y), x /= x)
            end function
            real(4) function fmax44(x, y) result (res)
                real(4), intent (in) :: x
                real(4), intent (in) :: y
                res = merge(y, merge(x, max(x, y), y /= y), x /= x)
            end function
            real(8) function fmax84(x, y) result(res)
                real(8), intent (in) :: x
                real(4), intent (in) :: y
                res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
            end function
            real(8) function fmax48(x, y) result(res)
                real(4), intent (in) :: x
                real(8), intent (in) :: y
                res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
            end function
            real(8) function fmin88(x, y) result (res)
                real(8), intent (in) :: x
                real(8), intent (in) :: y
                res = merge(y, merge(x, min(x, y), y /= y), x /= x)
            end function
            real(4) function fmin44(x, y) result (res)
                real(4), intent (in) :: x
                real(4), intent (in) :: y
                res = merge(y, merge(x, min(x, y), y /= y), x /= x)
            end function
            real(8) function fmin84(x, y) result(res)
                real(8), intent (in) :: x
                real(4), intent (in) :: y
                res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
            end function
            real(8) function fmin48(x, y) result(res)
                real(4), intent (in) :: x
                real(8), intent (in) :: y
                res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
            end function
        end module
        
        real(8) function code(x)
        use fmin_fmax_functions
            real(8), intent (in) :: x
            code = 0.0d0
        end function
        
        public static double code(double x) {
        	return 0.0;
        }
        
        def code(x):
        	return 0.0
        
        function code(x)
        	return 0.0
        end
        
        function tmp = code(x)
        	tmp = 0.0;
        end
        
        code[x_] := 0.0
        
        \begin{array}{l}
        
        \\
        0
        \end{array}
        
        Derivation
        1. Initial program 99.9%

          \[\left|\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right| \]
        2. Add Preprocessing
        3. Applied rewrites99.9%

          \[\leadsto \left|\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(0.2, {x}^{5}, \mathsf{fma}\left(0.6666666666666666, {x}^{3}, 2 \cdot x\right)\right), {\mathsf{PI}\left(\right)}^{-0.5}, \left({\mathsf{PI}\left(\right)}^{-0.5} \cdot 0.047619047619047616\right) \cdot {x}^{7}\right)}\right| \]
        4. Taylor expanded in x around 0

          \[\leadsto \left|\color{blue}{2 \cdot \left(x \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right)}\right| \]
        5. Step-by-step derivation
          1. associate-*r*N/A

            \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
          2. lower-*.f64N/A

            \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
          3. lower-*.f64N/A

            \[\leadsto \left|\color{blue}{\left(2 \cdot x\right)} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}\right| \]
          4. rem-square-sqrtN/A

            \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
          5. lower-sqrt.f64N/A

            \[\leadsto \left|\left(2 \cdot x\right) \cdot \color{blue}{\sqrt{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}}\right| \]
          6. rem-square-sqrtN/A

            \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
          7. lower-/.f64N/A

            \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
          8. lower-PI.f6474.2

            \[\leadsto \left|\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\color{blue}{\mathsf{PI}\left(\right)}}}\right| \]
        6. Applied rewrites74.2%

          \[\leadsto \left|\color{blue}{\left(2 \cdot x\right) \cdot \sqrt{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
        7. Step-by-step derivation
          1. Applied rewrites74.2%

            \[\leadsto \left|\left(x + x\right) \cdot \sqrt{\color{blue}{\frac{1}{\mathsf{PI}\left(\right)}}}\right| \]
          2. Applied rewrites4.3%

            \[\leadsto \color{blue}{\frac{0}{\sqrt{\mathsf{PI}\left(\right)}}} \]
          3. Final simplification4.3%

            \[\leadsto 0 \]
          4. Add Preprocessing

          Reproduce

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