2-ancestry mixing, positive discriminant

Percentage Accurate: 45.0% → 93.8%

Time: 12.6s

Alternatives: 9

Speedup: 1.4×

Specification

Math

\[\begin{array}{l} \\ \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} \end{array} \]

FPCore

(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))))))

C

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)));
}

Java

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)));
}

Julia

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

Wolfram

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]]]

Initial Program: 45.0% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \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} \end{array} \]

FPCore

(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))))))

C

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)));
}

Java

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)));
}

Julia

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

Wolfram

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]]]

Alternative 1: 93.8% accurate, 0.5× speedup

Math

\[\begin{array}{l} \\ \mathsf{fma}\left(\frac{\sqrt[3]{g \cdot -0.5}}{\sqrt[3]{a}}, {2}^{0.3333333333333333}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right) \end{array} \]

FPCore

(FPCore (g h a)
 :precision binary64
 (fma
  (/ (cbrt (* g -0.5)) (cbrt a))
  (pow 2.0 0.3333333333333333)
  (* (cbrt (* (/ h g) (/ h a))) (* (cbrt 0.5) (cbrt -0.5)))))

C

double code(double g, double h, double a) {
	return fma((cbrt((g * -0.5)) / cbrt(a)), pow(2.0, 0.3333333333333333), (cbrt(((h / g) * (h / a))) * (cbrt(0.5) * cbrt(-0.5))));
}

Julia

function code(g, h, a)
	return fma(Float64(cbrt(Float64(g * -0.5)) / cbrt(a)), (2.0 ^ 0.3333333333333333), Float64(cbrt(Float64(Float64(h / g) * Float64(h / a))) * Float64(cbrt(0.5) * cbrt(-0.5))))
end

Wolfram

code[g_, h_, a_] := N[(N[(N[Power[N[(g * -0.5), $MachinePrecision], 1/3], $MachinePrecision] / N[Power[a, 1/3], $MachinePrecision]), $MachinePrecision] * N[Power[2.0, 0.3333333333333333], $MachinePrecision] + N[(N[Power[N[(N[(h / g), $MachinePrecision] * N[(h / a), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * N[(N[Power[0.5, 1/3], $MachinePrecision] * N[Power[-0.5, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 43.7%

   \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
2. Add Preprocessing

3. Taylor expanded in h around 0

   \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{2}\right) + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)} \]
4. Step-by-step derivation

   1. associate-*r* N/A

      \[\leadsto \color{blue}{\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}\right) \cdot \sqrt[3]{2}} + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right) \]

   2. lower-fma.f64 N/A

      \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right)} \]

   3. lower-*.f64 N/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   4. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   5. lower-/.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   6. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   7. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \color{blue}{\sqrt[3]{2}}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   8. lower-*.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)}\right) \]

   9. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   10. unpow2 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{\color{blue}{h \cdot h}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   11. *-commutative N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h \cdot h}{\color{blue}{g \cdot a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   12. times-frac N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   13. lower-*.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   14. lower-/.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g}} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   15. lower-/.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \color{blue}{\frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   16. *-commutative N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

   17. lower-*.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

   18. lower-cbrt.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\color{blue}{\sqrt[3]{\frac{1}{2}}} \cdot \sqrt[3]{\frac{-1}{2}}\right)\right) \]

   19. lower-cbrt.f64 71.9

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \color{blue}{\sqrt[3]{-0.5}}\right)\right) \]

5. Applied rewrites 71.9%

   \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right)} \]
6. Step-by-step derivation

7. Applied rewrites 89.7%

   \[\leadsto \mathsf{fma}\left(\frac{\sqrt[3]{g}}{\sqrt[3]{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right) \]
8. Step-by-step derivation

9. Applied rewrites 90.4%

   \[\leadsto \mathsf{fma}\left(\frac{\sqrt[3]{g}}{\sqrt[3]{a}} \cdot \sqrt[3]{-0.5}, {2}^{\color{blue}{0.3333333333333333}}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right) \]
10. Step-by-step derivation

11. Applied rewrites 90.6%

    \[\leadsto \mathsf{fma}\left(\frac{\sqrt[3]{g \cdot -0.5}}{\sqrt[3]{a}}, {\color{blue}{2}}^{0.3333333333333333}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right) \]
12. Add Preprocessing

Alternative 2: 93.3% accurate, 0.5× speedup

Math

\[\begin{array}{l} \\ \mathsf{fma}\left(\frac{\sqrt[3]{g \cdot -0.5}}{\sqrt[3]{a}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right) \end{array} \]

FPCore

(FPCore (g h a)
 :precision binary64
 (fma
  (/ (cbrt (* g -0.5)) (cbrt a))
  (cbrt 2.0)
  (* (cbrt (* (/ h g) (/ h a))) (* (cbrt 0.5) (cbrt -0.5)))))

C

double code(double g, double h, double a) {
	return fma((cbrt((g * -0.5)) / cbrt(a)), cbrt(2.0), (cbrt(((h / g) * (h / a))) * (cbrt(0.5) * cbrt(-0.5))));
}

Julia

function code(g, h, a)
	return fma(Float64(cbrt(Float64(g * -0.5)) / cbrt(a)), cbrt(2.0), Float64(cbrt(Float64(Float64(h / g) * Float64(h / a))) * Float64(cbrt(0.5) * cbrt(-0.5))))
end

Wolfram

code[g_, h_, a_] := N[(N[(N[Power[N[(g * -0.5), $MachinePrecision], 1/3], $MachinePrecision] / N[Power[a, 1/3], $MachinePrecision]), $MachinePrecision] * N[Power[2.0, 1/3], $MachinePrecision] + N[(N[Power[N[(N[(h / g), $MachinePrecision] * N[(h / a), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * N[(N[Power[0.5, 1/3], $MachinePrecision] * N[Power[-0.5, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 43.7%

   \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
2. Add Preprocessing

3. Taylor expanded in h around 0

   \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{2}\right) + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)} \]
4. Step-by-step derivation

   1. associate-*r* N/A

      \[\leadsto \color{blue}{\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}\right) \cdot \sqrt[3]{2}} + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right) \]

   2. lower-fma.f64 N/A

      \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right)} \]

   3. lower-*.f64 N/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   4. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   5. lower-/.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   6. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   7. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \color{blue}{\sqrt[3]{2}}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   8. lower-*.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)}\right) \]

   9. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   10. unpow2 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{\color{blue}{h \cdot h}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   11. *-commutative N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h \cdot h}{\color{blue}{g \cdot a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   12. times-frac N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   13. lower-*.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   14. lower-/.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g}} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   15. lower-/.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \color{blue}{\frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   16. *-commutative N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

   17. lower-*.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

   18. lower-cbrt.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\color{blue}{\sqrt[3]{\frac{1}{2}}} \cdot \sqrt[3]{\frac{-1}{2}}\right)\right) \]

   19. lower-cbrt.f64 71.9

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \color{blue}{\sqrt[3]{-0.5}}\right)\right) \]

5. Applied rewrites 71.9%

   \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right)} \]
6. Step-by-step derivation

7. Applied rewrites 89.7%

   \[\leadsto \mathsf{fma}\left(\frac{\sqrt[3]{g}}{\sqrt[3]{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right) \]
8. Step-by-step derivation

9. Applied rewrites 89.9%

   \[\leadsto \mathsf{fma}\left(\frac{\sqrt[3]{g \cdot -0.5}}{\sqrt[3]{a}}, \sqrt[3]{\color{blue}{2}}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right) \]
10. Add Preprocessing

Alternative 3: 76.8% accurate, 0.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;g \leq 9.5 \cdot 10^{-198}:\\ \;\;\;\;\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \sqrt[3]{-1}, \sqrt[3]{\frac{h}{a} \cdot \frac{h}{g}} \cdot \sqrt[3]{-0.25}\right)\\ \mathbf{elif}\;g \leq 1.65 \cdot 10^{+148}:\\ \;\;\;\;\sqrt[3]{{\left(2 \cdot a\right)}^{-1} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot 2}}{\sqrt[3]{a \cdot 2}}\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\sqrt[3]{\frac{g}{a} \cdot -0.5}, {2}^{0.3333333333333333}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (g h a)
 :precision binary64
 (if (<= g 9.5e-198)
   (fma (cbrt (/ g a)) (cbrt -1.0) (* (cbrt (* (/ h a) (/ h g))) (cbrt -0.25)))
   (if (<= g 1.65e+148)
     (+
      (cbrt (* (pow (* 2.0 a) -1.0) (+ (- g) (sqrt (- (* g g) (* h h))))))
      (/ (- (cbrt (* g 2.0))) (cbrt (* a 2.0))))
     (fma
      (cbrt (* (/ g a) -0.5))
      (pow 2.0 0.3333333333333333)
      (* (cbrt (* (/ h g) (/ h a))) (* (cbrt 0.5) (cbrt -0.5)))))))

C

double code(double g, double h, double a) {
	double tmp;
	if (g <= 9.5e-198) {
		tmp = fma(cbrt((g / a)), cbrt(-1.0), (cbrt(((h / a) * (h / g))) * cbrt(-0.25)));
	} else if (g <= 1.65e+148) {
		tmp = cbrt((pow((2.0 * a), -1.0) * (-g + sqrt(((g * g) - (h * h)))))) + (-cbrt((g * 2.0)) / cbrt((a * 2.0)));
	} else {
		tmp = fma(cbrt(((g / a) * -0.5)), pow(2.0, 0.3333333333333333), (cbrt(((h / g) * (h / a))) * (cbrt(0.5) * cbrt(-0.5))));
	}
	return tmp;
}

Julia

function code(g, h, a)
	tmp = 0.0
	if (g <= 9.5e-198)
		tmp = fma(cbrt(Float64(g / a)), cbrt(-1.0), Float64(cbrt(Float64(Float64(h / a) * Float64(h / g))) * cbrt(-0.25)));
	elseif (g <= 1.65e+148)
		tmp = Float64(cbrt(Float64((Float64(2.0 * a) ^ -1.0) * Float64(Float64(-g) + sqrt(Float64(Float64(g * g) - Float64(h * h)))))) + Float64(Float64(-cbrt(Float64(g * 2.0))) / cbrt(Float64(a * 2.0))));
	else
		tmp = fma(cbrt(Float64(Float64(g / a) * -0.5)), (2.0 ^ 0.3333333333333333), Float64(cbrt(Float64(Float64(h / g) * Float64(h / a))) * Float64(cbrt(0.5) * cbrt(-0.5))));
	end
	return tmp
end

Wolfram

code[g_, h_, a_] := If[LessEqual[g, 9.5e-198], N[(N[Power[N[(g / a), $MachinePrecision], 1/3], $MachinePrecision] * N[Power[-1.0, 1/3], $MachinePrecision] + N[(N[Power[N[(N[(h / a), $MachinePrecision] * N[(h / g), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * N[Power[-0.25, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[g, 1.65e+148], N[(N[Power[N[(N[Power[N[(2.0 * a), $MachinePrecision], -1.0], $MachinePrecision] * N[((-g) + N[Sqrt[N[(N[(g * g), $MachinePrecision] - N[(h * h), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] + N[((-N[Power[N[(g * 2.0), $MachinePrecision], 1/3], $MachinePrecision]) / N[Power[N[(a * 2.0), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[N[(N[(g / a), $MachinePrecision] * -0.5), $MachinePrecision], 1/3], $MachinePrecision] * N[Power[2.0, 0.3333333333333333], $MachinePrecision] + N[(N[Power[N[(N[(h / g), $MachinePrecision] * N[(h / a), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * N[(N[Power[0.5, 1/3], $MachinePrecision] * N[Power[-0.5, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if g < 9.4999999999999997e-198

   1. Initial program 41.9%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing

   3. Taylor expanded in h around 0

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{2}\right) + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)} \]
   4. Step-by-step derivation

      1. associate-*r* N/A

         \[\leadsto \color{blue}{\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}\right) \cdot \sqrt[3]{2}} + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right) \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right)} \]

      3. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      4. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      5. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      6. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      7. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \color{blue}{\sqrt[3]{2}}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      8. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)}\right) \]

      9. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      10. unpow2 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{\color{blue}{h \cdot h}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      11. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h \cdot h}{\color{blue}{g \cdot a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      12. times-frac N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      13. lower-*.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      14. lower-/.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g}} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      15. lower-/.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \color{blue}{\frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      16. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

      17. lower-*.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

      18. lower-cbrt.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\color{blue}{\sqrt[3]{\frac{1}{2}}} \cdot \sqrt[3]{\frac{-1}{2}}\right)\right) \]

      19. lower-cbrt.f64 71.3

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \color{blue}{\sqrt[3]{-0.5}}\right)\right) \]

   5. Applied rewrites 71.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 72.0%

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \color{blue}{\sqrt[3]{-1}}, \left(-\left(-\sqrt[3]{\frac{h}{a} \cdot \frac{h}{g}}\right)\right) \cdot \sqrt[3]{-0.25}\right) \]

   if 9.4999999999999997e-198 < g < 1.65000000000000005e148

   1. Initial program 82.6%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      2. lift-*.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      3. lift-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a}} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

      4. associate-*l/ N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}{2 \cdot a}}} \]

      5. cbrt-div N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}}{\sqrt[3]{2 \cdot a}}} \]

      6. *-lft-identity N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\sqrt[3]{\color{blue}{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}}{\sqrt[3]{2 \cdot a}} \]

      7. lower-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}{\sqrt[3]{2 \cdot a}}} \]

   4. Applied rewrites 95.7%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{\left(g - h\right) \cdot \left(h + g\right)}}}{\sqrt[3]{a \cdot 2}}} \]

   5. Taylor expanded in g around -inf

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{-1 \cdot \sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}}}{\sqrt[3]{a \cdot 2}} \]
   6. Step-by-step derivation

      1. mul-1-neg N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{\mathsf{neg}\left(\sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}\right)}}{\sqrt[3]{a \cdot 2}} \]

      2. lower-neg.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{-\sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}}}{\sqrt[3]{a \cdot 2}} \]

      3. unpow2 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot \left(1 - \color{blue}{\sqrt{-1} \cdot \sqrt{-1}}\right)}}{\sqrt[3]{a \cdot 2}} \]

      4. rem-square-sqrt N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot \left(1 - \color{blue}{-1}\right)}}{\sqrt[3]{a \cdot 2}} \]

      5. metadata-eval N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot \color{blue}{2}}}{\sqrt[3]{a \cdot 2}} \]

      6. lower-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\color{blue}{\sqrt[3]{g \cdot 2}}}{\sqrt[3]{a \cdot 2}} \]

      7. lower-*.f64 98.1

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{\color{blue}{g \cdot 2}}}{\sqrt[3]{a \cdot 2}} \]

   7. Applied rewrites 98.1%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{-\sqrt[3]{g \cdot 2}}}{\sqrt[3]{a \cdot 2}} \]

   if 1.65000000000000005e148 < g

   1. Initial program 3.8%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing

   3. Taylor expanded in h around 0

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{2}\right) + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)} \]
   4. Step-by-step derivation

      1. associate-*r* N/A

         \[\leadsto \color{blue}{\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}\right) \cdot \sqrt[3]{2}} + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right) \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right)} \]

      3. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      4. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      5. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      6. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      7. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \color{blue}{\sqrt[3]{2}}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      8. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)}\right) \]

      9. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      10. unpow2 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{\color{blue}{h \cdot h}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      11. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h \cdot h}{\color{blue}{g \cdot a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      12. times-frac N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      13. lower-*.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      14. lower-/.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g}} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      15. lower-/.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \color{blue}{\frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      16. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

      17. lower-*.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

      18. lower-cbrt.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\color{blue}{\sqrt[3]{\frac{1}{2}}} \cdot \sqrt[3]{\frac{-1}{2}}\right)\right) \]

      19. lower-cbrt.f64 62.2

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \color{blue}{\sqrt[3]{-0.5}}\right)\right) \]

   5. Applied rewrites 62.2%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 86.9%

      \[\leadsto \mathsf{fma}\left(\frac{\sqrt[3]{g}}{\sqrt[3]{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 87.6%

      \[\leadsto \mathsf{fma}\left(\frac{\sqrt[3]{g}}{\sqrt[3]{a}} \cdot \sqrt[3]{-0.5}, {2}^{\color{blue}{0.3333333333333333}}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right) \]
   10. Step-by-step derivation

   11. Applied rewrites 62.7%

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a} \cdot -0.5}, {\color{blue}{2}}^{0.3333333333333333}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 76.2%

   \[\leadsto \begin{array}{l} \mathbf{if}\;g \leq 9.5 \cdot 10^{-198}:\\ \;\;\;\;\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \sqrt[3]{-1}, \sqrt[3]{\frac{h}{a} \cdot \frac{h}{g}} \cdot \sqrt[3]{-0.25}\right)\\ \mathbf{elif}\;g \leq 1.65 \cdot 10^{+148}:\\ \;\;\;\;\sqrt[3]{{\left(2 \cdot a\right)}^{-1} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot 2}}{\sqrt[3]{a \cdot 2}}\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\sqrt[3]{\frac{g}{a} \cdot -0.5}, {2}^{0.3333333333333333}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 4: 76.6% accurate, 0.6× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \sqrt[3]{\frac{g}{a}}\\ \mathbf{if}\;g \leq 9.5 \cdot 10^{-198}:\\ \;\;\;\;\mathsf{fma}\left(t\_0, \sqrt[3]{-1}, \sqrt[3]{\frac{h}{a} \cdot \frac{h}{g}} \cdot \sqrt[3]{-0.25}\right)\\ \mathbf{elif}\;g \leq 1.65 \cdot 10^{+148}:\\ \;\;\;\;\sqrt[3]{{\left(2 \cdot a\right)}^{-1} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot 2}}{\sqrt[3]{a \cdot 2}}\\ \mathbf{else}:\\ \;\;\;\;t\_0 \cdot \sqrt[3]{-1}\\ \end{array} \end{array} \]

FPCore

(FPCore (g h a)
 :precision binary64
 (let* ((t_0 (cbrt (/ g a))))
   (if (<= g 9.5e-198)
     (fma t_0 (cbrt -1.0) (* (cbrt (* (/ h a) (/ h g))) (cbrt -0.25)))
     (if (<= g 1.65e+148)
       (+
        (cbrt (* (pow (* 2.0 a) -1.0) (+ (- g) (sqrt (- (* g g) (* h h))))))
        (/ (- (cbrt (* g 2.0))) (cbrt (* a 2.0))))
       (* t_0 (cbrt -1.0))))))

C

double code(double g, double h, double a) {
	double t_0 = cbrt((g / a));
	double tmp;
	if (g <= 9.5e-198) {
		tmp = fma(t_0, cbrt(-1.0), (cbrt(((h / a) * (h / g))) * cbrt(-0.25)));
	} else if (g <= 1.65e+148) {
		tmp = cbrt((pow((2.0 * a), -1.0) * (-g + sqrt(((g * g) - (h * h)))))) + (-cbrt((g * 2.0)) / cbrt((a * 2.0)));
	} else {
		tmp = t_0 * cbrt(-1.0);
	}
	return tmp;
}

Julia

function code(g, h, a)
	t_0 = cbrt(Float64(g / a))
	tmp = 0.0
	if (g <= 9.5e-198)
		tmp = fma(t_0, cbrt(-1.0), Float64(cbrt(Float64(Float64(h / a) * Float64(h / g))) * cbrt(-0.25)));
	elseif (g <= 1.65e+148)
		tmp = Float64(cbrt(Float64((Float64(2.0 * a) ^ -1.0) * Float64(Float64(-g) + sqrt(Float64(Float64(g * g) - Float64(h * h)))))) + Float64(Float64(-cbrt(Float64(g * 2.0))) / cbrt(Float64(a * 2.0))));
	else
		tmp = Float64(t_0 * cbrt(-1.0));
	end
	return tmp
end

Wolfram

code[g_, h_, a_] := Block[{t$95$0 = N[Power[N[(g / a), $MachinePrecision], 1/3], $MachinePrecision]}, If[LessEqual[g, 9.5e-198], N[(t$95$0 * N[Power[-1.0, 1/3], $MachinePrecision] + N[(N[Power[N[(N[(h / a), $MachinePrecision] * N[(h / g), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * N[Power[-0.25, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[g, 1.65e+148], N[(N[Power[N[(N[Power[N[(2.0 * a), $MachinePrecision], -1.0], $MachinePrecision] * N[((-g) + N[Sqrt[N[(N[(g * g), $MachinePrecision] - N[(h * h), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] + N[((-N[Power[N[(g * 2.0), $MachinePrecision], 1/3], $MachinePrecision]) / N[Power[N[(a * 2.0), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[Power[-1.0, 1/3], $MachinePrecision]), $MachinePrecision]]]]

Derivation

1. Split input into 3 regimes

2. if g < 9.4999999999999997e-198

   1. Initial program 41.9%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing

   3. Taylor expanded in h around 0

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{2}\right) + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)} \]
   4. Step-by-step derivation

      1. associate-*r* N/A

         \[\leadsto \color{blue}{\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}\right) \cdot \sqrt[3]{2}} + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right) \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right)} \]

      3. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      4. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      5. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      6. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      7. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \color{blue}{\sqrt[3]{2}}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      8. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)}\right) \]

      9. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      10. unpow2 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{\color{blue}{h \cdot h}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      11. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h \cdot h}{\color{blue}{g \cdot a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      12. times-frac N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      13. lower-*.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      14. lower-/.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g}} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      15. lower-/.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \color{blue}{\frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      16. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

      17. lower-*.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

      18. lower-cbrt.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\color{blue}{\sqrt[3]{\frac{1}{2}}} \cdot \sqrt[3]{\frac{-1}{2}}\right)\right) \]

      19. lower-cbrt.f64 71.3

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \color{blue}{\sqrt[3]{-0.5}}\right)\right) \]

   5. Applied rewrites 71.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 72.0%

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \color{blue}{\sqrt[3]{-1}}, \left(-\left(-\sqrt[3]{\frac{h}{a} \cdot \frac{h}{g}}\right)\right) \cdot \sqrt[3]{-0.25}\right) \]

   if 9.4999999999999997e-198 < g < 1.65000000000000005e148

   1. Initial program 82.6%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      2. lift-*.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      3. lift-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a}} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

      4. associate-*l/ N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}{2 \cdot a}}} \]

      5. cbrt-div N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}}{\sqrt[3]{2 \cdot a}}} \]

      6. *-lft-identity N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\sqrt[3]{\color{blue}{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}}{\sqrt[3]{2 \cdot a}} \]

      7. lower-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}{\sqrt[3]{2 \cdot a}}} \]

   4. Applied rewrites 95.7%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{\left(g - h\right) \cdot \left(h + g\right)}}}{\sqrt[3]{a \cdot 2}}} \]

   5. Taylor expanded in g around -inf

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{-1 \cdot \sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}}}{\sqrt[3]{a \cdot 2}} \]
   6. Step-by-step derivation

      1. mul-1-neg N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{\mathsf{neg}\left(\sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}\right)}}{\sqrt[3]{a \cdot 2}} \]

      2. lower-neg.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{-\sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}}}{\sqrt[3]{a \cdot 2}} \]

      3. unpow2 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot \left(1 - \color{blue}{\sqrt{-1} \cdot \sqrt{-1}}\right)}}{\sqrt[3]{a \cdot 2}} \]

      4. rem-square-sqrt N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot \left(1 - \color{blue}{-1}\right)}}{\sqrt[3]{a \cdot 2}} \]

      5. metadata-eval N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot \color{blue}{2}}}{\sqrt[3]{a \cdot 2}} \]

      6. lower-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\color{blue}{\sqrt[3]{g \cdot 2}}}{\sqrt[3]{a \cdot 2}} \]

      7. lower-*.f64 98.1

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{\color{blue}{g \cdot 2}}}{\sqrt[3]{a \cdot 2}} \]

   7. Applied rewrites 98.1%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{-\sqrt[3]{g \cdot 2}}}{\sqrt[3]{a \cdot 2}} \]

   if 1.65000000000000005e148 < g

   1. Initial program 3.8%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      2. lift-*.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      3. lift-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a}} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

      4. associate-*l/ N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}{2 \cdot a}}} \]

      5. cbrt-div N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}}{\sqrt[3]{2 \cdot a}}} \]

      6. *-lft-identity N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\sqrt[3]{\color{blue}{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}}{\sqrt[3]{2 \cdot a}} \]

      7. lower-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}{\sqrt[3]{2 \cdot a}}} \]

   4. Applied rewrites 3.8%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{\left(g - h\right) \cdot \left(h + g\right)}}}{\sqrt[3]{a \cdot 2}}} \]

   5. Taylor expanded in g around inf

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]
   6. Step-by-step derivation

      1. lower-*.f64 N/A

         \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]

      2. lower-cbrt.f64 N/A

         \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{-1} \]

      3. lower-/.f64 N/A

         \[\leadsto \sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{-1} \]

      4. lower-cbrt.f64 62.6

         \[\leadsto \sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{-1}} \]

   7. Applied rewrites 62.6%

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]
3. Recombined 3 regimes into one program.

4. Final simplification 76.2%

   \[\leadsto \begin{array}{l} \mathbf{if}\;g \leq 9.5 \cdot 10^{-198}:\\ \;\;\;\;\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \sqrt[3]{-1}, \sqrt[3]{\frac{h}{a} \cdot \frac{h}{g}} \cdot \sqrt[3]{-0.25}\right)\\ \mathbf{elif}\;g \leq 1.65 \cdot 10^{+148}:\\ \;\;\;\;\sqrt[3]{{\left(2 \cdot a\right)}^{-1} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot 2}}{\sqrt[3]{a \cdot 2}}\\ \mathbf{else}:\\ \;\;\;\;\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}\\ \end{array} \]
5. Add Preprocessing

Alternative 5: 76.6% accurate, 0.6× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \sqrt[3]{\frac{g}{a}}\\ \mathbf{if}\;g \leq 9.5 \cdot 10^{-198}:\\ \;\;\;\;\mathsf{fma}\left(t\_0, \sqrt[3]{-1}, \sqrt[3]{-0.25 \cdot \left(\frac{h}{a} \cdot \frac{h}{g}\right)}\right)\\ \mathbf{elif}\;g \leq 1.65 \cdot 10^{+148}:\\ \;\;\;\;\sqrt[3]{{\left(2 \cdot a\right)}^{-1} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot 2}}{\sqrt[3]{a \cdot 2}}\\ \mathbf{else}:\\ \;\;\;\;t\_0 \cdot \sqrt[3]{-1}\\ \end{array} \end{array} \]

FPCore

(FPCore (g h a)
 :precision binary64
 (let* ((t_0 (cbrt (/ g a))))
   (if (<= g 9.5e-198)
     (fma t_0 (cbrt -1.0) (cbrt (* -0.25 (* (/ h a) (/ h g)))))
     (if (<= g 1.65e+148)
       (+
        (cbrt (* (pow (* 2.0 a) -1.0) (+ (- g) (sqrt (- (* g g) (* h h))))))
        (/ (- (cbrt (* g 2.0))) (cbrt (* a 2.0))))
       (* t_0 (cbrt -1.0))))))

C

double code(double g, double h, double a) {
	double t_0 = cbrt((g / a));
	double tmp;
	if (g <= 9.5e-198) {
		tmp = fma(t_0, cbrt(-1.0), cbrt((-0.25 * ((h / a) * (h / g)))));
	} else if (g <= 1.65e+148) {
		tmp = cbrt((pow((2.0 * a), -1.0) * (-g + sqrt(((g * g) - (h * h)))))) + (-cbrt((g * 2.0)) / cbrt((a * 2.0)));
	} else {
		tmp = t_0 * cbrt(-1.0);
	}
	return tmp;
}

Julia

function code(g, h, a)
	t_0 = cbrt(Float64(g / a))
	tmp = 0.0
	if (g <= 9.5e-198)
		tmp = fma(t_0, cbrt(-1.0), cbrt(Float64(-0.25 * Float64(Float64(h / a) * Float64(h / g)))));
	elseif (g <= 1.65e+148)
		tmp = Float64(cbrt(Float64((Float64(2.0 * a) ^ -1.0) * Float64(Float64(-g) + sqrt(Float64(Float64(g * g) - Float64(h * h)))))) + Float64(Float64(-cbrt(Float64(g * 2.0))) / cbrt(Float64(a * 2.0))));
	else
		tmp = Float64(t_0 * cbrt(-1.0));
	end
	return tmp
end

Wolfram

code[g_, h_, a_] := Block[{t$95$0 = N[Power[N[(g / a), $MachinePrecision], 1/3], $MachinePrecision]}, If[LessEqual[g, 9.5e-198], N[(t$95$0 * N[Power[-1.0, 1/3], $MachinePrecision] + N[Power[N[(-0.25 * N[(N[(h / a), $MachinePrecision] * N[(h / g), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], If[LessEqual[g, 1.65e+148], N[(N[Power[N[(N[Power[N[(2.0 * a), $MachinePrecision], -1.0], $MachinePrecision] * N[((-g) + N[Sqrt[N[(N[(g * g), $MachinePrecision] - N[(h * h), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] + N[((-N[Power[N[(g * 2.0), $MachinePrecision], 1/3], $MachinePrecision]) / N[Power[N[(a * 2.0), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[Power[-1.0, 1/3], $MachinePrecision]), $MachinePrecision]]]]

Derivation

1. Split input into 3 regimes

2. if g < 9.4999999999999997e-198

   1. Initial program 41.9%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing

   3. Taylor expanded in h around 0

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{2}\right) + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)} \]
   4. Step-by-step derivation

      1. associate-*r* N/A

         \[\leadsto \color{blue}{\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}\right) \cdot \sqrt[3]{2}} + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right) \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right)} \]

      3. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      4. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      5. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      6. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      7. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \color{blue}{\sqrt[3]{2}}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      8. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)}\right) \]

      9. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      10. unpow2 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{\color{blue}{h \cdot h}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      11. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h \cdot h}{\color{blue}{g \cdot a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      12. times-frac N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      13. lower-*.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      14. lower-/.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g}} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      15. lower-/.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \color{blue}{\frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      16. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

      17. lower-*.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

      18. lower-cbrt.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\color{blue}{\sqrt[3]{\frac{1}{2}}} \cdot \sqrt[3]{\frac{-1}{2}}\right)\right) \]

      19. lower-cbrt.f64 71.3

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \color{blue}{\sqrt[3]{-0.5}}\right)\right) \]

   5. Applied rewrites 71.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 72.0%

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \color{blue}{\sqrt[3]{-1}}, \sqrt[3]{-0.25 \cdot \left(\frac{h}{a} \cdot \frac{h}{g}\right)}\right) \]

   if 9.4999999999999997e-198 < g < 1.65000000000000005e148

   1. Initial program 82.6%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      2. lift-*.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      3. lift-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a}} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

      4. associate-*l/ N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}{2 \cdot a}}} \]

      5. cbrt-div N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}}{\sqrt[3]{2 \cdot a}}} \]

      6. *-lft-identity N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\sqrt[3]{\color{blue}{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}}{\sqrt[3]{2 \cdot a}} \]

      7. lower-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}{\sqrt[3]{2 \cdot a}}} \]

   4. Applied rewrites 95.7%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{\left(g - h\right) \cdot \left(h + g\right)}}}{\sqrt[3]{a \cdot 2}}} \]

   5. Taylor expanded in g around -inf

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{-1 \cdot \sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}}}{\sqrt[3]{a \cdot 2}} \]
   6. Step-by-step derivation

      1. mul-1-neg N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{\mathsf{neg}\left(\sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}\right)}}{\sqrt[3]{a \cdot 2}} \]

      2. lower-neg.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{-\sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}}}{\sqrt[3]{a \cdot 2}} \]

      3. unpow2 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot \left(1 - \color{blue}{\sqrt{-1} \cdot \sqrt{-1}}\right)}}{\sqrt[3]{a \cdot 2}} \]

      4. rem-square-sqrt N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot \left(1 - \color{blue}{-1}\right)}}{\sqrt[3]{a \cdot 2}} \]

      5. metadata-eval N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot \color{blue}{2}}}{\sqrt[3]{a \cdot 2}} \]

      6. lower-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\color{blue}{\sqrt[3]{g \cdot 2}}}{\sqrt[3]{a \cdot 2}} \]

      7. lower-*.f64 98.1

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{\color{blue}{g \cdot 2}}}{\sqrt[3]{a \cdot 2}} \]

   7. Applied rewrites 98.1%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\color{blue}{-\sqrt[3]{g \cdot 2}}}{\sqrt[3]{a \cdot 2}} \]

   if 1.65000000000000005e148 < g

   1. Initial program 3.8%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      2. lift-*.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      3. lift-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a}} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

      4. associate-*l/ N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}{2 \cdot a}}} \]

      5. cbrt-div N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}}{\sqrt[3]{2 \cdot a}}} \]

      6. *-lft-identity N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\sqrt[3]{\color{blue}{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}}{\sqrt[3]{2 \cdot a}} \]

      7. lower-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}{\sqrt[3]{2 \cdot a}}} \]

   4. Applied rewrites 3.8%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{\left(g - h\right) \cdot \left(h + g\right)}}}{\sqrt[3]{a \cdot 2}}} \]

   5. Taylor expanded in g around inf

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]
   6. Step-by-step derivation

      1. lower-*.f64 N/A

         \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]

      2. lower-cbrt.f64 N/A

         \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{-1} \]

      3. lower-/.f64 N/A

         \[\leadsto \sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{-1} \]

      4. lower-cbrt.f64 62.6

         \[\leadsto \sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{-1}} \]

   7. Applied rewrites 62.6%

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]
3. Recombined 3 regimes into one program.

4. Final simplification 76.2%

   \[\leadsto \begin{array}{l} \mathbf{if}\;g \leq 9.5 \cdot 10^{-198}:\\ \;\;\;\;\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \sqrt[3]{-1}, \sqrt[3]{-0.25 \cdot \left(\frac{h}{a} \cdot \frac{h}{g}\right)}\right)\\ \mathbf{elif}\;g \leq 1.65 \cdot 10^{+148}:\\ \;\;\;\;\sqrt[3]{{\left(2 \cdot a\right)}^{-1} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{-\sqrt[3]{g \cdot 2}}{\sqrt[3]{a \cdot 2}}\\ \mathbf{else}:\\ \;\;\;\;\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}\\ \end{array} \]
5. Add Preprocessing

Alternative 6: 76.5% accurate, 0.6× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \sqrt[3]{\frac{g}{a}}\\ \mathbf{if}\;g \leq 9.5 \cdot 10^{-198}:\\ \;\;\;\;\mathsf{fma}\left(t\_0, \sqrt[3]{-1}, \sqrt[3]{-0.25 \cdot \left(\frac{h}{a} \cdot \frac{h}{g}\right)}\right)\\ \mathbf{elif}\;g \leq 3.4 \cdot 10^{+145}:\\ \;\;\;\;\sqrt[3]{{\left(2 \cdot a\right)}^{-1} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{0.5}{a}} \cdot \left(-\sqrt[3]{g \cdot 2}\right)\\ \mathbf{else}:\\ \;\;\;\;t\_0 \cdot \sqrt[3]{-1}\\ \end{array} \end{array} \]

FPCore

(FPCore (g h a)
 :precision binary64
 (let* ((t_0 (cbrt (/ g a))))
   (if (<= g 9.5e-198)
     (fma t_0 (cbrt -1.0) (cbrt (* -0.25 (* (/ h a) (/ h g)))))
     (if (<= g 3.4e+145)
       (+
        (cbrt (* (pow (* 2.0 a) -1.0) (+ (- g) (sqrt (- (* g g) (* h h))))))
        (* (cbrt (/ 0.5 a)) (- (cbrt (* g 2.0)))))
       (* t_0 (cbrt -1.0))))))

C

double code(double g, double h, double a) {
	double t_0 = cbrt((g / a));
	double tmp;
	if (g <= 9.5e-198) {
		tmp = fma(t_0, cbrt(-1.0), cbrt((-0.25 * ((h / a) * (h / g)))));
	} else if (g <= 3.4e+145) {
		tmp = cbrt((pow((2.0 * a), -1.0) * (-g + sqrt(((g * g) - (h * h)))))) + (cbrt((0.5 / a)) * -cbrt((g * 2.0)));
	} else {
		tmp = t_0 * cbrt(-1.0);
	}
	return tmp;
}

Julia

function code(g, h, a)
	t_0 = cbrt(Float64(g / a))
	tmp = 0.0
	if (g <= 9.5e-198)
		tmp = fma(t_0, cbrt(-1.0), cbrt(Float64(-0.25 * Float64(Float64(h / a) * Float64(h / g)))));
	elseif (g <= 3.4e+145)
		tmp = Float64(cbrt(Float64((Float64(2.0 * a) ^ -1.0) * Float64(Float64(-g) + sqrt(Float64(Float64(g * g) - Float64(h * h)))))) + Float64(cbrt(Float64(0.5 / a)) * Float64(-cbrt(Float64(g * 2.0)))));
	else
		tmp = Float64(t_0 * cbrt(-1.0));
	end
	return tmp
end

Wolfram

code[g_, h_, a_] := Block[{t$95$0 = N[Power[N[(g / a), $MachinePrecision], 1/3], $MachinePrecision]}, If[LessEqual[g, 9.5e-198], N[(t$95$0 * N[Power[-1.0, 1/3], $MachinePrecision] + N[Power[N[(-0.25 * N[(N[(h / a), $MachinePrecision] * N[(h / g), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], If[LessEqual[g, 3.4e+145], N[(N[Power[N[(N[Power[N[(2.0 * a), $MachinePrecision], -1.0], $MachinePrecision] * N[((-g) + N[Sqrt[N[(N[(g * g), $MachinePrecision] - N[(h * h), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] + N[(N[Power[N[(0.5 / a), $MachinePrecision], 1/3], $MachinePrecision] * (-N[Power[N[(g * 2.0), $MachinePrecision], 1/3], $MachinePrecision])), $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[Power[-1.0, 1/3], $MachinePrecision]), $MachinePrecision]]]]

Derivation

1. Split input into 3 regimes

2. if g < 9.4999999999999997e-198

   1. Initial program 41.9%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing

   3. Taylor expanded in h around 0

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{2}\right) + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)} \]
   4. Step-by-step derivation

      1. associate-*r* N/A

         \[\leadsto \color{blue}{\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}\right) \cdot \sqrt[3]{2}} + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right) \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right)} \]

      3. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      4. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      5. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      6. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      7. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \color{blue}{\sqrt[3]{2}}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      8. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)}\right) \]

      9. lower-cbrt.f64 N/A

         \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      10. unpow2 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{\color{blue}{h \cdot h}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      11. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h \cdot h}{\color{blue}{g \cdot a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      12. times-frac N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      13. lower-*.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      14. lower-/.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g}} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      15. lower-/.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \color{blue}{\frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

      16. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

      17. lower-*.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

      18. lower-cbrt.f64 N/A

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\color{blue}{\sqrt[3]{\frac{1}{2}}} \cdot \sqrt[3]{\frac{-1}{2}}\right)\right) \]

      19. lower-cbrt.f64 71.3

          \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \color{blue}{\sqrt[3]{-0.5}}\right)\right) \]

   5. Applied rewrites 71.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 72.0%

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \color{blue}{\sqrt[3]{-1}}, \sqrt[3]{-0.25 \cdot \left(\frac{h}{a} \cdot \frac{h}{g}\right)}\right) \]

   if 9.4999999999999997e-198 < g < 3.3999999999999999e145

   1. Initial program 82.0%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      2. lift-*.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      3. cbrt-prod N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}} \]

      4. pow1/3 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{{\left(\frac{1}{2 \cdot a}\right)}^{\frac{1}{3}}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}} \]

      5. lower-*.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{{\left(\frac{1}{2 \cdot a}\right)}^{\frac{1}{3}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}} \]

      6. pow1/3 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a}}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}} \]

      7. lower-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a}}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}} \]

      8. lift-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a}}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}} \]

      9. lift-*.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{\color{blue}{2 \cdot a}}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}} \]

      10. associate-/r* N/A

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{\frac{1}{2}}{a}}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}} \]

      11. metadata-eval N/A

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\color{blue}{\frac{1}{2}}}{a}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}} \]

      12. lower-/.f64 N/A

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{\frac{1}{2}}{a}}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}} \]

      13. lower-cbrt.f64 95.4

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{0.5}{a}} \cdot \color{blue}{\sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}} \]

      14. lift--.f64 N/A

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{\color{blue}{g \cdot g - h \cdot h}}} \]

      15. lift-*.f64 N/A

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{\color{blue}{g \cdot g} - h \cdot h}} \]

      16. lift-*.f64 N/A

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - \color{blue}{h \cdot h}}} \]

      17. difference-of-squares N/A

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{\color{blue}{\left(g + h\right) \cdot \left(g - h\right)}}} \]

      18. *-commutative N/A

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{\color{blue}{\left(g - h\right) \cdot \left(g + h\right)}}} \]

      19. lower-*.f64 N/A

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{\color{blue}{\left(g - h\right) \cdot \left(g + h\right)}}} \]

      20. lower--.f64 N/A

          \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{\color{blue}{\left(g - h\right)} \cdot \left(g + h\right)}} \]

   4. Applied rewrites 95.4%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{0.5}{a}} \cdot \sqrt[3]{\left(-g\right) - \sqrt{\left(g - h\right) \cdot \left(h + g\right)}}} \]

   5. Taylor expanded in g around -inf

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \color{blue}{\left(-1 \cdot \sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}\right)} \]
   6. Step-by-step derivation

      1. mul-1-neg N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \color{blue}{\left(\mathsf{neg}\left(\sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}\right)\right)} \]

      2. lower-neg.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \color{blue}{\left(-\sqrt[3]{g \cdot \left(1 - {\left(\sqrt{-1}\right)}^{2}\right)}\right)} \]

      3. unpow2 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \left(-\sqrt[3]{g \cdot \left(1 - \color{blue}{\sqrt{-1} \cdot \sqrt{-1}}\right)}\right) \]

      4. rem-square-sqrt N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \left(-\sqrt[3]{g \cdot \left(1 - \color{blue}{-1}\right)}\right) \]

      5. metadata-eval N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \left(-\sqrt[3]{g \cdot \color{blue}{2}}\right) \]

      6. lower-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{\frac{1}{2}}{a}} \cdot \left(-\color{blue}{\sqrt[3]{g \cdot 2}}\right) \]

      7. lower-*.f64 97.9

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{0.5}{a}} \cdot \left(-\sqrt[3]{\color{blue}{g \cdot 2}}\right) \]

   7. Applied rewrites 97.9%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{0.5}{a}} \cdot \color{blue}{\left(-\sqrt[3]{g \cdot 2}\right)} \]

   if 3.3999999999999999e145 < g

   1. Initial program 7.3%

      \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-cbrt.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      2. lift-*.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

      3. lift-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a}} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

      4. associate-*l/ N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}{2 \cdot a}}} \]

      5. cbrt-div N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}}{\sqrt[3]{2 \cdot a}}} \]

      6. *-lft-identity N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\sqrt[3]{\color{blue}{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}}{\sqrt[3]{2 \cdot a}} \]

      7. lower-/.f64 N/A

         \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}{\sqrt[3]{2 \cdot a}}} \]

   4. Applied rewrites 7.3%

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{\left(g - h\right) \cdot \left(h + g\right)}}}{\sqrt[3]{a \cdot 2}}} \]

   5. Taylor expanded in g around inf

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]
   6. Step-by-step derivation

      1. lower-*.f64 N/A

         \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]

      2. lower-cbrt.f64 N/A

         \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{-1} \]

      3. lower-/.f64 N/A

         \[\leadsto \sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{-1} \]

      4. lower-cbrt.f64 64.0

         \[\leadsto \sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{-1}} \]

   7. Applied rewrites 64.0%

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]
3. Recombined 3 regimes into one program.

4. Final simplification 76.2%

   \[\leadsto \begin{array}{l} \mathbf{if}\;g \leq 9.5 \cdot 10^{-198}:\\ \;\;\;\;\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \sqrt[3]{-1}, \sqrt[3]{-0.25 \cdot \left(\frac{h}{a} \cdot \frac{h}{g}\right)}\right)\\ \mathbf{elif}\;g \leq 3.4 \cdot 10^{+145}:\\ \;\;\;\;\sqrt[3]{{\left(2 \cdot a\right)}^{-1} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{0.5}{a}} \cdot \left(-\sqrt[3]{g \cdot 2}\right)\\ \mathbf{else}:\\ \;\;\;\;\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}\\ \end{array} \]
5. Add Preprocessing

Alternative 7: 75.4% accurate, 0.9× speedup

Math

\[\begin{array}{l} \\ \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \sqrt[3]{-1}, \sqrt[3]{-0.25 \cdot \left(\frac{h}{a} \cdot \frac{h}{g}\right)}\right) \end{array} \]

FPCore

(FPCore (g h a)
 :precision binary64
 (fma (cbrt (/ g a)) (cbrt -1.0) (cbrt (* -0.25 (* (/ h a) (/ h g))))))

C

double code(double g, double h, double a) {
	return fma(cbrt((g / a)), cbrt(-1.0), cbrt((-0.25 * ((h / a) * (h / g)))));
}

Julia

function code(g, h, a)
	return fma(cbrt(Float64(g / a)), cbrt(-1.0), cbrt(Float64(-0.25 * Float64(Float64(h / a) * Float64(h / g)))))
end

Wolfram

code[g_, h_, a_] := N[(N[Power[N[(g / a), $MachinePrecision], 1/3], $MachinePrecision] * N[Power[-1.0, 1/3], $MachinePrecision] + N[Power[N[(-0.25 * N[(N[(h / a), $MachinePrecision] * N[(h / g), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 43.7%

   \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
2. Add Preprocessing

3. Taylor expanded in h around 0

   \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{2}\right) + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)} \]
4. Step-by-step derivation

   1. associate-*r* N/A

      \[\leadsto \color{blue}{\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}\right) \cdot \sqrt[3]{2}} + \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right) \]

   2. lower-fma.f64 N/A

      \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right)} \]

   3. lower-*.f64 N/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   4. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   5. lower-/.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   6. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{\frac{-1}{2}}}, \sqrt[3]{2}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   7. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \color{blue}{\sqrt[3]{2}}, \sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   8. lower-*.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)}\right) \]

   9. lower-cbrt.f64 N/A

      \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \color{blue}{\sqrt[3]{\frac{{h}^{2}}{a \cdot g}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   10. unpow2 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{\color{blue}{h \cdot h}}{a \cdot g}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   11. *-commutative N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h \cdot h}{\color{blue}{g \cdot a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   12. times-frac N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   13. lower-*.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g} \cdot \frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   14. lower-/.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\color{blue}{\frac{h}{g}} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   15. lower-/.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \color{blue}{\frac{h}{a}}} \cdot \left(\sqrt[3]{\frac{-1}{2}} \cdot \sqrt[3]{\frac{1}{2}}\right)\right) \]

   16. *-commutative N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

   17. lower-*.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \color{blue}{\left(\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{-1}{2}}\right)}\right) \]

   18. lower-cbrt.f64 N/A

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{\frac{-1}{2}}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\color{blue}{\sqrt[3]{\frac{1}{2}}} \cdot \sqrt[3]{\frac{-1}{2}}\right)\right) \]

   19. lower-cbrt.f64 71.9

       \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \color{blue}{\sqrt[3]{-0.5}}\right)\right) \]

5. Applied rewrites 71.9%

   \[\leadsto \color{blue}{\mathsf{fma}\left(\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-0.5}, \sqrt[3]{2}, \sqrt[3]{\frac{h}{g} \cdot \frac{h}{a}} \cdot \left(\sqrt[3]{0.5} \cdot \sqrt[3]{-0.5}\right)\right)} \]
6. Step-by-step derivation

7. Applied rewrites 72.6%

   \[\leadsto \mathsf{fma}\left(\sqrt[3]{\frac{g}{a}}, \color{blue}{\sqrt[3]{-1}}, \sqrt[3]{-0.25 \cdot \left(\frac{h}{a} \cdot \frac{h}{g}\right)}\right) \]
8. Add Preprocessing

Alternative 8: 73.8% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1} \end{array} \]

FPCore

(FPCore (g h a) :precision binary64 (* (cbrt (/ g a)) (cbrt -1.0)))

C

double code(double g, double h, double a) {
	return cbrt((g / a)) * cbrt(-1.0);
}

Java

public static double code(double g, double h, double a) {
	return Math.cbrt((g / a)) * Math.cbrt(-1.0);
}

Julia

function code(g, h, a)
	return Float64(cbrt(Float64(g / a)) * cbrt(-1.0))
end

Wolfram

code[g_, h_, a_] := N[(N[Power[N[(g / a), $MachinePrecision], 1/3], $MachinePrecision] * N[Power[-1.0, 1/3], $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 43.7%

   \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
2. Add Preprocessing
3. Step-by-step derivation

   1. lift-cbrt.f64 N/A

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

   2. lift-*.f64 N/A

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}} \]

   3. lift-/.f64 N/A

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1}{2 \cdot a}} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

   4. associate-*l/ N/A

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\color{blue}{\frac{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}{2 \cdot a}}} \]

   5. cbrt-div N/A

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{1 \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)}}{\sqrt[3]{2 \cdot a}}} \]

   6. *-lft-identity N/A

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \frac{\sqrt[3]{\color{blue}{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}}{\sqrt[3]{2 \cdot a}} \]

   7. lower-/.f64 N/A

      \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{g \cdot g - h \cdot h}}}{\sqrt[3]{2 \cdot a}}} \]

4. Applied rewrites 46.9%

   \[\leadsto \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \color{blue}{\frac{\sqrt[3]{\left(-g\right) - \sqrt{\left(g - h\right) \cdot \left(h + g\right)}}}{\sqrt[3]{a \cdot 2}}} \]

5. Taylor expanded in g around inf

   \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]
6. Step-by-step derivation

   1. lower-*.f64 N/A

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]

   2. lower-cbrt.f64 N/A

      \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}}} \cdot \sqrt[3]{-1} \]

   3. lower-/.f64 N/A

      \[\leadsto \sqrt[3]{\color{blue}{\frac{g}{a}}} \cdot \sqrt[3]{-1} \]

   4. lower-cbrt.f64 71.5

      \[\leadsto \sqrt[3]{\frac{g}{a}} \cdot \color{blue}{\sqrt[3]{-1}} \]

7. Applied rewrites 71.5%

   \[\leadsto \color{blue}{\sqrt[3]{\frac{g}{a}} \cdot \sqrt[3]{-1}} \]
8. Add Preprocessing

Alternative 9: 0.0% accurate, 2.3× speedup

Math

\[\begin{array}{l} \\ -\sqrt[3]{\frac{\sqrt{-1} \cdot h}{a} \cdot -0.5} \end{array} \]

FPCore

(FPCore (g h a)
 :precision binary64
 (- (cbrt (* (/ (* (sqrt -1.0) h) a) -0.5))))

C

double code(double g, double h, double a) {
	return -cbrt((((sqrt(-1.0) * h) / a) * -0.5));
}

Java

public static double code(double g, double h, double a) {
	return -Math.cbrt((((Math.sqrt(-1.0) * h) / a) * -0.5));
}

Julia

function code(g, h, a)
	return Float64(-cbrt(Float64(Float64(Float64(sqrt(-1.0) * h) / a) * -0.5)))
end

Wolfram

code[g_, h_, a_] := (-N[Power[N[(N[(N[(N[Sqrt[-1.0], $MachinePrecision] * h), $MachinePrecision] / a), $MachinePrecision] * -0.5), $MachinePrecision], 1/3], $MachinePrecision])

Derivation

1. Initial program 43.7%

   \[\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]
2. Add Preprocessing
3. Step-by-step derivation

   1. lift-cbrt.f64 N/A

      \[\leadsto \color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)}} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

   2. pow1/3 N/A

      \[\leadsto \color{blue}{{\left(\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)\right)}^{\frac{1}{3}}} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

   3. sqr-pow N/A

      \[\leadsto \color{blue}{{\left(\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)\right)}^{\left(\frac{\frac{1}{3}}{2}\right)} \cdot {\left(\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)\right)}^{\left(\frac{\frac{1}{3}}{2}\right)}} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

   4. fabs-sqr N/A

      \[\leadsto \color{blue}{\left|{\left(\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)\right)}^{\left(\frac{\frac{1}{3}}{2}\right)} \cdot {\left(\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)\right)}^{\left(\frac{\frac{1}{3}}{2}\right)}\right|} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

   5. sqr-pow N/A

      \[\leadsto \left|\color{blue}{{\left(\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)\right)}^{\frac{1}{3}}}\right| + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

   6. pow1/3 N/A

      \[\leadsto \left|\color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)}}\right| + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

   7. lift-cbrt.f64 N/A

      \[\leadsto \left|\color{blue}{\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)}}\right| + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

   8. lower-fabs.f64 34.8

      \[\leadsto \color{blue}{\left|\sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)}\right|} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

   9. lift-*.f64 N/A

      \[\leadsto \left|\sqrt[3]{\color{blue}{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right)}}\right| + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

   10. *-commutative N/A

       \[\leadsto \left|\sqrt[3]{\color{blue}{\left(\left(-g\right) + \sqrt{g \cdot g - h \cdot h}\right) \cdot \frac{1}{2 \cdot a}}}\right| + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

4. Applied rewrites 10.6%

   \[\leadsto \color{blue}{\left|\sqrt[3]{\mathsf{fma}\left(\sqrt{g - h}, \sqrt{h + g}, -g\right) \cdot \frac{0.5}{a}}\right|} + \sqrt[3]{\frac{1}{2 \cdot a} \cdot \left(\left(-g\right) - \sqrt{g \cdot g - h \cdot h}\right)} \]

5. Taylor expanded in h around -inf

   \[\leadsto \color{blue}{-1 \cdot \left(\sqrt[3]{\frac{h \cdot \sqrt{-1}}{a}} \cdot \sqrt[3]{\frac{-1}{2}}\right)} \]
6. Step-by-step derivation

   1. mul-1-neg N/A

      \[\leadsto \color{blue}{\mathsf{neg}\left(\sqrt[3]{\frac{h \cdot \sqrt{-1}}{a}} \cdot \sqrt[3]{\frac{-1}{2}}\right)} \]

   2. lower-neg.f64 N/A

      \[\leadsto \color{blue}{-\sqrt[3]{\frac{h \cdot \sqrt{-1}}{a}} \cdot \sqrt[3]{\frac{-1}{2}}} \]

   3. lower-*.f64 N/A

      \[\leadsto -\color{blue}{\sqrt[3]{\frac{h \cdot \sqrt{-1}}{a}} \cdot \sqrt[3]{\frac{-1}{2}}} \]

   4. lower-cbrt.f64 N/A

      \[\leadsto -\color{blue}{\sqrt[3]{\frac{h \cdot \sqrt{-1}}{a}}} \cdot \sqrt[3]{\frac{-1}{2}} \]

   5. lower-/.f64 N/A

      \[\leadsto -\sqrt[3]{\color{blue}{\frac{h \cdot \sqrt{-1}}{a}}} \cdot \sqrt[3]{\frac{-1}{2}} \]

   6. lower-*.f64 N/A

      \[\leadsto -\sqrt[3]{\frac{\color{blue}{h \cdot \sqrt{-1}}}{a}} \cdot \sqrt[3]{\frac{-1}{2}} \]

   7. lower-sqrt.f64 N/A

      \[\leadsto -\sqrt[3]{\frac{h \cdot \color{blue}{\sqrt{-1}}}{a}} \cdot \sqrt[3]{\frac{-1}{2}} \]

   8. lower-cbrt.f64 0.0

      \[\leadsto -\sqrt[3]{\frac{h \cdot \sqrt{-1}}{a}} \cdot \color{blue}{\sqrt[3]{-0.5}} \]

7. Applied rewrites 0.0%

   \[\leadsto \color{blue}{-\sqrt[3]{\frac{h \cdot \sqrt{-1}}{a}} \cdot \sqrt[3]{-0.5}} \]
8. Step-by-step derivation

9. Applied rewrites 0.0%

   \[\leadsto \color{blue}{-\sqrt[3]{\frac{\sqrt{-1} \cdot h}{a} \cdot -0.5}} \]
10. Add Preprocessing

Reproduce

herbie shell --seed 2024346 
(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))))))))