
(FPCore (g h) :precision binary64 (* 2.0 (cos (+ (/ (* 2.0 (PI)) 3.0) (/ (acos (/ (- g) h)) 3.0)))))
\begin{array}{l}
\\
2 \cdot \cos \left(\frac{2 \cdot \mathsf{PI}\left(\right)}{3} + \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{3}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 4 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (g h) :precision binary64 (* 2.0 (cos (+ (/ (* 2.0 (PI)) 3.0) (/ (acos (/ (- g) h)) 3.0)))))
\begin{array}{l}
\\
2 \cdot \cos \left(\frac{2 \cdot \mathsf{PI}\left(\right)}{3} + \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{3}\right)
\end{array}
(FPCore (g h)
:precision binary64
(let* ((t_0 (acos (/ (- g) h)))
(t_1
(cos
(*
(fma
(* 0.4444444444444444 (PI))
(PI)
(* 0.1111111111111111 (pow t_0 2.0)))
(/ -3.0 (fma -2.0 (PI) t_0))))))
(*
(*
(pow t_1 -1.0)
(* t_1 (cos (fma t_0 0.3333333333333333 (* (PI) 0.6666666666666666)))))
2.0)))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(\frac{-g}{h}\right)\\
t_1 := \cos \left(\mathsf{fma}\left(0.4444444444444444 \cdot \mathsf{PI}\left(\right), \mathsf{PI}\left(\right), 0.1111111111111111 \cdot {t\_0}^{2}\right) \cdot \frac{-3}{\mathsf{fma}\left(-2, \mathsf{PI}\left(\right), t\_0\right)}\right)\\
\left({t\_1}^{-1} \cdot \left(t\_1 \cdot \cos \left(\mathsf{fma}\left(t\_0, 0.3333333333333333, \mathsf{PI}\left(\right) \cdot 0.6666666666666666\right)\right)\right)\right) \cdot 2
\end{array}
\end{array}
Initial program 98.5%
lift-+.f64N/A
lift-/.f64N/A
div-invN/A
lift-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
metadata-evalN/A
metadata-eval98.5
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lower-*.f64N/A
metadata-eval98.5
Applied rewrites98.5%
Applied rewrites100.0%
Final simplification100.0%
(FPCore (g h)
:precision binary64
(let* ((t_0 (acos (/ (- g) h))) (t_1 (* -0.6666666666666666 (PI))))
(*
(- 2.0)
(*
(/
(cos (fma 0.3333333333333333 t_0 t_1))
(/ -1.0 (cos (fma 0.3333333333333333 t_0 (* (PI) 0.6666666666666666)))))
(pow (cos (fma t_0 0.3333333333333333 t_1)) -1.0)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(\frac{-g}{h}\right)\\
t_1 := -0.6666666666666666 \cdot \mathsf{PI}\left(\right)\\
\left(-2\right) \cdot \left(\frac{\cos \left(\mathsf{fma}\left(0.3333333333333333, t\_0, t\_1\right)\right)}{\frac{-1}{\cos \left(\mathsf{fma}\left(0.3333333333333333, t\_0, \mathsf{PI}\left(\right) \cdot 0.6666666666666666\right)\right)}} \cdot {\cos \left(\mathsf{fma}\left(t\_0, 0.3333333333333333, t\_1\right)\right)}^{-1}\right)
\end{array}
\end{array}
Initial program 98.5%
Applied rewrites98.5%
Applied rewrites98.5%
lift-pow.f64N/A
unpow-1N/A
frac-2negN/A
metadata-evalN/A
lift-/.f64N/A
/-rgt-identityN/A
lower-/.f64N/A
Applied rewrites98.5%
Applied rewrites100.0%
Final simplification100.0%
(FPCore (g h)
:precision binary64
(let* ((t_0 (acos (/ (- g) h)))
(t_1 (cos (fma t_0 0.3333333333333333 (* -0.6666666666666666 (PI))))))
(*
(/ (cos (fma t_0 0.3333333333333333 (* (PI) 0.6666666666666666))) t_1)
(* t_1 2.0))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(\frac{-g}{h}\right)\\
t_1 := \cos \left(\mathsf{fma}\left(t\_0, 0.3333333333333333, -0.6666666666666666 \cdot \mathsf{PI}\left(\right)\right)\right)\\
\frac{\cos \left(\mathsf{fma}\left(t\_0, 0.3333333333333333, \mathsf{PI}\left(\right) \cdot 0.6666666666666666\right)\right)}{t\_1} \cdot \left(t\_1 \cdot 2\right)
\end{array}
\end{array}
Initial program 98.5%
lift-+.f64N/A
lift-/.f64N/A
div-invN/A
lift-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
metadata-evalN/A
metadata-eval98.5
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lower-*.f64N/A
metadata-eval98.5
Applied rewrites98.5%
Applied rewrites99.9%
Final simplification99.9%
(FPCore (g h) :precision binary64 (* (cos (fma (PI) 0.6666666666666666 (* 0.3333333333333333 (acos (/ (- g) h))))) 2.0))
\begin{array}{l}
\\
\cos \left(\mathsf{fma}\left(\mathsf{PI}\left(\right), 0.6666666666666666, 0.3333333333333333 \cdot \cos^{-1} \left(\frac{-g}{h}\right)\right)\right) \cdot 2
\end{array}
Initial program 98.5%
lift-+.f64N/A
lift-/.f64N/A
div-invN/A
lift-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
metadata-evalN/A
metadata-eval98.5
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lower-*.f64N/A
metadata-eval98.5
Applied rewrites98.5%
Final simplification98.5%
herbie shell --seed 2024332
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
:precision binary64
(* 2.0 (cos (+ (/ (* 2.0 (PI)) 3.0) (/ (acos (/ (- g) h)) 3.0)))))