
(FPCore (a b c) :precision binary64 (let* ((t_0 (sqrt (- (* b b) (* (* 4.0 a) c))))) (if (>= b 0.0) (/ (* 2.0 c) (- (- b) t_0)) (/ (+ (- b) t_0) (* 2.0 a)))))
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 = (2.0 * c) / (-b - t_0);
} else {
tmp = (-b + t_0) / (2.0 * a);
}
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 = (2.0d0 * c) / (-b - t_0)
else
tmp = (-b + t_0) / (2.0d0 * a)
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 = (2.0 * c) / (-b - t_0);
} else {
tmp = (-b + t_0) / (2.0 * a);
}
return tmp;
}
def code(a, b, c): t_0 = math.sqrt(((b * b) - ((4.0 * a) * c))) tmp = 0 if b >= 0.0: tmp = (2.0 * c) / (-b - t_0) else: tmp = (-b + t_0) / (2.0 * a) 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(2.0 * c) / Float64(Float64(-b) - t_0)); else tmp = Float64(Float64(Float64(-b) + t_0) / Float64(2.0 * a)); 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 = (2.0 * c) / (-b - t_0); else tmp = (-b + t_0) / (2.0 * a); 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[(2.0 * c), $MachinePrecision] / N[((-b) - t$95$0), $MachinePrecision]), $MachinePrecision], N[(N[((-b) + t$95$0), $MachinePrecision] / N[(2.0 * a), $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{2 \cdot c}{\left(-b\right) - t\_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + t\_0}{2 \cdot a}\\
\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) (/ (* 2.0 c) (- (- b) t_0)) (/ (+ (- b) t_0) (* 2.0 a)))))
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 = (2.0 * c) / (-b - t_0);
} else {
tmp = (-b + t_0) / (2.0 * a);
}
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 = (2.0d0 * c) / (-b - t_0)
else
tmp = (-b + t_0) / (2.0d0 * a)
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 = (2.0 * c) / (-b - t_0);
} else {
tmp = (-b + t_0) / (2.0 * a);
}
return tmp;
}
def code(a, b, c): t_0 = math.sqrt(((b * b) - ((4.0 * a) * c))) tmp = 0 if b >= 0.0: tmp = (2.0 * c) / (-b - t_0) else: tmp = (-b + t_0) / (2.0 * a) 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(2.0 * c) / Float64(Float64(-b) - t_0)); else tmp = Float64(Float64(Float64(-b) + t_0) / Float64(2.0 * a)); 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 = (2.0 * c) / (-b - t_0); else tmp = (-b + t_0) / (2.0 * a); 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[(2.0 * c), $MachinePrecision] / N[((-b) - t$95$0), $MachinePrecision]), $MachinePrecision], N[(N[((-b) + t$95$0), $MachinePrecision] / N[(2.0 * a), $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{2 \cdot c}{\left(-b\right) - t\_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + t\_0}{2 \cdot a}\\
\end{array}
\end{array}
(FPCore (a b c)
:precision binary64
(if (<= b -1e+142)
(if (>= b 0.0)
(/ (* (- 2.0) c) (+ b (sqrt (- (* b b) (* (* 4.0 a) c)))))
(/
(/
(* (- (* 4.0 (pow (/ (* (/ c b) a) b) 2.0)) 4.0) (- b))
(- (* (* (/ -2.0 b) a) (/ c b)) 2.0))
(* 2.0 a)))
(if (<= b 6e+34)
(if (>= b 0.0)
(/ (* -2.0 c) (+ (sqrt (fma (* c a) -4.0 (* b b))) b))
(* (- (/ (sqrt (fma -4.0 (* c a) (* b b))) a) (/ b a)) 0.5))
(/ (* -2.0 c) (* -2.0 (- (* a (/ c b)) b))))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= -1e+142) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = (-2.0 * c) / (b + sqrt(((b * b) - ((4.0 * a) * c))));
} else {
tmp_2 = ((((4.0 * pow((((c / b) * a) / b), 2.0)) - 4.0) * -b) / ((((-2.0 / b) * a) * (c / b)) - 2.0)) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 6e+34) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-2.0 * c) / (sqrt(fma((c * a), -4.0, (b * b))) + b);
} else {
tmp_3 = ((sqrt(fma(-4.0, (c * a), (b * b))) / a) - (b / a)) * 0.5;
}
tmp_1 = tmp_3;
} else {
tmp_1 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b));
}
return tmp_1;
}
function code(a, b, c) tmp_1 = 0.0 if (b <= -1e+142) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(Float64(-2.0) * c) / Float64(b + sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c))))); else tmp_2 = Float64(Float64(Float64(Float64(Float64(4.0 * (Float64(Float64(Float64(c / b) * a) / b) ^ 2.0)) - 4.0) * Float64(-b)) / Float64(Float64(Float64(Float64(-2.0 / b) * a) * Float64(c / b)) - 2.0)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= 6e+34) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-2.0 * c) / Float64(sqrt(fma(Float64(c * a), -4.0, Float64(b * b))) + b)); else tmp_3 = Float64(Float64(Float64(sqrt(fma(-4.0, Float64(c * a), Float64(b * b))) / a) - Float64(b / a)) * 0.5); end tmp_1 = tmp_3; else tmp_1 = Float64(Float64(-2.0 * c) / Float64(-2.0 * Float64(Float64(a * Float64(c / b)) - b))); end return tmp_1 end
code[a_, b_, c_] := If[LessEqual[b, -1e+142], If[GreaterEqual[b, 0.0], N[(N[((-2.0) * c), $MachinePrecision] / N[(b + N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(4.0 * N[Power[N[(N[(N[(c / b), $MachinePrecision] * a), $MachinePrecision] / b), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] - 4.0), $MachinePrecision] * (-b)), $MachinePrecision] / N[(N[(N[(N[(-2.0 / b), $MachinePrecision] * a), $MachinePrecision] * N[(c / b), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 6e+34], If[GreaterEqual[b, 0.0], N[(N[(-2.0 * c), $MachinePrecision] / N[(N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0 + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[Sqrt[N[(-4.0 * N[(c * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / a), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision] * 0.5), $MachinePrecision]], N[(N[(-2.0 * c), $MachinePrecision] / N[(-2.0 * N[(N[(a * N[(c / b), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1 \cdot 10^{+142}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-2\right) \cdot c}{b + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\left(4 \cdot {\left(\frac{\frac{c}{b} \cdot a}{b}\right)}^{2} - 4\right) \cdot \left(-b\right)}{\left(\frac{-2}{b} \cdot a\right) \cdot \frac{c}{b} - 2}}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 6 \cdot 10^{+34}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-2 \cdot c}{\sqrt{\mathsf{fma}\left(c \cdot a, -4, b \cdot b\right)} + b}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\sqrt{\mathsf{fma}\left(-4, c \cdot a, b \cdot b\right)}}{a} - \frac{b}{a}\right) \cdot 0.5\\
\end{array}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{-2 \cdot \left(a \cdot \frac{c}{b} - b\right)}\\
\end{array}
\end{array}
if b < -1.00000000000000005e142Initial program 28.1%
Taylor expanded in b around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r/N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6496.0
Applied rewrites96.0%
Applied rewrites96.0%
if -1.00000000000000005e142 < b < 6.00000000000000037e34Initial program 85.3%
Taylor expanded in a around 0
Applied rewrites85.3%
Applied rewrites85.3%
if 6.00000000000000037e34 < b Initial program 55.7%
Applied rewrites55.7%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites55.7%
Taylor expanded in a around 0
Applied rewrites93.3%
Taylor expanded in a around 0
Applied rewrites93.7%
Final simplification89.9%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (* a (/ c b))))
(if (<= b -6.6e+141)
(if (>= b 0.0)
(/ (* (- 2.0) c) (+ b (sqrt (- (* b b) (* (* 4.0 a) c)))))
(/ (* (fma (/ -2.0 b) t_0 2.0) (- b)) (* 2.0 a)))
(if (<= b 6e+34)
(if (>= b 0.0)
(/ (* -2.0 c) (+ (sqrt (fma (* c a) -4.0 (* b b))) b))
(* (- (/ (sqrt (fma -4.0 (* c a) (* b b))) a) (/ b a)) 0.5))
(/ (* -2.0 c) (* -2.0 (- t_0 b)))))))
double code(double a, double b, double c) {
double t_0 = a * (c / b);
double tmp_1;
if (b <= -6.6e+141) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = (-2.0 * c) / (b + sqrt(((b * b) - ((4.0 * a) * c))));
} else {
tmp_2 = (fma((-2.0 / b), t_0, 2.0) * -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 6e+34) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-2.0 * c) / (sqrt(fma((c * a), -4.0, (b * b))) + b);
} else {
tmp_3 = ((sqrt(fma(-4.0, (c * a), (b * b))) / a) - (b / a)) * 0.5;
}
tmp_1 = tmp_3;
} else {
tmp_1 = (-2.0 * c) / (-2.0 * (t_0 - b));
}
return tmp_1;
}
function code(a, b, c) t_0 = Float64(a * Float64(c / b)) tmp_1 = 0.0 if (b <= -6.6e+141) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(Float64(-2.0) * c) / Float64(b + sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c))))); else tmp_2 = Float64(Float64(fma(Float64(-2.0 / b), t_0, 2.0) * Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= 6e+34) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-2.0 * c) / Float64(sqrt(fma(Float64(c * a), -4.0, Float64(b * b))) + b)); else tmp_3 = Float64(Float64(Float64(sqrt(fma(-4.0, Float64(c * a), Float64(b * b))) / a) - Float64(b / a)) * 0.5); end tmp_1 = tmp_3; else tmp_1 = Float64(Float64(-2.0 * c) / Float64(-2.0 * Float64(t_0 - b))); end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[(a * N[(c / b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -6.6e+141], If[GreaterEqual[b, 0.0], N[(N[((-2.0) * c), $MachinePrecision] / N[(b + N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(-2.0 / b), $MachinePrecision] * t$95$0 + 2.0), $MachinePrecision] * (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 6e+34], If[GreaterEqual[b, 0.0], N[(N[(-2.0 * c), $MachinePrecision] / N[(N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0 + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[Sqrt[N[(-4.0 * N[(c * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / a), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision] * 0.5), $MachinePrecision]], N[(N[(-2.0 * c), $MachinePrecision] / N[(-2.0 * N[(t$95$0 - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot \frac{c}{b}\\
\mathbf{if}\;b \leq -6.6 \cdot 10^{+141}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-2\right) \cdot c}{b + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{-2}{b}, t\_0, 2\right) \cdot \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 6 \cdot 10^{+34}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-2 \cdot c}{\sqrt{\mathsf{fma}\left(c \cdot a, -4, b \cdot b\right)} + b}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\sqrt{\mathsf{fma}\left(-4, c \cdot a, b \cdot b\right)}}{a} - \frac{b}{a}\right) \cdot 0.5\\
\end{array}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{-2 \cdot \left(t\_0 - b\right)}\\
\end{array}
\end{array}
if b < -6.5999999999999993e141Initial program 28.1%
Taylor expanded in b around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r/N/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6496.0
Applied rewrites96.0%
if -6.5999999999999993e141 < b < 6.00000000000000037e34Initial program 85.3%
Taylor expanded in a around 0
Applied rewrites85.3%
Applied rewrites85.3%
if 6.00000000000000037e34 < b Initial program 55.7%
Applied rewrites55.7%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites55.7%
Taylor expanded in a around 0
Applied rewrites93.3%
Taylor expanded in a around 0
Applied rewrites93.7%
Final simplification89.9%
(FPCore (a b c)
:precision binary64
(if (<= b -6.2e+82)
(if (>= b 0.0) (/ (- b) a) (/ (+ (- b) (- b)) (* 2.0 a)))
(if (<= b 6e+34)
(if (>= b 0.0)
(/ (* -2.0 c) (+ (sqrt (fma (* c a) -4.0 (* b b))) b))
(* (- (/ (sqrt (fma -4.0 (* c a) (* b b))) a) (/ b a)) 0.5))
(/ (* -2.0 c) (* -2.0 (- (* a (/ c b)) b))))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= -6.2e+82) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 6e+34) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-2.0 * c) / (sqrt(fma((c * a), -4.0, (b * b))) + b);
} else {
tmp_3 = ((sqrt(fma(-4.0, (c * a), (b * b))) / a) - (b / a)) * 0.5;
}
tmp_1 = tmp_3;
} else {
tmp_1 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b));
}
return tmp_1;
}
function code(a, b, c) tmp_1 = 0.0 if (b <= -6.2e+82) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(-b) / a); else tmp_2 = Float64(Float64(Float64(-b) + Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= 6e+34) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-2.0 * c) / Float64(sqrt(fma(Float64(c * a), -4.0, Float64(b * b))) + b)); else tmp_3 = Float64(Float64(Float64(sqrt(fma(-4.0, Float64(c * a), Float64(b * b))) / a) - Float64(b / a)) * 0.5); end tmp_1 = tmp_3; else tmp_1 = Float64(Float64(-2.0 * c) / Float64(-2.0 * Float64(Float64(a * Float64(c / b)) - b))); end return tmp_1 end
code[a_, b_, c_] := If[LessEqual[b, -6.2e+82], If[GreaterEqual[b, 0.0], N[((-b) / a), $MachinePrecision], N[(N[((-b) + (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 6e+34], If[GreaterEqual[b, 0.0], N[(N[(-2.0 * c), $MachinePrecision] / N[(N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0 + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[Sqrt[N[(-4.0 * N[(c * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / a), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision] * 0.5), $MachinePrecision]], N[(N[(-2.0 * c), $MachinePrecision] / N[(-2.0 * N[(N[(a * N[(c / b), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -6.2 \cdot 10^{+82}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 6 \cdot 10^{+34}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-2 \cdot c}{\sqrt{\mathsf{fma}\left(c \cdot a, -4, b \cdot b\right)} + b}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\sqrt{\mathsf{fma}\left(-4, c \cdot a, b \cdot b\right)}}{a} - \frac{b}{a}\right) \cdot 0.5\\
\end{array}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{-2 \cdot \left(a \cdot \frac{c}{b} - b\right)}\\
\end{array}
\end{array}
if b < -6.20000000000000065e82Initial program 44.4%
Taylor expanded in a around 0
distribute-lft-out--N/A
lower-*.f64N/A
lower--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6444.4
Applied rewrites44.4%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6495.2
Applied rewrites95.2%
Taylor expanded in b around -inf
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6495.2
Applied rewrites95.2%
if -6.20000000000000065e82 < b < 6.00000000000000037e34Initial program 83.9%
Taylor expanded in a around 0
Applied rewrites83.9%
Applied rewrites83.9%
if 6.00000000000000037e34 < b Initial program 55.7%
Applied rewrites55.7%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites55.7%
Taylor expanded in a around 0
Applied rewrites93.3%
Taylor expanded in a around 0
Applied rewrites93.7%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (sqrt (fma (* c a) -4.0 (* b b)))))
(if (<= b -6.2e+82)
(if (>= b 0.0) (/ (- b) a) (/ (+ (- b) (- b)) (* 2.0 a)))
(if (<= b 6e+34)
(if (>= b 0.0) (/ (* -2.0 c) (+ t_0 b)) (* (/ (- t_0 b) a) 0.5))
(/ (* -2.0 c) (* -2.0 (- (* a (/ c b)) b)))))))
double code(double a, double b, double c) {
double t_0 = sqrt(fma((c * a), -4.0, (b * b)));
double tmp_1;
if (b <= -6.2e+82) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 6e+34) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-2.0 * c) / (t_0 + b);
} else {
tmp_3 = ((t_0 - b) / a) * 0.5;
}
tmp_1 = tmp_3;
} else {
tmp_1 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b));
}
return tmp_1;
}
function code(a, b, c) t_0 = sqrt(fma(Float64(c * a), -4.0, Float64(b * b))) tmp_1 = 0.0 if (b <= -6.2e+82) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(-b) / a); else tmp_2 = Float64(Float64(Float64(-b) + Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= 6e+34) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-2.0 * c) / Float64(t_0 + b)); else tmp_3 = Float64(Float64(Float64(t_0 - b) / a) * 0.5); end tmp_1 = tmp_3; else tmp_1 = Float64(Float64(-2.0 * c) / Float64(-2.0 * Float64(Float64(a * Float64(c / b)) - b))); end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0 + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[b, -6.2e+82], If[GreaterEqual[b, 0.0], N[((-b) / a), $MachinePrecision], N[(N[((-b) + (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 6e+34], If[GreaterEqual[b, 0.0], N[(N[(-2.0 * c), $MachinePrecision] / N[(t$95$0 + b), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t$95$0 - b), $MachinePrecision] / a), $MachinePrecision] * 0.5), $MachinePrecision]], N[(N[(-2.0 * c), $MachinePrecision] / N[(-2.0 * N[(N[(a * N[(c / b), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{\mathsf{fma}\left(c \cdot a, -4, b \cdot b\right)}\\
\mathbf{if}\;b \leq -6.2 \cdot 10^{+82}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 6 \cdot 10^{+34}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-2 \cdot c}{t\_0 + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_0 - b}{a} \cdot 0.5\\
\end{array}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{-2 \cdot \left(a \cdot \frac{c}{b} - b\right)}\\
\end{array}
\end{array}
if b < -6.20000000000000065e82Initial program 44.4%
Taylor expanded in a around 0
distribute-lft-out--N/A
lower-*.f64N/A
lower--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6444.4
Applied rewrites44.4%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6495.2
Applied rewrites95.2%
Taylor expanded in b around -inf
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6495.2
Applied rewrites95.2%
if -6.20000000000000065e82 < b < 6.00000000000000037e34Initial program 83.9%
Taylor expanded in a around 0
Applied rewrites83.9%
if 6.00000000000000037e34 < b Initial program 55.7%
Applied rewrites55.7%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites55.7%
Taylor expanded in a around 0
Applied rewrites93.3%
Taylor expanded in a around 0
Applied rewrites93.7%
(FPCore (a b c)
:precision binary64
(if (<= b -6.2e+82)
(if (>= b 0.0) (/ (- b) a) (/ (+ (- b) (- b)) (* 2.0 a)))
(if (<= b 6e+34)
(if (>= b 0.0)
(* c (/ -2.0 (+ (sqrt (fma -4.0 (* c a) (* b b))) b)))
(* (/ (- (sqrt (fma (* c a) -4.0 (* b b))) b) a) 0.5))
(/ (* -2.0 c) (* -2.0 (- (* a (/ c b)) b))))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= -6.2e+82) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 6e+34) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = c * (-2.0 / (sqrt(fma(-4.0, (c * a), (b * b))) + b));
} else {
tmp_3 = ((sqrt(fma((c * a), -4.0, (b * b))) - b) / a) * 0.5;
}
tmp_1 = tmp_3;
} else {
tmp_1 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b));
}
return tmp_1;
}
function code(a, b, c) tmp_1 = 0.0 if (b <= -6.2e+82) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(-b) / a); else tmp_2 = Float64(Float64(Float64(-b) + Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= 6e+34) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(c * Float64(-2.0 / Float64(sqrt(fma(-4.0, Float64(c * a), Float64(b * b))) + b))); else tmp_3 = Float64(Float64(Float64(sqrt(fma(Float64(c * a), -4.0, Float64(b * b))) - b) / a) * 0.5); end tmp_1 = tmp_3; else tmp_1 = Float64(Float64(-2.0 * c) / Float64(-2.0 * Float64(Float64(a * Float64(c / b)) - b))); end return tmp_1 end
code[a_, b_, c_] := If[LessEqual[b, -6.2e+82], If[GreaterEqual[b, 0.0], N[((-b) / a), $MachinePrecision], N[(N[((-b) + (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 6e+34], If[GreaterEqual[b, 0.0], N[(c * N[(-2.0 / N[(N[Sqrt[N[(-4.0 * N[(c * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0 + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] / a), $MachinePrecision] * 0.5), $MachinePrecision]], N[(N[(-2.0 * c), $MachinePrecision] / N[(-2.0 * N[(N[(a * N[(c / b), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -6.2 \cdot 10^{+82}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 6 \cdot 10^{+34}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;c \cdot \frac{-2}{\sqrt{\mathsf{fma}\left(-4, c \cdot a, b \cdot b\right)} + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{\mathsf{fma}\left(c \cdot a, -4, b \cdot b\right)} - b}{a} \cdot 0.5\\
\end{array}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{-2 \cdot \left(a \cdot \frac{c}{b} - b\right)}\\
\end{array}
\end{array}
if b < -6.20000000000000065e82Initial program 44.4%
Taylor expanded in a around 0
distribute-lft-out--N/A
lower-*.f64N/A
lower--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6444.4
Applied rewrites44.4%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6495.2
Applied rewrites95.2%
Taylor expanded in b around -inf
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6495.2
Applied rewrites95.2%
if -6.20000000000000065e82 < b < 6.00000000000000037e34Initial program 83.9%
Taylor expanded in a around 0
Applied rewrites83.9%
Applied rewrites83.8%
if 6.00000000000000037e34 < b Initial program 55.7%
Applied rewrites55.7%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites55.7%
Taylor expanded in a around 0
Applied rewrites93.3%
Taylor expanded in a around 0
Applied rewrites93.7%
(FPCore (a b c)
:precision binary64
(if (<= b -4.1e+82)
(if (>= b 0.0) (/ (- b) a) (/ (+ (- b) (- b)) (* 2.0 a)))
(if (<= b -2.9e-265)
(* (/ 0.5 a) (- (sqrt (fma (* -4.0 c) a (* b b))) b))
(if (<= b 6e+34)
(/ (* -2.0 c) (+ (sqrt (fma (* c a) -4.0 (* b b))) b))
(/ (* -2.0 c) (* -2.0 (- (* a (/ c b)) b)))))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= -4.1e+82) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= -2.9e-265) {
tmp_1 = (0.5 / a) * (sqrt(fma((-4.0 * c), a, (b * b))) - b);
} else if (b <= 6e+34) {
tmp_1 = (-2.0 * c) / (sqrt(fma((c * a), -4.0, (b * b))) + b);
} else {
tmp_1 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b));
}
return tmp_1;
}
function code(a, b, c) tmp_1 = 0.0 if (b <= -4.1e+82) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(-b) / a); else tmp_2 = Float64(Float64(Float64(-b) + Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= -2.9e-265) tmp_1 = Float64(Float64(0.5 / a) * Float64(sqrt(fma(Float64(-4.0 * c), a, Float64(b * b))) - b)); elseif (b <= 6e+34) tmp_1 = Float64(Float64(-2.0 * c) / Float64(sqrt(fma(Float64(c * a), -4.0, Float64(b * b))) + b)); else tmp_1 = Float64(Float64(-2.0 * c) / Float64(-2.0 * Float64(Float64(a * Float64(c / b)) - b))); end return tmp_1 end
code[a_, b_, c_] := If[LessEqual[b, -4.1e+82], If[GreaterEqual[b, 0.0], N[((-b) / a), $MachinePrecision], N[(N[((-b) + (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, -2.9e-265], N[(N[(0.5 / a), $MachinePrecision] * N[(N[Sqrt[N[(N[(-4.0 * c), $MachinePrecision] * a + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 6e+34], N[(N[(-2.0 * c), $MachinePrecision] / N[(N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0 + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision], N[(N[(-2.0 * c), $MachinePrecision] / N[(-2.0 * N[(N[(a * N[(c / b), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -4.1 \cdot 10^{+82}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq -2.9 \cdot 10^{-265}:\\
\;\;\;\;\frac{0.5}{a} \cdot \left(\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} - b\right)\\
\mathbf{elif}\;b \leq 6 \cdot 10^{+34}:\\
\;\;\;\;\frac{-2 \cdot c}{\sqrt{\mathsf{fma}\left(c \cdot a, -4, b \cdot b\right)} + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{-2 \cdot \left(a \cdot \frac{c}{b} - b\right)}\\
\end{array}
\end{array}
if b < -4.09999999999999995e82Initial program 44.4%
Taylor expanded in a around 0
distribute-lft-out--N/A
lower-*.f64N/A
lower--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6444.4
Applied rewrites44.4%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6495.2
Applied rewrites95.2%
Taylor expanded in b around -inf
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6495.2
Applied rewrites95.2%
if -4.09999999999999995e82 < b < -2.89999999999999975e-265Initial program 83.9%
Applied rewrites29.6%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites43.1%
Applied rewrites42.6%
Taylor expanded in a around 0
Applied rewrites83.7%
if -2.89999999999999975e-265 < b < 6.00000000000000037e34Initial program 84.0%
Applied rewrites83.8%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites84.1%
if 6.00000000000000037e34 < b Initial program 55.7%
Applied rewrites55.7%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites55.7%
Taylor expanded in a around 0
Applied rewrites93.3%
Taylor expanded in a around 0
Applied rewrites93.7%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (sqrt (fma (* -4.0 c) a (* b b)))))
(if (<= b -4.1e+82)
(if (>= b 0.0) (/ (- b) a) (/ (+ (- b) (- b)) (* 2.0 a)))
(if (<= b -2.35e-266)
(* (/ 0.5 a) (- t_0 b))
(if (<= b 6e+34)
(* c (/ -2.0 (+ t_0 b)))
(/ (* -2.0 c) (* -2.0 (- (* a (/ c b)) b))))))))
double code(double a, double b, double c) {
double t_0 = sqrt(fma((-4.0 * c), a, (b * b)));
double tmp_1;
if (b <= -4.1e+82) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= -2.35e-266) {
tmp_1 = (0.5 / a) * (t_0 - b);
} else if (b <= 6e+34) {
tmp_1 = c * (-2.0 / (t_0 + b));
} else {
tmp_1 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b));
}
return tmp_1;
}
function code(a, b, c) t_0 = sqrt(fma(Float64(-4.0 * c), a, Float64(b * b))) tmp_1 = 0.0 if (b <= -4.1e+82) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(-b) / a); else tmp_2 = Float64(Float64(Float64(-b) + Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= -2.35e-266) tmp_1 = Float64(Float64(0.5 / a) * Float64(t_0 - b)); elseif (b <= 6e+34) tmp_1 = Float64(c * Float64(-2.0 / Float64(t_0 + b))); else tmp_1 = Float64(Float64(-2.0 * c) / Float64(-2.0 * Float64(Float64(a * Float64(c / b)) - b))); end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[Sqrt[N[(N[(-4.0 * c), $MachinePrecision] * a + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[b, -4.1e+82], If[GreaterEqual[b, 0.0], N[((-b) / a), $MachinePrecision], N[(N[((-b) + (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, -2.35e-266], N[(N[(0.5 / a), $MachinePrecision] * N[(t$95$0 - b), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 6e+34], N[(c * N[(-2.0 / N[(t$95$0 + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(-2.0 * c), $MachinePrecision] / N[(-2.0 * N[(N[(a * N[(c / b), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)}\\
\mathbf{if}\;b \leq -4.1 \cdot 10^{+82}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq -2.35 \cdot 10^{-266}:\\
\;\;\;\;\frac{0.5}{a} \cdot \left(t\_0 - b\right)\\
\mathbf{elif}\;b \leq 6 \cdot 10^{+34}:\\
\;\;\;\;c \cdot \frac{-2}{t\_0 + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{-2 \cdot \left(a \cdot \frac{c}{b} - b\right)}\\
\end{array}
\end{array}
if b < -4.09999999999999995e82Initial program 44.4%
Taylor expanded in a around 0
distribute-lft-out--N/A
lower-*.f64N/A
lower--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6444.4
Applied rewrites44.4%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6495.2
Applied rewrites95.2%
Taylor expanded in b around -inf
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6495.2
Applied rewrites95.2%
if -4.09999999999999995e82 < b < -2.35000000000000014e-266Initial program 83.9%
Applied rewrites29.6%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites43.1%
Applied rewrites42.6%
Taylor expanded in a around 0
Applied rewrites83.7%
if -2.35000000000000014e-266 < b < 6.00000000000000037e34Initial program 84.0%
Applied rewrites83.8%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites84.1%
Applied rewrites83.8%
if 6.00000000000000037e34 < b Initial program 55.7%
Applied rewrites55.7%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites55.7%
Taylor expanded in a around 0
Applied rewrites93.3%
Taylor expanded in a around 0
Applied rewrites93.7%
(FPCore (a b c)
:precision binary64
(if (<= b -1.12e-11)
(if (>= b 0.0) (/ (- b) a) (/ (+ (- b) (- b)) (* 2.0 a)))
(if (<= b 6e+34)
(* c (/ -2.0 (+ (sqrt (fma (* -4.0 c) a (* b b))) b)))
(/ (* -2.0 c) (* -2.0 (- (* a (/ c b)) b))))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= -1.12e-11) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 6e+34) {
tmp_1 = c * (-2.0 / (sqrt(fma((-4.0 * c), a, (b * b))) + b));
} else {
tmp_1 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b));
}
return tmp_1;
}
function code(a, b, c) tmp_1 = 0.0 if (b <= -1.12e-11) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(-b) / a); else tmp_2 = Float64(Float64(Float64(-b) + Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= 6e+34) tmp_1 = Float64(c * Float64(-2.0 / Float64(sqrt(fma(Float64(-4.0 * c), a, Float64(b * b))) + b))); else tmp_1 = Float64(Float64(-2.0 * c) / Float64(-2.0 * Float64(Float64(a * Float64(c / b)) - b))); end return tmp_1 end
code[a_, b_, c_] := If[LessEqual[b, -1.12e-11], If[GreaterEqual[b, 0.0], N[((-b) / a), $MachinePrecision], N[(N[((-b) + (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 6e+34], N[(c * N[(-2.0 / N[(N[Sqrt[N[(N[(-4.0 * c), $MachinePrecision] * a + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(-2.0 * c), $MachinePrecision] / N[(-2.0 * N[(N[(a * N[(c / b), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.12 \cdot 10^{-11}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 6 \cdot 10^{+34}:\\
\;\;\;\;c \cdot \frac{-2}{\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{-2 \cdot \left(a \cdot \frac{c}{b} - b\right)}\\
\end{array}
\end{array}
if b < -1.1200000000000001e-11Initial program 55.2%
Taylor expanded in a around 0
distribute-lft-out--N/A
lower-*.f64N/A
lower--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6455.2
Applied rewrites55.2%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6491.3
Applied rewrites91.3%
Taylor expanded in b around -inf
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6491.3
Applied rewrites91.3%
if -1.1200000000000001e-11 < b < 6.00000000000000037e34Initial program 82.4%
Applied rewrites73.0%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites77.6%
Applied rewrites77.4%
if 6.00000000000000037e34 < b Initial program 55.7%
Applied rewrites55.7%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites55.7%
Taylor expanded in a around 0
Applied rewrites93.3%
Taylor expanded in a around 0
Applied rewrites93.7%
(FPCore (a b c)
:precision binary64
(if (<= b -1.12e-11)
(if (>= b 0.0) (/ (- b) a) (/ (+ (- b) (- b)) (* 2.0 a)))
(if (<= b 9.5e-64)
(/ (* -2.0 c) (+ (sqrt (* (* a c) -4.0)) b))
(/ (* -2.0 c) (* -2.0 (- (* (/ a b) c) b))))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= -1.12e-11) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 9.5e-64) {
tmp_1 = (-2.0 * c) / (sqrt(((a * c) * -4.0)) + b);
} else {
tmp_1 = (-2.0 * c) / (-2.0 * (((a / b) * c) - b));
}
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) :: tmp
real(8) :: tmp_1
real(8) :: tmp_2
if (b <= (-1.12d-11)) then
if (b >= 0.0d0) then
tmp_2 = -b / a
else
tmp_2 = (-b + -b) / (2.0d0 * a)
end if
tmp_1 = tmp_2
else if (b <= 9.5d-64) then
tmp_1 = ((-2.0d0) * c) / (sqrt(((a * c) * (-4.0d0))) + b)
else
tmp_1 = ((-2.0d0) * c) / ((-2.0d0) * (((a / b) * c) - b))
end if
code = tmp_1
end function
public static double code(double a, double b, double c) {
double tmp_1;
if (b <= -1.12e-11) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 9.5e-64) {
tmp_1 = (-2.0 * c) / (Math.sqrt(((a * c) * -4.0)) + b);
} else {
tmp_1 = (-2.0 * c) / (-2.0 * (((a / b) * c) - b));
}
return tmp_1;
}
def code(a, b, c): tmp_1 = 0 if b <= -1.12e-11: tmp_2 = 0 if b >= 0.0: tmp_2 = -b / a else: tmp_2 = (-b + -b) / (2.0 * a) tmp_1 = tmp_2 elif b <= 9.5e-64: tmp_1 = (-2.0 * c) / (math.sqrt(((a * c) * -4.0)) + b) else: tmp_1 = (-2.0 * c) / (-2.0 * (((a / b) * c) - b)) return tmp_1
function code(a, b, c) tmp_1 = 0.0 if (b <= -1.12e-11) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(-b) / a); else tmp_2 = Float64(Float64(Float64(-b) + Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= 9.5e-64) tmp_1 = Float64(Float64(-2.0 * c) / Float64(sqrt(Float64(Float64(a * c) * -4.0)) + b)); else tmp_1 = Float64(Float64(-2.0 * c) / Float64(-2.0 * Float64(Float64(Float64(a / b) * c) - b))); end return tmp_1 end
function tmp_4 = code(a, b, c) tmp_2 = 0.0; if (b <= -1.12e-11) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = -b / a; else tmp_3 = (-b + -b) / (2.0 * a); end tmp_2 = tmp_3; elseif (b <= 9.5e-64) tmp_2 = (-2.0 * c) / (sqrt(((a * c) * -4.0)) + b); else tmp_2 = (-2.0 * c) / (-2.0 * (((a / b) * c) - b)); end tmp_4 = tmp_2; end
code[a_, b_, c_] := If[LessEqual[b, -1.12e-11], If[GreaterEqual[b, 0.0], N[((-b) / a), $MachinePrecision], N[(N[((-b) + (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 9.5e-64], N[(N[(-2.0 * c), $MachinePrecision] / N[(N[Sqrt[N[(N[(a * c), $MachinePrecision] * -4.0), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision], N[(N[(-2.0 * c), $MachinePrecision] / N[(-2.0 * N[(N[(N[(a / b), $MachinePrecision] * c), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.12 \cdot 10^{-11}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 9.5 \cdot 10^{-64}:\\
\;\;\;\;\frac{-2 \cdot c}{\sqrt{\left(a \cdot c\right) \cdot -4} + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{-2 \cdot \left(\frac{a}{b} \cdot c - b\right)}\\
\end{array}
\end{array}
if b < -1.1200000000000001e-11Initial program 55.2%
Taylor expanded in a around 0
distribute-lft-out--N/A
lower-*.f64N/A
lower--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6455.2
Applied rewrites55.2%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6491.3
Applied rewrites91.3%
Taylor expanded in b around -inf
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6491.3
Applied rewrites91.3%
if -1.1200000000000001e-11 < b < 9.50000000000000043e-64Initial program 78.1%
Applied rewrites66.5%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites72.2%
Taylor expanded in a around inf
Applied rewrites62.5%
if 9.50000000000000043e-64 < b Initial program 64.3%
Applied rewrites64.3%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites64.3%
Taylor expanded in a around 0
Applied rewrites90.9%
Applied rewrites90.9%
(FPCore (a b c)
:precision binary64
(if (<= b -1.12e-11)
(if (>= b 0.0) (/ (- b) a) (/ (+ (- b) (- b)) (* 2.0 a)))
(if (<= b 9.5e-64)
(/ (* -2.0 c) (+ (sqrt (* (* a c) -4.0)) b))
(/ (* -2.0 c) (* -2.0 (- (* a (/ c b)) b))))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= -1.12e-11) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 9.5e-64) {
tmp_1 = (-2.0 * c) / (sqrt(((a * c) * -4.0)) + b);
} else {
tmp_1 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b));
}
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) :: tmp
real(8) :: tmp_1
real(8) :: tmp_2
if (b <= (-1.12d-11)) then
if (b >= 0.0d0) then
tmp_2 = -b / a
else
tmp_2 = (-b + -b) / (2.0d0 * a)
end if
tmp_1 = tmp_2
else if (b <= 9.5d-64) then
tmp_1 = ((-2.0d0) * c) / (sqrt(((a * c) * (-4.0d0))) + b)
else
tmp_1 = ((-2.0d0) * c) / ((-2.0d0) * ((a * (c / b)) - b))
end if
code = tmp_1
end function
public static double code(double a, double b, double c) {
double tmp_1;
if (b <= -1.12e-11) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 9.5e-64) {
tmp_1 = (-2.0 * c) / (Math.sqrt(((a * c) * -4.0)) + b);
} else {
tmp_1 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b));
}
return tmp_1;
}
def code(a, b, c): tmp_1 = 0 if b <= -1.12e-11: tmp_2 = 0 if b >= 0.0: tmp_2 = -b / a else: tmp_2 = (-b + -b) / (2.0 * a) tmp_1 = tmp_2 elif b <= 9.5e-64: tmp_1 = (-2.0 * c) / (math.sqrt(((a * c) * -4.0)) + b) else: tmp_1 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b)) return tmp_1
function code(a, b, c) tmp_1 = 0.0 if (b <= -1.12e-11) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(-b) / a); else tmp_2 = Float64(Float64(Float64(-b) + Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= 9.5e-64) tmp_1 = Float64(Float64(-2.0 * c) / Float64(sqrt(Float64(Float64(a * c) * -4.0)) + b)); else tmp_1 = Float64(Float64(-2.0 * c) / Float64(-2.0 * Float64(Float64(a * Float64(c / b)) - b))); end return tmp_1 end
function tmp_4 = code(a, b, c) tmp_2 = 0.0; if (b <= -1.12e-11) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = -b / a; else tmp_3 = (-b + -b) / (2.0 * a); end tmp_2 = tmp_3; elseif (b <= 9.5e-64) tmp_2 = (-2.0 * c) / (sqrt(((a * c) * -4.0)) + b); else tmp_2 = (-2.0 * c) / (-2.0 * ((a * (c / b)) - b)); end tmp_4 = tmp_2; end
code[a_, b_, c_] := If[LessEqual[b, -1.12e-11], If[GreaterEqual[b, 0.0], N[((-b) / a), $MachinePrecision], N[(N[((-b) + (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 9.5e-64], N[(N[(-2.0 * c), $MachinePrecision] / N[(N[Sqrt[N[(N[(a * c), $MachinePrecision] * -4.0), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision], N[(N[(-2.0 * c), $MachinePrecision] / N[(-2.0 * N[(N[(a * N[(c / b), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.12 \cdot 10^{-11}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 9.5 \cdot 10^{-64}:\\
\;\;\;\;\frac{-2 \cdot c}{\sqrt{\left(a \cdot c\right) \cdot -4} + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{-2 \cdot \left(a \cdot \frac{c}{b} - b\right)}\\
\end{array}
\end{array}
if b < -1.1200000000000001e-11Initial program 55.2%
Taylor expanded in a around 0
distribute-lft-out--N/A
lower-*.f64N/A
lower--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6455.2
Applied rewrites55.2%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6491.3
Applied rewrites91.3%
Taylor expanded in b around -inf
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6491.3
Applied rewrites91.3%
if -1.1200000000000001e-11 < b < 9.50000000000000043e-64Initial program 78.1%
Applied rewrites66.5%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites72.2%
Taylor expanded in a around inf
Applied rewrites62.5%
if 9.50000000000000043e-64 < b Initial program 64.3%
Applied rewrites64.3%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites64.3%
Taylor expanded in a around 0
Applied rewrites90.4%
Taylor expanded in a around 0
Applied rewrites90.9%
(FPCore (a b c)
:precision binary64
(if (<= b -1.12e-11)
(if (>= b 0.0) (/ (- b) a) (/ (+ (- b) (- b)) (* 2.0 a)))
(if (<= b 9.5e-64)
(/ (* -2.0 c) (+ (sqrt (* (* a c) -4.0)) b))
(/ (* -2.0 c) (* 2.0 b)))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= -1.12e-11) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 9.5e-64) {
tmp_1 = (-2.0 * c) / (sqrt(((a * c) * -4.0)) + b);
} else {
tmp_1 = (-2.0 * c) / (2.0 * b);
}
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) :: tmp
real(8) :: tmp_1
real(8) :: tmp_2
if (b <= (-1.12d-11)) then
if (b >= 0.0d0) then
tmp_2 = -b / a
else
tmp_2 = (-b + -b) / (2.0d0 * a)
end if
tmp_1 = tmp_2
else if (b <= 9.5d-64) then
tmp_1 = ((-2.0d0) * c) / (sqrt(((a * c) * (-4.0d0))) + b)
else
tmp_1 = ((-2.0d0) * c) / (2.0d0 * b)
end if
code = tmp_1
end function
public static double code(double a, double b, double c) {
double tmp_1;
if (b <= -1.12e-11) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else if (b <= 9.5e-64) {
tmp_1 = (-2.0 * c) / (Math.sqrt(((a * c) * -4.0)) + b);
} else {
tmp_1 = (-2.0 * c) / (2.0 * b);
}
return tmp_1;
}
def code(a, b, c): tmp_1 = 0 if b <= -1.12e-11: tmp_2 = 0 if b >= 0.0: tmp_2 = -b / a else: tmp_2 = (-b + -b) / (2.0 * a) tmp_1 = tmp_2 elif b <= 9.5e-64: tmp_1 = (-2.0 * c) / (math.sqrt(((a * c) * -4.0)) + b) else: tmp_1 = (-2.0 * c) / (2.0 * b) return tmp_1
function code(a, b, c) tmp_1 = 0.0 if (b <= -1.12e-11) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(-b) / a); else tmp_2 = Float64(Float64(Float64(-b) + Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; elseif (b <= 9.5e-64) tmp_1 = Float64(Float64(-2.0 * c) / Float64(sqrt(Float64(Float64(a * c) * -4.0)) + b)); else tmp_1 = Float64(Float64(-2.0 * c) / Float64(2.0 * b)); end return tmp_1 end
function tmp_4 = code(a, b, c) tmp_2 = 0.0; if (b <= -1.12e-11) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = -b / a; else tmp_3 = (-b + -b) / (2.0 * a); end tmp_2 = tmp_3; elseif (b <= 9.5e-64) tmp_2 = (-2.0 * c) / (sqrt(((a * c) * -4.0)) + b); else tmp_2 = (-2.0 * c) / (2.0 * b); end tmp_4 = tmp_2; end
code[a_, b_, c_] := If[LessEqual[b, -1.12e-11], If[GreaterEqual[b, 0.0], N[((-b) / a), $MachinePrecision], N[(N[((-b) + (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 9.5e-64], N[(N[(-2.0 * c), $MachinePrecision] / N[(N[Sqrt[N[(N[(a * c), $MachinePrecision] * -4.0), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision], N[(N[(-2.0 * c), $MachinePrecision] / N[(2.0 * b), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.12 \cdot 10^{-11}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 9.5 \cdot 10^{-64}:\\
\;\;\;\;\frac{-2 \cdot c}{\sqrt{\left(a \cdot c\right) \cdot -4} + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{2 \cdot b}\\
\end{array}
\end{array}
if b < -1.1200000000000001e-11Initial program 55.2%
Taylor expanded in a around 0
distribute-lft-out--N/A
lower-*.f64N/A
lower--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6455.2
Applied rewrites55.2%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6491.3
Applied rewrites91.3%
Taylor expanded in b around -inf
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6491.3
Applied rewrites91.3%
if -1.1200000000000001e-11 < b < 9.50000000000000043e-64Initial program 78.1%
Applied rewrites66.5%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites72.2%
Taylor expanded in a around inf
Applied rewrites62.5%
if 9.50000000000000043e-64 < b Initial program 64.3%
Applied rewrites64.3%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites64.3%
Taylor expanded in a around 0
Applied rewrites90.4%
(FPCore (a b c) :precision binary64 (if (<= b 2.3e-308) (if (>= b 0.0) (/ (- b) a) (/ (+ (- b) (- b)) (* 2.0 a))) (/ (* -2.0 c) (* 2.0 b))))
double code(double a, double b, double c) {
double tmp_1;
if (b <= 2.3e-308) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else {
tmp_1 = (-2.0 * c) / (2.0 * b);
}
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) :: tmp
real(8) :: tmp_1
real(8) :: tmp_2
if (b <= 2.3d-308) then
if (b >= 0.0d0) then
tmp_2 = -b / a
else
tmp_2 = (-b + -b) / (2.0d0 * a)
end if
tmp_1 = tmp_2
else
tmp_1 = ((-2.0d0) * c) / (2.0d0 * b)
end if
code = tmp_1
end function
public static double code(double a, double b, double c) {
double tmp_1;
if (b <= 2.3e-308) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = -b / a;
} else {
tmp_2 = (-b + -b) / (2.0 * a);
}
tmp_1 = tmp_2;
} else {
tmp_1 = (-2.0 * c) / (2.0 * b);
}
return tmp_1;
}
def code(a, b, c): tmp_1 = 0 if b <= 2.3e-308: tmp_2 = 0 if b >= 0.0: tmp_2 = -b / a else: tmp_2 = (-b + -b) / (2.0 * a) tmp_1 = tmp_2 else: tmp_1 = (-2.0 * c) / (2.0 * b) return tmp_1
function code(a, b, c) tmp_1 = 0.0 if (b <= 2.3e-308) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = Float64(Float64(-b) / a); else tmp_2 = Float64(Float64(Float64(-b) + Float64(-b)) / Float64(2.0 * a)); end tmp_1 = tmp_2; else tmp_1 = Float64(Float64(-2.0 * c) / Float64(2.0 * b)); end return tmp_1 end
function tmp_4 = code(a, b, c) tmp_2 = 0.0; if (b <= 2.3e-308) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = -b / a; else tmp_3 = (-b + -b) / (2.0 * a); end tmp_2 = tmp_3; else tmp_2 = (-2.0 * c) / (2.0 * b); end tmp_4 = tmp_2; end
code[a_, b_, c_] := If[LessEqual[b, 2.3e-308], If[GreaterEqual[b, 0.0], N[((-b) / a), $MachinePrecision], N[(N[((-b) + (-b)), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], N[(N[(-2.0 * c), $MachinePrecision] / N[(2.0 * b), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 2.3 \cdot 10^{-308}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-b}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) + \left(-b\right)}{2 \cdot a}\\
\end{array}\\
\mathbf{else}:\\
\;\;\;\;\frac{-2 \cdot c}{2 \cdot b}\\
\end{array}
\end{array}
if b < 2.2999999999999999e-308Initial program 61.0%
Taylor expanded in a around 0
distribute-lft-out--N/A
lower-*.f64N/A
lower--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6461.0
Applied rewrites61.0%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6472.4
Applied rewrites72.4%
Taylor expanded in b around -inf
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6472.4
Applied rewrites72.4%
if 2.2999999999999999e-308 < b Initial program 68.9%
Applied rewrites68.9%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites68.9%
Taylor expanded in a around 0
Applied rewrites69.2%
(FPCore (a b c) :precision binary64 (/ (* -2.0 c) (* 2.0 b)))
double code(double a, double b, double c) {
return (-2.0 * c) / (2.0 * b);
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = ((-2.0d0) * c) / (2.0d0 * b)
end function
public static double code(double a, double b, double c) {
return (-2.0 * c) / (2.0 * b);
}
def code(a, b, c): return (-2.0 * c) / (2.0 * b)
function code(a, b, c) return Float64(Float64(-2.0 * c) / Float64(2.0 * b)) end
function tmp = code(a, b, c) tmp = (-2.0 * c) / (2.0 * b); end
code[a_, b_, c_] := N[(N[(-2.0 * c), $MachinePrecision] / N[(2.0 * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-2 \cdot c}{2 \cdot b}
\end{array}
Initial program 65.1%
Applied rewrites43.1%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites47.0%
Taylor expanded in a around 0
Applied rewrites37.2%
(FPCore (a b c) :precision binary64 (* c (/ -2.0 (* 2.0 b))))
double code(double a, double b, double c) {
return c * (-2.0 / (2.0 * b));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = c * ((-2.0d0) / (2.0d0 * b))
end function
public static double code(double a, double b, double c) {
return c * (-2.0 / (2.0 * b));
}
def code(a, b, c): return c * (-2.0 / (2.0 * b))
function code(a, b, c) return Float64(c * Float64(-2.0 / Float64(2.0 * b))) end
function tmp = code(a, b, c) tmp = c * (-2.0 / (2.0 * b)); end
code[a_, b_, c_] := N[(c * N[(-2.0 / N[(2.0 * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
c \cdot \frac{-2}{2 \cdot b}
\end{array}
Initial program 65.1%
Applied rewrites43.1%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
if-sameN/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites47.0%
Taylor expanded in a around 0
Applied rewrites37.2%
Applied rewrites37.1%
herbie shell --seed 2024326
(FPCore (a b c)
:name "jeff quadratic root 2"
:precision binary64
(if (>= b 0.0) (/ (* 2.0 c) (- (- b) (sqrt (- (* b b) (* (* 4.0 a) c))))) (/ (+ (- b) (sqrt (- (* b b) (* (* 4.0 a) c)))) (* 2.0 a))))