
(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 12 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
(let* ((t_0 (/ (- b) a)))
(if (<= b -2e+151)
(if (>= b 0.0) t_0 t_0)
(if (<= b 7.6e-208)
(if (>= b 0.0)
(* (/ -2.0 (- (sqrt (fma (* -4.0 c) a (* b b))) b)) c)
(/ (- (sqrt (fma (* c a) -4.0 (* b b))) b) (* 2.0 a)))
(if (<= b 1.5e+154)
(if (>= b 0.0)
(/ (* (- c) 2.0) (+ (sqrt (fma b b (* (* -4.0 a) c))) b))
(/ (- (- b) b) (* 2.0 a)))
(if (>= b 0.0) (/ (- c) b) (/ c b)))))))
double code(double a, double b, double c) {
double t_0 = -b / a;
double tmp_1;
if (b <= -2e+151) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_0;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b <= 7.6e-208) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-2.0 / (sqrt(fma((-4.0 * c), a, (b * b))) - b)) * c;
} else {
tmp_3 = (sqrt(fma((c * a), -4.0, (b * b))) - b) / (2.0 * a);
}
tmp_1 = tmp_3;
} else if (b <= 1.5e+154) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (-c * 2.0) / (sqrt(fma(b, b, ((-4.0 * a) * c))) + b);
} else {
tmp_4 = (-b - b) / (2.0 * a);
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = -c / b;
} else {
tmp_1 = c / b;
}
return tmp_1;
}
function code(a, b, c) t_0 = Float64(Float64(-b) / a) tmp_1 = 0.0 if (b <= -2e+151) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = t_0; else tmp_2 = t_0; end tmp_1 = tmp_2; elseif (b <= 7.6e-208) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-2.0 / Float64(sqrt(fma(Float64(-4.0 * c), a, Float64(b * b))) - b)) * c); else tmp_3 = Float64(Float64(sqrt(fma(Float64(c * a), -4.0, Float64(b * b))) - b) / Float64(2.0 * a)); end tmp_1 = tmp_3; elseif (b <= 1.5e+154) tmp_4 = 0.0 if (b >= 0.0) tmp_4 = Float64(Float64(Float64(-c) * 2.0) / Float64(sqrt(fma(b, b, Float64(Float64(-4.0 * a) * c))) + b)); else tmp_4 = Float64(Float64(Float64(-b) - b) / Float64(2.0 * a)); end tmp_1 = tmp_4; elseif (b >= 0.0) tmp_1 = Float64(Float64(-c) / b); else tmp_1 = Float64(c / b); end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[((-b) / a), $MachinePrecision]}, If[LessEqual[b, -2e+151], If[GreaterEqual[b, 0.0], t$95$0, t$95$0], If[LessEqual[b, 7.6e-208], If[GreaterEqual[b, 0.0], N[(N[(-2.0 / N[(N[Sqrt[N[(N[(-4.0 * c), $MachinePrecision] * a + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision] * c), $MachinePrecision], N[(N[(N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0 + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 1.5e+154], If[GreaterEqual[b, 0.0], N[(N[((-c) * 2.0), $MachinePrecision] / N[(N[Sqrt[N[(b * b + N[(N[(-4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision], N[(N[((-b) - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[((-c) / b), $MachinePrecision], N[(c / b), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-b}{a}\\
\mathbf{if}\;b \leq -2 \cdot 10^{+151}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}\\
\mathbf{elif}\;b \leq 7.6 \cdot 10^{-208}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-2}{\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} - b} \cdot c\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{\mathsf{fma}\left(c \cdot a, -4, b \cdot b\right)} - b}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 1.5 \cdot 10^{+154}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-c\right) \cdot 2}{\sqrt{\mathsf{fma}\left(b, b, \left(-4 \cdot a\right) \cdot c\right)} + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) - b}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{-c}{b}\\
\mathbf{else}:\\
\;\;\;\;\frac{c}{b}\\
\end{array}
\end{array}
if b < -2.00000000000000003e151Initial program 40.3%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6440.3
Applied rewrites40.3%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.2
Applied rewrites96.2%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.2
Applied rewrites96.2%
if -2.00000000000000003e151 < b < 7.60000000000000023e-208Initial program 91.4%
Applied rewrites91.4%
lift-+.f64N/A
+-commutativeN/A
lift-neg.f64N/A
unsub-negN/A
lower--.f6491.4
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f64N/A
+-commutativeN/A
lift-fma.f6491.4
Applied rewrites91.4%
if 7.60000000000000023e-208 < b < 1.50000000000000013e154Initial program 82.0%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6482.0
Applied rewrites82.0%
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
cancel-sign-sub-invN/A
lift-*.f64N/A
metadata-evalN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f64N/A
lower-fma.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
lift-*.f64N/A
lower-*.f6483.6
Applied rewrites83.6%
if 1.50000000000000013e154 < b Initial program 32.9%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in b around -inf
associate-*r*N/A
mul-1-negN/A
lower-*.f64N/A
lower-neg.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in c around inf
Applied rewrites100.0%
Final simplification92.1%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- b) a)) (t_1 (sqrt (fma (* c a) -4.0 (* b b)))))
(if (<= b -2e+151)
(if (>= b 0.0) t_0 t_0)
(if (<= b -5e-304)
(if (>= b 0.0)
(* (/ -2.0 (- (sqrt (fma (* -4.0 c) a (* b b))) b)) c)
(/ (- t_1 b) (* 2.0 a)))
(if (<= b 8e+59)
(if (>= b 0.0) (* (/ 2.0 (+ t_1 b)) (- c)) (/ (- (- b) b) (* 2.0 a)))
(if (>= b 0.0) (/ (- c) b) (/ c b)))))))
double code(double a, double b, double c) {
double t_0 = -b / a;
double t_1 = sqrt(fma((c * a), -4.0, (b * b)));
double tmp_1;
if (b <= -2e+151) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_0;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b <= -5e-304) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-2.0 / (sqrt(fma((-4.0 * c), a, (b * b))) - b)) * c;
} else {
tmp_3 = (t_1 - b) / (2.0 * a);
}
tmp_1 = tmp_3;
} else if (b <= 8e+59) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (2.0 / (t_1 + b)) * -c;
} else {
tmp_4 = (-b - b) / (2.0 * a);
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = -c / b;
} else {
tmp_1 = c / b;
}
return tmp_1;
}
function code(a, b, c) t_0 = Float64(Float64(-b) / a) t_1 = sqrt(fma(Float64(c * a), -4.0, Float64(b * b))) tmp_1 = 0.0 if (b <= -2e+151) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = t_0; else tmp_2 = t_0; end tmp_1 = tmp_2; elseif (b <= -5e-304) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-2.0 / Float64(sqrt(fma(Float64(-4.0 * c), a, Float64(b * b))) - b)) * c); else tmp_3 = Float64(Float64(t_1 - b) / Float64(2.0 * a)); end tmp_1 = tmp_3; elseif (b <= 8e+59) tmp_4 = 0.0 if (b >= 0.0) tmp_4 = Float64(Float64(2.0 / Float64(t_1 + b)) * Float64(-c)); else tmp_4 = Float64(Float64(Float64(-b) - b) / Float64(2.0 * a)); end tmp_1 = tmp_4; elseif (b >= 0.0) tmp_1 = Float64(Float64(-c) / b); else tmp_1 = Float64(c / b); end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[((-b) / a), $MachinePrecision]}, Block[{t$95$1 = N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0 + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[b, -2e+151], If[GreaterEqual[b, 0.0], t$95$0, t$95$0], If[LessEqual[b, -5e-304], If[GreaterEqual[b, 0.0], N[(N[(-2.0 / N[(N[Sqrt[N[(N[(-4.0 * c), $MachinePrecision] * a + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision] * c), $MachinePrecision], N[(N[(t$95$1 - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 8e+59], If[GreaterEqual[b, 0.0], N[(N[(2.0 / N[(t$95$1 + b), $MachinePrecision]), $MachinePrecision] * (-c)), $MachinePrecision], N[(N[((-b) - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[((-c) / b), $MachinePrecision], N[(c / b), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-b}{a}\\
t_1 := \sqrt{\mathsf{fma}\left(c \cdot a, -4, b \cdot b\right)}\\
\mathbf{if}\;b \leq -2 \cdot 10^{+151}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}\\
\mathbf{elif}\;b \leq -5 \cdot 10^{-304}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-2}{\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} - b} \cdot c\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_1 - b}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 8 \cdot 10^{+59}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{2}{t\_1 + b} \cdot \left(-c\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) - b}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{-c}{b}\\
\mathbf{else}:\\
\;\;\;\;\frac{c}{b}\\
\end{array}
\end{array}
if b < -2.00000000000000003e151Initial program 40.3%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6440.3
Applied rewrites40.3%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.2
Applied rewrites96.2%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.2
Applied rewrites96.2%
if -2.00000000000000003e151 < b < -4.99999999999999965e-304Initial program 90.6%
Applied rewrites90.6%
lift-+.f64N/A
+-commutativeN/A
lift-neg.f64N/A
unsub-negN/A
lower--.f6490.6
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f64N/A
+-commutativeN/A
lift-fma.f6490.6
Applied rewrites90.6%
if -4.99999999999999965e-304 < b < 7.99999999999999977e59Initial program 83.1%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6483.1
Applied rewrites83.1%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6483.1
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f64N/A
+-commutativeN/A
Applied rewrites83.1%
if 7.99999999999999977e59 < b Initial program 50.1%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6496.1
Applied rewrites96.1%
Taylor expanded in b around -inf
associate-*r*N/A
mul-1-negN/A
lower-*.f64N/A
lower-neg.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6496.1
Applied rewrites96.1%
Taylor expanded in c around inf
Applied rewrites96.1%
Final simplification91.8%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (sqrt (- (* b b) (* (* 4.0 a) c)))) (t_1 (/ (- b) a)))
(if (<= b -2e+151)
(if (>= b 0.0) t_1 t_1)
(if (<= b 8e+59)
(if (>= b 0.0) (/ (* (- c) 2.0) (+ t_0 b)) (/ (- t_0 b) (* 2.0 a)))
(if (>= b 0.0) (/ (- c) b) (/ c b))))))
double code(double a, double b, double c) {
double t_0 = sqrt(((b * b) - ((4.0 * a) * c)));
double t_1 = -b / a;
double tmp_1;
if (b <= -2e+151) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_1;
} else {
tmp_2 = t_1;
}
tmp_1 = tmp_2;
} else if (b <= 8e+59) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-c * 2.0) / (t_0 + b);
} else {
tmp_3 = (t_0 - b) / (2.0 * a);
}
tmp_1 = tmp_3;
} else if (b >= 0.0) {
tmp_1 = -c / b;
} else {
tmp_1 = 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) :: 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 * a) * c)))
t_1 = -b / a
if (b <= (-2d+151)) then
if (b >= 0.0d0) then
tmp_2 = t_1
else
tmp_2 = t_1
end if
tmp_1 = tmp_2
else if (b <= 8d+59) then
if (b >= 0.0d0) then
tmp_3 = (-c * 2.0d0) / (t_0 + b)
else
tmp_3 = (t_0 - b) / (2.0d0 * a)
end if
tmp_1 = tmp_3
else if (b >= 0.0d0) then
tmp_1 = -c / b
else
tmp_1 = c / b
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 * a) * c)));
double t_1 = -b / a;
double tmp_1;
if (b <= -2e+151) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_1;
} else {
tmp_2 = t_1;
}
tmp_1 = tmp_2;
} else if (b <= 8e+59) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-c * 2.0) / (t_0 + b);
} else {
tmp_3 = (t_0 - b) / (2.0 * a);
}
tmp_1 = tmp_3;
} else if (b >= 0.0) {
tmp_1 = -c / b;
} else {
tmp_1 = c / b;
}
return tmp_1;
}
def code(a, b, c): t_0 = math.sqrt(((b * b) - ((4.0 * a) * c))) t_1 = -b / a tmp_1 = 0 if b <= -2e+151: tmp_2 = 0 if b >= 0.0: tmp_2 = t_1 else: tmp_2 = t_1 tmp_1 = tmp_2 elif b <= 8e+59: tmp_3 = 0 if b >= 0.0: tmp_3 = (-c * 2.0) / (t_0 + b) else: tmp_3 = (t_0 - b) / (2.0 * a) tmp_1 = tmp_3 elif b >= 0.0: tmp_1 = -c / b else: tmp_1 = c / b return tmp_1
function code(a, b, c) t_0 = sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c))) t_1 = Float64(Float64(-b) / a) tmp_1 = 0.0 if (b <= -2e+151) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = t_1; else tmp_2 = t_1; end tmp_1 = tmp_2; elseif (b <= 8e+59) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(Float64(-c) * 2.0) / Float64(t_0 + b)); else tmp_3 = Float64(Float64(t_0 - b) / Float64(2.0 * a)); end tmp_1 = tmp_3; elseif (b >= 0.0) tmp_1 = Float64(Float64(-c) / b); else tmp_1 = Float64(c / b); end return tmp_1 end
function tmp_5 = code(a, b, c) t_0 = sqrt(((b * b) - ((4.0 * a) * c))); t_1 = -b / a; tmp_2 = 0.0; if (b <= -2e+151) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = t_1; else tmp_3 = t_1; end tmp_2 = tmp_3; elseif (b <= 8e+59) tmp_4 = 0.0; if (b >= 0.0) tmp_4 = (-c * 2.0) / (t_0 + b); else tmp_4 = (t_0 - b) / (2.0 * a); end tmp_2 = tmp_4; elseif (b >= 0.0) tmp_2 = -c / b; else tmp_2 = c / b; end tmp_5 = tmp_2; 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]}, Block[{t$95$1 = N[((-b) / a), $MachinePrecision]}, If[LessEqual[b, -2e+151], If[GreaterEqual[b, 0.0], t$95$1, t$95$1], If[LessEqual[b, 8e+59], If[GreaterEqual[b, 0.0], N[(N[((-c) * 2.0), $MachinePrecision] / N[(t$95$0 + b), $MachinePrecision]), $MachinePrecision], N[(N[(t$95$0 - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[((-c) / b), $MachinePrecision], N[(c / b), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}\\
t_1 := \frac{-b}{a}\\
\mathbf{if}\;b \leq -2 \cdot 10^{+151}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}\\
\mathbf{elif}\;b \leq 8 \cdot 10^{+59}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-c\right) \cdot 2}{t\_0 + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_0 - b}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{-c}{b}\\
\mathbf{else}:\\
\;\;\;\;\frac{c}{b}\\
\end{array}
\end{array}
if b < -2.00000000000000003e151Initial program 40.3%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6440.3
Applied rewrites40.3%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.2
Applied rewrites96.2%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.2
Applied rewrites96.2%
if -2.00000000000000003e151 < b < 7.99999999999999977e59Initial program 88.2%
if 7.99999999999999977e59 < b Initial program 50.1%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6496.1
Applied rewrites96.1%
Taylor expanded in b around -inf
associate-*r*N/A
mul-1-negN/A
lower-*.f64N/A
lower-neg.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6496.1
Applied rewrites96.1%
Taylor expanded in c around inf
Applied rewrites96.1%
Final simplification91.8%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- b) a)))
(if (<= b -1.45e+113)
(if (>= b 0.0) t_0 t_0)
(if (<= b 8e+59)
(if (>= b 0.0)
(/ (* (- c) 2.0) (+ (sqrt (- (* b b) (* (* 4.0 a) c))) b))
(* (- b (sqrt (fma (* -4.0 c) a (* b b)))) (/ -0.5 a)))
(if (>= b 0.0) (/ (- c) b) (/ c b))))))
double code(double a, double b, double c) {
double t_0 = -b / a;
double tmp_1;
if (b <= -1.45e+113) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_0;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b <= 8e+59) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-c * 2.0) / (sqrt(((b * b) - ((4.0 * a) * c))) + b);
} else {
tmp_3 = (b - sqrt(fma((-4.0 * c), a, (b * b)))) * (-0.5 / a);
}
tmp_1 = tmp_3;
} else if (b >= 0.0) {
tmp_1 = -c / b;
} else {
tmp_1 = c / b;
}
return tmp_1;
}
function code(a, b, c) t_0 = Float64(Float64(-b) / a) tmp_1 = 0.0 if (b <= -1.45e+113) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = t_0; else tmp_2 = t_0; end tmp_1 = tmp_2; elseif (b <= 8e+59) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(Float64(-c) * 2.0) / Float64(sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c))) + b)); else tmp_3 = Float64(Float64(b - sqrt(fma(Float64(-4.0 * c), a, Float64(b * b)))) * Float64(-0.5 / a)); end tmp_1 = tmp_3; elseif (b >= 0.0) tmp_1 = Float64(Float64(-c) / b); else tmp_1 = Float64(c / b); end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[((-b) / a), $MachinePrecision]}, If[LessEqual[b, -1.45e+113], If[GreaterEqual[b, 0.0], t$95$0, t$95$0], If[LessEqual[b, 8e+59], If[GreaterEqual[b, 0.0], N[(N[((-c) * 2.0), $MachinePrecision] / N[(N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision], N[(N[(b - N[Sqrt[N[(N[(-4.0 * c), $MachinePrecision] * a + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(-0.5 / a), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[((-c) / b), $MachinePrecision], N[(c / b), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-b}{a}\\
\mathbf{if}\;b \leq -1.45 \cdot 10^{+113}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}\\
\mathbf{elif}\;b \leq 8 \cdot 10^{+59}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-c\right) \cdot 2}{\sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c} + b}\\
\mathbf{else}:\\
\;\;\;\;\left(b - \sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)}\right) \cdot \frac{-0.5}{a}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{-c}{b}\\
\mathbf{else}:\\
\;\;\;\;\frac{c}{b}\\
\end{array}
\end{array}
if b < -1.44999999999999992e113Initial program 49.7%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6449.7
Applied rewrites49.7%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.8
Applied rewrites96.8%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.8
Applied rewrites96.8%
if -1.44999999999999992e113 < b < 7.99999999999999977e59Initial program 87.4%
lift-/.f64N/A
frac-2negN/A
distribute-frac-negN/A
neg-sub0N/A
lower--.f64N/A
div-invN/A
lower-*.f64N/A
Applied rewrites87.2%
if 7.99999999999999977e59 < b Initial program 50.1%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6496.1
Applied rewrites96.1%
Taylor expanded in b around -inf
associate-*r*N/A
mul-1-negN/A
lower-*.f64N/A
lower-neg.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6496.1
Applied rewrites96.1%
Taylor expanded in c around inf
Applied rewrites96.1%
Final simplification91.7%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- (- b) b) (* 2.0 a))) (t_1 (/ (- c) b)))
(if (<= b -6e-20)
(if (>= b 0.0) t_1 (- (/ c b) (/ b a)))
(if (<= b -2e-310)
(if (>= b 0.0) t_1 (/ (- (sqrt (* (* c a) -4.0)) b) (* 2.0 a)))
(if (<= b 6.2e-27)
(if (>= b 0.0) (* (/ 2.0 (+ (sqrt (* (* -4.0 a) c)) b)) (- c)) t_0)
(if (>= b 0.0) (/ 1.0 (- (/ a b) (/ b c))) t_0))))))
double code(double a, double b, double c) {
double t_0 = (-b - b) / (2.0 * a);
double t_1 = -c / b;
double tmp_1;
if (b <= -6e-20) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_1;
} else {
tmp_2 = (c / b) - (b / a);
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = t_1;
} else {
tmp_3 = (sqrt(((c * a) * -4.0)) - b) / (2.0 * a);
}
tmp_1 = tmp_3;
} else if (b <= 6.2e-27) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (2.0 / (sqrt(((-4.0 * a) * c)) + b)) * -c;
} else {
tmp_4 = t_0;
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = 1.0 / ((a / b) - (b / c));
} 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) :: t_1
real(8) :: tmp
real(8) :: tmp_1
real(8) :: tmp_2
real(8) :: tmp_3
real(8) :: tmp_4
t_0 = (-b - b) / (2.0d0 * a)
t_1 = -c / b
if (b <= (-6d-20)) then
if (b >= 0.0d0) then
tmp_2 = t_1
else
tmp_2 = (c / b) - (b / a)
end if
tmp_1 = tmp_2
else if (b <= (-2d-310)) then
if (b >= 0.0d0) then
tmp_3 = t_1
else
tmp_3 = (sqrt(((c * a) * (-4.0d0))) - b) / (2.0d0 * a)
end if
tmp_1 = tmp_3
else if (b <= 6.2d-27) then
if (b >= 0.0d0) then
tmp_4 = (2.0d0 / (sqrt((((-4.0d0) * a) * c)) + b)) * -c
else
tmp_4 = t_0
end if
tmp_1 = tmp_4
else if (b >= 0.0d0) then
tmp_1 = 1.0d0 / ((a / b) - (b / c))
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 = (-b - b) / (2.0 * a);
double t_1 = -c / b;
double tmp_1;
if (b <= -6e-20) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_1;
} else {
tmp_2 = (c / b) - (b / a);
}
tmp_1 = tmp_2;
} else if (b <= -2e-310) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = t_1;
} else {
tmp_3 = (Math.sqrt(((c * a) * -4.0)) - b) / (2.0 * a);
}
tmp_1 = tmp_3;
} else if (b <= 6.2e-27) {
double tmp_4;
if (b >= 0.0) {
tmp_4 = (2.0 / (Math.sqrt(((-4.0 * a) * c)) + b)) * -c;
} else {
tmp_4 = t_0;
}
tmp_1 = tmp_4;
} else if (b >= 0.0) {
tmp_1 = 1.0 / ((a / b) - (b / c));
} else {
tmp_1 = t_0;
}
return tmp_1;
}
def code(a, b, c): t_0 = (-b - b) / (2.0 * a) t_1 = -c / b tmp_1 = 0 if b <= -6e-20: tmp_2 = 0 if b >= 0.0: tmp_2 = t_1 else: tmp_2 = (c / b) - (b / a) tmp_1 = tmp_2 elif b <= -2e-310: tmp_3 = 0 if b >= 0.0: tmp_3 = t_1 else: tmp_3 = (math.sqrt(((c * a) * -4.0)) - b) / (2.0 * a) tmp_1 = tmp_3 elif b <= 6.2e-27: tmp_4 = 0 if b >= 0.0: tmp_4 = (2.0 / (math.sqrt(((-4.0 * a) * c)) + b)) * -c else: tmp_4 = t_0 tmp_1 = tmp_4 elif b >= 0.0: tmp_1 = 1.0 / ((a / b) - (b / c)) else: tmp_1 = t_0 return tmp_1
function code(a, b, c) t_0 = Float64(Float64(Float64(-b) - b) / Float64(2.0 * a)) t_1 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -6e-20) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = t_1; else tmp_2 = Float64(Float64(c / b) - Float64(b / a)); end tmp_1 = tmp_2; elseif (b <= -2e-310) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = t_1; else tmp_3 = Float64(Float64(sqrt(Float64(Float64(c * a) * -4.0)) - b) / Float64(2.0 * a)); end tmp_1 = tmp_3; elseif (b <= 6.2e-27) tmp_4 = 0.0 if (b >= 0.0) tmp_4 = Float64(Float64(2.0 / Float64(sqrt(Float64(Float64(-4.0 * a) * c)) + b)) * Float64(-c)); else tmp_4 = t_0; end tmp_1 = tmp_4; elseif (b >= 0.0) tmp_1 = Float64(1.0 / Float64(Float64(a / b) - Float64(b / c))); else tmp_1 = t_0; end return tmp_1 end
function tmp_6 = code(a, b, c) t_0 = (-b - b) / (2.0 * a); t_1 = -c / b; tmp_2 = 0.0; if (b <= -6e-20) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = t_1; else tmp_3 = (c / b) - (b / a); end tmp_2 = tmp_3; elseif (b <= -2e-310) tmp_4 = 0.0; if (b >= 0.0) tmp_4 = t_1; else tmp_4 = (sqrt(((c * a) * -4.0)) - b) / (2.0 * a); end tmp_2 = tmp_4; elseif (b <= 6.2e-27) tmp_5 = 0.0; if (b >= 0.0) tmp_5 = (2.0 / (sqrt(((-4.0 * a) * c)) + b)) * -c; else tmp_5 = t_0; end tmp_2 = tmp_5; elseif (b >= 0.0) tmp_2 = 1.0 / ((a / b) - (b / c)); else tmp_2 = t_0; end tmp_6 = tmp_2; end
code[a_, b_, c_] := Block[{t$95$0 = N[(N[((-b) - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -6e-20], If[GreaterEqual[b, 0.0], t$95$1, N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, -2e-310], If[GreaterEqual[b, 0.0], t$95$1, N[(N[(N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[LessEqual[b, 6.2e-27], If[GreaterEqual[b, 0.0], N[(N[(2.0 / N[(N[Sqrt[N[(N[(-4.0 * a), $MachinePrecision] * c), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] * (-c)), $MachinePrecision], t$95$0], If[GreaterEqual[b, 0.0], N[(1.0 / N[(N[(a / b), $MachinePrecision] - N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\left(-b\right) - b}{2 \cdot a}\\
t_1 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -6 \cdot 10^{-20}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\end{array}\\
\mathbf{elif}\;b \leq -2 \cdot 10^{-310}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{\left(c \cdot a\right) \cdot -4} - b}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \leq 6.2 \cdot 10^{-27}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{2}{\sqrt{\left(-4 \cdot a\right) \cdot c} + b} \cdot \left(-c\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{1}{\frac{a}{b} - \frac{b}{c}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if b < -6.00000000000000057e-20Initial program 66.1%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6466.1
Applied rewrites66.1%
Taylor expanded in b around -inf
associate-*r*N/A
mul-1-negN/A
lower-*.f64N/A
lower-neg.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6491.0
Applied rewrites91.0%
Taylor expanded in c around 0
Applied rewrites91.3%
if -6.00000000000000057e-20 < b < -1.999999999999994e-310Initial program 88.0%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6488.0
Applied rewrites88.0%
Taylor expanded in c around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6471.9
Applied rewrites71.9%
if -1.999999999999994e-310 < b < 6.1999999999999997e-27Initial program 82.1%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6482.1
Applied rewrites82.1%
Taylor expanded in c around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6469.3
Applied rewrites69.3%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6469.4
Applied rewrites69.4%
if 6.1999999999999997e-27 < b Initial program 54.0%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6454.0
Applied rewrites54.0%
Taylor expanded in c around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6429.2
Applied rewrites29.2%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6429.2
Applied rewrites29.2%
Taylor expanded in a around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f6492.7
Applied rewrites92.7%
Final simplification84.8%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- b) a)))
(if (<= b -2e+151)
(if (>= b 0.0) t_0 t_0)
(if (<= b 1.5e+154)
(if (>= b 0.0)
(/ (* (- c) 2.0) (+ (sqrt (fma b b (* (* -4.0 c) a))) b))
(/ (- (sqrt (- (* b b) (* (* 4.0 a) c))) b) (* 2.0 a)))
(if (>= b 0.0) (/ (- c) b) (/ c b))))))
double code(double a, double b, double c) {
double t_0 = -b / a;
double tmp_1;
if (b <= -2e+151) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_0;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b <= 1.5e+154) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-c * 2.0) / (sqrt(fma(b, b, ((-4.0 * c) * a))) + b);
} else {
tmp_3 = (sqrt(((b * b) - ((4.0 * a) * c))) - b) / (2.0 * a);
}
tmp_1 = tmp_3;
} else if (b >= 0.0) {
tmp_1 = -c / b;
} else {
tmp_1 = c / b;
}
return tmp_1;
}
function code(a, b, c) t_0 = Float64(Float64(-b) / a) tmp_1 = 0.0 if (b <= -2e+151) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = t_0; else tmp_2 = t_0; end tmp_1 = tmp_2; elseif (b <= 1.5e+154) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(Float64(-c) * 2.0) / Float64(sqrt(fma(b, b, Float64(Float64(-4.0 * c) * a))) + b)); else tmp_3 = Float64(Float64(sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c))) - b) / Float64(2.0 * a)); end tmp_1 = tmp_3; elseif (b >= 0.0) tmp_1 = Float64(Float64(-c) / b); else tmp_1 = Float64(c / b); end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[((-b) / a), $MachinePrecision]}, If[LessEqual[b, -2e+151], If[GreaterEqual[b, 0.0], t$95$0, t$95$0], If[LessEqual[b, 1.5e+154], If[GreaterEqual[b, 0.0], N[(N[((-c) * 2.0), $MachinePrecision] / N[(N[Sqrt[N[(b * b + N[(N[(-4.0 * c), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision], N[(N[(N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[((-c) / b), $MachinePrecision], N[(c / b), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-b}{a}\\
\mathbf{if}\;b \leq -2 \cdot 10^{+151}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}\\
\mathbf{elif}\;b \leq 1.5 \cdot 10^{+154}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{\left(-c\right) \cdot 2}{\sqrt{\mathsf{fma}\left(b, b, \left(-4 \cdot c\right) \cdot a\right)} + b}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c} - b}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{-c}{b}\\
\mathbf{else}:\\
\;\;\;\;\frac{c}{b}\\
\end{array}
\end{array}
if b < -2.00000000000000003e151Initial program 40.3%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6440.3
Applied rewrites40.3%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.2
Applied rewrites96.2%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.2
Applied rewrites96.2%
if -2.00000000000000003e151 < b < 1.50000000000000013e154Initial program 88.0%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
lower-fma.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
distribute-lft-neg-inN/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
metadata-eval88.6
Applied rewrites88.6%
if 1.50000000000000013e154 < b Initial program 32.9%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in b around -inf
associate-*r*N/A
mul-1-negN/A
lower-*.f64N/A
lower-neg.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in c around inf
Applied rewrites100.0%
Final simplification92.1%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- b) a)))
(if (<= b -2e+151)
(if (>= b 0.0) t_0 t_0)
(if (<= b 6.2e-27)
(if (>= b 0.0)
(* (/ -2.0 (- (sqrt (fma (* -4.0 c) a (* b b))) b)) c)
(/ (- (sqrt (fma (* c a) -4.0 (* b b))) b) (* 2.0 a)))
(if (>= b 0.0)
(/ 1.0 (- (/ a b) (/ b c)))
(/ (- (- b) b) (* 2.0 a)))))))
double code(double a, double b, double c) {
double t_0 = -b / a;
double tmp_1;
if (b <= -2e+151) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_0;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b <= 6.2e-27) {
double tmp_3;
if (b >= 0.0) {
tmp_3 = (-2.0 / (sqrt(fma((-4.0 * c), a, (b * b))) - b)) * c;
} else {
tmp_3 = (sqrt(fma((c * a), -4.0, (b * b))) - b) / (2.0 * a);
}
tmp_1 = tmp_3;
} else if (b >= 0.0) {
tmp_1 = 1.0 / ((a / b) - (b / c));
} else {
tmp_1 = (-b - b) / (2.0 * a);
}
return tmp_1;
}
function code(a, b, c) t_0 = Float64(Float64(-b) / a) tmp_1 = 0.0 if (b <= -2e+151) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = t_0; else tmp_2 = t_0; end tmp_1 = tmp_2; elseif (b <= 6.2e-27) tmp_3 = 0.0 if (b >= 0.0) tmp_3 = Float64(Float64(-2.0 / Float64(sqrt(fma(Float64(-4.0 * c), a, Float64(b * b))) - b)) * c); else tmp_3 = Float64(Float64(sqrt(fma(Float64(c * a), -4.0, Float64(b * b))) - b) / Float64(2.0 * a)); end tmp_1 = tmp_3; elseif (b >= 0.0) tmp_1 = Float64(1.0 / Float64(Float64(a / b) - Float64(b / c))); else tmp_1 = Float64(Float64(Float64(-b) - b) / Float64(2.0 * a)); end return tmp_1 end
code[a_, b_, c_] := Block[{t$95$0 = N[((-b) / a), $MachinePrecision]}, If[LessEqual[b, -2e+151], If[GreaterEqual[b, 0.0], t$95$0, t$95$0], If[LessEqual[b, 6.2e-27], If[GreaterEqual[b, 0.0], N[(N[(-2.0 / N[(N[Sqrt[N[(N[(-4.0 * c), $MachinePrecision] * a + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision] * c), $MachinePrecision], N[(N[(N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0 + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], N[(1.0 / N[(N[(a / b), $MachinePrecision] - N[(b / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[((-b) - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-b}{a}\\
\mathbf{if}\;b \leq -2 \cdot 10^{+151}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}\\
\mathbf{elif}\;b \leq 6.2 \cdot 10^{-27}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-2}{\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} - b} \cdot c\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{\mathsf{fma}\left(c \cdot a, -4, b \cdot b\right)} - b}{2 \cdot a}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{1}{\frac{a}{b} - \frac{b}{c}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-b\right) - b}{2 \cdot a}\\
\end{array}
\end{array}
if b < -2.00000000000000003e151Initial program 40.3%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6440.3
Applied rewrites40.3%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.2
Applied rewrites96.2%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6496.2
Applied rewrites96.2%
if -2.00000000000000003e151 < b < 6.1999999999999997e-27Initial program 88.2%
Applied rewrites83.4%
lift-+.f64N/A
+-commutativeN/A
lift-neg.f64N/A
unsub-negN/A
lower--.f6483.4
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f64N/A
+-commutativeN/A
lift-fma.f6483.4
Applied rewrites83.4%
if 6.1999999999999997e-27 < b Initial program 54.0%
Taylor expanded in b around -inf
mul-1-negN/A
lower-neg.f6454.0
Applied rewrites54.0%
Taylor expanded in c around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6429.2
Applied rewrites29.2%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6429.2
Applied rewrites29.2%
Taylor expanded in a around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f6492.7
Applied rewrites92.7%
Final simplification88.6%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- c) b)))
(if (<= b -6e-20)
(if (>= b 0.0) t_0 (- (/ c b) (/ b a)))
(if (>= b 0.0) t_0 (/ (- (sqrt (* (* c a) -4.0)) b) (* 2.0 a))))))
double code(double a, double b, double c) {
double t_0 = -c / b;
double tmp_1;
if (b <= -6e-20) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_0;
} else {
tmp_2 = (c / b) - (b / a);
}
tmp_1 = tmp_2;
} else if (b >= 0.0) {
tmp_1 = t_0;
} else {
tmp_1 = (sqrt(((c * a) * -4.0)) - b) / (2.0 * a);
}
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
t_0 = -c / b
if (b <= (-6d-20)) then
if (b >= 0.0d0) then
tmp_2 = t_0
else
tmp_2 = (c / b) - (b / a)
end if
tmp_1 = tmp_2
else if (b >= 0.0d0) then
tmp_1 = t_0
else
tmp_1 = (sqrt(((c * a) * (-4.0d0))) - b) / (2.0d0 * a)
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 <= -6e-20) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_0;
} else {
tmp_2 = (c / b) - (b / a);
}
tmp_1 = tmp_2;
} else if (b >= 0.0) {
tmp_1 = t_0;
} else {
tmp_1 = (Math.sqrt(((c * a) * -4.0)) - b) / (2.0 * a);
}
return tmp_1;
}
def code(a, b, c): t_0 = -c / b tmp_1 = 0 if b <= -6e-20: tmp_2 = 0 if b >= 0.0: tmp_2 = t_0 else: tmp_2 = (c / b) - (b / a) tmp_1 = tmp_2 elif b >= 0.0: tmp_1 = t_0 else: tmp_1 = (math.sqrt(((c * a) * -4.0)) - b) / (2.0 * a) return tmp_1
function code(a, b, c) t_0 = Float64(Float64(-c) / b) tmp_1 = 0.0 if (b <= -6e-20) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = t_0; else tmp_2 = Float64(Float64(c / b) - Float64(b / a)); end tmp_1 = tmp_2; elseif (b >= 0.0) tmp_1 = t_0; else tmp_1 = Float64(Float64(sqrt(Float64(Float64(c * a) * -4.0)) - b) / Float64(2.0 * a)); end return tmp_1 end
function tmp_4 = code(a, b, c) t_0 = -c / b; tmp_2 = 0.0; if (b <= -6e-20) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = t_0; else tmp_3 = (c / b) - (b / a); end tmp_2 = tmp_3; elseif (b >= 0.0) tmp_2 = t_0; else tmp_2 = (sqrt(((c * a) * -4.0)) - b) / (2.0 * a); end tmp_4 = tmp_2; end
code[a_, b_, c_] := Block[{t$95$0 = N[((-c) / b), $MachinePrecision]}, If[LessEqual[b, -6e-20], If[GreaterEqual[b, 0.0], t$95$0, N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision]], If[GreaterEqual[b, 0.0], t$95$0, N[(N[(N[Sqrt[N[(N[(c * a), $MachinePrecision] * -4.0), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-c}{b}\\
\mathbf{if}\;b \leq -6 \cdot 10^{-20}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{\left(c \cdot a\right) \cdot -4} - b}{2 \cdot a}\\
\end{array}
\end{array}
if b < -6.00000000000000057e-20Initial program 66.1%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6466.1
Applied rewrites66.1%
Taylor expanded in b around -inf
associate-*r*N/A
mul-1-negN/A
lower-*.f64N/A
lower-neg.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6491.0
Applied rewrites91.0%
Taylor expanded in c around 0
Applied rewrites91.3%
if -6.00000000000000057e-20 < b Initial program 70.6%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6476.5
Applied rewrites76.5%
Taylor expanded in c around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6471.7
Applied rewrites71.7%
Final simplification78.8%
(FPCore (a b c) :precision binary64 (if (>= b 0.0) (/ (- c) b) (- (/ c b) (/ b a))))
double code(double a, double b, double c) {
double tmp;
if (b >= 0.0) {
tmp = -c / b;
} else {
tmp = (c / b) - (b / 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) :: tmp
if (b >= 0.0d0) then
tmp = -c / b
else
tmp = (c / b) - (b / a)
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) - (b / a);
}
return tmp;
}
def code(a, b, c): tmp = 0 if b >= 0.0: tmp = -c / b else: tmp = (c / b) - (b / a) return tmp
function code(a, b, c) tmp = 0.0 if (b >= 0.0) tmp = Float64(Float64(-c) / b); else tmp = Float64(Float64(c / b) - Float64(b / a)); 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) - (b / a); end tmp_2 = tmp; end
code[a_, b_, c_] := If[GreaterEqual[b, 0.0], N[((-c) / b), $MachinePrecision], N[(N[(c / b), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-c}{b}\\
\mathbf{else}:\\
\;\;\;\;\frac{c}{b} - \frac{b}{a}\\
\end{array}
\end{array}
Initial program 68.9%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6472.7
Applied rewrites72.7%
Taylor expanded in b around -inf
associate-*r*N/A
mul-1-negN/A
lower-*.f64N/A
lower-neg.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6469.0
Applied rewrites69.0%
Taylor expanded in c around 0
Applied rewrites69.3%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (/ (- b) a)))
(if (<= b -2.15e-303)
(if (>= b 0.0) t_0 t_0)
(if (>= b 0.0) (/ (- c) b) (/ c b)))))
double code(double a, double b, double c) {
double t_0 = -b / a;
double tmp_1;
if (b <= -2.15e-303) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_0;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b >= 0.0) {
tmp_1 = -c / b;
} else {
tmp_1 = 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) :: t_0
real(8) :: tmp
real(8) :: tmp_1
real(8) :: tmp_2
t_0 = -b / a
if (b <= (-2.15d-303)) then
if (b >= 0.0d0) then
tmp_2 = t_0
else
tmp_2 = t_0
end if
tmp_1 = tmp_2
else if (b >= 0.0d0) then
tmp_1 = -c / b
else
tmp_1 = c / b
end if
code = tmp_1
end function
public static double code(double a, double b, double c) {
double t_0 = -b / a;
double tmp_1;
if (b <= -2.15e-303) {
double tmp_2;
if (b >= 0.0) {
tmp_2 = t_0;
} else {
tmp_2 = t_0;
}
tmp_1 = tmp_2;
} else if (b >= 0.0) {
tmp_1 = -c / b;
} else {
tmp_1 = c / b;
}
return tmp_1;
}
def code(a, b, c): t_0 = -b / a tmp_1 = 0 if b <= -2.15e-303: tmp_2 = 0 if b >= 0.0: tmp_2 = t_0 else: tmp_2 = t_0 tmp_1 = tmp_2 elif b >= 0.0: tmp_1 = -c / b else: tmp_1 = c / b return tmp_1
function code(a, b, c) t_0 = Float64(Float64(-b) / a) tmp_1 = 0.0 if (b <= -2.15e-303) tmp_2 = 0.0 if (b >= 0.0) tmp_2 = t_0; else tmp_2 = t_0; end tmp_1 = tmp_2; elseif (b >= 0.0) tmp_1 = Float64(Float64(-c) / b); else tmp_1 = Float64(c / b); end return tmp_1 end
function tmp_4 = code(a, b, c) t_0 = -b / a; tmp_2 = 0.0; if (b <= -2.15e-303) tmp_3 = 0.0; if (b >= 0.0) tmp_3 = t_0; else tmp_3 = t_0; end tmp_2 = tmp_3; elseif (b >= 0.0) tmp_2 = -c / b; else tmp_2 = c / b; end tmp_4 = tmp_2; end
code[a_, b_, c_] := Block[{t$95$0 = N[((-b) / a), $MachinePrecision]}, If[LessEqual[b, -2.15e-303], If[GreaterEqual[b, 0.0], t$95$0, t$95$0], If[GreaterEqual[b, 0.0], N[((-c) / b), $MachinePrecision], N[(c / b), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-b}{a}\\
\mathbf{if}\;b \leq -2.15 \cdot 10^{-303}:\\
\;\;\;\;\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}\\
\mathbf{elif}\;b \geq 0:\\
\;\;\;\;\frac{-c}{b}\\
\mathbf{else}:\\
\;\;\;\;\frac{c}{b}\\
\end{array}
\end{array}
if b < -2.14999999999999991e-303Initial program 73.6%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6473.6
Applied rewrites73.6%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6467.1
Applied rewrites67.1%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6467.1
Applied rewrites67.1%
if -2.14999999999999991e-303 < b Initial program 63.1%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6471.6
Applied rewrites71.6%
Taylor expanded in b around -inf
associate-*r*N/A
mul-1-negN/A
lower-*.f64N/A
lower-neg.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6471.6
Applied rewrites71.6%
Taylor expanded in c around inf
Applied rewrites71.6%
(FPCore (a b c) :precision binary64 (if (>= b 0.0) (/ (- c) b) (/ (- b) a)))
double code(double a, double b, double c) {
double tmp;
if (b >= 0.0) {
tmp = -c / b;
} else {
tmp = -b / 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) :: tmp
if (b >= 0.0d0) then
tmp = -c / b
else
tmp = -b / a
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 = -b / a;
}
return tmp;
}
def code(a, b, c): tmp = 0 if b >= 0.0: tmp = -c / b else: tmp = -b / a return tmp
function code(a, b, c) tmp = 0.0 if (b >= 0.0) tmp = Float64(Float64(-c) / b); else tmp = Float64(Float64(-b) / a); end return tmp end
function tmp_2 = code(a, b, c) tmp = 0.0; if (b >= 0.0) tmp = -c / b; else tmp = -b / a; end tmp_2 = tmp; end
code[a_, b_, c_] := If[GreaterEqual[b, 0.0], N[((-c) / b), $MachinePrecision], N[((-b) / a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \geq 0:\\
\;\;\;\;\frac{-c}{b}\\
\mathbf{else}:\\
\;\;\;\;\frac{-b}{a}\\
\end{array}
\end{array}
Initial program 68.9%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6472.7
Applied rewrites72.7%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6469.1
Applied rewrites69.1%
(FPCore (a b c) :precision binary64 (let* ((t_0 (/ (- b) a))) (if (>= b 0.0) t_0 t_0)))
double code(double a, double b, double c) {
double t_0 = -b / a;
double tmp;
if (b >= 0.0) {
tmp = t_0;
} else {
tmp = 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 = -b / a
if (b >= 0.0d0) then
tmp = t_0
else
tmp = t_0
end if
code = tmp
end function
public static double code(double a, double b, double c) {
double t_0 = -b / a;
double tmp;
if (b >= 0.0) {
tmp = t_0;
} else {
tmp = t_0;
}
return tmp;
}
def code(a, b, c): t_0 = -b / a tmp = 0 if b >= 0.0: tmp = t_0 else: tmp = t_0 return tmp
function code(a, b, c) t_0 = Float64(Float64(-b) / a) tmp = 0.0 if (b >= 0.0) tmp = t_0; else tmp = t_0; end return tmp end
function tmp_2 = code(a, b, c) t_0 = -b / a; tmp = 0.0; if (b >= 0.0) tmp = t_0; else tmp = t_0; end tmp_2 = tmp; end
code[a_, b_, c_] := Block[{t$95$0 = N[((-b) / a), $MachinePrecision]}, If[GreaterEqual[b, 0.0], t$95$0, t$95$0]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-b}{a}\\
\mathbf{if}\;b \geq 0:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
Initial program 68.9%
Taylor expanded in c around 0
associate-*r/N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f6472.7
Applied rewrites72.7%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6469.1
Applied rewrites69.1%
Taylor expanded in b around -inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6438.4
Applied rewrites38.4%
herbie shell --seed 2024277
(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))))