
(FPCore (g h) :precision binary64 (* 2.0 (cos (+ (/ (* 2.0 PI) 3.0) (/ (acos (/ (- g) h)) 3.0)))))
double code(double g, double h) {
return 2.0 * cos((((2.0 * ((double) M_PI)) / 3.0) + (acos((-g / h)) / 3.0)));
}
public static double code(double g, double h) {
return 2.0 * Math.cos((((2.0 * Math.PI) / 3.0) + (Math.acos((-g / h)) / 3.0)));
}
def code(g, h): return 2.0 * math.cos((((2.0 * math.pi) / 3.0) + (math.acos((-g / h)) / 3.0)))
function code(g, h) return Float64(2.0 * cos(Float64(Float64(Float64(2.0 * pi) / 3.0) + Float64(acos(Float64(Float64(-g) / h)) / 3.0)))) end
function tmp = code(g, h) tmp = 2.0 * cos((((2.0 * pi) / 3.0) + (acos((-g / h)) / 3.0))); end
code[g_, h_] := N[(2.0 * N[Cos[N[(N[(N[(2.0 * Pi), $MachinePrecision] / 3.0), $MachinePrecision] + N[(N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \cos \left(\frac{2 \cdot \pi}{3} + \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{3}\right)
\end{array}
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (g h) :precision binary64 (* 2.0 (cos (+ (/ (* 2.0 PI) 3.0) (/ (acos (/ (- g) h)) 3.0)))))
double code(double g, double h) {
return 2.0 * cos((((2.0 * ((double) M_PI)) / 3.0) + (acos((-g / h)) / 3.0)));
}
public static double code(double g, double h) {
return 2.0 * Math.cos((((2.0 * Math.PI) / 3.0) + (Math.acos((-g / h)) / 3.0)));
}
def code(g, h): return 2.0 * math.cos((((2.0 * math.pi) / 3.0) + (math.acos((-g / h)) / 3.0)))
function code(g, h) return Float64(2.0 * cos(Float64(Float64(Float64(2.0 * pi) / 3.0) + Float64(acos(Float64(Float64(-g) / h)) / 3.0)))) end
function tmp = code(g, h) tmp = 2.0 * cos((((2.0 * pi) / 3.0) + (acos((-g / h)) / 3.0))); end
code[g_, h_] := N[(2.0 * N[Cos[N[(N[(N[(2.0 * Pi), $MachinePrecision] / 3.0), $MachinePrecision] + N[(N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \cos \left(\frac{2 \cdot \pi}{3} + \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{3}\right)
\end{array}
(FPCore (g h)
:precision binary64
(*
(sin
(fma
(sqrt PI)
(/ (sqrt PI) 2.0)
(fma (acos (/ (- g) h)) -0.3333333333333333 (* -0.6666666666666666 PI))))
2.0))
double code(double g, double h) {
return sin(fma(sqrt(((double) M_PI)), (sqrt(((double) M_PI)) / 2.0), fma(acos((-g / h)), -0.3333333333333333, (-0.6666666666666666 * ((double) M_PI))))) * 2.0;
}
function code(g, h) return Float64(sin(fma(sqrt(pi), Float64(sqrt(pi) / 2.0), fma(acos(Float64(Float64(-g) / h)), -0.3333333333333333, Float64(-0.6666666666666666 * pi)))) * 2.0) end
code[g_, h_] := N[(N[Sin[N[(N[Sqrt[Pi], $MachinePrecision] * N[(N[Sqrt[Pi], $MachinePrecision] / 2.0), $MachinePrecision] + N[(N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision] * -0.3333333333333333 + N[(-0.6666666666666666 * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 2.0), $MachinePrecision]
\begin{array}{l}
\\
\sin \left(\mathsf{fma}\left(\sqrt{\pi}, \frac{\sqrt{\pi}}{2}, \mathsf{fma}\left(\cos^{-1} \left(\frac{-g}{h}\right), -0.3333333333333333, -0.6666666666666666 \cdot \pi\right)\right)\right) \cdot 2
\end{array}
Initial program 98.5%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.5
lift-+.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-add-revN/A
lower-/.f64N/A
lift-PI.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lift-PI.f6498.5
Applied rewrites98.5%
lift-/.f64N/A
lift-PI.f64N/A
lift-fma.f64N/A
div-addN/A
*-commutativeN/A
frac-addN/A
*-commutativeN/A
metadata-evalN/A
div-addN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
lift-*.f64N/A
lift-PI.f64N/A
Applied rewrites98.4%
lift-cos.f64N/A
lift-*.f64N/A
lift-PI.f64N/A
lift-fma.f64N/A
lift-acos.f64N/A
lift-/.f64N/A
lift-neg.f64N/A
cos-neg-revN/A
sin-+PI/2-revN/A
lower-sin.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
lift-/.f64N/A
lift-PI.f64N/A
Applied rewrites99.9%
lift-fma.f64N/A
lift-PI.f64N/A
lift-fma.f64N/A
lift-acos.f64N/A
lift-/.f64N/A
lift-neg.f64N/A
lift-PI.f64N/A
lift-/.f64N/A
+-commutativeN/A
add-sqr-sqrtN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites100.0%
(FPCore (g h)
:precision binary64
(*
(sin
(fma
(acos (/ (- g) h))
-0.3333333333333333
(fma -0.6666666666666666 PI (* 0.5 PI))))
2.0))
double code(double g, double h) {
return sin(fma(acos((-g / h)), -0.3333333333333333, fma(-0.6666666666666666, ((double) M_PI), (0.5 * ((double) M_PI))))) * 2.0;
}
function code(g, h) return Float64(sin(fma(acos(Float64(Float64(-g) / h)), -0.3333333333333333, fma(-0.6666666666666666, pi, Float64(0.5 * pi)))) * 2.0) end
code[g_, h_] := N[(N[Sin[N[(N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision] * -0.3333333333333333 + N[(-0.6666666666666666 * Pi + N[(0.5 * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 2.0), $MachinePrecision]
\begin{array}{l}
\\
\sin \left(\mathsf{fma}\left(\cos^{-1} \left(\frac{-g}{h}\right), -0.3333333333333333, \mathsf{fma}\left(-0.6666666666666666, \pi, 0.5 \cdot \pi\right)\right)\right) \cdot 2
\end{array}
Initial program 98.5%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.5
lift-+.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-add-revN/A
lower-/.f64N/A
lift-PI.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lift-PI.f6498.5
Applied rewrites98.5%
lift-/.f64N/A
lift-PI.f64N/A
lift-fma.f64N/A
div-addN/A
*-commutativeN/A
frac-addN/A
*-commutativeN/A
metadata-evalN/A
div-addN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
lift-*.f64N/A
lift-PI.f64N/A
Applied rewrites98.4%
lift-cos.f64N/A
lift-*.f64N/A
lift-PI.f64N/A
lift-fma.f64N/A
lift-acos.f64N/A
lift-/.f64N/A
lift-neg.f64N/A
cos-neg-revN/A
sin-+PI/2-revN/A
lower-sin.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
lift-/.f64N/A
lift-PI.f64N/A
Applied rewrites99.9%
Taylor expanded in g around 0
distribute-lft-inN/A
mul-1-negN/A
distribute-frac-negN/A
lift-neg.f64N/A
lift-/.f64N/A
lift-acos.f64N/A
lift-*.f64N/A
associate-*r*N/A
metadata-evalN/A
associate-+l+N/A
lift-*.f64N/A
*-commutativeN/A
Applied rewrites99.9%
(FPCore (g h)
:precision binary64
(*
(cos
(fma
(* 0.1111111111111111 PI)
6.0
(* 0.3333333333333333 (acos (/ (- g) h)))))
2.0))
double code(double g, double h) {
return cos(fma((0.1111111111111111 * ((double) M_PI)), 6.0, (0.3333333333333333 * acos((-g / h))))) * 2.0;
}
function code(g, h) return Float64(cos(fma(Float64(0.1111111111111111 * pi), 6.0, Float64(0.3333333333333333 * acos(Float64(Float64(-g) / h))))) * 2.0) end
code[g_, h_] := N[(N[Cos[N[(N[(0.1111111111111111 * Pi), $MachinePrecision] * 6.0 + N[(0.3333333333333333 * N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 2.0), $MachinePrecision]
\begin{array}{l}
\\
\cos \left(\mathsf{fma}\left(0.1111111111111111 \cdot \pi, 6, 0.3333333333333333 \cdot \cos^{-1} \left(\frac{-g}{h}\right)\right)\right) \cdot 2
\end{array}
Initial program 98.5%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.5
lift-+.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-add-revN/A
lower-/.f64N/A
lift-PI.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lift-PI.f6498.5
Applied rewrites98.5%
lift-/.f64N/A
lift-PI.f64N/A
lift-fma.f64N/A
div-addN/A
*-commutativeN/A
frac-addN/A
*-commutativeN/A
metadata-evalN/A
div-addN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
lift-*.f64N/A
lift-PI.f64N/A
Applied rewrites98.5%
(FPCore (g h) :precision binary64 (* (cos (/ (fma PI 2.0 (acos (/ (- g) h))) 3.0)) 2.0))
double code(double g, double h) {
return cos((fma(((double) M_PI), 2.0, acos((-g / h))) / 3.0)) * 2.0;
}
function code(g, h) return Float64(cos(Float64(fma(pi, 2.0, acos(Float64(Float64(-g) / h))) / 3.0)) * 2.0) end
code[g_, h_] := N[(N[Cos[N[(N[(Pi * 2.0 + N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / 3.0), $MachinePrecision]], $MachinePrecision] * 2.0), $MachinePrecision]
\begin{array}{l}
\\
\cos \left(\frac{\mathsf{fma}\left(\pi, 2, \cos^{-1} \left(\frac{-g}{h}\right)\right)}{3}\right) \cdot 2
\end{array}
Initial program 98.5%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.5
lift-+.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-add-revN/A
lower-/.f64N/A
lift-PI.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lift-PI.f6498.5
Applied rewrites98.5%
(FPCore (g h) :precision binary64 (* (cos (fma PI 0.6666666666666666 (* 0.3333333333333333 (acos (/ (- g) h))))) 2.0))
double code(double g, double h) {
return cos(fma(((double) M_PI), 0.6666666666666666, (0.3333333333333333 * acos((-g / h))))) * 2.0;
}
function code(g, h) return Float64(cos(fma(pi, 0.6666666666666666, Float64(0.3333333333333333 * acos(Float64(Float64(-g) / h))))) * 2.0) end
code[g_, h_] := N[(N[Cos[N[(Pi * 0.6666666666666666 + N[(0.3333333333333333 * N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 2.0), $MachinePrecision]
\begin{array}{l}
\\
\cos \left(\mathsf{fma}\left(\pi, 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
*-commutativeN/A
lower-*.f6498.5
lift-+.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-add-revN/A
lower-/.f64N/A
lift-PI.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lift-PI.f6498.5
Applied rewrites98.5%
lift-/.f64N/A
lift-PI.f64N/A
lift-fma.f64N/A
div-addN/A
*-commutativeN/A
frac-addN/A
*-commutativeN/A
metadata-evalN/A
div-addN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
lift-*.f64N/A
lift-PI.f64N/A
Applied rewrites98.5%
(FPCore (g h) :precision binary64 (* 2.0 (cos (fma 0.3333333333333333 (acos (/ (- g) h)) (* 0.6666666666666666 PI)))))
double code(double g, double h) {
return 2.0 * cos(fma(0.3333333333333333, acos((-g / h)), (0.6666666666666666 * ((double) M_PI))));
}
function code(g, h) return Float64(2.0 * cos(fma(0.3333333333333333, acos(Float64(Float64(-g) / h)), Float64(0.6666666666666666 * pi)))) end
code[g_, h_] := N[(2.0 * N[Cos[N[(0.3333333333333333 * N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision] + N[(0.6666666666666666 * Pi), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \cos \left(\mathsf{fma}\left(0.3333333333333333, \cos^{-1} \left(\frac{-g}{h}\right), 0.6666666666666666 \cdot \pi\right)\right)
\end{array}
Initial program 98.5%
Taylor expanded in g around 0
mul-1-negN/A
distribute-frac-negN/A
lift-neg.f64N/A
lift-/.f64N/A
lift-acos.f64N/A
lower-fma.f64N/A
lower-*.f64N/A
lift-PI.f6498.4
Applied rewrites98.4%
(FPCore (g h) :precision binary64 (* (sin (fma 0.3333333333333333 (acos (/ (- g) h)) (* 1.1666666666666667 PI))) 2.0))
double code(double g, double h) {
return sin(fma(0.3333333333333333, acos((-g / h)), (1.1666666666666667 * ((double) M_PI)))) * 2.0;
}
function code(g, h) return Float64(sin(fma(0.3333333333333333, acos(Float64(Float64(-g) / h)), Float64(1.1666666666666667 * pi))) * 2.0) end
code[g_, h_] := N[(N[Sin[N[(0.3333333333333333 * N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision] + N[(1.1666666666666667 * Pi), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 2.0), $MachinePrecision]
\begin{array}{l}
\\
\sin \left(\mathsf{fma}\left(0.3333333333333333, \cos^{-1} \left(\frac{-g}{h}\right), 1.1666666666666667 \cdot \pi\right)\right) \cdot 2
\end{array}
Initial program 98.5%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.5
lift-+.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-add-revN/A
lower-/.f64N/A
lift-PI.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lift-PI.f6498.5
Applied rewrites98.5%
Applied rewrites97.5%
Taylor expanded in g around 0
mul-1-negN/A
distribute-frac-negN/A
lift-neg.f64N/A
lift-/.f64N/A
lift-acos.f64N/A
lower-fma.f64N/A
distribute-rgt-outN/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
lift-PI.f6497.6
Applied rewrites97.6%
herbie shell --seed 2025130
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
:precision binary64
(* 2.0 (cos (+ (/ (* 2.0 PI) 3.0) (/ (acos (/ (- g) h)) 3.0)))))