
(FPCore (g h a) :precision binary64 (let* ((t_0 (/ 1.0 (* 2.0 a))) (t_1 (sqrt (- (* g g) (* h h))))) (+ (cbrt (* t_0 (+ (- g) t_1))) (cbrt (* t_0 (- (- g) t_1))))))
double code(double g, double h, double a) {
double t_0 = 1.0 / (2.0 * a);
double t_1 = sqrt(((g * g) - (h * h)));
return cbrt((t_0 * (-g + t_1))) + cbrt((t_0 * (-g - t_1)));
}
public static double code(double g, double h, double a) {
double t_0 = 1.0 / (2.0 * a);
double t_1 = Math.sqrt(((g * g) - (h * h)));
return Math.cbrt((t_0 * (-g + t_1))) + Math.cbrt((t_0 * (-g - t_1)));
}
function code(g, h, a) t_0 = Float64(1.0 / Float64(2.0 * a)) t_1 = sqrt(Float64(Float64(g * g) - Float64(h * h))) return Float64(cbrt(Float64(t_0 * Float64(Float64(-g) + t_1))) + cbrt(Float64(t_0 * Float64(Float64(-g) - t_1)))) end
code[g_, h_, a_] := Block[{t$95$0 = N[(1.0 / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sqrt[N[(N[(g * g), $MachinePrecision] - N[(h * h), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, N[(N[Power[N[(t$95$0 * N[((-g) + t$95$1), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] + N[Power[N[(t$95$0 * N[((-g) - t$95$1), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
t_0 := \frac{1}{2 \cdot a}\\
t_1 := \sqrt{g \cdot g - h \cdot h}\\
\sqrt[3]{t\_0 \cdot \left(\left(-g\right) + t\_1\right)} + \sqrt[3]{t\_0 \cdot \left(\left(-g\right) - t\_1\right)}
\end{array}
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (g h a) :precision binary64 (let* ((t_0 (/ 1.0 (* 2.0 a))) (t_1 (sqrt (- (* g g) (* h h))))) (+ (cbrt (* t_0 (+ (- g) t_1))) (cbrt (* t_0 (- (- g) t_1))))))
double code(double g, double h, double a) {
double t_0 = 1.0 / (2.0 * a);
double t_1 = sqrt(((g * g) - (h * h)));
return cbrt((t_0 * (-g + t_1))) + cbrt((t_0 * (-g - t_1)));
}
public static double code(double g, double h, double a) {
double t_0 = 1.0 / (2.0 * a);
double t_1 = Math.sqrt(((g * g) - (h * h)));
return Math.cbrt((t_0 * (-g + t_1))) + Math.cbrt((t_0 * (-g - t_1)));
}
function code(g, h, a) t_0 = Float64(1.0 / Float64(2.0 * a)) t_1 = sqrt(Float64(Float64(g * g) - Float64(h * h))) return Float64(cbrt(Float64(t_0 * Float64(Float64(-g) + t_1))) + cbrt(Float64(t_0 * Float64(Float64(-g) - t_1)))) end
code[g_, h_, a_] := Block[{t$95$0 = N[(1.0 / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sqrt[N[(N[(g * g), $MachinePrecision] - N[(h * h), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, N[(N[Power[N[(t$95$0 * N[((-g) + t$95$1), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] + N[Power[N[(t$95$0 * N[((-g) - t$95$1), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
t_0 := \frac{1}{2 \cdot a}\\
t_1 := \sqrt{g \cdot g - h \cdot h}\\
\sqrt[3]{t\_0 \cdot \left(\left(-g\right) + t\_1\right)} + \sqrt[3]{t\_0 \cdot \left(\left(-g\right) - t\_1\right)}
\end{array}
(FPCore (g h a) :precision binary64 (* (cbrt (/ 0.5 a)) (cbrt (* -2.0 g))))
double code(double g, double h, double a) {
return cbrt((0.5 / a)) * cbrt((-2.0 * g));
}
public static double code(double g, double h, double a) {
return Math.cbrt((0.5 / a)) * Math.cbrt((-2.0 * g));
}
function code(g, h, a) return Float64(cbrt(Float64(0.5 / a)) * cbrt(Float64(-2.0 * g))) end
code[g_, h_, a_] := N[(N[Power[N[(0.5 / a), $MachinePrecision], 1/3], $MachinePrecision] * N[Power[N[(-2.0 * g), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]
\sqrt[3]{\frac{0.5}{a}} \cdot \sqrt[3]{-2 \cdot g}
Initial program 44.8%
Applied rewrites50.5%
Taylor expanded in g around inf
lower-*.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f6495.0
Applied rewrites95.0%
lift-/.f64N/A
mult-flip-revN/A
metadata-evalN/A
lift-cbrt.f64N/A
lift-+.f64N/A
count-2-revN/A
lift-*.f64N/A
cbrt-undivN/A
lift-/.f64N/A
lower-*.f64N/A
Applied rewrites95.7%
Applied rewrites95.7%
(FPCore (g h a) :precision binary64 (* (cbrt g) (cbrt (/ -1.0 a))))
double code(double g, double h, double a) {
return cbrt(g) * cbrt((-1.0 / a));
}
public static double code(double g, double h, double a) {
return Math.cbrt(g) * Math.cbrt((-1.0 / a));
}
function code(g, h, a) return Float64(cbrt(g) * cbrt(Float64(-1.0 / a))) end
code[g_, h_, a_] := N[(N[Power[g, 1/3], $MachinePrecision] * N[Power[N[(-1.0 / a), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]
\sqrt[3]{g} \cdot \sqrt[3]{\frac{-1}{a}}
Initial program 44.8%
Taylor expanded in g around inf
lower-/.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f6495.1
Applied rewrites95.1%
lift-/.f64N/A
div-flipN/A
lower-unsound-/.f64N/A
lower-unsound-/.f6495.1
lift-*.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
lift-*.f64N/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
cbrt-unprodN/A
metadata-evalN/A
pow1/3N/A
pow-prod-downN/A
*-commutativeN/A
mul-1-negN/A
lift-neg.f64N/A
pow1/3N/A
lift-neg.f64N/A
mul-1-negN/A
*-commutativeN/A
*-commutativeN/A
mul-1-negN/A
cbrt-neg-revN/A
lift-cbrt.f64N/A
lower-neg.f6495.7
Applied rewrites95.7%
lift-/.f64N/A
lift-/.f64N/A
div-flip-revN/A
lift-neg.f64N/A
distribute-neg-fracN/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
cbrt-divN/A
lift-/.f64N/A
cbrt-neg-revN/A
pow1/3N/A
lift-/.f64N/A
mult-flipN/A
lift-/.f64N/A
distribute-rgt-neg-inN/A
unpow-prod-downN/A
lower-unsound-pow.f32N/A
lower-pow.f32N/A
pow1/3N/A
lift-cbrt.f64N/A
lower-unsound-*.f64N/A
lower-unsound-pow.f64N/A
lower-neg.f6445.0
Applied rewrites45.0%
lift-pow.f64N/A
unpow1/3N/A
lower-cbrt.f6495.7
lift-neg.f64N/A
lift-/.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
lower-/.f6495.7
Applied rewrites95.7%
(FPCore (g h a) :precision binary64 (/ (cbrt g) (- (cbrt a))))
double code(double g, double h, double a) {
return cbrt(g) / -cbrt(a);
}
public static double code(double g, double h, double a) {
return Math.cbrt(g) / -Math.cbrt(a);
}
function code(g, h, a) return Float64(cbrt(g) / Float64(-cbrt(a))) end
code[g_, h_, a_] := N[(N[Power[g, 1/3], $MachinePrecision] / (-N[Power[a, 1/3], $MachinePrecision])), $MachinePrecision]
\frac{\sqrt[3]{g}}{-\sqrt[3]{a}}
Initial program 44.8%
Taylor expanded in g around inf
lower-/.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f6495.1
Applied rewrites95.1%
lift-/.f64N/A
div-flipN/A
lower-unsound-/.f64N/A
lower-unsound-/.f6495.1
lift-*.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
lift-*.f64N/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
cbrt-unprodN/A
metadata-evalN/A
pow1/3N/A
pow-prod-downN/A
*-commutativeN/A
mul-1-negN/A
lift-neg.f64N/A
pow1/3N/A
lift-neg.f64N/A
mul-1-negN/A
*-commutativeN/A
*-commutativeN/A
mul-1-negN/A
cbrt-neg-revN/A
lift-cbrt.f64N/A
lower-neg.f6495.7
Applied rewrites95.7%
lift-/.f64N/A
lift-/.f64N/A
div-flip-revN/A
lift-neg.f64N/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
lower-/.f64N/A
lower-neg.f6495.7
Applied rewrites95.7%
(FPCore (g h a) :precision binary64 (/ 1.0 (cbrt (* (/ -1.0 g) a))))
double code(double g, double h, double a) {
return 1.0 / cbrt(((-1.0 / g) * a));
}
public static double code(double g, double h, double a) {
return 1.0 / Math.cbrt(((-1.0 / g) * a));
}
function code(g, h, a) return Float64(1.0 / cbrt(Float64(Float64(-1.0 / g) * a))) end
code[g_, h_, a_] := N[(1.0 / N[Power[N[(N[(-1.0 / g), $MachinePrecision] * a), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]
\frac{1}{\sqrt[3]{\frac{-1}{g} \cdot a}}
Initial program 44.8%
Taylor expanded in g around inf
lower-/.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f6495.1
Applied rewrites95.1%
lift-/.f64N/A
div-flipN/A
lower-unsound-/.f64N/A
lower-unsound-/.f6495.1
lift-*.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
lift-*.f64N/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
cbrt-unprodN/A
metadata-evalN/A
pow1/3N/A
pow-prod-downN/A
*-commutativeN/A
mul-1-negN/A
lift-neg.f64N/A
pow1/3N/A
lift-neg.f64N/A
mul-1-negN/A
*-commutativeN/A
*-commutativeN/A
mul-1-negN/A
cbrt-neg-revN/A
lift-cbrt.f64N/A
lower-neg.f6495.7
Applied rewrites95.7%
lift-/.f64N/A
mult-flipN/A
*-commutativeN/A
metadata-evalN/A
lift-neg.f64N/A
lift-cbrt.f64N/A
cbrt-neg-revN/A
lift-neg.f64N/A
cbrt-undivN/A
lift-cbrt.f64N/A
cbrt-unprodN/A
lower-cbrt.f64N/A
lower-*.f64N/A
metadata-evalN/A
lift-neg.f64N/A
frac-2neg-revN/A
lower-/.f6473.8
Applied rewrites73.8%
(FPCore (g h a) :precision binary64 (/ -1.0 (cbrt (/ a g))))
double code(double g, double h, double a) {
return -1.0 / cbrt((a / g));
}
public static double code(double g, double h, double a) {
return -1.0 / Math.cbrt((a / g));
}
function code(g, h, a) return Float64(-1.0 / cbrt(Float64(a / g))) end
code[g_, h_, a_] := N[(-1.0 / N[Power[N[(a / g), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]
\frac{-1}{\sqrt[3]{\frac{a}{g}}}
Initial program 44.8%
Taylor expanded in g around inf
lower-/.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f6495.1
Applied rewrites95.1%
lift-/.f64N/A
div-flipN/A
lower-unsound-/.f64N/A
lower-unsound-/.f6495.1
lift-*.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
lift-*.f64N/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
cbrt-unprodN/A
metadata-evalN/A
pow1/3N/A
pow-prod-downN/A
*-commutativeN/A
mul-1-negN/A
lift-neg.f64N/A
pow1/3N/A
lift-neg.f64N/A
mul-1-negN/A
*-commutativeN/A
*-commutativeN/A
mul-1-negN/A
cbrt-neg-revN/A
lift-cbrt.f64N/A
lower-neg.f6495.7
Applied rewrites95.7%
lift-/.f64N/A
frac-2negN/A
metadata-evalN/A
lower-/.f64N/A
lift-/.f64N/A
distribute-neg-fracN/A
lift-neg.f64N/A
frac-2negN/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
cbrt-undivN/A
lower-cbrt.f64N/A
lower-/.f6473.8
Applied rewrites73.8%
(FPCore (g h a) :precision binary64 (cbrt (* (/ -1.0 a) g)))
double code(double g, double h, double a) {
return cbrt(((-1.0 / a) * g));
}
public static double code(double g, double h, double a) {
return Math.cbrt(((-1.0 / a) * g));
}
function code(g, h, a) return cbrt(Float64(Float64(-1.0 / a) * g)) end
code[g_, h_, a_] := N[Power[N[(N[(-1.0 / a), $MachinePrecision] * g), $MachinePrecision], 1/3], $MachinePrecision]
\sqrt[3]{\frac{-1}{a} \cdot g}
Initial program 44.8%
Taylor expanded in g around inf
lower-/.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f6495.1
Applied rewrites95.1%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lift-*.f64N/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
cbrt-unprodN/A
metadata-evalN/A
lift-cbrt.f64N/A
cbrt-undivN/A
lift-cbrt.f64N/A
cbrt-unprodN/A
lower-cbrt.f64N/A
lower-*.f64N/A
lower-/.f6473.0
Applied rewrites73.0%
(FPCore (g h a) :precision binary64 (- (cbrt (/ g a))))
double code(double g, double h, double a) {
return -cbrt((g / a));
}
public static double code(double g, double h, double a) {
return -Math.cbrt((g / a));
}
function code(g, h, a) return Float64(-cbrt(Float64(g / a))) end
code[g_, h_, a_] := (-N[Power[N[(g / a), $MachinePrecision], 1/3], $MachinePrecision])
-\sqrt[3]{\frac{g}{a}}
Initial program 44.8%
Taylor expanded in g around inf
lower-/.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f64N/A
lower-cbrt.f6495.1
Applied rewrites95.1%
lift-/.f64N/A
frac-2negN/A
distribute-frac-negN/A
lift-*.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
lift-*.f64N/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
cbrt-unprodN/A
metadata-evalN/A
pow1/3N/A
pow-prod-downN/A
*-commutativeN/A
mul-1-negN/A
lift-neg.f64N/A
pow1/3N/A
lift-neg.f64N/A
mul-1-negN/A
*-commutativeN/A
*-commutativeN/A
mul-1-negN/A
cbrt-neg-revN/A
lift-cbrt.f64N/A
Applied rewrites73.1%
herbie shell --seed 2025180
(FPCore (g h a)
:name "2-ancestry mixing, positive discriminant"
:precision binary64
(+ (cbrt (* (/ 1.0 (* 2.0 a)) (+ (- g) (sqrt (- (* g g) (* h h)))))) (cbrt (* (/ 1.0 (* 2.0 a)) (- (- g) (sqrt (- (* g g) (* h h))))))))