
(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}
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)))))
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
(let* ((t_0 (acos (/ (- g) h)))
(t_1 (* 0.5 (sqrt 3.0)))
(t_2 (/ (sqrt 3.0) 2.0))
(t_3 (* (- 0.25 (* t_2 t_2)) (cos (* t_0 -0.3333333333333333))))
(t_4 (* 0.3333333333333333 t_0))
(t_5 (sin t_4)))
(*
(/
(- (pow t_3 3.0) (pow (* 0.75 (- 0.5 (* (cos (* t_4 2.0)) 0.5))) 1.5))
(+ (* (fma t_5 t_1 t_3) (* t_1 t_5)) (pow t_3 2.0)))
2.0)))
double code(double g, double h) {
double t_0 = acos((-g / h));
double t_1 = 0.5 * sqrt(3.0);
double t_2 = sqrt(3.0) / 2.0;
double t_3 = (0.25 - (t_2 * t_2)) * cos((t_0 * -0.3333333333333333));
double t_4 = 0.3333333333333333 * t_0;
double t_5 = sin(t_4);
return ((pow(t_3, 3.0) - pow((0.75 * (0.5 - (cos((t_4 * 2.0)) * 0.5))), 1.5)) / ((fma(t_5, t_1, t_3) * (t_1 * t_5)) + pow(t_3, 2.0))) * 2.0;
}
function code(g, h) t_0 = acos(Float64(Float64(-g) / h)) t_1 = Float64(0.5 * sqrt(3.0)) t_2 = Float64(sqrt(3.0) / 2.0) t_3 = Float64(Float64(0.25 - Float64(t_2 * t_2)) * cos(Float64(t_0 * -0.3333333333333333))) t_4 = Float64(0.3333333333333333 * t_0) t_5 = sin(t_4) return Float64(Float64(Float64((t_3 ^ 3.0) - (Float64(0.75 * Float64(0.5 - Float64(cos(Float64(t_4 * 2.0)) * 0.5))) ^ 1.5)) / Float64(Float64(fma(t_5, t_1, t_3) * Float64(t_1 * t_5)) + (t_3 ^ 2.0))) * 2.0) end
code[g_, h_] := Block[{t$95$0 = N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(0.5 * N[Sqrt[3.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Sqrt[3.0], $MachinePrecision] / 2.0), $MachinePrecision]}, Block[{t$95$3 = N[(N[(0.25 - N[(t$95$2 * t$95$2), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(t$95$0 * -0.3333333333333333), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(0.3333333333333333 * t$95$0), $MachinePrecision]}, Block[{t$95$5 = N[Sin[t$95$4], $MachinePrecision]}, N[(N[(N[(N[Power[t$95$3, 3.0], $MachinePrecision] - N[Power[N[(0.75 * N[(0.5 - N[(N[Cos[N[(t$95$4 * 2.0), $MachinePrecision]], $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.5], $MachinePrecision]), $MachinePrecision] / N[(N[(N[(t$95$5 * t$95$1 + t$95$3), $MachinePrecision] * N[(t$95$1 * t$95$5), $MachinePrecision]), $MachinePrecision] + N[Power[t$95$3, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(\frac{-g}{h}\right)\\
t_1 := 0.5 \cdot \sqrt{3}\\
t_2 := \frac{\sqrt{3}}{2}\\
t_3 := \left(0.25 - t\_2 \cdot t\_2\right) \cdot \cos \left(t\_0 \cdot -0.3333333333333333\right)\\
t_4 := 0.3333333333333333 \cdot t\_0\\
t_5 := \sin t\_4\\
\frac{{t\_3}^{3} - {\left(0.75 \cdot \left(0.5 - \cos \left(t\_4 \cdot 2\right) \cdot 0.5\right)\right)}^{1.5}}{\mathsf{fma}\left(t\_5, t\_1, t\_3\right) \cdot \left(t\_1 \cdot t\_5\right) + {t\_3}^{2}} \cdot 2
\end{array}
\end{array}
Initial program 98.5%
Applied rewrites99.9%
lift-pow.f64N/A
sqr-powN/A
pow-prod-downN/A
lower-pow.f64N/A
Applied rewrites99.9%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
div-invN/A
metadata-evalN/A
associate-*l*N/A
metadata-evalN/A
*-rgt-identity99.9
Applied rewrites99.9%
lift-*.f64N/A
lift-/.f64N/A
div-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identity99.9
Applied rewrites99.9%
Final simplification99.9%
(FPCore (g h) :precision binary64 (* (cos (/ (fma (acos (/ (- g) h)) 3.0 (* PI 6.0)) 9.0)) 2.0))
double code(double g, double h) {
return cos((fma(acos((-g / h)), 3.0, (((double) M_PI) * 6.0)) / 9.0)) * 2.0;
}
function code(g, h) return Float64(cos(Float64(fma(acos(Float64(Float64(-g) / h)), 3.0, Float64(pi * 6.0)) / 9.0)) * 2.0) end
code[g_, h_] := N[(N[Cos[N[(N[(N[ArcCos[N[((-g) / h), $MachinePrecision]], $MachinePrecision] * 3.0 + N[(Pi * 6.0), $MachinePrecision]), $MachinePrecision] / 9.0), $MachinePrecision]], $MachinePrecision] * 2.0), $MachinePrecision]
\begin{array}{l}
\\
\cos \left(\frac{\mathsf{fma}\left(\cos^{-1} \left(\frac{-g}{h}\right), 3, \pi \cdot 6\right)}{9}\right) \cdot 2
\end{array}
Initial program 98.5%
lift-+.f64N/A
lift-/.f64N/A
lift-/.f64N/A
frac-addN/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
metadata-evalN/A
metadata-evalN/A
lower-*.f64N/A
metadata-evalN/A
metadata-eval98.5
Applied rewrites98.5%
Final simplification98.5%
(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
lift-/.f64N/A
lift-/.f64N/A
frac-addN/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
metadata-evalN/A
metadata-evalN/A
lower-*.f64N/A
metadata-evalN/A
metadata-eval98.5
Applied rewrites98.5%
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lift-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-eval98.5
Applied rewrites98.5%
lift-fma.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6498.5
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
metadata-evalN/A
lift-*.f6498.5
Applied rewrites98.5%
Final simplification98.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
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 2024237
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
:precision binary64
(* 2.0 (cos (+ (/ (* 2.0 PI) 3.0) (/ (acos (/ (- g) h)) 3.0)))))