
(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 6 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 (* -0.3333333333333333 (asin (/ g h)))))
(+
(sin
(+
t_0
(*
(pow (* (* PI PI) (* PI (* PI (* PI PI)))) 0.16666666666666666)
-0.6666666666666666)))
(/ (cos (+ t_0 (* PI -0.6666666666666666))) (sqrt 0.3333333333333333)))))
double code(double g, double h) {
double t_0 = -0.3333333333333333 * asin((g / h));
return sin((t_0 + (pow(((((double) M_PI) * ((double) M_PI)) * (((double) M_PI) * (((double) M_PI) * (((double) M_PI) * ((double) M_PI))))), 0.16666666666666666) * -0.6666666666666666))) + (cos((t_0 + (((double) M_PI) * -0.6666666666666666))) / sqrt(0.3333333333333333));
}
public static double code(double g, double h) {
double t_0 = -0.3333333333333333 * Math.asin((g / h));
return Math.sin((t_0 + (Math.pow(((Math.PI * Math.PI) * (Math.PI * (Math.PI * (Math.PI * Math.PI)))), 0.16666666666666666) * -0.6666666666666666))) + (Math.cos((t_0 + (Math.PI * -0.6666666666666666))) / Math.sqrt(0.3333333333333333));
}
def code(g, h): t_0 = -0.3333333333333333 * math.asin((g / h)) return math.sin((t_0 + (math.pow(((math.pi * math.pi) * (math.pi * (math.pi * (math.pi * math.pi)))), 0.16666666666666666) * -0.6666666666666666))) + (math.cos((t_0 + (math.pi * -0.6666666666666666))) / math.sqrt(0.3333333333333333))
function code(g, h) t_0 = Float64(-0.3333333333333333 * asin(Float64(g / h))) return Float64(sin(Float64(t_0 + Float64((Float64(Float64(pi * pi) * Float64(pi * Float64(pi * Float64(pi * pi)))) ^ 0.16666666666666666) * -0.6666666666666666))) + Float64(cos(Float64(t_0 + Float64(pi * -0.6666666666666666))) / sqrt(0.3333333333333333))) end
function tmp = code(g, h) t_0 = -0.3333333333333333 * asin((g / h)); tmp = sin((t_0 + ((((pi * pi) * (pi * (pi * (pi * pi)))) ^ 0.16666666666666666) * -0.6666666666666666))) + (cos((t_0 + (pi * -0.6666666666666666))) / sqrt(0.3333333333333333)); end
code[g_, h_] := Block[{t$95$0 = N[(-0.3333333333333333 * N[ArcSin[N[(g / h), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(N[Sin[N[(t$95$0 + N[(N[Power[N[(N[(Pi * Pi), $MachinePrecision] * N[(Pi * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.16666666666666666], $MachinePrecision] * -0.6666666666666666), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + N[(N[Cos[N[(t$95$0 + N[(Pi * -0.6666666666666666), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[Sqrt[0.3333333333333333], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -0.3333333333333333 \cdot \sin^{-1} \left(\frac{g}{h}\right)\\
\sin \left(t\_0 + {\left(\left(\pi \cdot \pi\right) \cdot \left(\pi \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right)}^{0.16666666666666666} \cdot -0.6666666666666666\right) + \frac{\cos \left(t\_0 + \pi \cdot -0.6666666666666666\right)}{\sqrt{0.3333333333333333}}
\end{array}
\end{array}
Initial program 98.4%
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
acos-lowering-acos.f64N/A
distribute-frac-negN/A
distribute-neg-frac2N/A
neg-mul-1N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6498.4%
Simplified98.4%
Applied egg-rr98.4%
clear-numN/A
inv-powN/A
div-invN/A
unpow-prod-downN/A
metadata-evalN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
metadata-eval100.0%
Applied egg-rr100.0%
Taylor expanded in g around 0
distribute-lft-outN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified98.5%
add-cbrt-cubeN/A
associate-*r*N/A
pow1/3N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
pow-prod-upN/A
pow-prod-downN/A
pow-lowering-pow.f64N/A
Applied egg-rr100.0%
(FPCore (g h)
:precision binary64
(let* ((t_0 (asin (/ g h))))
(+
(sin (+ (* -0.3333333333333333 t_0) (* PI -0.6666666666666666)))
(*
(pow 0.3333333333333333 -0.5)
(cos (+ (* PI -0.6666666666666666) (/ t_0 -3.0)))))))
double code(double g, double h) {
double t_0 = asin((g / h));
return sin(((-0.3333333333333333 * t_0) + (((double) M_PI) * -0.6666666666666666))) + (pow(0.3333333333333333, -0.5) * cos(((((double) M_PI) * -0.6666666666666666) + (t_0 / -3.0))));
}
public static double code(double g, double h) {
double t_0 = Math.asin((g / h));
return Math.sin(((-0.3333333333333333 * t_0) + (Math.PI * -0.6666666666666666))) + (Math.pow(0.3333333333333333, -0.5) * Math.cos(((Math.PI * -0.6666666666666666) + (t_0 / -3.0))));
}
def code(g, h): t_0 = math.asin((g / h)) return math.sin(((-0.3333333333333333 * t_0) + (math.pi * -0.6666666666666666))) + (math.pow(0.3333333333333333, -0.5) * math.cos(((math.pi * -0.6666666666666666) + (t_0 / -3.0))))
function code(g, h) t_0 = asin(Float64(g / h)) return Float64(sin(Float64(Float64(-0.3333333333333333 * t_0) + Float64(pi * -0.6666666666666666))) + Float64((0.3333333333333333 ^ -0.5) * cos(Float64(Float64(pi * -0.6666666666666666) + Float64(t_0 / -3.0))))) end
function tmp = code(g, h) t_0 = asin((g / h)); tmp = sin(((-0.3333333333333333 * t_0) + (pi * -0.6666666666666666))) + ((0.3333333333333333 ^ -0.5) * cos(((pi * -0.6666666666666666) + (t_0 / -3.0)))); end
code[g_, h_] := Block[{t$95$0 = N[ArcSin[N[(g / h), $MachinePrecision]], $MachinePrecision]}, N[(N[Sin[N[(N[(-0.3333333333333333 * t$95$0), $MachinePrecision] + N[(Pi * -0.6666666666666666), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + N[(N[Power[0.3333333333333333, -0.5], $MachinePrecision] * N[Cos[N[(N[(Pi * -0.6666666666666666), $MachinePrecision] + N[(t$95$0 / -3.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin^{-1} \left(\frac{g}{h}\right)\\
\sin \left(-0.3333333333333333 \cdot t\_0 + \pi \cdot -0.6666666666666666\right) + {0.3333333333333333}^{-0.5} \cdot \cos \left(\pi \cdot -0.6666666666666666 + \frac{t\_0}{-3}\right)
\end{array}
\end{array}
Initial program 98.4%
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
acos-lowering-acos.f64N/A
distribute-frac-negN/A
distribute-neg-frac2N/A
neg-mul-1N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6498.4%
Simplified98.4%
Applied egg-rr98.4%
clear-numN/A
inv-powN/A
div-invN/A
unpow-prod-downN/A
metadata-evalN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
metadata-eval100.0%
Applied egg-rr100.0%
Taylor expanded in g around 0
distribute-lft-outN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Simplified98.5%
clear-numN/A
associate-/r/N/A
*-commutativeN/A
metadata-evalN/A
div-invN/A
*-commutativeN/A
metadata-evalN/A
cancel-sign-sub-invN/A
*-commutativeN/A
*-lowering-*.f64N/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
cos-lowering-cos.f64N/A
sub-negN/A
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (g h)
:precision binary64
(*
2.0
(cos
(fma
(cbrt (* PI (* PI PI)))
0.6666666666666666
(/ (acos (- 0.0 (/ g h))) 3.0)))))
double code(double g, double h) {
return 2.0 * cos(fma(cbrt((((double) M_PI) * (((double) M_PI) * ((double) M_PI)))), 0.6666666666666666, (acos((0.0 - (g / h))) / 3.0)));
}
function code(g, h) return Float64(2.0 * cos(fma(cbrt(Float64(pi * Float64(pi * pi))), 0.6666666666666666, Float64(acos(Float64(0.0 - Float64(g / h))) / 3.0)))) end
code[g_, h_] := N[(2.0 * N[Cos[N[(N[Power[N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * 0.6666666666666666 + N[(N[ArcCos[N[(0.0 - N[(g / h), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \cos \left(\mathsf{fma}\left(\sqrt[3]{\pi \cdot \left(\pi \cdot \pi\right)}, 0.6666666666666666, \frac{\cos^{-1} \left(0 - \frac{g}{h}\right)}{3}\right)\right)
\end{array}
Initial program 98.4%
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
acos-lowering-acos.f64N/A
distribute-frac-negN/A
distribute-neg-frac2N/A
neg-mul-1N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6498.4%
Simplified98.4%
fma-defineN/A
fma-lowering-fma.f64N/A
PI-lowering-PI.f64N/A
/-lowering-/.f64N/A
Applied egg-rr98.4%
div-invN/A
clear-numN/A
mul-1-negN/A
neg-lowering-neg.f64N/A
/-lowering-/.f6498.4%
Applied egg-rr98.4%
add-cbrt-cubeN/A
associate-*r*N/A
cbrt-lowering-cbrt.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f6498.5%
Applied egg-rr98.5%
Final simplification98.5%
(FPCore (g h) :precision binary64 (* 2.0 (cos (fma PI 0.6666666666666666 (/ 1.0 (/ 3.0 (acos (- 0.0 (/ g h)))))))))
double code(double g, double h) {
return 2.0 * cos(fma(((double) M_PI), 0.6666666666666666, (1.0 / (3.0 / acos((0.0 - (g / h)))))));
}
function code(g, h) return Float64(2.0 * cos(fma(pi, 0.6666666666666666, Float64(1.0 / Float64(3.0 / acos(Float64(0.0 - Float64(g / h)))))))) end
code[g_, h_] := N[(2.0 * N[Cos[N[(Pi * 0.6666666666666666 + N[(1.0 / N[(3.0 / N[ArcCos[N[(0.0 - N[(g / h), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \cos \left(\mathsf{fma}\left(\pi, 0.6666666666666666, \frac{1}{\frac{3}{\cos^{-1} \left(0 - \frac{g}{h}\right)}}\right)\right)
\end{array}
Initial program 98.4%
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
acos-lowering-acos.f64N/A
distribute-frac-negN/A
distribute-neg-frac2N/A
neg-mul-1N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6498.4%
Simplified98.4%
fma-defineN/A
fma-lowering-fma.f64N/A
PI-lowering-PI.f64N/A
/-lowering-/.f64N/A
Applied egg-rr98.4%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
acos-lowering-acos.f64N/A
div-invN/A
clear-numN/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
/-lowering-/.f6498.5%
Applied egg-rr98.5%
(FPCore (g h) :precision binary64 (* 2.0 (cos (- (+ (* PI 0.6666666666666666) (/ PI 3.0)) (/ (acos (/ g h)) 3.0)))))
double code(double g, double h) {
return 2.0 * cos((((((double) M_PI) * 0.6666666666666666) + (((double) M_PI) / 3.0)) - (acos((g / h)) / 3.0)));
}
public static double code(double g, double h) {
return 2.0 * Math.cos((((Math.PI * 0.6666666666666666) + (Math.PI / 3.0)) - (Math.acos((g / h)) / 3.0)));
}
def code(g, h): return 2.0 * math.cos((((math.pi * 0.6666666666666666) + (math.pi / 3.0)) - (math.acos((g / h)) / 3.0)))
function code(g, h) return Float64(2.0 * cos(Float64(Float64(Float64(pi * 0.6666666666666666) + Float64(pi / 3.0)) - Float64(acos(Float64(g / h)) / 3.0)))) end
function tmp = code(g, h) tmp = 2.0 * cos((((pi * 0.6666666666666666) + (pi / 3.0)) - (acos((g / h)) / 3.0))); end
code[g_, h_] := N[(2.0 * N[Cos[N[(N[(N[(Pi * 0.6666666666666666), $MachinePrecision] + N[(Pi / 3.0), $MachinePrecision]), $MachinePrecision] - N[(N[ArcCos[N[(g / h), $MachinePrecision]], $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \cos \left(\left(\pi \cdot 0.6666666666666666 + \frac{\pi}{3}\right) - \frac{\cos^{-1} \left(\frac{g}{h}\right)}{3}\right)
\end{array}
Initial program 98.4%
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
acos-lowering-acos.f64N/A
distribute-frac-negN/A
distribute-neg-frac2N/A
neg-mul-1N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6498.4%
Simplified98.4%
fma-defineN/A
fma-lowering-fma.f64N/A
PI-lowering-PI.f64N/A
/-lowering-/.f64N/A
Applied egg-rr98.4%
div-invN/A
clear-numN/A
mul-1-negN/A
neg-lowering-neg.f64N/A
/-lowering-/.f6498.4%
Applied egg-rr98.4%
acos-negN/A
div-subN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
/-lowering-/.f64N/A
PI-lowering-PI.f64N/A
/-lowering-/.f64N/A
acos-lowering-acos.f64N/A
/-lowering-/.f6498.5%
Applied egg-rr98.5%
(FPCore (g h) :precision binary64 (* 2.0 (cos (+ (* PI 0.6666666666666666) (/ (acos (/ (/ g -1.0) h)) 3.0)))))
double code(double g, double h) {
return 2.0 * cos(((((double) M_PI) * 0.6666666666666666) + (acos(((g / -1.0) / h)) / 3.0)));
}
public static double code(double g, double h) {
return 2.0 * Math.cos(((Math.PI * 0.6666666666666666) + (Math.acos(((g / -1.0) / h)) / 3.0)));
}
def code(g, h): return 2.0 * math.cos(((math.pi * 0.6666666666666666) + (math.acos(((g / -1.0) / h)) / 3.0)))
function code(g, h) return Float64(2.0 * cos(Float64(Float64(pi * 0.6666666666666666) + Float64(acos(Float64(Float64(g / -1.0) / h)) / 3.0)))) end
function tmp = code(g, h) tmp = 2.0 * cos(((pi * 0.6666666666666666) + (acos(((g / -1.0) / h)) / 3.0))); end
code[g_, h_] := N[(2.0 * N[Cos[N[(N[(Pi * 0.6666666666666666), $MachinePrecision] + N[(N[ArcCos[N[(N[(g / -1.0), $MachinePrecision] / h), $MachinePrecision]], $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \cos \left(\pi \cdot 0.6666666666666666 + \frac{\cos^{-1} \left(\frac{\frac{g}{-1}}{h}\right)}{3}\right)
\end{array}
Initial program 98.4%
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
acos-lowering-acos.f64N/A
distribute-frac-negN/A
distribute-neg-frac2N/A
neg-mul-1N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6498.4%
Simplified98.4%
herbie shell --seed 2024164
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
:precision binary64
(* 2.0 (cos (+ (/ (* 2.0 PI) 3.0) (/ (acos (/ (- g) h)) 3.0)))))