Ian Simplification

Percentage Accurate: 6.5% → 8.0%
Time: 41.5s
Alternatives: 4
Speedup: 1.1×

Specification

?
\[\begin{array}{l} \\ \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right) \end{array} \]
(FPCore (x)
 :precision binary64
 (- (/ (PI) 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))
\begin{array}{l}

\\
\frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)
\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 4 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: 6.5% accurate, 1.0× speedup?

\[\begin{array}{l} \\ \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right) \end{array} \]
(FPCore (x)
 :precision binary64
 (- (/ (PI) 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))
\begin{array}{l}

\\
\frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)
\end{array}

Alternative 1: 8.0% accurate, 0.4× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \mathsf{PI}\left(\right) \cdot 0.5\\ t_1 := \sqrt{\mathsf{fma}\left(-0.5, x, 0.5\right)}\\ \frac{\mathsf{fma}\left({\left(t\_0 - \cos^{-1} t\_1\right)}^{2}, -4, \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot 0.25\right)}{\mathsf{fma}\left(\sin^{-1} t\_1, 2, t\_0\right)} \end{array} \end{array} \]
(FPCore (x)
 :precision binary64
 (let* ((t_0 (* (PI) 0.5)) (t_1 (sqrt (fma -0.5 x 0.5))))
   (/
    (fma (pow (- t_0 (acos t_1)) 2.0) -4.0 (* (* (PI) (PI)) 0.25))
    (fma (asin t_1) 2.0 t_0))))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \mathsf{PI}\left(\right) \cdot 0.5\\
t_1 := \sqrt{\mathsf{fma}\left(-0.5, x, 0.5\right)}\\
\frac{\mathsf{fma}\left({\left(t\_0 - \cos^{-1} t\_1\right)}^{2}, -4, \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot 0.25\right)}{\mathsf{fma}\left(\sin^{-1} t\_1, 2, t\_0\right)}
\end{array}
\end{array}
Derivation
  1. Initial program 6.9%

    \[\frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right) \]
  2. Add Preprocessing
  3. Applied rewrites6.9%

    \[\leadsto \color{blue}{\left(-\mathsf{fma}\left({\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)}^{2}, -4, 0.25 \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right), 2, 0.5 \cdot \mathsf{PI}\left(\right)\right)}} \]
  4. Step-by-step derivation
    1. lift-asin.f64N/A

      \[\leadsto \left(-\mathsf{fma}\left({\color{blue}{\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right)}}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    2. asin-acosN/A

      \[\leadsto \left(-\mathsf{fma}\left({\color{blue}{\left(\frac{\mathsf{PI}\left(\right)}{2} - \cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right)\right)}}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    3. lift-PI.f64N/A

      \[\leadsto \left(-\mathsf{fma}\left({\left(\frac{\color{blue}{\mathsf{PI}\left(\right)}}{2} - \cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right)\right)}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    4. lift-/.f64N/A

      \[\leadsto \left(-\mathsf{fma}\left({\left(\color{blue}{\frac{\mathsf{PI}\left(\right)}{2}} - \cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right)\right)}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    5. lift-acos.f64N/A

      \[\leadsto \left(-\mathsf{fma}\left({\left(\frac{\mathsf{PI}\left(\right)}{2} - \color{blue}{\cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right)}\right)}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    6. unsub-negN/A

      \[\leadsto \left(-\mathsf{fma}\left({\color{blue}{\left(\frac{\mathsf{PI}\left(\right)}{2} + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right)\right)\right)\right)}}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    7. lift-/.f64N/A

      \[\leadsto \left(-\mathsf{fma}\left({\left(\color{blue}{\frac{\mathsf{PI}\left(\right)}{2}} + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right)\right)\right)\right)}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    8. div-invN/A

      \[\leadsto \left(-\mathsf{fma}\left({\left(\color{blue}{\mathsf{PI}\left(\right) \cdot \frac{1}{2}} + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right)\right)\right)\right)}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    9. metadata-evalN/A

      \[\leadsto \left(-\mathsf{fma}\left({\left(\mathsf{PI}\left(\right) \cdot \color{blue}{\frac{1}{2}} + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right)\right)\right)\right)}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    10. lift-neg.f64N/A

      \[\leadsto \left(-\mathsf{fma}\left({\left(\mathsf{PI}\left(\right) \cdot \frac{1}{2} + \color{blue}{\left(-\cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right)\right)}\right)}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    11. lift-fma.f648.5

      \[\leadsto \left(-\mathsf{fma}\left({\color{blue}{\left(\mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, -\cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)\right)\right)}}^{2}, -4, 0.25 \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right), 2, 0.5 \cdot \mathsf{PI}\left(\right)\right)} \]
    12. lift-fma.f64N/A

      \[\leadsto \left(-\mathsf{fma}\left({\left(\mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{x \cdot \frac{-1}{2} + \frac{1}{2}}}\right)\right)\right)}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    13. *-commutativeN/A

      \[\leadsto \left(-\mathsf{fma}\left({\left(\mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\frac{-1}{2} \cdot x} + \frac{1}{2}}\right)\right)\right)}^{2}, -4, \frac{1}{4} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, \frac{-1}{2}, \frac{1}{2}\right)}\right), 2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)} \]
    14. lower-fma.f648.5

      \[\leadsto \left(-\mathsf{fma}\left({\left(\mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, -\cos^{-1} \left(\sqrt{\color{blue}{\mathsf{fma}\left(-0.5, x, 0.5\right)}}\right)\right)\right)}^{2}, -4, 0.25 \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right), 2, 0.5 \cdot \mathsf{PI}\left(\right)\right)} \]
  5. Applied rewrites8.5%

    \[\leadsto \left(-\mathsf{fma}\left({\color{blue}{\left(\mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, -\cos^{-1} \left(\sqrt{\mathsf{fma}\left(-0.5, x, 0.5\right)}\right)\right)\right)}}^{2}, -4, 0.25 \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\right) \cdot \frac{1}{-\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right), 2, 0.5 \cdot \mathsf{PI}\left(\right)\right)} \]
  6. Taylor expanded in x around 0

    \[\leadsto \color{blue}{\frac{-4 \cdot {\left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) - \cos^{-1} \left(\sqrt{\frac{1}{2} + \frac{-1}{2} \cdot x}\right)\right)}^{2} + \frac{1}{4} \cdot {\mathsf{PI}\left(\right)}^{2}}{\frac{1}{2} \cdot \mathsf{PI}\left(\right) + 2 \cdot \sin^{-1} \left(\sqrt{\frac{1}{2} + \frac{-1}{2} \cdot x}\right)}} \]
  7. Applied rewrites8.5%

    \[\leadsto \color{blue}{\frac{\mathsf{fma}\left({\left(\mathsf{PI}\left(\right) \cdot 0.5 - \cos^{-1} \left(\sqrt{\mathsf{fma}\left(-0.5, x, 0.5\right)}\right)\right)}^{2}, -4, \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot 0.25\right)}{\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(-0.5, x, 0.5\right)}\right), 2, \mathsf{PI}\left(\right) \cdot 0.5\right)}} \]
  8. Add Preprocessing

Alternative 2: 8.0% accurate, 0.4× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)\\ \frac{\left(\mathsf{PI}\left(\right) \cdot 0.25\right) \cdot \mathsf{PI}\left(\right) - {t\_0}^{2} \cdot 4}{\mathsf{fma}\left(\mathsf{PI}\left(\right), -0.5, -2 \cdot t\_0\right)} \end{array} \end{array} \]
(FPCore (x)
 :precision binary64
 (let* ((t_0 (acos (sqrt (fma x -0.5 0.5)))))
   (/
    (- (* (* (PI) 0.25) (PI)) (* (pow t_0 2.0) 4.0))
    (fma (PI) -0.5 (* -2.0 t_0)))))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)\\
\frac{\left(\mathsf{PI}\left(\right) \cdot 0.25\right) \cdot \mathsf{PI}\left(\right) - {t\_0}^{2} \cdot 4}{\mathsf{fma}\left(\mathsf{PI}\left(\right), -0.5, -2 \cdot t\_0\right)}
\end{array}
\end{array}
Derivation
  1. Initial program 6.9%

    \[\frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right) \]
  2. Add Preprocessing
  3. Step-by-step derivation
    1. lift-asin.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \color{blue}{\sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)} \]
    2. asin-acosN/A

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

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

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \left(\color{blue}{\frac{\mathsf{PI}\left(\right)}{2}} - \cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)\right) \]
    5. sub-negN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \color{blue}{\left(\frac{\mathsf{PI}\left(\right)}{2} + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)\right)\right)\right)} \]
    6. lift-/.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \left(\color{blue}{\frac{\mathsf{PI}\left(\right)}{2}} + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)\right)\right)\right) \]
    7. div-invN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \left(\color{blue}{\mathsf{PI}\left(\right) \cdot \frac{1}{2}} + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)\right)\right)\right) \]
    8. metadata-evalN/A

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

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \color{blue}{\mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, \mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)\right)\right)} \]
    10. lower-neg.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, \color{blue}{-\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)}\right) \]
    11. lower-acos.f648.5

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, -\color{blue}{\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)}\right) \]
    12. lift-/.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\frac{1 - x}{2}}}\right)\right) \]
    13. lift--.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\frac{\color{blue}{1 - x}}{2}}\right)\right) \]
    14. div-subN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\frac{1}{2} - \frac{x}{2}}}\right)\right) \]
    15. metadata-evalN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\frac{1}{2}} - \frac{x}{2}}\right)\right) \]
    16. sub-negN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\frac{1}{2} + \left(\mathsf{neg}\left(\frac{x}{2}\right)\right)}}\right)\right) \]
    17. +-commutativeN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\left(\mathsf{neg}\left(\frac{x}{2}\right)\right) + \frac{1}{2}}}\right)\right) \]
    18. div-invN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\left(\mathsf{neg}\left(\color{blue}{x \cdot \frac{1}{2}}\right)\right) + \frac{1}{2}}\right)\right) \]
    19. metadata-evalN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\left(\mathsf{neg}\left(x \cdot \color{blue}{\frac{1}{2}}\right)\right) + \frac{1}{2}}\right)\right) \]
    20. distribute-rgt-neg-inN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{x \cdot \left(\mathsf{neg}\left(\frac{1}{2}\right)\right)} + \frac{1}{2}}\right)\right) \]
    21. metadata-evalN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{x \cdot \color{blue}{\frac{-1}{2}} + \frac{1}{2}}\right)\right) \]
    22. metadata-evalN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{x \cdot \color{blue}{\frac{1}{-2}} + \frac{1}{2}}\right)\right) \]
    23. metadata-evalN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{x \cdot \frac{1}{\color{blue}{\mathsf{neg}\left(2\right)}} + \frac{1}{2}}\right)\right) \]
    24. lower-fma.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\mathsf{fma}\left(x, \frac{1}{\mathsf{neg}\left(2\right)}, \frac{1}{2}\right)}}\right)\right) \]
  4. Applied rewrites8.5%

    \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \color{blue}{\mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, -\cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)\right)} \]
  5. Taylor expanded in x around 0

    \[\leadsto \color{blue}{\frac{1}{2} \cdot \mathsf{PI}\left(\right) - 2 \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) - \cos^{-1} \left(\sqrt{\frac{1}{2} + \frac{-1}{2} \cdot x}\right)\right)} \]
  6. Step-by-step derivation
    1. cancel-sign-sub-invN/A

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

      \[\leadsto \frac{1}{2} \cdot \mathsf{PI}\left(\right) + \color{blue}{-2} \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) - \cos^{-1} \left(\sqrt{\frac{1}{2} + \frac{-1}{2} \cdot x}\right)\right) \]
    3. sub-negN/A

      \[\leadsto \frac{1}{2} \cdot \mathsf{PI}\left(\right) + -2 \cdot \color{blue}{\left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} + \frac{-1}{2} \cdot x}\right)\right)\right)\right)} \]
    4. metadata-evalN/A

      \[\leadsto \frac{1}{2} \cdot \mathsf{PI}\left(\right) + -2 \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} + \color{blue}{\left(\mathsf{neg}\left(\frac{1}{2}\right)\right)} \cdot x}\right)\right)\right)\right) \]
    5. cancel-sign-sub-invN/A

      \[\leadsto \frac{1}{2} \cdot \mathsf{PI}\left(\right) + -2 \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\color{blue}{\frac{1}{2} - \frac{1}{2} \cdot x}}\right)\right)\right)\right) \]
    6. distribute-lft-inN/A

      \[\leadsto \frac{1}{2} \cdot \mathsf{PI}\left(\right) + \color{blue}{\left(-2 \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right)\right) + -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right)\right)} \]
    7. associate-+r+N/A

      \[\leadsto \color{blue}{\left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) + -2 \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)\right) + -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right)} \]
    8. distribute-rgt1-inN/A

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

      \[\leadsto \color{blue}{-1} \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right)\right) + -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right) \]
    10. associate-*r*N/A

      \[\leadsto \color{blue}{\left(-1 \cdot \frac{1}{2}\right) \cdot \mathsf{PI}\left(\right)} + -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right) \]
    11. metadata-evalN/A

      \[\leadsto \color{blue}{\frac{-1}{2}} \cdot \mathsf{PI}\left(\right) + -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right) \]
    12. lower-fma.f64N/A

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{2}, \mathsf{PI}\left(\right), -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right)\right)} \]
  7. Applied rewrites8.5%

    \[\leadsto \color{blue}{\mathsf{fma}\left(-0.5, \mathsf{PI}\left(\right), \cos^{-1} \left(\sqrt{\mathsf{fma}\left(-0.5, x, 0.5\right)}\right) \cdot 2\right)} \]
  8. Applied rewrites8.5%

    \[\leadsto \frac{\left(\mathsf{PI}\left(\right) \cdot 0.25\right) \cdot \mathsf{PI}\left(\right) - {\cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)}^{2} \cdot 4}{\color{blue}{\mathsf{fma}\left(\mathsf{PI}\left(\right), -0.5, -2 \cdot \cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)\right)}} \]
  9. Add Preprocessing

Alternative 3: 8.0% accurate, 1.1× speedup?

\[\begin{array}{l} \\ \mathsf{fma}\left(-0.5, \mathsf{PI}\left(\right), \cos^{-1} \left(\sqrt{\mathsf{fma}\left(-0.5, x, 0.5\right)}\right) \cdot 2\right) \end{array} \]
(FPCore (x)
 :precision binary64
 (fma -0.5 (PI) (* (acos (sqrt (fma -0.5 x 0.5))) 2.0)))
\begin{array}{l}

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

    \[\frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right) \]
  2. Add Preprocessing
  3. Step-by-step derivation
    1. lift-asin.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \color{blue}{\sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)} \]
    2. asin-acosN/A

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

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

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \left(\color{blue}{\frac{\mathsf{PI}\left(\right)}{2}} - \cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)\right) \]
    5. sub-negN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \color{blue}{\left(\frac{\mathsf{PI}\left(\right)}{2} + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)\right)\right)\right)} \]
    6. lift-/.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \left(\color{blue}{\frac{\mathsf{PI}\left(\right)}{2}} + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)\right)\right)\right) \]
    7. div-invN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \left(\color{blue}{\mathsf{PI}\left(\right) \cdot \frac{1}{2}} + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)\right)\right)\right) \]
    8. metadata-evalN/A

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

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \color{blue}{\mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, \mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)\right)\right)} \]
    10. lower-neg.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, \color{blue}{-\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)}\right) \]
    11. lower-acos.f648.5

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, -\color{blue}{\cos^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)}\right) \]
    12. lift-/.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\frac{1 - x}{2}}}\right)\right) \]
    13. lift--.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\frac{\color{blue}{1 - x}}{2}}\right)\right) \]
    14. div-subN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\frac{1}{2} - \frac{x}{2}}}\right)\right) \]
    15. metadata-evalN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\frac{1}{2}} - \frac{x}{2}}\right)\right) \]
    16. sub-negN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\frac{1}{2} + \left(\mathsf{neg}\left(\frac{x}{2}\right)\right)}}\right)\right) \]
    17. +-commutativeN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\left(\mathsf{neg}\left(\frac{x}{2}\right)\right) + \frac{1}{2}}}\right)\right) \]
    18. div-invN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\left(\mathsf{neg}\left(\color{blue}{x \cdot \frac{1}{2}}\right)\right) + \frac{1}{2}}\right)\right) \]
    19. metadata-evalN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\left(\mathsf{neg}\left(x \cdot \color{blue}{\frac{1}{2}}\right)\right) + \frac{1}{2}}\right)\right) \]
    20. distribute-rgt-neg-inN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{x \cdot \left(\mathsf{neg}\left(\frac{1}{2}\right)\right)} + \frac{1}{2}}\right)\right) \]
    21. metadata-evalN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{x \cdot \color{blue}{\frac{-1}{2}} + \frac{1}{2}}\right)\right) \]
    22. metadata-evalN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{x \cdot \color{blue}{\frac{1}{-2}} + \frac{1}{2}}\right)\right) \]
    23. metadata-evalN/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{x \cdot \frac{1}{\color{blue}{\mathsf{neg}\left(2\right)}} + \frac{1}{2}}\right)\right) \]
    24. lower-fma.f64N/A

      \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \mathsf{fma}\left(\mathsf{PI}\left(\right), \frac{1}{2}, -\cos^{-1} \left(\sqrt{\color{blue}{\mathsf{fma}\left(x, \frac{1}{\mathsf{neg}\left(2\right)}, \frac{1}{2}\right)}}\right)\right) \]
  4. Applied rewrites8.5%

    \[\leadsto \frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \color{blue}{\mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, -\cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)\right)} \]
  5. Taylor expanded in x around 0

    \[\leadsto \color{blue}{\frac{1}{2} \cdot \mathsf{PI}\left(\right) - 2 \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) - \cos^{-1} \left(\sqrt{\frac{1}{2} + \frac{-1}{2} \cdot x}\right)\right)} \]
  6. Step-by-step derivation
    1. cancel-sign-sub-invN/A

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

      \[\leadsto \frac{1}{2} \cdot \mathsf{PI}\left(\right) + \color{blue}{-2} \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) - \cos^{-1} \left(\sqrt{\frac{1}{2} + \frac{-1}{2} \cdot x}\right)\right) \]
    3. sub-negN/A

      \[\leadsto \frac{1}{2} \cdot \mathsf{PI}\left(\right) + -2 \cdot \color{blue}{\left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} + \frac{-1}{2} \cdot x}\right)\right)\right)\right)} \]
    4. metadata-evalN/A

      \[\leadsto \frac{1}{2} \cdot \mathsf{PI}\left(\right) + -2 \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} + \color{blue}{\left(\mathsf{neg}\left(\frac{1}{2}\right)\right)} \cdot x}\right)\right)\right)\right) \]
    5. cancel-sign-sub-invN/A

      \[\leadsto \frac{1}{2} \cdot \mathsf{PI}\left(\right) + -2 \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) + \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\color{blue}{\frac{1}{2} - \frac{1}{2} \cdot x}}\right)\right)\right)\right) \]
    6. distribute-lft-inN/A

      \[\leadsto \frac{1}{2} \cdot \mathsf{PI}\left(\right) + \color{blue}{\left(-2 \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right)\right) + -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right)\right)} \]
    7. associate-+r+N/A

      \[\leadsto \color{blue}{\left(\frac{1}{2} \cdot \mathsf{PI}\left(\right) + -2 \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right)\right)\right) + -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right)} \]
    8. distribute-rgt1-inN/A

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

      \[\leadsto \color{blue}{-1} \cdot \left(\frac{1}{2} \cdot \mathsf{PI}\left(\right)\right) + -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right) \]
    10. associate-*r*N/A

      \[\leadsto \color{blue}{\left(-1 \cdot \frac{1}{2}\right) \cdot \mathsf{PI}\left(\right)} + -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right) \]
    11. metadata-evalN/A

      \[\leadsto \color{blue}{\frac{-1}{2}} \cdot \mathsf{PI}\left(\right) + -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right) \]
    12. lower-fma.f64N/A

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{2}, \mathsf{PI}\left(\right), -2 \cdot \left(\mathsf{neg}\left(\cos^{-1} \left(\sqrt{\frac{1}{2} - \frac{1}{2} \cdot x}\right)\right)\right)\right)} \]
  7. Applied rewrites8.5%

    \[\leadsto \color{blue}{\mathsf{fma}\left(-0.5, \mathsf{PI}\left(\right), \cos^{-1} \left(\sqrt{\mathsf{fma}\left(-0.5, x, 0.5\right)}\right) \cdot 2\right)} \]
  8. Add Preprocessing

Alternative 4: 4.1% accurate, 1.2× speedup?

\[\begin{array}{l} \\ \mathsf{fma}\left(\sin^{-1} \left(\sqrt{0.5}\right), -2, 0.5 \cdot \mathsf{PI}\left(\right)\right) \end{array} \]
(FPCore (x) :precision binary64 (fma (asin (sqrt 0.5)) -2.0 (* 0.5 (PI))))
\begin{array}{l}

\\
\mathsf{fma}\left(\sin^{-1} \left(\sqrt{0.5}\right), -2, 0.5 \cdot \mathsf{PI}\left(\right)\right)
\end{array}
Derivation
  1. Initial program 6.9%

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

      \[\leadsto \color{blue}{\frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)} \]
    2. sub-negN/A

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

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

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

      \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right) \cdot 2}\right)\right) + \frac{\mathsf{PI}\left(\right)}{2} \]
    6. distribute-rgt-neg-inN/A

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

      \[\leadsto \color{blue}{\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right), \mathsf{neg}\left(2\right), \frac{\mathsf{PI}\left(\right)}{2}\right)} \]
  4. Applied rewrites6.9%

    \[\leadsto \color{blue}{\mathsf{fma}\left(\sin^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right), -2, 0.5 \cdot \mathsf{PI}\left(\right)\right)} \]
  5. Taylor expanded in x around 0

    \[\leadsto \mathsf{fma}\left(\sin^{-1} \left(\sqrt{\color{blue}{\frac{1}{2}}}\right), -2, \frac{1}{2} \cdot \mathsf{PI}\left(\right)\right) \]
  6. Step-by-step derivation
    1. Applied rewrites4.2%

      \[\leadsto \mathsf{fma}\left(\sin^{-1} \left(\sqrt{\color{blue}{0.5}}\right), -2, 0.5 \cdot \mathsf{PI}\left(\right)\right) \]
    2. Add Preprocessing

    Developer Target 1: 100.0% accurate, 1.4× speedup?

    \[\begin{array}{l} \\ \sin^{-1} x \end{array} \]
    (FPCore (x) :precision binary64 (asin x))
    double code(double x) {
    	return asin(x);
    }
    
    real(8) function code(x)
        real(8), intent (in) :: x
        code = asin(x)
    end function
    
    public static double code(double x) {
    	return Math.asin(x);
    }
    
    def code(x):
    	return math.asin(x)
    
    function code(x)
    	return asin(x)
    end
    
    function tmp = code(x)
    	tmp = asin(x);
    end
    
    code[x_] := N[ArcSin[x], $MachinePrecision]
    
    \begin{array}{l}
    
    \\
    \sin^{-1} x
    \end{array}
    

    Reproduce

    ?
    herbie shell --seed 2024313 
    (FPCore (x)
      :name "Ian Simplification"
      :precision binary64
    
      :alt
      (! :herbie-platform default (asin x))
    
      (- (/ (PI) 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))