
(FPCore (a b c) :precision binary64 (let* ((t_0 (sqrt (- (* b b) (* (* 4.0 a) c))))) (if (>= b 0.0) (/ (- (- b) t_0) (* 2.0 a)) (/ (* 2.0 c) (+ (- b) t_0)))))
double code(double a, double b, double c) {
double t_0 = sqrt(((b * b) - ((4.0 * a) * c)));
double tmp;
if (b >= 0.0) {
tmp = (-b - t_0) / (2.0 * a);
} else {
tmp = (2.0 * c) / (-b + t_0);
}
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) :: t_0
real(8) :: tmp
t_0 = sqrt(((b * b) - ((4.0d0 * a) * c)))
if (b >= 0.0d0) then
tmp = (-b - t_0) / (2.0d0 * a)
else
tmp = (2.0d0 * c) / (-b + t_0)
end if
code = tmp
end function
public static double code(double a, double b, double c) {
double t_0 = Math.sqrt(((b * b) - ((4.0 * a) * c)));
double tmp;
if (b >= 0.0) {
tmp = (-b - t_0) / (2.0 * a);
} else {
tmp = (2.0 * c) / (-b + t_0);
}
return tmp;
}
def code(a, b, c): t_0 = math.sqrt(((b * b) - ((4.0 * a) * c))) tmp = 0 if b >= 0.0: tmp = (-b - t_0) / (2.0 * a) else: tmp = (2.0 * c) / (-b + t_0) return tmp
function code(a, b, c) t_0 = sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c))) tmp = 0.0 if (b >= 0.0) tmp = Float64(Float64(Float64(-b) - t_0) / Float64(2.0 * a)); else tmp = Float64(Float64(2.0 * c) / Float64(Float64(-b) + t_0)); end return tmp end
function tmp_2 = code(a, b, c) t_0 = sqrt(((b * b) - ((4.0 * a) * c))); tmp = 0.0; if (b >= 0.0) tmp = (-b - t_0) / (2.0 * a); else tmp = (2.0 * c) / (-b + t_0); end tmp_2 = tmp; end
code[a_, b_, c_] := Block[{t$95$0 = N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[GreaterEqual[b, 0.0], N[(N[((-b) - t$95$0), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * c), $MachinePrecision] / N[((-b) + t$95$0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}\\
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-b\right) - t_0}{2 \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot c}{\left(-b\right) + t_0}\\
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 14 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b c) :precision binary64 (let* ((t_0 (sqrt (- (* b b) (* (* 4.0 a) c))))) (if (>= b 0.0) (/ (- (- b) t_0) (* 2.0 a)) (/ (* 2.0 c) (+ (- b) t_0)))))
double code(double a, double b, double c) {
double t_0 = sqrt(((b * b) - ((4.0 * a) * c)));
double tmp;
if (b >= 0.0) {
tmp = (-b - t_0) / (2.0 * a);
} else {
tmp = (2.0 * c) / (-b + t_0);
}
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) :: t_0
real(8) :: tmp
t_0 = sqrt(((b * b) - ((4.0d0 * a) * c)))
if (b >= 0.0d0) then
tmp = (-b - t_0) / (2.0d0 * a)
else
tmp = (2.0d0 * c) / (-b + t_0)
end if
code = tmp
end function
public static double code(double a, double b, double c) {
double t_0 = Math.sqrt(((b * b) - ((4.0 * a) * c)));
double tmp;
if (b >= 0.0) {
tmp = (-b - t_0) / (2.0 * a);
} else {
tmp = (2.0 * c) / (-b + t_0);
}
return tmp;
}
def code(a, b, c): t_0 = math.sqrt(((b * b) - ((4.0 * a) * c))) tmp = 0 if b >= 0.0: tmp = (-b - t_0) / (2.0 * a) else: tmp = (2.0 * c) / (-b + t_0) return tmp
function code(a, b, c) t_0 = sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c))) tmp = 0.0 if (b >= 0.0) tmp = Float64(Float64(Float64(-b) - t_0) / Float64(2.0 * a)); else tmp = Float64(Float64(2.0 * c) / Float64(Float64(-b) + t_0)); end return tmp end
function tmp_2 = code(a, b, c) t_0 = sqrt(((b * b) - ((4.0 * a) * c))); tmp = 0.0; if (b >= 0.0) tmp = (-b - t_0) / (2.0 * a); else tmp = (2.0 * c) / (-b + t_0); end tmp_2 = tmp; end
code[a_, b_, c_] := Block[{t$95$0 = N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[GreaterEqual[b, 0.0], N[(N[((-b) - t$95$0), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * c), $MachinePrecision] / N[((-b) + t$95$0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}\\
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-b\right) - t_0}{2 \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot c}{\left(-b\right) + t_0}\\
\end{array}
\end{array}
(FPCore (a b c)
:precision binary64
(let* ((t_0 (sqrt (- (* b b) (* c (* 4.0 a))))) (t_1 (/ (- c) b)))
(if (<= b -3e+136)
(if (>= b 0.0) (/ c b) t_1)
(if (<= b 8.2e+116)
(if (>= b 0.0) (/ (- (- b) t_0) (* a 2.0)) (/ (* c 2.0) (- t_0 b)))
(if (>= b 0.0) (- (/ b a)) t_1)))))
double code(double a, double b, double c) {
double t_0 = sqrt(((b * b) - (c * (4.0 * a))));
double t_1 = -c / b;
double tmp_1;
if (b <= -3e+136) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_1;
}
tmp_1 = tmp_2;
} else if (b <= 8.2e+116) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-b - t_0) / (a * 2.0);
} else {
tmp_3 = (c * 2.0) / (t_0 - b);
}
tmp_1 = tmp_3;
} else if (b >= 0.0) {
tmp_1 = -(b / a);
} else {
tmp_1 = t_1;
}
return tmp_1;
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
real(8) :: tmp_1
real(8) :: tmp_2
real(8) :: tmp_3
t_0 = sqrt(((b * b) - (c * (4.0d0 * a))))
t_1 = -c / b
if (b <= (-3d+136)) then
if (b >= 0.0d0) then
tmp_2 = c / b
else
tmp_2 = t_1
end if
tmp_1 = tmp_2
else if (b <= 8.2d+116) then
if (b >= 0.0d0) then
tmp_3 = (-b - t_0) / (a * 2.0d0)
else
tmp_3 = (c * 2.0d0) / (t_0 - b)
end if
tmp_1 = tmp_3
else if (b >= 0.0d0) then
tmp_1 = -(b / a)
else
tmp_1 = t_1
end if
code = tmp_1
end function
public static double code(double a, double b, double c) {
double t_0 = Math.sqrt(((b * b) - (c * (4.0 * a))));
double t_1 = -c / b;
double tmp_1;
if (b <= -3e+136) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_1;
}
tmp_1 = tmp_2;
} else if (b <= 8.2e+116) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-b - t_0) / (a * 2.0);
} else {
tmp_3 = (c * 2.0) / (t_0 - b);
}
tmp_1 = tmp_3;
} else if (b >= 0.0) {
tmp_1 = -(b / a);
} else {
tmp_1 = t_1;
}
return tmp_1;
}
def code(a, b, c): t_0 = math.sqrt(((b * b) - (c * (4.0 * a)))) t_1 = -c / b tmp_1 = 0 if b <= -3e+136: tmp_2 = 0 if b >= 0.0: tmp_2 = c / b else: tmp_2 = t_1 tmp_1 = tmp_2 elif b <= 8.2e+116: tmp_3 = 0 if b >= 0.0: tmp_3 = (-b - t_0) / (a * 2.0) else: tmp_3 = (c * 2.0) / (t_0 - b) tmp_1 = tmp_3 elif b >= 0.0: tmp_1 = -(b / a) else: tmp_1 = t_1 return tmp_1
function code(a, b, c) t_0 = sqrt(Float64(Float64(b * b) - Float64(c * Float64(4.0 * a)))) t_1 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -3e+136) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(c / b); else tmp_2 = t_1; end tmp_1 = tmp_2; elseif (b <= 8.2e+116) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(Float64(-b) - t_0) / Float64(a * 2.0)); else tmp_3 = Float64(Float64(c * 2.0) / Float64(t_0 - b)); end tmp_1 = tmp_3; elseif (b >= 0.0) tmp_1 = Float64(-Float64(b / a)); else tmp_1 = t_1; end return tmp_1 end
function tmp_5 = code(a, b, c) t_0 = sqrt(((b * b) - (c * (4.0 * a)))); t_1 = -c / b; tmp_2 = 0.0; if (b <= -3e+136) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = c / b; else tmp_3 = t_1; end tmp_2 = tmp_3; elseif (b <= 8.2e+116) tmp_4 = 0.0; if (b >= 0.0) tmp_4 = (-b - t_0) / (a * 2.0); else tmp_4 = (c * 2.0) / (t_0 - b); end tmp_2 = tmp_4; elseif (b >= 0.0) tmp_2 = -(b / a); else tmp_2 = t_1; end tmp_5 = tmp_2; end
code[a_, b_, c_] := Block[{t$95$0 = N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(c * N[(4.0 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -3e+136], If[GreaterEqual[b, 0.0], N[(c / b), $MachinePrecision], t$95$1], If[LessEqual[b, 8.2e+116], If[GreaterEqual[b, 0.0], N[(N[((-b) - t$95$0), $MachinePrecision] / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(c * 2.0), $MachinePrecision] / N[(t$95$0 - b), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], (-N[(b / a), $MachinePrecision]), t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{b \cdot b - c \cdot \left(4 \cdot a\right)}\\
t_1 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -3 \cdot 10^{+136}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}\\
\mathbf{elif}\;b \leq 8.2 \cdot 10^{+116}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-b\right) - t_0}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{c \cdot 2}{t_0 - b}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;-\frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if b < -2.99999999999999979e136Initial program 43.9%
Taylor expanded in b around -inf 95.6%
associate-*r/95.6%
neg-mul-195.6%
Simplified95.6%
Taylor expanded in b around inf 95.6%
Taylor expanded in b around 0 95.6%
if -2.99999999999999979e136 < b < 8.1999999999999996e116Initial program 84.3%
if 8.1999999999999996e116 < b Initial program 43.7%
Taylor expanded in b around -inf 43.7%
associate-*r/43.7%
neg-mul-143.7%
Simplified43.7%
Taylor expanded in b around inf 94.7%
Taylor expanded in b around inf 98.3%
associate-*r/98.3%
mul-1-neg98.3%
Simplified98.3%
Final simplification89.3%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (sqrt (- (* b b) (* 4.0 (* c a))))) (t_1 (/ (- c) b)))
(if (<= b -1.3e+85)
(if (>= b 0.0) (/ c b) t_1)
(if (<= b 1e+117)
(if (>= b 0.0) (/ (- (- b) t_0) (* a 2.0)) (/ 2.0 (/ (- t_0 b) c)))
(if (>= b 0.0) (- (/ b a)) t_1)))))
double code(double a, double b, double c) {
double t_0 = sqrt(((b * b) - (4.0 * (c * a))));
double t_1 = -c / b;
double tmp_1;
if (b <= -1.3e+85) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_1;
}
tmp_1 = tmp_2;
} else if (b <= 1e+117) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-b - t_0) / (a * 2.0);
} else {
tmp_3 = 2.0 / ((t_0 - b) / c);
}
tmp_1 = tmp_3;
} else if (b >= 0.0) {
tmp_1 = -(b / a);
} else {
tmp_1 = t_1;
}
return tmp_1;
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
real(8) :: tmp_1
real(8) :: tmp_2
real(8) :: tmp_3
t_0 = sqrt(((b * b) - (4.0d0 * (c * a))))
t_1 = -c / b
if (b <= (-1.3d+85)) then
if (b >= 0.0d0) then
tmp_2 = c / b
else
tmp_2 = t_1
end if
tmp_1 = tmp_2
else if (b <= 1d+117) then
if (b >= 0.0d0) then
tmp_3 = (-b - t_0) / (a * 2.0d0)
else
tmp_3 = 2.0d0 / ((t_0 - b) / c)
end if
tmp_1 = tmp_3
else if (b >= 0.0d0) then
tmp_1 = -(b / a)
else
tmp_1 = t_1
end if
code = tmp_1
end function
public static double code(double a, double b, double c) {
double t_0 = Math.sqrt(((b * b) - (4.0 * (c * a))));
double t_1 = -c / b;
double tmp_1;
if (b <= -1.3e+85) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_1;
}
tmp_1 = tmp_2;
} else if (b <= 1e+117) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-b - t_0) / (a * 2.0);
} else {
tmp_3 = 2.0 / ((t_0 - b) / c);
}
tmp_1 = tmp_3;
} else if (b >= 0.0) {
tmp_1 = -(b / a);
} else {
tmp_1 = t_1;
}
return tmp_1;
}
def code(a, b, c): t_0 = math.sqrt(((b * b) - (4.0 * (c * a)))) t_1 = -c / b tmp_1 = 0 if b <= -1.3e+85: tmp_2 = 0 if b >= 0.0: tmp_2 = c / b else: tmp_2 = t_1 tmp_1 = tmp_2 elif b <= 1e+117: tmp_3 = 0 if b >= 0.0: tmp_3 = (-b - t_0) / (a * 2.0) else: tmp_3 = 2.0 / ((t_0 - b) / c) tmp_1 = tmp_3 elif b >= 0.0: tmp_1 = -(b / a) else: tmp_1 = t_1 return tmp_1
function code(a, b, c) t_0 = sqrt(Float64(Float64(b * b) - Float64(4.0 * Float64(c * a)))) t_1 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -1.3e+85) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(c / b); else tmp_2 = t_1; end tmp_1 = tmp_2; elseif (b <= 1e+117) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(Float64(-b) - t_0) / Float64(a * 2.0)); else tmp_3 = Float64(2.0 / Float64(Float64(t_0 - b) / c)); end tmp_1 = tmp_3; elseif (b >= 0.0) tmp_1 = Float64(-Float64(b / a)); else tmp_1 = t_1; end return tmp_1 end
function tmp_5 = code(a, b, c) t_0 = sqrt(((b * b) - (4.0 * (c * a)))); t_1 = -c / b; tmp_2 = 0.0; if (b <= -1.3e+85) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = c / b; else tmp_3 = t_1; end tmp_2 = tmp_3; elseif (b <= 1e+117) tmp_4 = 0.0; if (b >= 0.0) tmp_4 = (-b - t_0) / (a * 2.0); else tmp_4 = 2.0 / ((t_0 - b) / c); end tmp_2 = tmp_4; elseif (b >= 0.0) tmp_2 = -(b / a); else tmp_2 = t_1; end tmp_5 = tmp_2; end
code[a_, b_, c_] := Block[{t$95$0 = N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(4.0 * N[(c * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -1.3e+85], If[GreaterEqual[b, 0.0], N[(c / b), $MachinePrecision], t$95$1], If[LessEqual[b, 1e+117], If[GreaterEqual[b, 0.0], N[(N[((-b) - t$95$0), $MachinePrecision] / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t$95$0 - b), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], (-N[(b / a), $MachinePrecision]), t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{b \cdot b - 4 \cdot \left(c \cdot a\right)}\\
t_1 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -1.3 \cdot 10^{+85}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}\\
\mathbf{elif}\;b \leq 10^{+117}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-b\right) - t_0}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{t_0 - b}{c}}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;-\frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if b < -1.30000000000000005e85Initial program 51.6%
Taylor expanded in b around -inf 96.2%
associate-*r/96.2%
neg-mul-196.2%
Simplified96.2%
Taylor expanded in b around inf 96.2%
Taylor expanded in b around 0 96.2%
if -1.30000000000000005e85 < b < 1.00000000000000005e117Initial program 83.6%
Simplified83.5%
if 1.00000000000000005e117 < b Initial program 43.7%
Taylor expanded in b around -inf 43.7%
associate-*r/43.7%
neg-mul-143.7%
Simplified43.7%
Taylor expanded in b around inf 94.7%
Taylor expanded in b around inf 98.3%
associate-*r/98.3%
mul-1-neg98.3%
Simplified98.3%
Final simplification89.3%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- c) b)))
(if (<= b -1.95e+77)
(if (>= b 0.0) (/ c b) t_0)
(if (<= b -2e-310)
(if (>= b 0.0)
(/ (- (+ b b)) (* a 2.0))
(/ 2.0 (/ (- (sqrt (- (* b b) (* 4.0 (* c a)))) b) c)))
(if (<= b 5.5e-39)
(if (>= b 0.0)
(/ (- (- b) (sqrt (* c (* a -4.0)))) (* a 2.0))
(/ b a))
(if (>= b 0.0) (- (/ c b) (/ b a)) t_0))))))
double code(double a, double b, double c) {
double t_0 = -c / b;
double tmp_1;
if (b <= -1.95e+77) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = -(b + b) / (a * 2.0);
} else {
tmp_3 = 2.0 / ((sqrt(((b * b) - (4.0 * (c * a)))) - b) / c);
}
tmp_1 = tmp_3;
} else if (b <= 5.5e-39) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (-b - sqrt((c * (a * -4.0)))) / (a * 2.0);
} else {
tmp_4 = b / a;
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = t_0;
}
return tmp_1;
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_0
real(8) :: tmp
real(8) :: tmp_1
real(8) :: tmp_2
real(8) :: tmp_3
real(8) :: tmp_4
t_0 = -c / b
if (b <= (-1.95d+77)) then
if (b >= 0.0d0) then
tmp_2 = c / b
else
tmp_2 = t_0
end if
tmp_1 = tmp_2
else if (b <= (-2d-310)) then
if (b >= 0.0d0) then
tmp_3 = -(b + b) / (a * 2.0d0)
else
tmp_3 = 2.0d0 / ((sqrt(((b * b) - (4.0d0 * (c * a)))) - b) / c)
end if
tmp_1 = tmp_3
else if (b <= 5.5d-39) then
if (b >= 0.0d0) then
tmp_4 = (-b - sqrt((c * (a * (-4.0d0))))) / (a * 2.0d0)
else
tmp_4 = b / a
end if
tmp_1 = tmp_4
else if (b >= 0.0d0) then
tmp_1 = (c / b) - (b / a)
else
tmp_1 = t_0
end if
code = tmp_1
end function
public static double code(double a, double b, double c) {
double t_0 = -c / b;
double tmp_1;
if (b <= -1.95e+77) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = -(b + b) / (a * 2.0);
} else {
tmp_3 = 2.0 / ((Math.sqrt(((b * b) - (4.0 * (c * a)))) - b) / c);
}
tmp_1 = tmp_3;
} else if (b <= 5.5e-39) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (-b - Math.sqrt((c * (a * -4.0)))) / (a * 2.0);
} else {
tmp_4 = b / a;
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = t_0;
}
return tmp_1;
}
def code(a, b, c): t_0 = -c / b tmp_1 = 0 if b <= -1.95e+77: tmp_2 = 0 if b >= 0.0: tmp_2 = c / b else: tmp_2 = t_0 tmp_1 = tmp_2 elif b <= -2e-310: tmp_3 = 0 if b >= 0.0: tmp_3 = -(b + b) / (a * 2.0) else: tmp_3 = 2.0 / ((math.sqrt(((b * b) - (4.0 * (c * a)))) - b) / c) tmp_1 = tmp_3 elif b <= 5.5e-39: tmp_4 = 0 if b >= 0.0: tmp_4 = (-b - math.sqrt((c * (a * -4.0)))) / (a * 2.0) else: tmp_4 = b / a tmp_1 = tmp_4 elif b >= 0.0: tmp_1 = (c / b) - (b / a) else: tmp_1 = t_0 return tmp_1
function code(a, b, c) t_0 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -1.95e+77) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(c / b); else tmp_2 = t_0; end tmp_1 = tmp_2; elseif (b <= -2e-310) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-Float64(b + b)) / Float64(a * 2.0)); else tmp_3 = Float64(2.0 / Float64(Float64(sqrt(Float64(Float64(b * b) - Float64(4.0 * Float64(c * a)))) - b) / c)); end tmp_1 = tmp_3; elseif (b <= 5.5e-39) tmp_4 = 0.0 if (b >= 0.0) tmp_4 = Float64(Float64(Float64(-b) - sqrt(Float64(c * Float64(a * -4.0)))) / Float64(a * 2.0)); else tmp_4 = Float64(b / a); end tmp_1 = tmp_4; elseif (b >= 0.0) tmp_1 = Float64(Float64(c / b) - Float64(b / a)); else tmp_1 = t_0; end return tmp_1 end
function tmp_6 = code(a, b, c) t_0 = -c / b; tmp_2 = 0.0; if (b <= -1.95e+77) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = c / b; else tmp_3 = t_0; end tmp_2 = tmp_3; elseif (b <= -2e-310) tmp_4 = 0.0; if (b >= 0.0) tmp_4 = -(b + b) / (a * 2.0); else tmp_4 = 2.0 / ((sqrt(((b * b) - (4.0 * (c * a)))) - b) / c); end tmp_2 = tmp_4; elseif (b <= 5.5e-39) tmp_5 = 0.0; if (b >= 0.0) tmp_5 = (-b - sqrt((c * (a * -4.0)))) / (a * 2.0); else tmp_5 = b / a; end tmp_2 = tmp_5; elseif (b >= 0.0) tmp_2 = (c / b) - (b / a); else tmp_2 = t_0; end tmp_6 = tmp_2; end
code[a_, b_, c_] := Block[{t$95$0 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -1.95e+77], If[GreaterEqual[b, 0.0], N[(c / b), $MachinePrecision], t$95$0], If[LessEqual[b, -2e-310], If[GreaterEqual[b, 0.0], N[((-N[(b + b), $MachinePrecision]) / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(4.0 * N[(c * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 5.5e-39], If[GreaterEqual[b, 0.0], N[(N[((-b) - N[Sqrt[N[(c * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(b / a), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -1.95 \cdot 10^{+77}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}\\
\mathbf{elif}\;b \leq -2 \cdot 10^{-310}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-\left(b + b\right)}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\sqrt{b \cdot b - 4 \cdot \left(c \cdot a\right)} - b}{c}}\\
\end{array}\\
\mathbf{elif}\;b \leq 5.5 \cdot 10^{-39}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-b\right) - \sqrt{c \cdot \left(a \cdot -4\right)}}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{b}{a}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
if b < -1.9499999999999999e77Initial program 51.6%
Taylor expanded in b around -inf 96.2%
associate-*r/96.2%
neg-mul-196.2%
Simplified96.2%
Taylor expanded in b around inf 96.2%
Taylor expanded in b around 0 96.2%
if -1.9499999999999999e77 < b < -1.999999999999994e-310Initial program 81.9%
Simplified81.7%
Taylor expanded in b around inf 81.7%
if -1.999999999999994e-310 < b < 5.50000000000000018e-39Initial program 82.3%
Simplified82.3%
Taylor expanded in b around -inf 82.3%
+-commutative82.3%
fma-def82.3%
Simplified82.3%
Taylor expanded in b around 0 73.1%
associate-*r*73.1%
metadata-eval73.1%
distribute-lft-neg-in73.1%
*-commutative73.1%
distribute-lft-neg-in73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in a around inf 73.1%
if 5.50000000000000018e-39 < b Initial program 62.1%
Taylor expanded in b around -inf 62.1%
associate-*r/62.1%
neg-mul-162.1%
Simplified62.1%
Taylor expanded in b around inf 86.2%
Taylor expanded in b around 0 88.4%
+-commutative88.4%
mul-1-neg88.4%
unsub-neg88.4%
Simplified88.4%
Final simplification85.9%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (sqrt (* c (* a -4.0)))) (t_1 (/ (- c) b)))
(if (<= b -2.15e-44)
(if (>= b 0.0) (/ c b) t_1)
(if (<= b -2e-310)
(if (>= b 0.0) (/ (- (+ b b)) (* a 2.0)) (/ 2.0 (/ (- t_0 b) c)))
(if (<= b 1.1e-33)
(if (>= b 0.0)
(/ (- (- b) t_0) (* a 2.0))
(/ 2.0 (fma 2.0 (/ a b) (* -2.0 (/ b c)))))
(if (>= b 0.0) (- (/ c b) (/ b a)) t_1))))))
double code(double a, double b, double c) {
double t_0 = sqrt((c * (a * -4.0)));
double t_1 = -c / b;
double tmp_1;
if (b <= -2.15e-44) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_1;
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = -(b + b) / (a * 2.0);
} else {
tmp_3 = 2.0 / ((t_0 - b) / c);
}
tmp_1 = tmp_3;
} else if (b <= 1.1e-33) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (-b - t_0) / (a * 2.0);
} else {
tmp_4 = 2.0 / fma(2.0, (a / b), (-2.0 * (b / c)));
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = t_1;
}
return tmp_1;
}
function code(a, b, c) t_0 = sqrt(Float64(c * Float64(a * -4.0))) t_1 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -2.15e-44) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(c / b); else tmp_2 = t_1; end tmp_1 = tmp_2; elseif (b <= -2e-310) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-Float64(b + b)) / Float64(a * 2.0)); else tmp_3 = Float64(2.0 / Float64(Float64(t_0 - b) / c)); end tmp_1 = tmp_3; elseif (b <= 1.1e-33) tmp_4 = 0.0 if (b >= 0.0) tmp_4 = Float64(Float64(Float64(-b) - t_0) / Float64(a * 2.0)); else tmp_4 = Float64(2.0 / fma(2.0, Float64(a / b), Float64(-2.0 * Float64(b / c)))); end tmp_1 = tmp_4; elseif (b >= 0.0) tmp_1 = Float64(Float64(c / b) - Float64(b / a)); else tmp_1 = t_1; end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[Sqrt[N[(c * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -2.15e-44], If[GreaterEqual[b, 0.0], N[(c / b), $MachinePrecision], t$95$1], If[LessEqual[b, -2e-310], If[GreaterEqual[b, 0.0], N[((-N[(b + b), $MachinePrecision]) / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t$95$0 - b), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 1.1e-33], If[GreaterEqual[b, 0.0], N[(N[((-b) - t$95$0), $MachinePrecision] / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(2.0 * N[(a / b), $MachinePrecision] + N[(-2.0 * N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{c \cdot \left(a \cdot -4\right)}\\
t_1 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -2.15 \cdot 10^{-44}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}\\
\mathbf{elif}\;b \leq -2 \cdot 10^{-310}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-\left(b + b\right)}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{t_0 - b}{c}}\\
\end{array}\\
\mathbf{elif}\;b \leq 1.1 \cdot 10^{-33}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-b\right) - t_0}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(2, \frac{a}{b}, -2 \cdot \frac{b}{c}\right)}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if b < -2.15000000000000007e-44Initial program 62.7%
Taylor expanded in b around -inf 87.2%
associate-*r/87.2%
neg-mul-187.2%
Simplified87.2%
Taylor expanded in b around inf 87.2%
Taylor expanded in b around 0 87.2%
if -2.15000000000000007e-44 < b < -1.999999999999994e-310Initial program 79.3%
Simplified79.2%
Taylor expanded in b around inf 79.2%
Taylor expanded in b around 0 68.5%
associate-*r*23.6%
metadata-eval23.6%
distribute-lft-neg-in23.6%
*-commutative23.6%
distribute-lft-neg-in23.6%
metadata-eval23.6%
Simplified68.5%
if -1.999999999999994e-310 < b < 1.10000000000000003e-33Initial program 82.3%
Simplified82.3%
Taylor expanded in b around -inf 82.3%
+-commutative82.3%
fma-def82.3%
Simplified82.3%
Taylor expanded in b around 0 73.1%
associate-*r*73.1%
metadata-eval73.1%
distribute-lft-neg-in73.1%
*-commutative73.1%
distribute-lft-neg-in73.1%
metadata-eval73.1%
Simplified73.1%
if 1.10000000000000003e-33 < b Initial program 62.1%
Taylor expanded in b around -inf 62.1%
associate-*r/62.1%
neg-mul-162.1%
Simplified62.1%
Taylor expanded in b around inf 86.2%
Taylor expanded in b around 0 88.4%
+-commutative88.4%
mul-1-neg88.4%
unsub-neg88.4%
Simplified88.4%
Final simplification82.0%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (sqrt (* c (* a -4.0)))) (t_1 (/ (- c) b)))
(if (<= b -6.5e-43)
(if (>= b 0.0) (/ c b) t_1)
(if (<= b -2e-310)
(if (>= b 0.0) (/ (- (+ b b)) (* a 2.0)) (/ 2.0 (/ (- t_0 b) c)))
(if (<= b 1.15e-37)
(if (>= b 0.0) (/ (- (- b) t_0) (* a 2.0)) (/ b a))
(if (>= b 0.0) (- (/ c b) (/ b a)) t_1))))))
double code(double a, double b, double c) {
double t_0 = sqrt((c * (a * -4.0)));
double t_1 = -c / b;
double tmp_1;
if (b <= -6.5e-43) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_1;
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = -(b + b) / (a * 2.0);
} else {
tmp_3 = 2.0 / ((t_0 - b) / c);
}
tmp_1 = tmp_3;
} else if (b <= 1.15e-37) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (-b - t_0) / (a * 2.0);
} else {
tmp_4 = b / a;
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = t_1;
}
return tmp_1;
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
real(8) :: tmp_1
real(8) :: tmp_2
real(8) :: tmp_3
real(8) :: tmp_4
t_0 = sqrt((c * (a * (-4.0d0))))
t_1 = -c / b
if (b <= (-6.5d-43)) then
if (b >= 0.0d0) then
tmp_2 = c / b
else
tmp_2 = t_1
end if
tmp_1 = tmp_2
else if (b <= (-2d-310)) then
if (b >= 0.0d0) then
tmp_3 = -(b + b) / (a * 2.0d0)
else
tmp_3 = 2.0d0 / ((t_0 - b) / c)
end if
tmp_1 = tmp_3
else if (b <= 1.15d-37) then
if (b >= 0.0d0) then
tmp_4 = (-b - t_0) / (a * 2.0d0)
else
tmp_4 = b / a
end if
tmp_1 = tmp_4
else if (b >= 0.0d0) then
tmp_1 = (c / b) - (b / a)
else
tmp_1 = t_1
end if
code = tmp_1
end function
public static double code(double a, double b, double c) {
double t_0 = Math.sqrt((c * (a * -4.0)));
double t_1 = -c / b;
double tmp_1;
if (b <= -6.5e-43) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_1;
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = -(b + b) / (a * 2.0);
} else {
tmp_3 = 2.0 / ((t_0 - b) / c);
}
tmp_1 = tmp_3;
} else if (b <= 1.15e-37) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (-b - t_0) / (a * 2.0);
} else {
tmp_4 = b / a;
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = t_1;
}
return tmp_1;
}
def code(a, b, c): t_0 = math.sqrt((c * (a * -4.0))) t_1 = -c / b tmp_1 = 0 if b <= -6.5e-43: tmp_2 = 0 if b >= 0.0: tmp_2 = c / b else: tmp_2 = t_1 tmp_1 = tmp_2 elif b <= -2e-310: tmp_3 = 0 if b >= 0.0: tmp_3 = -(b + b) / (a * 2.0) else: tmp_3 = 2.0 / ((t_0 - b) / c) tmp_1 = tmp_3 elif b <= 1.15e-37: tmp_4 = 0 if b >= 0.0: tmp_4 = (-b - t_0) / (a * 2.0) else: tmp_4 = b / a tmp_1 = tmp_4 elif b >= 0.0: tmp_1 = (c / b) - (b / a) else: tmp_1 = t_1 return tmp_1
function code(a, b, c) t_0 = sqrt(Float64(c * Float64(a * -4.0))) t_1 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -6.5e-43) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(c / b); else tmp_2 = t_1; end tmp_1 = tmp_2; elseif (b <= -2e-310) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-Float64(b + b)) / Float64(a * 2.0)); else tmp_3 = Float64(2.0 / Float64(Float64(t_0 - b) / c)); end tmp_1 = tmp_3; elseif (b <= 1.15e-37) tmp_4 = 0.0 if (b >= 0.0) tmp_4 = Float64(Float64(Float64(-b) - t_0) / Float64(a * 2.0)); else tmp_4 = Float64(b / a); end tmp_1 = tmp_4; elseif (b >= 0.0) tmp_1 = Float64(Float64(c / b) - Float64(b / a)); else tmp_1 = t_1; end return tmp_1 end
function tmp_6 = code(a, b, c) t_0 = sqrt((c * (a * -4.0))); t_1 = -c / b; tmp_2 = 0.0; if (b <= -6.5e-43) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = c / b; else tmp_3 = t_1; end tmp_2 = tmp_3; elseif (b <= -2e-310) tmp_4 = 0.0; if (b >= 0.0) tmp_4 = -(b + b) / (a * 2.0); else tmp_4 = 2.0 / ((t_0 - b) / c); end tmp_2 = tmp_4; elseif (b <= 1.15e-37) tmp_5 = 0.0; if (b >= 0.0) tmp_5 = (-b - t_0) / (a * 2.0); else tmp_5 = b / a; end tmp_2 = tmp_5; elseif (b >= 0.0) tmp_2 = (c / b) - (b / a); else tmp_2 = t_1; end tmp_6 = tmp_2; end
code[a_, b_, c_] := Block[{t$95$0 = N[Sqrt[N[(c * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -6.5e-43], If[GreaterEqual[b, 0.0], N[(c / b), $MachinePrecision], t$95$1], If[LessEqual[b, -2e-310], If[GreaterEqual[b, 0.0], N[((-N[(b + b), $MachinePrecision]) / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t$95$0 - b), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 1.15e-37], If[GreaterEqual[b, 0.0], N[(N[((-b) - t$95$0), $MachinePrecision] / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(b / a), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{c \cdot \left(a \cdot -4\right)}\\
t_1 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -6.5 \cdot 10^{-43}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}\\
\mathbf{elif}\;b \leq -2 \cdot 10^{-310}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-\left(b + b\right)}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{t_0 - b}{c}}\\
\end{array}\\
\mathbf{elif}\;b \leq 1.15 \cdot 10^{-37}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-b\right) - t_0}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{b}{a}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if b < -6.50000000000000001e-43Initial program 62.7%
Taylor expanded in b around -inf 87.2%
associate-*r/87.2%
neg-mul-187.2%
Simplified87.2%
Taylor expanded in b around inf 87.2%
Taylor expanded in b around 0 87.2%
if -6.50000000000000001e-43 < b < -1.999999999999994e-310Initial program 79.3%
Simplified79.2%
Taylor expanded in b around inf 79.2%
Taylor expanded in b around 0 68.5%
associate-*r*23.6%
metadata-eval23.6%
distribute-lft-neg-in23.6%
*-commutative23.6%
distribute-lft-neg-in23.6%
metadata-eval23.6%
Simplified68.5%
if -1.999999999999994e-310 < b < 1.15e-37Initial program 82.3%
Simplified82.3%
Taylor expanded in b around -inf 82.3%
+-commutative82.3%
fma-def82.3%
Simplified82.3%
Taylor expanded in b around 0 73.1%
associate-*r*73.1%
metadata-eval73.1%
distribute-lft-neg-in73.1%
*-commutative73.1%
distribute-lft-neg-in73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in a around inf 73.1%
if 1.15e-37 < b Initial program 62.1%
Taylor expanded in b around -inf 62.1%
associate-*r/62.1%
neg-mul-162.1%
Simplified62.1%
Taylor expanded in b around inf 86.2%
Taylor expanded in b around 0 88.4%
+-commutative88.4%
mul-1-neg88.4%
unsub-neg88.4%
Simplified88.4%
Final simplification82.0%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (sqrt (* a (* c -4.0)))) (t_1 (/ (- c) b)))
(if (<= b -2.35e-92)
(if (>= b 0.0) (/ c b) t_1)
(if (<= b -2e-310)
(if (>= b 0.0) (/ (- (+ b b)) (* a 2.0)) (* c (/ 2.0 (+ b t_0))))
(if (<= b 4.2e-37)
(if (>= b 0.0)
(/ (* 0.5 (- b t_0)) a)
(/ 2.0 (fma 2.0 (/ a b) (* -2.0 (/ b c)))))
(if (>= b 0.0) (- (/ c b) (/ b a)) t_1))))))
double code(double a, double b, double c) {
double t_0 = sqrt((a * (c * -4.0)));
double t_1 = -c / b;
double tmp_1;
if (b <= -2.35e-92) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_1;
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = -(b + b) / (a * 2.0);
} else {
tmp_3 = c * (2.0 / (b + t_0));
}
tmp_1 = tmp_3;
} else if (b <= 4.2e-37) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (0.5 * (b - t_0)) / a;
} else {
tmp_4 = 2.0 / fma(2.0, (a / b), (-2.0 * (b / c)));
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = t_1;
}
return tmp_1;
}
function code(a, b, c) t_0 = sqrt(Float64(a * Float64(c * -4.0))) t_1 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -2.35e-92) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(c / b); else tmp_2 = t_1; end tmp_1 = tmp_2; elseif (b <= -2e-310) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-Float64(b + b)) / Float64(a * 2.0)); else tmp_3 = Float64(c * Float64(2.0 / Float64(b + t_0))); end tmp_1 = tmp_3; elseif (b <= 4.2e-37) tmp_4 = 0.0 if (b >= 0.0) tmp_4 = Float64(Float64(0.5 * Float64(b - t_0)) / a); else tmp_4 = Float64(2.0 / fma(2.0, Float64(a / b), Float64(-2.0 * Float64(b / c)))); end tmp_1 = tmp_4; elseif (b >= 0.0) tmp_1 = Float64(Float64(c / b) - Float64(b / a)); else tmp_1 = t_1; end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[Sqrt[N[(a * N[(c * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -2.35e-92], If[GreaterEqual[b, 0.0], N[(c / b), $MachinePrecision], t$95$1], If[LessEqual[b, -2e-310], If[GreaterEqual[b, 0.0], N[((-N[(b + b), $MachinePrecision]) / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(c * N[(2.0 / N[(b + t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 4.2e-37], If[GreaterEqual[b, 0.0], N[(N[(0.5 * N[(b - t$95$0), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], N[(2.0 / N[(2.0 * N[(a / b), $MachinePrecision] + N[(-2.0 * N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{a \cdot \left(c \cdot -4\right)}\\
t_1 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -2.35 \cdot 10^{-92}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}\\
\mathbf{elif}\;b \leq -2 \cdot 10^{-310}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-\left(b + b\right)}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;c \cdot \frac{2}{b + t_0}\\
\end{array}\\
\mathbf{elif}\;b \leq 4.2 \cdot 10^{-37}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{0.5 \cdot \left(b - t_0\right)}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(2, \frac{a}{b}, -2 \cdot \frac{b}{c}\right)}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if b < -2.34999999999999996e-92Initial program 67.4%
Taylor expanded in b around -inf 82.6%
associate-*r/82.6%
neg-mul-182.6%
Simplified82.6%
Taylor expanded in b around inf 82.6%
Taylor expanded in b around 0 82.6%
if -2.34999999999999996e-92 < b < -1.999999999999994e-310Initial program 73.5%
Simplified73.3%
Taylor expanded in b around inf 73.3%
Taylor expanded in b around 0 73.3%
associate-*r*14.8%
metadata-eval14.8%
distribute-lft-neg-in14.8%
*-commutative14.8%
distribute-lft-neg-in14.8%
metadata-eval14.8%
Simplified73.3%
associate-/r/73.2%
add-sqr-sqrt73.2%
sqrt-unprod71.5%
sqr-neg71.5%
sqrt-unprod0.0%
add-sqr-sqrt71.3%
associate-*r*71.3%
Applied egg-rr71.3%
if -1.999999999999994e-310 < b < 4.2000000000000002e-37Initial program 82.3%
Simplified82.3%
Taylor expanded in b around -inf 82.3%
+-commutative82.3%
fma-def82.3%
Simplified82.3%
Taylor expanded in b around 0 73.1%
associate-*r*73.1%
metadata-eval73.1%
distribute-lft-neg-in73.1%
*-commutative73.1%
distribute-lft-neg-in73.1%
metadata-eval73.1%
Simplified73.1%
div-sub73.1%
sub-neg73.1%
*-un-lft-identity73.1%
*-commutative73.1%
times-frac73.1%
metadata-eval73.1%
add-sqr-sqrt0.0%
sqrt-unprod69.8%
sqr-neg69.8%
sqrt-unprod69.6%
add-sqr-sqrt69.6%
*-un-lft-identity69.6%
*-commutative69.6%
times-frac69.6%
metadata-eval69.6%
associate-*r*69.6%
*-commutative69.6%
Applied egg-rr69.6%
sub-neg69.6%
distribute-lft-out--69.6%
div-sub69.6%
associate-*r/69.6%
Simplified69.6%
if 4.2000000000000002e-37 < b Initial program 62.1%
Taylor expanded in b around -inf 62.1%
associate-*r/62.1%
neg-mul-162.1%
Simplified62.1%
Taylor expanded in b around inf 86.2%
Taylor expanded in b around 0 88.4%
+-commutative88.4%
mul-1-neg88.4%
unsub-neg88.4%
Simplified88.4%
Final simplification81.2%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- c) b)))
(if (<= b -3.9e-94)
(if (>= b 0.0) (/ c b) t_0)
(if (<= b -2e-310)
(if (>= b 0.0)
(/ (- (+ b b)) (* a 2.0))
(/ 2.0 (/ (+ b (sqrt (* c (* a -4.0)))) c)))
(if (<= b 4.4e-33)
(if (>= b 0.0)
(/ (* 0.5 (- b (sqrt (* a (* c -4.0))))) a)
(/ 2.0 (fma 2.0 (/ a b) (* -2.0 (/ b c)))))
(if (>= b 0.0) (- (/ c b) (/ b a)) t_0))))))
double code(double a, double b, double c) {
double t_0 = -c / b;
double tmp_1;
if (b <= -3.9e-94) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = -(b + b) / (a * 2.0);
} else {
tmp_3 = 2.0 / ((b + sqrt((c * (a * -4.0)))) / c);
}
tmp_1 = tmp_3;
} else if (b <= 4.4e-33) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (0.5 * (b - sqrt((a * (c * -4.0))))) / a;
} else {
tmp_4 = 2.0 / fma(2.0, (a / b), (-2.0 * (b / c)));
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = t_0;
}
return tmp_1;
}
function code(a, b, c) t_0 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -3.9e-94) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(c / b); else tmp_2 = t_0; end tmp_1 = tmp_2; elseif (b <= -2e-310) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-Float64(b + b)) / Float64(a * 2.0)); else tmp_3 = Float64(2.0 / Float64(Float64(b + sqrt(Float64(c * Float64(a * -4.0)))) / c)); end tmp_1 = tmp_3; elseif (b <= 4.4e-33) tmp_4 = 0.0 if (b >= 0.0) tmp_4 = Float64(Float64(0.5 * Float64(b - sqrt(Float64(a * Float64(c * -4.0))))) / a); else tmp_4 = Float64(2.0 / fma(2.0, Float64(a / b), Float64(-2.0 * Float64(b / c)))); end tmp_1 = tmp_4; elseif (b >= 0.0) tmp_1 = Float64(Float64(c / b) - Float64(b / a)); else tmp_1 = t_0; end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -3.9e-94], If[GreaterEqual[b, 0.0], N[(c / b), $MachinePrecision], t$95$0], If[LessEqual[b, -2e-310], If[GreaterEqual[b, 0.0], N[((-N[(b + b), $MachinePrecision]) / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(b + N[Sqrt[N[(c * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 4.4e-33], If[GreaterEqual[b, 0.0], N[(N[(0.5 * N[(b - N[Sqrt[N[(a * N[(c * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], N[(2.0 / N[(2.0 * N[(a / b), $MachinePrecision] + N[(-2.0 * N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -3.9 \cdot 10^{-94}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}\\
\mathbf{elif}\;b \leq -2 \cdot 10^{-310}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-\left(b + b\right)}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{b + \sqrt{c \cdot \left(a \cdot -4\right)}}{c}}\\
\end{array}\\
\mathbf{elif}\;b \leq 4.4 \cdot 10^{-33}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{0.5 \cdot \left(b - \sqrt{a \cdot \left(c \cdot -4\right)}\right)}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(2, \frac{a}{b}, -2 \cdot \frac{b}{c}\right)}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
if b < -3.9000000000000002e-94Initial program 67.4%
Taylor expanded in b around -inf 82.6%
associate-*r/82.6%
neg-mul-182.6%
Simplified82.6%
Taylor expanded in b around inf 82.6%
Taylor expanded in b around 0 82.6%
if -3.9000000000000002e-94 < b < -1.999999999999994e-310Initial program 73.5%
Simplified73.3%
Taylor expanded in b around inf 73.3%
Taylor expanded in b around 0 73.3%
associate-*r*14.8%
metadata-eval14.8%
distribute-lft-neg-in14.8%
*-commutative14.8%
distribute-lft-neg-in14.8%
metadata-eval14.8%
Simplified73.3%
expm1-log1p-u50.6%
expm1-udef26.0%
add-sqr-sqrt26.0%
sqrt-unprod24.8%
sqr-neg24.8%
sqrt-unprod0.0%
add-sqr-sqrt24.6%
associate-*r*24.6%
Applied egg-rr24.6%
expm1-def48.5%
expm1-log1p71.3%
associate-*l*71.3%
*-commutative71.3%
Simplified71.3%
if -1.999999999999994e-310 < b < 4.40000000000000011e-33Initial program 82.3%
Simplified82.3%
Taylor expanded in b around -inf 82.3%
+-commutative82.3%
fma-def82.3%
Simplified82.3%
Taylor expanded in b around 0 73.1%
associate-*r*73.1%
metadata-eval73.1%
distribute-lft-neg-in73.1%
*-commutative73.1%
distribute-lft-neg-in73.1%
metadata-eval73.1%
Simplified73.1%
div-sub73.1%
sub-neg73.1%
*-un-lft-identity73.1%
*-commutative73.1%
times-frac73.1%
metadata-eval73.1%
add-sqr-sqrt0.0%
sqrt-unprod69.8%
sqr-neg69.8%
sqrt-unprod69.6%
add-sqr-sqrt69.6%
*-un-lft-identity69.6%
*-commutative69.6%
times-frac69.6%
metadata-eval69.6%
associate-*r*69.6%
*-commutative69.6%
Applied egg-rr69.6%
sub-neg69.6%
distribute-lft-out--69.6%
div-sub69.6%
associate-*r/69.6%
Simplified69.6%
if 4.40000000000000011e-33 < b Initial program 62.1%
Taylor expanded in b around -inf 62.1%
associate-*r/62.1%
neg-mul-162.1%
Simplified62.1%
Taylor expanded in b around inf 86.2%
Taylor expanded in b around 0 88.4%
+-commutative88.4%
mul-1-neg88.4%
unsub-neg88.4%
Simplified88.4%
Final simplification81.2%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- c) b)))
(if (<= b -1.2e-94)
(if (>= b 0.0) (/ c b) t_0)
(if (<= b -2e-310)
(if (>= b 0.0)
(/ (- (+ b b)) (* a 2.0))
(/ (* c 2.0) (+ b (sqrt (* c (* a -4.0))))))
(if (<= b 4.2e-39)
(if (>= b 0.0)
(/ (* 0.5 (- b (sqrt (* a (* c -4.0))))) a)
(/ 2.0 (fma 2.0 (/ a b) (* -2.0 (/ b c)))))
(if (>= b 0.0) (- (/ c b) (/ b a)) t_0))))))
double code(double a, double b, double c) {
double t_0 = -c / b;
double tmp_1;
if (b <= -1.2e-94) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = -(b + b) / (a * 2.0);
} else {
tmp_3 = (c * 2.0) / (b + sqrt((c * (a * -4.0))));
}
tmp_1 = tmp_3;
} else if (b <= 4.2e-39) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (0.5 * (b - sqrt((a * (c * -4.0))))) / a;
} else {
tmp_4 = 2.0 / fma(2.0, (a / b), (-2.0 * (b / c)));
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = t_0;
}
return tmp_1;
}
function code(a, b, c) t_0 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -1.2e-94) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(c / b); else tmp_2 = t_0; end tmp_1 = tmp_2; elseif (b <= -2e-310) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-Float64(b + b)) / Float64(a * 2.0)); else tmp_3 = Float64(Float64(c * 2.0) / Float64(b + sqrt(Float64(c * Float64(a * -4.0))))); end tmp_1 = tmp_3; elseif (b <= 4.2e-39) tmp_4 = 0.0 if (b >= 0.0) tmp_4 = Float64(Float64(0.5 * Float64(b - sqrt(Float64(a * Float64(c * -4.0))))) / a); else tmp_4 = Float64(2.0 / fma(2.0, Float64(a / b), Float64(-2.0 * Float64(b / c)))); end tmp_1 = tmp_4; elseif (b >= 0.0) tmp_1 = Float64(Float64(c / b) - Float64(b / a)); else tmp_1 = t_0; end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -1.2e-94], If[GreaterEqual[b, 0.0], N[(c / b), $MachinePrecision], t$95$0], If[LessEqual[b, -2e-310], If[GreaterEqual[b, 0.0], N[((-N[(b + b), $MachinePrecision]) / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(c * 2.0), $MachinePrecision] / N[(b + N[Sqrt[N[(c * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 4.2e-39], If[GreaterEqual[b, 0.0], N[(N[(0.5 * N[(b - N[Sqrt[N[(a * N[(c * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], N[(2.0 / N[(2.0 * N[(a / b), $MachinePrecision] + N[(-2.0 * N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -1.2 \cdot 10^{-94}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}\\
\mathbf{elif}\;b \leq -2 \cdot 10^{-310}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-\left(b + b\right)}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{c \cdot 2}{b + \sqrt{c \cdot \left(a \cdot -4\right)}}\\
\end{array}\\
\mathbf{elif}\;b \leq 4.2 \cdot 10^{-39}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{0.5 \cdot \left(b - \sqrt{a \cdot \left(c \cdot -4\right)}\right)}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(2, \frac{a}{b}, -2 \cdot \frac{b}{c}\right)}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
if b < -1.2e-94Initial program 67.4%
Taylor expanded in b around -inf 82.6%
associate-*r/82.6%
neg-mul-182.6%
Simplified82.6%
Taylor expanded in b around inf 82.6%
Taylor expanded in b around 0 82.6%
if -1.2e-94 < b < -1.999999999999994e-310Initial program 73.5%
Simplified73.3%
Taylor expanded in b around inf 73.3%
Taylor expanded in b around 0 73.3%
associate-*r*14.8%
metadata-eval14.8%
distribute-lft-neg-in14.8%
*-commutative14.8%
distribute-lft-neg-in14.8%
metadata-eval14.8%
Simplified73.3%
expm1-log1p-u51.6%
expm1-udef10.1%
div-inv10.1%
clear-num10.1%
add-sqr-sqrt10.1%
sqrt-unprod9.2%
sqr-neg9.2%
sqrt-unprod0.0%
add-sqr-sqrt9.0%
associate-*r*9.0%
Applied egg-rr9.0%
expm1-def50.1%
expm1-log1p71.5%
associate-*r/71.5%
associate-*l*71.5%
*-commutative71.5%
Simplified71.5%
if -1.999999999999994e-310 < b < 4.19999999999999987e-39Initial program 82.3%
Simplified82.3%
Taylor expanded in b around -inf 82.3%
+-commutative82.3%
fma-def82.3%
Simplified82.3%
Taylor expanded in b around 0 73.1%
associate-*r*73.1%
metadata-eval73.1%
distribute-lft-neg-in73.1%
*-commutative73.1%
distribute-lft-neg-in73.1%
metadata-eval73.1%
Simplified73.1%
div-sub73.1%
sub-neg73.1%
*-un-lft-identity73.1%
*-commutative73.1%
times-frac73.1%
metadata-eval73.1%
add-sqr-sqrt0.0%
sqrt-unprod69.8%
sqr-neg69.8%
sqrt-unprod69.6%
add-sqr-sqrt69.6%
*-un-lft-identity69.6%
*-commutative69.6%
times-frac69.6%
metadata-eval69.6%
associate-*r*69.6%
*-commutative69.6%
Applied egg-rr69.6%
sub-neg69.6%
distribute-lft-out--69.6%
div-sub69.6%
associate-*r/69.6%
Simplified69.6%
if 4.19999999999999987e-39 < b Initial program 62.1%
Taylor expanded in b around -inf 62.1%
associate-*r/62.1%
neg-mul-162.1%
Simplified62.1%
Taylor expanded in b around inf 86.2%
Taylor expanded in b around 0 88.4%
+-commutative88.4%
mul-1-neg88.4%
unsub-neg88.4%
Simplified88.4%
Final simplification81.3%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- c) b)))
(if (<= b -5.4e-48)
(if (>= b 0.0) (/ c b) t_0)
(if (<= b -2e-310)
(if (>= b 0.0)
(/ (- (+ b b)) (* a 2.0))
(/ 2.0 (/ (- (sqrt (* c (* a -4.0))) b) c)))
(if (<= b 8.5e-39)
(if (>= b 0.0)
(/ (* 0.5 (- b (sqrt (* a (* c -4.0))))) a)
(/ 2.0 (fma 2.0 (/ a b) (* -2.0 (/ b c)))))
(if (>= b 0.0) (- (/ c b) (/ b a)) t_0))))))
double code(double a, double b, double c) {
double t_0 = -c / b;
double tmp_1;
if (b <= -5.4e-48) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = c / b;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = -(b + b) / (a * 2.0);
} else {
tmp_3 = 2.0 / ((sqrt((c * (a * -4.0))) - b) / c);
}
tmp_1 = tmp_3;
} else if (b <= 8.5e-39) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (0.5 * (b - sqrt((a * (c * -4.0))))) / a;
} else {
tmp_4 = 2.0 / fma(2.0, (a / b), (-2.0 * (b / c)));
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = t_0;
}
return tmp_1;
}
function code(a, b, c) t_0 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -5.4e-48) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(c / b); else tmp_2 = t_0; end tmp_1 = tmp_2; elseif (b <= -2e-310) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-Float64(b + b)) / Float64(a * 2.0)); else tmp_3 = Float64(2.0 / Float64(Float64(sqrt(Float64(c * Float64(a * -4.0))) - b) / c)); end tmp_1 = tmp_3; elseif (b <= 8.5e-39) tmp_4 = 0.0 if (b >= 0.0) tmp_4 = Float64(Float64(0.5 * Float64(b - sqrt(Float64(a * Float64(c * -4.0))))) / a); else tmp_4 = Float64(2.0 / fma(2.0, Float64(a / b), Float64(-2.0 * Float64(b / c)))); end tmp_1 = tmp_4; elseif (b >= 0.0) tmp_1 = Float64(Float64(c / b) - Float64(b / a)); else tmp_1 = t_0; end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -5.4e-48], If[GreaterEqual[b, 0.0], N[(c / b), $MachinePrecision], t$95$0], If[LessEqual[b, -2e-310], If[GreaterEqual[b, 0.0], N[((-N[(b + b), $MachinePrecision]) / N[(a * 2.0), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Sqrt[N[(c * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 8.5e-39], If[GreaterEqual[b, 0.0], N[(N[(0.5 * N[(b - N[Sqrt[N[(a * N[(c * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], N[(2.0 / N[(2.0 * N[(a / b), $MachinePrecision] + N[(-2.0 * N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -5.4 \cdot 10^{-48}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}\\
\mathbf{elif}\;b \leq -2 \cdot 10^{-310}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-\left(b + b\right)}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\sqrt{c \cdot \left(a \cdot -4\right)} - b}{c}}\\
\end{array}\\
\mathbf{elif}\;b \leq 8.5 \cdot 10^{-39}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{0.5 \cdot \left(b - \sqrt{a \cdot \left(c \cdot -4\right)}\right)}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(2, \frac{a}{b}, -2 \cdot \frac{b}{c}\right)}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
if b < -5.40000000000000023e-48Initial program 62.7%
Taylor expanded in b around -inf 87.2%
associate-*r/87.2%
neg-mul-187.2%
Simplified87.2%
Taylor expanded in b around inf 87.2%
Taylor expanded in b around 0 87.2%
if -5.40000000000000023e-48 < b < -1.999999999999994e-310Initial program 79.3%
Simplified79.2%
Taylor expanded in b around inf 79.2%
Taylor expanded in b around 0 68.5%
associate-*r*23.6%
metadata-eval23.6%
distribute-lft-neg-in23.6%
*-commutative23.6%
distribute-lft-neg-in23.6%
metadata-eval23.6%
Simplified68.5%
if -1.999999999999994e-310 < b < 8.5000000000000005e-39Initial program 82.3%
Simplified82.3%
Taylor expanded in b around -inf 82.3%
+-commutative82.3%
fma-def82.3%
Simplified82.3%
Taylor expanded in b around 0 73.1%
associate-*r*73.1%
metadata-eval73.1%
distribute-lft-neg-in73.1%
*-commutative73.1%
distribute-lft-neg-in73.1%
metadata-eval73.1%
Simplified73.1%
div-sub73.1%
sub-neg73.1%
*-un-lft-identity73.1%
*-commutative73.1%
times-frac73.1%
metadata-eval73.1%
add-sqr-sqrt0.0%
sqrt-unprod69.8%
sqr-neg69.8%
sqrt-unprod69.6%
add-sqr-sqrt69.6%
*-un-lft-identity69.6%
*-commutative69.6%
times-frac69.6%
metadata-eval69.6%
associate-*r*69.6%
*-commutative69.6%
Applied egg-rr69.6%
sub-neg69.6%
distribute-lft-out--69.6%
div-sub69.6%
associate-*r/69.6%
Simplified69.6%
if 8.5000000000000005e-39 < b Initial program 62.1%
Taylor expanded in b around -inf 62.1%
associate-*r/62.1%
neg-mul-162.1%
Simplified62.1%
Taylor expanded in b around inf 86.2%
Taylor expanded in b around 0 88.4%
+-commutative88.4%
mul-1-neg88.4%
unsub-neg88.4%
Simplified88.4%
Final simplification81.5%
(FPCore (a b c)
:precision binary64
(if (<= b 1.32e-33)
(if (>= b 0.0)
(* (- b (sqrt (* a (* c -4.0)))) (/ 0.5 a))
(/ 2.0 (fma 2.0 (/ a b) (* -2.0 (/ b c)))))
(if (>= b 0.0) (- (/ c b) (/ b a)) (/ (- c) b))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= 1.32e-33) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = (b - sqrt((a * (c * -4.0)))) * (0.5 / a);
} else {
tmp_2 = 2.0 / fma(2.0, (a / b), (-2.0 * (b / c)));
}
tmp_1 = tmp_2;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = -c / b;
}
return tmp_1;
}
function code(a, b, c) tmp_1 = 0.0 if (b <= 1.32e-33) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(b - sqrt(Float64(a * Float64(c * -4.0)))) * Float64(0.5 / a)); else tmp_2 = Float64(2.0 / fma(2.0, Float64(a / b), Float64(-2.0 * Float64(b / c)))); end tmp_1 = tmp_2; elseif (b >= 0.0) tmp_1 = Float64(Float64(c / b) - Float64(b / a)); else tmp_1 = Float64(Float64(-c) / b); end return tmp_1 end
code[a_, b_, c_] := If[LessEqual[b, 1.32e-33], If[GreaterEqual[b, 0.0], N[(N[(b - N[Sqrt[N[(a * N[(c * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(0.5 / a), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(2.0 * N[(a / b), $MachinePrecision] + N[(-2.0 * N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision], N[((-c) / b), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 1.32 \cdot 10^{-33}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\left(b - \sqrt{a \cdot \left(c \cdot -4\right)}\right) \cdot \frac{0.5}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(2, \frac{a}{b}, -2 \cdot \frac{b}{c}\right)}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{-c}{b}\\
\end{array}
\end{array}
if b < 1.31999999999999993e-33Initial program 72.3%
Simplified72.2%
Taylor expanded in b around -inf 66.1%
+-commutative66.1%
fma-def66.1%
Simplified66.1%
Taylor expanded in b around 0 64.0%
associate-*r*64.0%
metadata-eval64.0%
distribute-lft-neg-in64.0%
*-commutative64.0%
distribute-lft-neg-in64.0%
metadata-eval64.0%
Simplified64.0%
div-sub64.0%
sub-neg64.0%
*-un-lft-identity64.0%
*-commutative64.0%
times-frac64.0%
metadata-eval64.0%
add-sqr-sqrt47.1%
sqrt-unprod63.2%
sqr-neg63.2%
sqrt-unprod63.2%
add-sqr-sqrt63.2%
*-un-lft-identity63.2%
*-commutative63.2%
times-frac63.2%
metadata-eval63.2%
associate-*r*63.2%
*-commutative63.2%
Applied egg-rr63.2%
sub-neg63.2%
distribute-lft-out--63.2%
div-sub63.2%
associate-*r/63.2%
associate-*l/63.1%
Simplified63.1%
if 1.31999999999999993e-33 < b Initial program 62.1%
Taylor expanded in b around -inf 62.1%
associate-*r/62.1%
neg-mul-162.1%
Simplified62.1%
Taylor expanded in b around inf 86.2%
Taylor expanded in b around 0 88.4%
+-commutative88.4%
mul-1-neg88.4%
unsub-neg88.4%
Simplified88.4%
Final simplification72.6%
(FPCore (a b c)
:precision binary64
(if (<= b 8.5e-39)
(if (>= b 0.0)
(/ (* 0.5 (- b (sqrt (* a (* c -4.0))))) a)
(/ 2.0 (fma 2.0 (/ a b) (* -2.0 (/ b c)))))
(if (>= b 0.0) (- (/ c b) (/ b a)) (/ (- c) b))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= 8.5e-39) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = (0.5 * (b - sqrt((a * (c * -4.0))))) / a;
} else {
tmp_2 = 2.0 / fma(2.0, (a / b), (-2.0 * (b / c)));
}
tmp_1 = tmp_2;
} else if (b >= 0.0) {
tmp_1 = (c / b) - (b / a);
} else {
tmp_1 = -c / b;
}
return tmp_1;
}
function code(a, b, c) tmp_1 = 0.0 if (b <= 8.5e-39) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(0.5 * Float64(b - sqrt(Float64(a * Float64(c * -4.0))))) / a); else tmp_2 = Float64(2.0 / fma(2.0, Float64(a / b), Float64(-2.0 * Float64(b / c)))); end tmp_1 = tmp_2; elseif (b >= 0.0) tmp_1 = Float64(Float64(c / b) - Float64(b / a)); else tmp_1 = Float64(Float64(-c) / b); end return tmp_1 end
code[a_, b_, c_] := If[LessEqual[b, 8.5e-39], If[GreaterEqual[b, 0.0], N[(N[(0.5 * N[(b - N[Sqrt[N[(a * N[(c * -4.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], N[(2.0 / N[(2.0 * N[(a / b), $MachinePrecision] + N[(-2.0 * N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision], N[((-c) / b), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 8.5 \cdot 10^{-39}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{0.5 \cdot \left(b - \sqrt{a \cdot \left(c \cdot -4\right)}\right)}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(2, \frac{a}{b}, -2 \cdot \frac{b}{c}\right)}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{-c}{b}\\
\end{array}
\end{array}
if b < 8.5000000000000005e-39Initial program 72.3%
Simplified72.2%
Taylor expanded in b around -inf 66.1%
+-commutative66.1%
fma-def66.1%
Simplified66.1%
Taylor expanded in b around 0 64.0%
associate-*r*64.0%
metadata-eval64.0%
distribute-lft-neg-in64.0%
*-commutative64.0%
distribute-lft-neg-in64.0%
metadata-eval64.0%
Simplified64.0%
div-sub64.0%
sub-neg64.0%
*-un-lft-identity64.0%
*-commutative64.0%
times-frac64.0%
metadata-eval64.0%
add-sqr-sqrt47.1%
sqrt-unprod63.2%
sqr-neg63.2%
sqrt-unprod63.2%
add-sqr-sqrt63.2%
*-un-lft-identity63.2%
*-commutative63.2%
times-frac63.2%
metadata-eval63.2%
associate-*r*63.2%
*-commutative63.2%
Applied egg-rr63.2%
sub-neg63.2%
distribute-lft-out--63.2%
div-sub63.2%
associate-*r/63.2%
Simplified63.2%
if 8.5000000000000005e-39 < b Initial program 62.1%
Taylor expanded in b around -inf 62.1%
associate-*r/62.1%
neg-mul-162.1%
Simplified62.1%
Taylor expanded in b around inf 86.2%
Taylor expanded in b around 0 88.4%
+-commutative88.4%
mul-1-neg88.4%
unsub-neg88.4%
Simplified88.4%
Final simplification72.6%
(FPCore (a b c) :precision binary64 (if (>= b 0.0) (- (/ c b) (/ b a)) (/ (- c) b)))
double code(double a, double b, double c) {
double tmp;
if (b >= 0.0) {
tmp = (c / b) - (b / a);
} else {
tmp = -c / b;
}
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.0d0) then
tmp = (c / b) - (b / a)
else
tmp = -c / b
end if
code = tmp
end function
public static double code(double a, double b, double c) {
double tmp;
if (b >= 0.0) {
tmp = (c / b) - (b / a);
} else {
tmp = -c / b;
}
return tmp;
}
def code(a, b, c): tmp = 0 if b >= 0.0: tmp = (c / b) - (b / a) else: tmp = -c / b return tmp
function code(a, b, c) tmp = 0.0 if (b >= 0.0) tmp = Float64(Float64(c / b) - Float64(b / a)); else tmp = Float64(Float64(-c) / b); end return tmp end
function tmp_2 = code(a, b, c) tmp = 0.0; if (b >= 0.0) tmp = (c / b) - (b / a); else tmp = -c / b; end tmp_2 = tmp; end
code[a_, b_, c_] := If[GreaterEqual[b, 0.0], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision], N[((-c) / b), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{-c}{b}\\
\end{array}
\end{array}
Initial program 68.5%
Taylor expanded in b around -inf 64.8%
associate-*r/64.8%
neg-mul-164.8%
Simplified64.8%
Taylor expanded in b around inf 64.9%
Taylor expanded in b around 0 65.8%
+-commutative65.8%
mul-1-neg65.8%
unsub-neg65.8%
Simplified65.8%
Final simplification65.8%
(FPCore (a b c) :precision binary64 (if (>= b 0.0) (/ c b) (/ (- c) b)))
double code(double a, double b, double c) {
double tmp;
if (b >= 0.0) {
tmp = c / b;
} else {
tmp = -c / b;
}
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.0d0) then
tmp = c / b
else
tmp = -c / b
end if
code = tmp
end function
public static double code(double a, double b, double c) {
double tmp;
if (b >= 0.0) {
tmp = c / b;
} else {
tmp = -c / b;
}
return tmp;
}
def code(a, b, c): tmp = 0 if b >= 0.0: tmp = c / b else: tmp = -c / b return tmp
function code(a, b, c) tmp = 0.0 if (b >= 0.0) tmp = Float64(c / b); else tmp = Float64(Float64(-c) / b); end return tmp end
function tmp_2 = code(a, b, c) tmp = 0.0; if (b >= 0.0) tmp = c / b; else tmp = -c / b; end tmp_2 = tmp; end
code[a_, b_, c_] := If[GreaterEqual[b, 0.0], N[(c / b), $MachinePrecision], N[((-c) / b), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{c}{b}\\
\mathbf{else}:\\
\;\;\;\;\frac{-c}{b}\\
\end{array}
\end{array}
Initial program 68.5%
Taylor expanded in b around -inf 64.8%
associate-*r/64.8%
neg-mul-164.8%
Simplified64.8%
Taylor expanded in b around inf 64.9%
Taylor expanded in b around 0 31.4%
Final simplification31.4%
(FPCore (a b c) :precision binary64 (if (>= b 0.0) (- (/ b a)) (/ (- c) b)))
double code(double a, double b, double c) {
double tmp;
if (b >= 0.0) {
tmp = -(b / a);
} else {
tmp = -c / b;
}
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.0d0) then
tmp = -(b / a)
else
tmp = -c / b
end if
code = tmp
end function
public static double code(double a, double b, double c) {
double tmp;
if (b >= 0.0) {
tmp = -(b / a);
} else {
tmp = -c / b;
}
return tmp;
}
def code(a, b, c): tmp = 0 if b >= 0.0: tmp = -(b / a) else: tmp = -c / b return tmp
function code(a, b, c) tmp = 0.0 if (b >= 0.0) tmp = Float64(-Float64(b / a)); else tmp = Float64(Float64(-c) / b); end return tmp end
function tmp_2 = code(a, b, c) tmp = 0.0; if (b >= 0.0) tmp = -(b / a); else tmp = -c / b; end tmp_2 = tmp; end
code[a_, b_, c_] := If[GreaterEqual[b, 0.0], (-N[(b / a), $MachinePrecision]), N[((-c) / b), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;-\frac{b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{-c}{b}\\
\end{array}
\end{array}
Initial program 68.5%
Taylor expanded in b around -inf 64.8%
associate-*r/64.8%
neg-mul-164.8%
Simplified64.8%
Taylor expanded in b around inf 64.9%
Taylor expanded in b around inf 65.5%
associate-*r/65.5%
mul-1-neg65.5%
Simplified65.5%
Final simplification65.5%
herbie shell --seed 2023322
(FPCore (a b c)
:name "jeff quadratic root 1"
:precision binary64
(if (>= b 0.0) (/ (- (- b) (sqrt (- (* b b) (* (* 4.0 a) c)))) (* 2.0 a)) (/ (* 2.0 c) (+ (- b) (sqrt (- (* b b) (* (* 4.0 a) c)))))))