
(FPCore (a b c) :precision binary64 (/ (+ (- b) (sqrt (- (* b b) (* (* 4.0 a) c)))) (* 2.0 a)))
double code(double a, double b, double c) {
return (-b + sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a);
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = (-b + sqrt(((b * b) - ((4.0d0 * a) * c)))) / (2.0d0 * a)
end function
public static double code(double a, double b, double c) {
return (-b + Math.sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a);
}
def code(a, b, c): return (-b + math.sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a)
function code(a, b, c) return Float64(Float64(Float64(-b) + sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c)))) / Float64(2.0 * a)) end
function tmp = code(a, b, c) tmp = (-b + sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a); end
code[a_, b_, c_] := N[(N[((-b) + N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(-b\right) + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}{2 \cdot a}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 10 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b c) :precision binary64 (/ (+ (- b) (sqrt (- (* b b) (* (* 4.0 a) c)))) (* 2.0 a)))
double code(double a, double b, double c) {
return (-b + sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a);
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = (-b + sqrt(((b * b) - ((4.0d0 * a) * c)))) / (2.0d0 * a)
end function
public static double code(double a, double b, double c) {
return (-b + Math.sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a);
}
def code(a, b, c): return (-b + math.sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a)
function code(a, b, c) return Float64(Float64(Float64(-b) + sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c)))) / Float64(2.0 * a)) end
function tmp = code(a, b, c) tmp = (-b + sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a); end
code[a_, b_, c_] := N[(N[((-b) + N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(-b\right) + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}{2 \cdot a}
\end{array}
(FPCore (a b c) :precision binary64 (/ c (* a (/ (+ b (sqrt (+ (* b b) (* -4.0 (* c a))))) (/ a -0.5)))))
double code(double a, double b, double c) {
return c / (a * ((b + sqrt(((b * b) + (-4.0 * (c * a))))) / (a / -0.5)));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = c / (a * ((b + sqrt(((b * b) + ((-4.0d0) * (c * a))))) / (a / (-0.5d0))))
end function
public static double code(double a, double b, double c) {
return c / (a * ((b + Math.sqrt(((b * b) + (-4.0 * (c * a))))) / (a / -0.5)));
}
def code(a, b, c): return c / (a * ((b + math.sqrt(((b * b) + (-4.0 * (c * a))))) / (a / -0.5)))
function code(a, b, c) return Float64(c / Float64(a * Float64(Float64(b + sqrt(Float64(Float64(b * b) + Float64(-4.0 * Float64(c * a))))) / Float64(a / -0.5)))) end
function tmp = code(a, b, c) tmp = c / (a * ((b + sqrt(((b * b) + (-4.0 * (c * a))))) / (a / -0.5))); end
code[a_, b_, c_] := N[(c / N[(a * N[(N[(b + N[Sqrt[N[(N[(b * b), $MachinePrecision] + N[(-4.0 * N[(c * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(a / -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{c}{a \cdot \frac{b + \sqrt{b \cdot b + -4 \cdot \left(c \cdot a\right)}}{\frac{a}{-0.5}}}
\end{array}
Initial program 53.2%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified53.2%
div-subN/A
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr53.2%
Taylor expanded in b around 0
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f6499.1%
Simplified99.1%
clear-numN/A
distribute-frac-negN/A
remove-double-negN/A
distribute-frac-negN/A
distribute-frac-negN/A
neg-lowering-neg.f64N/A
/-lowering-/.f64N/A
distribute-frac-negN/A
remove-double-negN/A
/-lowering-/.f64N/A
*-commutativeN/A
clear-numN/A
div-invN/A
metadata-evalN/A
*-commutativeN/A
un-div-invN/A
Applied egg-rr99.3%
distribute-neg-frac2N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
Final simplification99.5%
(FPCore (a b c)
:precision binary64
(if (<= (/ (- (sqrt (- (* b b) (* c (* a 4.0)))) b) (* a 2.0)) -2.5)
(/ 0.5 (* a (/ 1.0 (- (sqrt (+ (* b b) (* a (* c -4.0)))) b))))
(/
1.0
(-
(*
a
(+
(/ 1.0 b)
(/ (+ (* c a) (/ (* (* 2.0 (* a a)) (* c c)) (* b b))) (* b (* b b)))))
(/ b c)))))
double code(double a, double b, double c) {
double tmp;
if (((sqrt(((b * b) - (c * (a * 4.0)))) - b) / (a * 2.0)) <= -2.5) {
tmp = 0.5 / (a * (1.0 / (sqrt(((b * b) + (a * (c * -4.0)))) - b)));
} else {
tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c));
}
return tmp;
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (((sqrt(((b * b) - (c * (a * 4.0d0)))) - b) / (a * 2.0d0)) <= (-2.5d0)) then
tmp = 0.5d0 / (a * (1.0d0 / (sqrt(((b * b) + (a * (c * (-4.0d0))))) - b)))
else
tmp = 1.0d0 / ((a * ((1.0d0 / b) + (((c * a) + (((2.0d0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c))
end if
code = tmp
end function
public static double code(double a, double b, double c) {
double tmp;
if (((Math.sqrt(((b * b) - (c * (a * 4.0)))) - b) / (a * 2.0)) <= -2.5) {
tmp = 0.5 / (a * (1.0 / (Math.sqrt(((b * b) + (a * (c * -4.0)))) - b)));
} else {
tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c));
}
return tmp;
}
def code(a, b, c): tmp = 0 if ((math.sqrt(((b * b) - (c * (a * 4.0)))) - b) / (a * 2.0)) <= -2.5: tmp = 0.5 / (a * (1.0 / (math.sqrt(((b * b) + (a * (c * -4.0)))) - b))) else: tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c)) return tmp
function code(a, b, c) tmp = 0.0 if (Float64(Float64(sqrt(Float64(Float64(b * b) - Float64(c * Float64(a * 4.0)))) - b) / Float64(a * 2.0)) <= -2.5) tmp = Float64(0.5 / Float64(a * Float64(1.0 / Float64(sqrt(Float64(Float64(b * b) + Float64(a * Float64(c * -4.0)))) - b)))); else tmp = Float64(1.0 / Float64(Float64(a * Float64(Float64(1.0 / b) + Float64(Float64(Float64(c * a) + Float64(Float64(Float64(2.0 * Float64(a * a)) * Float64(c * c)) / Float64(b * b))) / Float64(b * Float64(b * b))))) - Float64(b / c))); end return tmp end
function tmp_2 = code(a, b, c) tmp = 0.0; if (((sqrt(((b * b) - (c * (a * 4.0)))) - b) / (a * 2.0)) <= -2.5) tmp = 0.5 / (a * (1.0 / (sqrt(((b * b) + (a * (c * -4.0)))) - b))); else tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c)); end tmp_2 = tmp; end
code[a_, b_, c_] := If[LessEqual[N[(N[(N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(c * N[(a * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], -2.5], N[(0.5 / N[(a * N[(1.0 / N[(N[Sqrt[N[(N[(b * b), $MachinePrecision] + N[(a * N[(c * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(a * N[(N[(1.0 / b), $MachinePrecision] + N[(N[(N[(c * a), $MachinePrecision] + N[(N[(N[(2.0 * N[(a * a), $MachinePrecision]), $MachinePrecision] * N[(c * c), $MachinePrecision]), $MachinePrecision] / N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{\sqrt{b \cdot b - c \cdot \left(a \cdot 4\right)} - b}{a \cdot 2} \leq -2.5:\\
\;\;\;\;\frac{0.5}{a \cdot \frac{1}{\sqrt{b \cdot b + a \cdot \left(c \cdot -4\right)} - b}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{a \cdot \left(\frac{1}{b} + \frac{c \cdot a + \frac{\left(2 \cdot \left(a \cdot a\right)\right) \cdot \left(c \cdot c\right)}{b \cdot b}}{b \cdot \left(b \cdot b\right)}\right) - \frac{b}{c}}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 #s(literal 4 binary64) a) c)))) (*.f64 #s(literal 2 binary64) a)) < -2.5Initial program 88.2%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified88.2%
div-invN/A
flip--N/A
clear-numN/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
frac-timesN/A
metadata-evalN/A
/-lowering-/.f64N/A
Applied egg-rr88.2%
if -2.5 < (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 #s(literal 4 binary64) a) c)))) (*.f64 #s(literal 2 binary64) a)) Initial program 48.0%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified48.0%
div-subN/A
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr48.1%
Taylor expanded in a around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified94.8%
Taylor expanded in b around inf
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6494.8%
Simplified94.8%
Final simplification94.0%
(FPCore (a b c)
:precision binary64
(if (<= (/ (- (sqrt (- (* b b) (* c (* a 4.0)))) b) (* a 2.0)) -2.5)
(/ 0.5 (/ a (- (sqrt (+ (* b b) (* a (* c -4.0)))) b)))
(/
1.0
(-
(*
a
(+
(/ 1.0 b)
(/ (+ (* c a) (/ (* (* 2.0 (* a a)) (* c c)) (* b b))) (* b (* b b)))))
(/ b c)))))
double code(double a, double b, double c) {
double tmp;
if (((sqrt(((b * b) - (c * (a * 4.0)))) - b) / (a * 2.0)) <= -2.5) {
tmp = 0.5 / (a / (sqrt(((b * b) + (a * (c * -4.0)))) - b));
} else {
tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c));
}
return tmp;
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (((sqrt(((b * b) - (c * (a * 4.0d0)))) - b) / (a * 2.0d0)) <= (-2.5d0)) then
tmp = 0.5d0 / (a / (sqrt(((b * b) + (a * (c * (-4.0d0))))) - b))
else
tmp = 1.0d0 / ((a * ((1.0d0 / b) + (((c * a) + (((2.0d0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c))
end if
code = tmp
end function
public static double code(double a, double b, double c) {
double tmp;
if (((Math.sqrt(((b * b) - (c * (a * 4.0)))) - b) / (a * 2.0)) <= -2.5) {
tmp = 0.5 / (a / (Math.sqrt(((b * b) + (a * (c * -4.0)))) - b));
} else {
tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c));
}
return tmp;
}
def code(a, b, c): tmp = 0 if ((math.sqrt(((b * b) - (c * (a * 4.0)))) - b) / (a * 2.0)) <= -2.5: tmp = 0.5 / (a / (math.sqrt(((b * b) + (a * (c * -4.0)))) - b)) else: tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c)) return tmp
function code(a, b, c) tmp = 0.0 if (Float64(Float64(sqrt(Float64(Float64(b * b) - Float64(c * Float64(a * 4.0)))) - b) / Float64(a * 2.0)) <= -2.5) tmp = Float64(0.5 / Float64(a / Float64(sqrt(Float64(Float64(b * b) + Float64(a * Float64(c * -4.0)))) - b))); else tmp = Float64(1.0 / Float64(Float64(a * Float64(Float64(1.0 / b) + Float64(Float64(Float64(c * a) + Float64(Float64(Float64(2.0 * Float64(a * a)) * Float64(c * c)) / Float64(b * b))) / Float64(b * Float64(b * b))))) - Float64(b / c))); end return tmp end
function tmp_2 = code(a, b, c) tmp = 0.0; if (((sqrt(((b * b) - (c * (a * 4.0)))) - b) / (a * 2.0)) <= -2.5) tmp = 0.5 / (a / (sqrt(((b * b) + (a * (c * -4.0)))) - b)); else tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c)); end tmp_2 = tmp; end
code[a_, b_, c_] := If[LessEqual[N[(N[(N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(c * N[(a * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], -2.5], N[(0.5 / N[(a / N[(N[Sqrt[N[(N[(b * b), $MachinePrecision] + N[(a * N[(c * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(a * N[(N[(1.0 / b), $MachinePrecision] + N[(N[(N[(c * a), $MachinePrecision] + N[(N[(N[(2.0 * N[(a * a), $MachinePrecision]), $MachinePrecision] * N[(c * c), $MachinePrecision]), $MachinePrecision] / N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{\sqrt{b \cdot b - c \cdot \left(a \cdot 4\right)} - b}{a \cdot 2} \leq -2.5:\\
\;\;\;\;\frac{0.5}{\frac{a}{\sqrt{b \cdot b + a \cdot \left(c \cdot -4\right)} - b}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{a \cdot \left(\frac{1}{b} + \frac{c \cdot a + \frac{\left(2 \cdot \left(a \cdot a\right)\right) \cdot \left(c \cdot c\right)}{b \cdot b}}{b \cdot \left(b \cdot b\right)}\right) - \frac{b}{c}}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 #s(literal 4 binary64) a) c)))) (*.f64 #s(literal 2 binary64) a)) < -2.5Initial program 88.2%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified88.2%
associate-/r*N/A
clear-numN/A
associate-/l/N/A
associate-/r*N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6488.2%
Applied egg-rr88.2%
if -2.5 < (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 #s(literal 4 binary64) a) c)))) (*.f64 #s(literal 2 binary64) a)) Initial program 48.0%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified48.0%
div-subN/A
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr48.1%
Taylor expanded in a around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified94.8%
Taylor expanded in b around inf
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6494.8%
Simplified94.8%
Final simplification94.0%
(FPCore (a b c)
:precision binary64
(if (<= b 0.00061)
(* (- (sqrt (+ (* b b) (* a (* c -4.0)))) b) (/ 0.5 a))
(/
1.0
(-
(*
a
(+
(/ 1.0 b)
(/ (+ (* c a) (/ (* (* 2.0 (* a a)) (* c c)) (* b b))) (* b (* b b)))))
(/ b c)))))
double code(double a, double b, double c) {
double tmp;
if (b <= 0.00061) {
tmp = (sqrt(((b * b) + (a * (c * -4.0)))) - b) * (0.5 / a);
} else {
tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c));
}
return tmp;
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (b <= 0.00061d0) then
tmp = (sqrt(((b * b) + (a * (c * (-4.0d0))))) - b) * (0.5d0 / a)
else
tmp = 1.0d0 / ((a * ((1.0d0 / b) + (((c * a) + (((2.0d0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c))
end if
code = tmp
end function
public static double code(double a, double b, double c) {
double tmp;
if (b <= 0.00061) {
tmp = (Math.sqrt(((b * b) + (a * (c * -4.0)))) - b) * (0.5 / a);
} else {
tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c));
}
return tmp;
}
def code(a, b, c): tmp = 0 if b <= 0.00061: tmp = (math.sqrt(((b * b) + (a * (c * -4.0)))) - b) * (0.5 / a) else: tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c)) return tmp
function code(a, b, c) tmp = 0.0 if (b <= 0.00061) tmp = Float64(Float64(sqrt(Float64(Float64(b * b) + Float64(a * Float64(c * -4.0)))) - b) * Float64(0.5 / a)); else tmp = Float64(1.0 / Float64(Float64(a * Float64(Float64(1.0 / b) + Float64(Float64(Float64(c * a) + Float64(Float64(Float64(2.0 * Float64(a * a)) * Float64(c * c)) / Float64(b * b))) / Float64(b * Float64(b * b))))) - Float64(b / c))); end return tmp end
function tmp_2 = code(a, b, c) tmp = 0.0; if (b <= 0.00061) tmp = (sqrt(((b * b) + (a * (c * -4.0)))) - b) * (0.5 / a); else tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c)); end tmp_2 = tmp; end
code[a_, b_, c_] := If[LessEqual[b, 0.00061], N[(N[(N[Sqrt[N[(N[(b * b), $MachinePrecision] + N[(a * N[(c * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] * N[(0.5 / a), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(a * N[(N[(1.0 / b), $MachinePrecision] + N[(N[(N[(c * a), $MachinePrecision] + N[(N[(N[(2.0 * N[(a * a), $MachinePrecision]), $MachinePrecision] * N[(c * c), $MachinePrecision]), $MachinePrecision] / N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 0.00061:\\
\;\;\;\;\left(\sqrt{b \cdot b + a \cdot \left(c \cdot -4\right)} - b\right) \cdot \frac{0.5}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{a \cdot \left(\frac{1}{b} + \frac{c \cdot a + \frac{\left(2 \cdot \left(a \cdot a\right)\right) \cdot \left(c \cdot c\right)}{b \cdot b}}{b \cdot \left(b \cdot b\right)}\right) - \frac{b}{c}}\\
\end{array}
\end{array}
if b < 6.09999999999999974e-4Initial program 93.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified93.4%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6493.4%
Applied egg-rr93.4%
if 6.09999999999999974e-4 < b Initial program 51.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified51.4%
div-subN/A
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr51.4%
Taylor expanded in a around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified93.7%
Taylor expanded in b around inf
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6493.7%
Simplified93.7%
Final simplification93.6%
(FPCore (a b c) :precision binary64 (* (/ a -0.5) (/ (/ c a) (+ b (sqrt (+ (* b b) (* -4.0 (* c a))))))))
double code(double a, double b, double c) {
return (a / -0.5) * ((c / a) / (b + sqrt(((b * b) + (-4.0 * (c * a))))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = (a / (-0.5d0)) * ((c / a) / (b + sqrt(((b * b) + ((-4.0d0) * (c * a))))))
end function
public static double code(double a, double b, double c) {
return (a / -0.5) * ((c / a) / (b + Math.sqrt(((b * b) + (-4.0 * (c * a))))));
}
def code(a, b, c): return (a / -0.5) * ((c / a) / (b + math.sqrt(((b * b) + (-4.0 * (c * a))))))
function code(a, b, c) return Float64(Float64(a / -0.5) * Float64(Float64(c / a) / Float64(b + sqrt(Float64(Float64(b * b) + Float64(-4.0 * Float64(c * a))))))) end
function tmp = code(a, b, c) tmp = (a / -0.5) * ((c / a) / (b + sqrt(((b * b) + (-4.0 * (c * a)))))); end
code[a_, b_, c_] := N[(N[(a / -0.5), $MachinePrecision] * N[(N[(c / a), $MachinePrecision] / N[(b + N[Sqrt[N[(N[(b * b), $MachinePrecision] + N[(-4.0 * N[(c * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{a}{-0.5} \cdot \frac{\frac{c}{a}}{b + \sqrt{b \cdot b + -4 \cdot \left(c \cdot a\right)}}
\end{array}
Initial program 53.2%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified53.2%
div-subN/A
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr53.2%
Taylor expanded in b around 0
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f6499.1%
Simplified99.1%
clear-numN/A
distribute-frac-negN/A
remove-double-negN/A
distribute-frac-negN/A
distribute-frac-negN/A
neg-lowering-neg.f64N/A
/-lowering-/.f64N/A
distribute-frac-negN/A
remove-double-negN/A
/-lowering-/.f64N/A
*-commutativeN/A
clear-numN/A
div-invN/A
metadata-evalN/A
*-commutativeN/A
un-div-invN/A
Applied egg-rr99.3%
associate-/r/N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
Applied egg-rr99.3%
Final simplification99.3%
(FPCore (a b c)
:precision binary64
(/
1.0
(-
(*
a
(+
(/ 1.0 b)
(/ (+ (* c a) (/ (* (* 2.0 (* a a)) (* c c)) (* b b))) (* b (* b b)))))
(/ b c))))
double code(double a, double b, double c) {
return 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = 1.0d0 / ((a * ((1.0d0 / b) + (((c * a) + (((2.0d0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c))
end function
public static double code(double a, double b, double c) {
return 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c));
}
def code(a, b, c): return 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c))
function code(a, b, c) return Float64(1.0 / Float64(Float64(a * Float64(Float64(1.0 / b) + Float64(Float64(Float64(c * a) + Float64(Float64(Float64(2.0 * Float64(a * a)) * Float64(c * c)) / Float64(b * b))) / Float64(b * Float64(b * b))))) - Float64(b / c))) end
function tmp = code(a, b, c) tmp = 1.0 / ((a * ((1.0 / b) + (((c * a) + (((2.0 * (a * a)) * (c * c)) / (b * b))) / (b * (b * b))))) - (b / c)); end
code[a_, b_, c_] := N[(1.0 / N[(N[(a * N[(N[(1.0 / b), $MachinePrecision] + N[(N[(N[(c * a), $MachinePrecision] + N[(N[(N[(2.0 * N[(a * a), $MachinePrecision]), $MachinePrecision] * N[(c * c), $MachinePrecision]), $MachinePrecision] / N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{a \cdot \left(\frac{1}{b} + \frac{c \cdot a + \frac{\left(2 \cdot \left(a \cdot a\right)\right) \cdot \left(c \cdot c\right)}{b \cdot b}}{b \cdot \left(b \cdot b\right)}\right) - \frac{b}{c}}
\end{array}
Initial program 53.2%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified53.2%
div-subN/A
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr53.2%
Taylor expanded in a around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified92.0%
Taylor expanded in b around inf
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.0%
Simplified92.0%
Final simplification92.0%
(FPCore (a b c) :precision binary64 (/ 1.0 (- (/ (+ a (/ (* c (* a a)) (* b b))) b) (/ b c))))
double code(double a, double b, double c) {
return 1.0 / (((a + ((c * (a * a)) / (b * b))) / b) - (b / c));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = 1.0d0 / (((a + ((c * (a * a)) / (b * b))) / b) - (b / c))
end function
public static double code(double a, double b, double c) {
return 1.0 / (((a + ((c * (a * a)) / (b * b))) / b) - (b / c));
}
def code(a, b, c): return 1.0 / (((a + ((c * (a * a)) / (b * b))) / b) - (b / c))
function code(a, b, c) return Float64(1.0 / Float64(Float64(Float64(a + Float64(Float64(c * Float64(a * a)) / Float64(b * b))) / b) - Float64(b / c))) end
function tmp = code(a, b, c) tmp = 1.0 / (((a + ((c * (a * a)) / (b * b))) / b) - (b / c)); end
code[a_, b_, c_] := N[(1.0 / N[(N[(N[(a + N[(N[(c * N[(a * a), $MachinePrecision]), $MachinePrecision] / N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / b), $MachinePrecision] - N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{a + \frac{c \cdot \left(a \cdot a\right)}{b \cdot b}}{b} - \frac{b}{c}}
\end{array}
Initial program 53.2%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified53.2%
div-subN/A
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr53.2%
Taylor expanded in a around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified92.0%
Taylor expanded in b around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6489.3%
Simplified89.3%
Final simplification89.3%
(FPCore (a b c) :precision binary64 (/ 1.0 (- (/ a b) (/ b c))))
double code(double a, double b, double c) {
return 1.0 / ((a / b) - (b / c));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = 1.0d0 / ((a / b) - (b / c))
end function
public static double code(double a, double b, double c) {
return 1.0 / ((a / b) - (b / c));
}
def code(a, b, c): return 1.0 / ((a / b) - (b / c))
function code(a, b, c) return Float64(1.0 / Float64(Float64(a / b) - Float64(b / c))) end
function tmp = code(a, b, c) tmp = 1.0 / ((a / b) - (b / c)); end
code[a_, b_, c_] := N[(1.0 / N[(N[(a / b), $MachinePrecision] - N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{a}{b} - \frac{b}{c}}
\end{array}
Initial program 53.2%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified53.2%
div-subN/A
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr53.2%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f6484.2%
Simplified84.2%
Final simplification84.2%
(FPCore (a b c) :precision binary64 (- 0.0 (/ c b)))
double code(double a, double b, double c) {
return 0.0 - (c / b);
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = 0.0d0 - (c / b)
end function
public static double code(double a, double b, double c) {
return 0.0 - (c / b);
}
def code(a, b, c): return 0.0 - (c / b)
function code(a, b, c) return Float64(0.0 - Float64(c / b)) end
function tmp = code(a, b, c) tmp = 0.0 - (c / b); end
code[a_, b_, c_] := N[(0.0 - N[(c / b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0 - \frac{c}{b}
\end{array}
Initial program 53.2%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified53.2%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
/-lowering-/.f6466.5%
Simplified66.5%
Taylor expanded in c around 0
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f6466.5%
Simplified66.5%
Final simplification66.5%
(FPCore (a b c) :precision binary64 (/ b a))
double code(double a, double b, double c) {
return b / a;
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = b / a
end function
public static double code(double a, double b, double c) {
return b / a;
}
def code(a, b, c): return b / a
function code(a, b, c) return Float64(b / a) end
function tmp = code(a, b, c) tmp = b / a; end
code[a_, b_, c_] := N[(b / a), $MachinePrecision]
\begin{array}{l}
\\
\frac{b}{a}
\end{array}
Initial program 53.2%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
distribute-rgt-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
Simplified53.2%
div-subN/A
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr53.2%
Taylor expanded in b around 0
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f6499.1%
Simplified99.1%
Taylor expanded in c around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-lowering-neg.f64N/A
/-lowering-/.f6484.1%
Simplified84.1%
Taylor expanded in c around inf
/-lowering-/.f641.6%
Simplified1.6%
herbie shell --seed 2024145
(FPCore (a b c)
:name "Quadratic roots, narrow range"
:precision binary64
:pre (and (and (and (< 1.0536712127723509e-8 a) (< a 94906265.62425156)) (and (< 1.0536712127723509e-8 b) (< b 94906265.62425156))) (and (< 1.0536712127723509e-8 c) (< c 94906265.62425156)))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4.0 a) c)))) (* 2.0 a)))