
(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 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 (- (fma (sin (* 0.3333333333333333 (acos (/ g (- h))))) (sqrt 3.0) (cos (* 0.3333333333333333 (fma (sqrt PI) (sqrt PI) (- (acos (/ g h)))))))))
double code(double g, double h) {
return -fma(sin((0.3333333333333333 * acos((g / -h)))), sqrt(3.0), cos((0.3333333333333333 * fma(sqrt(((double) M_PI)), sqrt(((double) M_PI)), -acos((g / h))))));
}
function code(g, h) return Float64(-fma(sin(Float64(0.3333333333333333 * acos(Float64(g / Float64(-h))))), sqrt(3.0), cos(Float64(0.3333333333333333 * fma(sqrt(pi), sqrt(pi), Float64(-acos(Float64(g / h)))))))) end
code[g_, h_] := (-N[(N[Sin[N[(0.3333333333333333 * N[ArcCos[N[(g / (-h)), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Sqrt[3.0], $MachinePrecision] + N[Cos[N[(0.3333333333333333 * N[(N[Sqrt[Pi], $MachinePrecision] * N[Sqrt[Pi], $MachinePrecision] + (-N[ArcCos[N[(g / h), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision])
\begin{array}{l}
\\
-\mathsf{fma}\left(\sin \left(0.3333333333333333 \cdot \cos^{-1} \left(\frac{g}{-h}\right)\right), \sqrt{3}, \cos \left(0.3333333333333333 \cdot \mathsf{fma}\left(\sqrt{\pi}, \sqrt{\pi}, -\cos^{-1} \left(\frac{g}{h}\right)\right)\right)\right)
\end{array}
Initial program 98.5%
Applied egg-rr99.9%
Applied egg-rr99.9%
Taylor expanded in g around 0
distribute-lft-outN/A
associate-*r*N/A
metadata-evalN/A
mul-1-negN/A
lower-neg.f64N/A
+-commutativeN/A
lower-fma.f64N/A
Simplified100.0%
lift-/.f64N/A
acos-negN/A
lift-PI.f64N/A
sub-negN/A
lift-PI.f64N/A
add-sqr-sqrtN/A
lower-fma.f64N/A
lift-PI.f64N/A
lower-sqrt.f64N/A
lift-PI.f64N/A
lower-sqrt.f64N/A
acos-asinN/A
sub-negN/A
asin-negN/A
lift-/.f64N/A
distribute-frac-neg2N/A
lift-neg.f64N/A
lift-/.f64N/A
lower-neg.f64N/A
lift-/.f64N/A
lift-neg.f64N/A
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (g h)
:precision binary64
(let* ((t_0 (acos (/ g (- h)))))
(-
(fma
(sin (* (* t_0 -3.0) -0.1111111111111111))
(sqrt 3.0)
(cos (* 0.3333333333333333 t_0))))))
double code(double g, double h) {
double t_0 = acos((g / -h));
return -fma(sin(((t_0 * -3.0) * -0.1111111111111111)), sqrt(3.0), cos((0.3333333333333333 * t_0)));
}
function code(g, h) t_0 = acos(Float64(g / Float64(-h))) return Float64(-fma(sin(Float64(Float64(t_0 * -3.0) * -0.1111111111111111)), sqrt(3.0), cos(Float64(0.3333333333333333 * t_0)))) end
code[g_, h_] := Block[{t$95$0 = N[ArcCos[N[(g / (-h)), $MachinePrecision]], $MachinePrecision]}, (-N[(N[Sin[N[(N[(t$95$0 * -3.0), $MachinePrecision] * -0.1111111111111111), $MachinePrecision]], $MachinePrecision] * N[Sqrt[3.0], $MachinePrecision] + N[Cos[N[(0.3333333333333333 * t$95$0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision])]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(\frac{g}{-h}\right)\\
-\mathsf{fma}\left(\sin \left(\left(t\_0 \cdot -3\right) \cdot -0.1111111111111111\right), \sqrt{3}, \cos \left(0.3333333333333333 \cdot t\_0\right)\right)
\end{array}
\end{array}
Initial program 98.5%
Applied egg-rr99.9%
Applied egg-rr99.9%
Taylor expanded in g around 0
distribute-lft-outN/A
associate-*r*N/A
metadata-evalN/A
mul-1-negN/A
lower-neg.f64N/A
+-commutativeN/A
lower-fma.f64N/A
Simplified100.0%
lift-/.f64N/A
lift-neg.f64N/A
lift-acos.f64N/A
*-commutativeN/A
metadata-evalN/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64100.0
lift-neg.f64N/A
lift-/.f64N/A
distribute-frac-neg2N/A
lift-neg.f64N/A
lift-/.f64100.0
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (g h)
:precision binary64
(let* ((t_0 (acos (/ g (- h)))))
(-
(fma
(sin (* 0.3333333333333333 t_0))
(sqrt 3.0)
(cos (* (* t_0 -3.0) -0.1111111111111111))))))
double code(double g, double h) {
double t_0 = acos((g / -h));
return -fma(sin((0.3333333333333333 * t_0)), sqrt(3.0), cos(((t_0 * -3.0) * -0.1111111111111111)));
}
function code(g, h) t_0 = acos(Float64(g / Float64(-h))) return Float64(-fma(sin(Float64(0.3333333333333333 * t_0)), sqrt(3.0), cos(Float64(Float64(t_0 * -3.0) * -0.1111111111111111)))) end
code[g_, h_] := Block[{t$95$0 = N[ArcCos[N[(g / (-h)), $MachinePrecision]], $MachinePrecision]}, (-N[(N[Sin[N[(0.3333333333333333 * t$95$0), $MachinePrecision]], $MachinePrecision] * N[Sqrt[3.0], $MachinePrecision] + N[Cos[N[(N[(t$95$0 * -3.0), $MachinePrecision] * -0.1111111111111111), $MachinePrecision]], $MachinePrecision]), $MachinePrecision])]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(\frac{g}{-h}\right)\\
-\mathsf{fma}\left(\sin \left(0.3333333333333333 \cdot t\_0\right), \sqrt{3}, \cos \left(\left(t\_0 \cdot -3\right) \cdot -0.1111111111111111\right)\right)
\end{array}
\end{array}
Initial program 98.5%
Applied egg-rr99.9%
Applied egg-rr99.9%
Taylor expanded in g around 0
distribute-lft-outN/A
associate-*r*N/A
metadata-evalN/A
mul-1-negN/A
lower-neg.f64N/A
+-commutativeN/A
lower-fma.f64N/A
Simplified100.0%
lift-/.f64N/A
lift-neg.f64N/A
lift-acos.f64N/A
*-commutativeN/A
metadata-evalN/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64100.0
lift-neg.f64N/A
lift-/.f64N/A
distribute-frac-neg2N/A
lift-neg.f64N/A
lift-/.f64100.0
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (g h) :precision binary64 (let* ((t_0 (* 0.3333333333333333 (acos (/ g (- h)))))) (- (fma (sin t_0) (sqrt 3.0) (cos t_0)))))
double code(double g, double h) {
double t_0 = 0.3333333333333333 * acos((g / -h));
return -fma(sin(t_0), sqrt(3.0), cos(t_0));
}
function code(g, h) t_0 = Float64(0.3333333333333333 * acos(Float64(g / Float64(-h)))) return Float64(-fma(sin(t_0), sqrt(3.0), cos(t_0))) end
code[g_, h_] := Block[{t$95$0 = N[(0.3333333333333333 * N[ArcCos[N[(g / (-h)), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, (-N[(N[Sin[t$95$0], $MachinePrecision] * N[Sqrt[3.0], $MachinePrecision] + N[Cos[t$95$0], $MachinePrecision]), $MachinePrecision])]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.3333333333333333 \cdot \cos^{-1} \left(\frac{g}{-h}\right)\\
-\mathsf{fma}\left(\sin t\_0, \sqrt{3}, \cos t\_0\right)
\end{array}
\end{array}
Initial program 98.5%
Applied egg-rr99.9%
Applied egg-rr99.9%
Taylor expanded in g around 0
distribute-lft-outN/A
associate-*r*N/A
metadata-evalN/A
mul-1-negN/A
lower-neg.f64N/A
+-commutativeN/A
lower-fma.f64N/A
Simplified100.0%
Final simplification100.0%
(FPCore (g h)
:precision binary64
(*
2.0
(cos
(fma
(* PI -0.1111111111111111)
-6.0
(* (* (acos (/ g (- h))) -3.0) -0.1111111111111111)))))
double code(double g, double h) {
return 2.0 * cos(fma((((double) M_PI) * -0.1111111111111111), -6.0, ((acos((g / -h)) * -3.0) * -0.1111111111111111)));
}
function code(g, h) return Float64(2.0 * cos(fma(Float64(pi * -0.1111111111111111), -6.0, Float64(Float64(acos(Float64(g / Float64(-h))) * -3.0) * -0.1111111111111111)))) end
code[g_, h_] := N[(2.0 * N[Cos[N[(N[(Pi * -0.1111111111111111), $MachinePrecision] * -6.0 + N[(N[(N[ArcCos[N[(g / (-h)), $MachinePrecision]], $MachinePrecision] * -3.0), $MachinePrecision] * -0.1111111111111111), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \cos \left(\mathsf{fma}\left(\pi \cdot -0.1111111111111111, -6, \left(\cos^{-1} \left(\frac{g}{-h}\right) \cdot -3\right) \cdot -0.1111111111111111\right)\right)
\end{array}
Initial program 98.5%
lift-PI.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-neg.f64N/A
lift-/.f64N/A
lift-acos.f64N/A
lift-/.f64N/A
+-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
frac-2negN/A
frac-addN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
Applied egg-rr98.5%
lift-neg.f64N/A
lift-/.f64N/A
lift-acos.f64N/A
lift-PI.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-/r/N/A
lift-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
Applied egg-rr98.5%
lift-PI.f64N/A
lift-*.f64N/A
lift-neg.f64N/A
lift-/.f64N/A
lift-acos.f64N/A
associate-*l*N/A
metadata-evalN/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6498.5
Applied egg-rr98.5%
lift-/.f64N/A
lift-neg.f64N/A
lift-acos.f64N/A
metadata-evalN/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f6498.5
lift-neg.f64N/A
lift-/.f64N/A
distribute-frac-neg2N/A
lift-neg.f64N/A
lift-/.f6498.5
Applied egg-rr98.5%
Final simplification98.5%
(FPCore (g h)
:precision binary64
(*
2.0
(cos
(fma
(* PI -0.1111111111111111)
-6.0
(* 0.3333333333333333 (acos (/ g (- h))))))))
double code(double g, double h) {
return 2.0 * cos(fma((((double) M_PI) * -0.1111111111111111), -6.0, (0.3333333333333333 * acos((g / -h)))));
}
function code(g, h) return Float64(2.0 * cos(fma(Float64(pi * -0.1111111111111111), -6.0, Float64(0.3333333333333333 * acos(Float64(g / Float64(-h))))))) end
code[g_, h_] := N[(2.0 * N[Cos[N[(N[(Pi * -0.1111111111111111), $MachinePrecision] * -6.0 + N[(0.3333333333333333 * N[ArcCos[N[(g / (-h)), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \cos \left(\mathsf{fma}\left(\pi \cdot -0.1111111111111111, -6, 0.3333333333333333 \cdot \cos^{-1} \left(\frac{g}{-h}\right)\right)\right)
\end{array}
Initial program 98.5%
lift-PI.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-neg.f64N/A
lift-/.f64N/A
lift-acos.f64N/A
lift-/.f64N/A
+-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
frac-2negN/A
frac-addN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
Applied egg-rr98.5%
lift-neg.f64N/A
lift-/.f64N/A
lift-acos.f64N/A
lift-PI.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-/r/N/A
lift-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
Applied egg-rr98.5%
lift-PI.f64N/A
lift-*.f64N/A
lift-neg.f64N/A
lift-/.f64N/A
lift-acos.f64N/A
associate-*l*N/A
metadata-evalN/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6498.5
Applied egg-rr98.5%
Final simplification98.5%
(FPCore (g h) :precision binary64 (* 2.0 (cos (* 0.3333333333333333 (fma 2.0 PI (acos (/ g (- h))))))))
double code(double g, double h) {
return 2.0 * cos((0.3333333333333333 * fma(2.0, ((double) M_PI), acos((g / -h)))));
}
function code(g, h) return Float64(2.0 * cos(Float64(0.3333333333333333 * fma(2.0, pi, acos(Float64(g / Float64(-h))))))) end
code[g_, h_] := N[(2.0 * N[Cos[N[(0.3333333333333333 * N[(2.0 * Pi + N[ArcCos[N[(g / (-h)), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \cos \left(0.3333333333333333 \cdot \mathsf{fma}\left(2, \pi, \cos^{-1} \left(\frac{g}{-h}\right)\right)\right)
\end{array}
Initial program 98.5%
lift-PI.f64N/A
lift-*.f64N/A
frac-2negN/A
lift-neg.f64N/A
lift-/.f64N/A
lift-acos.f64N/A
frac-2negN/A
frac-2negN/A
div-invN/A
frac-2negN/A
div-invN/A
distribute-rgt-outN/A
lower-*.f64N/A
metadata-evalN/A
lift-*.f64N/A
lower-fma.f6498.5
lift-/.f64N/A
frac-2negN/A
Applied egg-rr98.5%
herbie shell --seed 2024207
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
:precision binary64
(* 2.0 (cos (+ (/ (* 2.0 PI) 3.0) (/ (acos (/ (- g) h)) 3.0)))))