
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (pow B 2.0) (* (* 4.0 A) C))))
(/
(-
(sqrt
(*
(* 2.0 (* t_0 F))
(+ (+ A C) (sqrt (+ (pow (- A C) 2.0) (pow B 2.0)))))))
t_0)))
double code(double A, double B, double C, double F) {
double t_0 = pow(B, 2.0) - ((4.0 * A) * C);
return -sqrt(((2.0 * (t_0 * F)) * ((A + C) + sqrt((pow((A - C), 2.0) + pow(B, 2.0)))))) / t_0;
}
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
t_0 = (b ** 2.0d0) - ((4.0d0 * a) * c)
code = -sqrt(((2.0d0 * (t_0 * f)) * ((a + c) + sqrt((((a - c) ** 2.0d0) + (b ** 2.0d0)))))) / t_0
end function
public static double code(double A, double B, double C, double F) {
double t_0 = Math.pow(B, 2.0) - ((4.0 * A) * C);
return -Math.sqrt(((2.0 * (t_0 * F)) * ((A + C) + Math.sqrt((Math.pow((A - C), 2.0) + Math.pow(B, 2.0)))))) / t_0;
}
def code(A, B, C, F): t_0 = math.pow(B, 2.0) - ((4.0 * A) * C) return -math.sqrt(((2.0 * (t_0 * F)) * ((A + C) + math.sqrt((math.pow((A - C), 2.0) + math.pow(B, 2.0)))))) / t_0
function code(A, B, C, F) t_0 = Float64((B ^ 2.0) - Float64(Float64(4.0 * A) * C)) return Float64(Float64(-sqrt(Float64(Float64(2.0 * Float64(t_0 * F)) * Float64(Float64(A + C) + sqrt(Float64((Float64(A - C) ^ 2.0) + (B ^ 2.0))))))) / t_0) end
function tmp = code(A, B, C, F) t_0 = (B ^ 2.0) - ((4.0 * A) * C); tmp = -sqrt(((2.0 * (t_0 * F)) * ((A + C) + sqrt((((A - C) ^ 2.0) + (B ^ 2.0)))))) / t_0; end
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[Power[B, 2.0], $MachinePrecision] - N[(N[(4.0 * A), $MachinePrecision] * C), $MachinePrecision]), $MachinePrecision]}, N[((-N[Sqrt[N[(N[(2.0 * N[(t$95$0 * F), $MachinePrecision]), $MachinePrecision] * N[(N[(A + C), $MachinePrecision] + N[Sqrt[N[(N[Power[N[(A - C), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[B, 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {B}^{2} - \left(4 \cdot A\right) \cdot C\\
\frac{-\sqrt{\left(2 \cdot \left(t_0 \cdot F\right)\right) \cdot \left(\left(A + C\right) + \sqrt{{\left(A - C\right)}^{2} + {B}^{2}}\right)}}{t_0}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 23 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (pow B 2.0) (* (* 4.0 A) C))))
(/
(-
(sqrt
(*
(* 2.0 (* t_0 F))
(+ (+ A C) (sqrt (+ (pow (- A C) 2.0) (pow B 2.0)))))))
t_0)))
double code(double A, double B, double C, double F) {
double t_0 = pow(B, 2.0) - ((4.0 * A) * C);
return -sqrt(((2.0 * (t_0 * F)) * ((A + C) + sqrt((pow((A - C), 2.0) + pow(B, 2.0)))))) / t_0;
}
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
t_0 = (b ** 2.0d0) - ((4.0d0 * a) * c)
code = -sqrt(((2.0d0 * (t_0 * f)) * ((a + c) + sqrt((((a - c) ** 2.0d0) + (b ** 2.0d0)))))) / t_0
end function
public static double code(double A, double B, double C, double F) {
double t_0 = Math.pow(B, 2.0) - ((4.0 * A) * C);
return -Math.sqrt(((2.0 * (t_0 * F)) * ((A + C) + Math.sqrt((Math.pow((A - C), 2.0) + Math.pow(B, 2.0)))))) / t_0;
}
def code(A, B, C, F): t_0 = math.pow(B, 2.0) - ((4.0 * A) * C) return -math.sqrt(((2.0 * (t_0 * F)) * ((A + C) + math.sqrt((math.pow((A - C), 2.0) + math.pow(B, 2.0)))))) / t_0
function code(A, B, C, F) t_0 = Float64((B ^ 2.0) - Float64(Float64(4.0 * A) * C)) return Float64(Float64(-sqrt(Float64(Float64(2.0 * Float64(t_0 * F)) * Float64(Float64(A + C) + sqrt(Float64((Float64(A - C) ^ 2.0) + (B ^ 2.0))))))) / t_0) end
function tmp = code(A, B, C, F) t_0 = (B ^ 2.0) - ((4.0 * A) * C); tmp = -sqrt(((2.0 * (t_0 * F)) * ((A + C) + sqrt((((A - C) ^ 2.0) + (B ^ 2.0)))))) / t_0; end
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[Power[B, 2.0], $MachinePrecision] - N[(N[(4.0 * A), $MachinePrecision] * C), $MachinePrecision]), $MachinePrecision]}, N[((-N[Sqrt[N[(N[(2.0 * N[(t$95$0 * F), $MachinePrecision]), $MachinePrecision] * N[(N[(A + C), $MachinePrecision] + N[Sqrt[N[(N[Power[N[(A - C), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[B, 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {B}^{2} - \left(4 \cdot A\right) \cdot C\\
\frac{-\sqrt{\left(2 \cdot \left(t_0 \cdot F\right)\right) \cdot \left(\left(A + C\right) + \sqrt{{\left(A - C\right)}^{2} + {B}^{2}}\right)}}{t_0}
\end{array}
\end{array}
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (fma B B (* -4.0 (* A C))))
(t_1 (- (pow B 2.0) (* (* 4.0 A) C)))
(t_2
(/
(-
(sqrt
(*
(* 2.0 (* t_1 F))
(+ (+ A C) (sqrt (+ (pow B 2.0) (pow (- A C) 2.0)))))))
t_1))
(t_3 (sqrt (+ C (+ A (hypot B (- A C))))))
(t_4 (fma B B (* C (* A -4.0)))))
(if (<= t_2 -1e-224)
(/
(* (* (sqrt 2.0) (* (sqrt F) (sqrt t_0))) (- t_3))
(- (* B B) (* 4.0 (* A C))))
(if (<= t_2 5e-102)
(/
(- (sqrt (* (* 2.0 (* F t_4)) (+ A (+ A (* -0.5 (/ (* B B) C)))))))
t_4)
(if (<= t_2 INFINITY)
(- (/ (sqrt (* 2.0 (* F t_0))) (/ t_0 t_3)))
(* (/ (sqrt 2.0) B) (- (sqrt (* F (+ C (hypot B C)))))))))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = fma(B, B, (-4.0 * (A * C)));
double t_1 = pow(B, 2.0) - ((4.0 * A) * C);
double t_2 = -sqrt(((2.0 * (t_1 * F)) * ((A + C) + sqrt((pow(B, 2.0) + pow((A - C), 2.0)))))) / t_1;
double t_3 = sqrt((C + (A + hypot(B, (A - C)))));
double t_4 = fma(B, B, (C * (A * -4.0)));
double tmp;
if (t_2 <= -1e-224) {
tmp = ((sqrt(2.0) * (sqrt(F) * sqrt(t_0))) * -t_3) / ((B * B) - (4.0 * (A * C)));
} else if (t_2 <= 5e-102) {
tmp = -sqrt(((2.0 * (F * t_4)) * (A + (A + (-0.5 * ((B * B) / C)))))) / t_4;
} else if (t_2 <= ((double) INFINITY)) {
tmp = -(sqrt((2.0 * (F * t_0))) / (t_0 / t_3));
} else {
tmp = (sqrt(2.0) / B) * -sqrt((F * (C + hypot(B, C))));
}
return tmp;
}
B = abs(B) function code(A, B, C, F) t_0 = fma(B, B, Float64(-4.0 * Float64(A * C))) t_1 = Float64((B ^ 2.0) - Float64(Float64(4.0 * A) * C)) t_2 = Float64(Float64(-sqrt(Float64(Float64(2.0 * Float64(t_1 * F)) * Float64(Float64(A + C) + sqrt(Float64((B ^ 2.0) + (Float64(A - C) ^ 2.0))))))) / t_1) t_3 = sqrt(Float64(C + Float64(A + hypot(B, Float64(A - C))))) t_4 = fma(B, B, Float64(C * Float64(A * -4.0))) tmp = 0.0 if (t_2 <= -1e-224) tmp = Float64(Float64(Float64(sqrt(2.0) * Float64(sqrt(F) * sqrt(t_0))) * Float64(-t_3)) / Float64(Float64(B * B) - Float64(4.0 * Float64(A * C)))); elseif (t_2 <= 5e-102) tmp = Float64(Float64(-sqrt(Float64(Float64(2.0 * Float64(F * t_4)) * Float64(A + Float64(A + Float64(-0.5 * Float64(Float64(B * B) / C))))))) / t_4); elseif (t_2 <= Inf) tmp = Float64(-Float64(sqrt(Float64(2.0 * Float64(F * t_0))) / Float64(t_0 / t_3))); else tmp = Float64(Float64(sqrt(2.0) / B) * Float64(-sqrt(Float64(F * Float64(C + hypot(B, C)))))); end return tmp end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(B * B + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Power[B, 2.0], $MachinePrecision] - N[(N[(4.0 * A), $MachinePrecision] * C), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[((-N[Sqrt[N[(N[(2.0 * N[(t$95$1 * F), $MachinePrecision]), $MachinePrecision] * N[(N[(A + C), $MachinePrecision] + N[Sqrt[N[(N[Power[B, 2.0], $MachinePrecision] + N[Power[N[(A - C), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[Sqrt[N[(C + N[(A + N[Sqrt[B ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$4 = N[(B * B + N[(C * N[(A * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -1e-224], N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] * N[(N[Sqrt[F], $MachinePrecision] * N[Sqrt[t$95$0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * (-t$95$3)), $MachinePrecision] / N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 5e-102], N[((-N[Sqrt[N[(N[(2.0 * N[(F * t$95$4), $MachinePrecision]), $MachinePrecision] * N[(A + N[(A + N[(-0.5 * N[(N[(B * B), $MachinePrecision] / C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$4), $MachinePrecision], If[LessEqual[t$95$2, Infinity], (-N[(N[Sqrt[N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(t$95$0 / t$95$3), $MachinePrecision]), $MachinePrecision]), N[(N[(N[Sqrt[2.0], $MachinePrecision] / B), $MachinePrecision] * (-N[Sqrt[N[(F * N[(C + N[Sqrt[B ^ 2 + C ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]]]]]]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(B, B, -4 \cdot \left(A \cdot C\right)\right)\\
t_1 := {B}^{2} - \left(4 \cdot A\right) \cdot C\\
t_2 := \frac{-\sqrt{\left(2 \cdot \left(t_1 \cdot F\right)\right) \cdot \left(\left(A + C\right) + \sqrt{{B}^{2} + {\left(A - C\right)}^{2}}\right)}}{t_1}\\
t_3 := \sqrt{C + \left(A + \mathsf{hypot}\left(B, A - C\right)\right)}\\
t_4 := \mathsf{fma}\left(B, B, C \cdot \left(A \cdot -4\right)\right)\\
\mathbf{if}\;t_2 \leq -1 \cdot 10^{-224}:\\
\;\;\;\;\frac{\left(\sqrt{2} \cdot \left(\sqrt{F} \cdot \sqrt{t_0}\right)\right) \cdot \left(-t_3\right)}{B \cdot B - 4 \cdot \left(A \cdot C\right)}\\
\mathbf{elif}\;t_2 \leq 5 \cdot 10^{-102}:\\
\;\;\;\;\frac{-\sqrt{\left(2 \cdot \left(F \cdot t_4\right)\right) \cdot \left(A + \left(A + -0.5 \cdot \frac{B \cdot B}{C}\right)\right)}}{t_4}\\
\mathbf{elif}\;t_2 \leq \infty:\\
\;\;\;\;-\frac{\sqrt{2 \cdot \left(F \cdot t_0\right)}}{\frac{t_0}{t_3}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2}}{B} \cdot \left(-\sqrt{F \cdot \left(C + \mathsf{hypot}\left(B, C\right)\right)}\right)\\
\end{array}
\end{array}
if (/.f64 (neg.f64 (sqrt.f64 (*.f64 (*.f64 2 (*.f64 (-.f64 (pow.f64 B 2) (*.f64 (*.f64 4 A) C)) F)) (+.f64 (+.f64 A C) (sqrt.f64 (+.f64 (pow.f64 (-.f64 A C) 2) (pow.f64 B 2))))))) (-.f64 (pow.f64 B 2) (*.f64 (*.f64 4 A) C))) < -1e-224Initial program 48.9%
associate-*l*48.9%
unpow248.9%
+-commutative48.9%
unpow248.9%
associate-*l*48.9%
unpow248.9%
Simplified48.9%
sqrt-prod54.1%
*-commutative54.1%
*-commutative54.1%
associate-+l+54.1%
unpow254.1%
hypot-udef72.3%
associate-+r+72.0%
+-commutative72.0%
associate-+r+72.6%
Applied egg-rr72.6%
sqrt-prod72.4%
cancel-sign-sub-inv72.4%
metadata-eval72.4%
Applied egg-rr72.4%
sqrt-prod82.3%
fma-def82.3%
Applied egg-rr82.3%
if -1e-224 < (/.f64 (neg.f64 (sqrt.f64 (*.f64 (*.f64 2 (*.f64 (-.f64 (pow.f64 B 2) (*.f64 (*.f64 4 A) C)) F)) (+.f64 (+.f64 A C) (sqrt.f64 (+.f64 (pow.f64 (-.f64 A C) 2) (pow.f64 B 2))))))) (-.f64 (pow.f64 B 2) (*.f64 (*.f64 4 A) C))) < 5.00000000000000026e-102Initial program 6.3%
Simplified10.5%
Taylor expanded in C around -inf 35.3%
unpow235.3%
Simplified35.3%
if 5.00000000000000026e-102 < (/.f64 (neg.f64 (sqrt.f64 (*.f64 (*.f64 2 (*.f64 (-.f64 (pow.f64 B 2) (*.f64 (*.f64 4 A) C)) F)) (+.f64 (+.f64 A C) (sqrt.f64 (+.f64 (pow.f64 (-.f64 A C) 2) (pow.f64 B 2))))))) (-.f64 (pow.f64 B 2) (*.f64 (*.f64 4 A) C))) < +inf.0Initial program 36.0%
associate-*l*36.0%
unpow236.0%
+-commutative36.0%
unpow236.0%
associate-*l*36.0%
unpow236.0%
Simplified36.0%
sqrt-prod51.7%
*-commutative51.7%
*-commutative51.7%
associate-+l+51.7%
unpow251.7%
hypot-udef78.3%
associate-+r+78.3%
+-commutative78.3%
associate-+r+78.3%
Applied egg-rr78.3%
sqrt-prod78.2%
cancel-sign-sub-inv78.2%
metadata-eval78.2%
Applied egg-rr78.2%
*-un-lft-identity78.2%
distribute-rgt-neg-in78.2%
sqrt-unprod78.3%
fma-def78.3%
cancel-sign-sub-inv78.3%
metadata-eval78.3%
*-commutative78.3%
fma-def78.3%
Applied egg-rr78.3%
*-lft-identity78.3%
distribute-rgt-neg-out78.3%
distribute-neg-frac78.3%
associate-/l*78.7%
Simplified78.7%
if +inf.0 < (/.f64 (neg.f64 (sqrt.f64 (*.f64 (*.f64 2 (*.f64 (-.f64 (pow.f64 B 2) (*.f64 (*.f64 4 A) C)) F)) (+.f64 (+.f64 A C) (sqrt.f64 (+.f64 (pow.f64 (-.f64 A C) 2) (pow.f64 B 2))))))) (-.f64 (pow.f64 B 2) (*.f64 (*.f64 4 A) C))) Initial program 0.0%
Simplified0.8%
Taylor expanded in A around 0 1.8%
mul-1-neg1.8%
distribute-rgt-neg-in1.8%
*-commutative1.8%
unpow21.8%
unpow21.8%
hypot-def16.6%
Simplified16.6%
Final simplification46.7%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (fma B B (* -4.0 (* A C)))))
(if (<= B 8e+87)
(-
(/
(sqrt (* 2.0 (* F t_0)))
(/ t_0 (sqrt (+ C (+ A (hypot B (- A C))))))))
(if (or (<= B 5.5e+209) (not (<= B 1.4e+295)))
(* (/ (sqrt 2.0) B) (- (sqrt (* F (+ C (hypot B C))))))
(* (sqrt 2.0) (- (sqrt (/ F B))))))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = fma(B, B, (-4.0 * (A * C)));
double tmp;
if (B <= 8e+87) {
tmp = -(sqrt((2.0 * (F * t_0))) / (t_0 / sqrt((C + (A + hypot(B, (A - C)))))));
} else if ((B <= 5.5e+209) || !(B <= 1.4e+295)) {
tmp = (sqrt(2.0) / B) * -sqrt((F * (C + hypot(B, C))));
} else {
tmp = sqrt(2.0) * -sqrt((F / B));
}
return tmp;
}
B = abs(B) function code(A, B, C, F) t_0 = fma(B, B, Float64(-4.0 * Float64(A * C))) tmp = 0.0 if (B <= 8e+87) tmp = Float64(-Float64(sqrt(Float64(2.0 * Float64(F * t_0))) / Float64(t_0 / sqrt(Float64(C + Float64(A + hypot(B, Float64(A - C)))))))); elseif ((B <= 5.5e+209) || !(B <= 1.4e+295)) tmp = Float64(Float64(sqrt(2.0) / B) * Float64(-sqrt(Float64(F * Float64(C + hypot(B, C)))))); else tmp = Float64(sqrt(2.0) * Float64(-sqrt(Float64(F / B)))); end return tmp end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(B * B + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[B, 8e+87], (-N[(N[Sqrt[N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(t$95$0 / N[Sqrt[N[(C + N[(A + N[Sqrt[B ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), If[Or[LessEqual[B, 5.5e+209], N[Not[LessEqual[B, 1.4e+295]], $MachinePrecision]], N[(N[(N[Sqrt[2.0], $MachinePrecision] / B), $MachinePrecision] * (-N[Sqrt[N[(F * N[(C + N[Sqrt[B ^ 2 + C ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision])), $MachinePrecision], N[(N[Sqrt[2.0], $MachinePrecision] * (-N[Sqrt[N[(F / B), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(B, B, -4 \cdot \left(A \cdot C\right)\right)\\
\mathbf{if}\;B \leq 8 \cdot 10^{+87}:\\
\;\;\;\;-\frac{\sqrt{2 \cdot \left(F \cdot t_0\right)}}{\frac{t_0}{\sqrt{C + \left(A + \mathsf{hypot}\left(B, A - C\right)\right)}}}\\
\mathbf{elif}\;B \leq 5.5 \cdot 10^{+209} \lor \neg \left(B \leq 1.4 \cdot 10^{+295}\right):\\
\;\;\;\;\frac{\sqrt{2}}{B} \cdot \left(-\sqrt{F \cdot \left(C + \mathsf{hypot}\left(B, C\right)\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2} \cdot \left(-\sqrt{\frac{F}{B}}\right)\\
\end{array}
\end{array}
if B < 7.9999999999999997e87Initial program 24.8%
associate-*l*24.8%
unpow224.8%
+-commutative24.8%
unpow224.8%
associate-*l*24.8%
unpow224.8%
Simplified24.8%
sqrt-prod27.8%
*-commutative27.8%
*-commutative27.8%
associate-+l+28.4%
unpow228.4%
hypot-udef38.0%
associate-+r+37.3%
+-commutative37.3%
associate-+r+37.9%
Applied egg-rr37.9%
sqrt-prod37.7%
cancel-sign-sub-inv37.7%
metadata-eval37.7%
Applied egg-rr37.7%
*-un-lft-identity37.7%
distribute-rgt-neg-in37.7%
sqrt-unprod37.9%
fma-def37.9%
cancel-sign-sub-inv37.9%
metadata-eval37.9%
*-commutative37.9%
fma-def37.9%
Applied egg-rr37.9%
*-lft-identity37.9%
distribute-rgt-neg-out37.9%
distribute-neg-frac37.9%
associate-/l*37.9%
Simplified37.9%
if 7.9999999999999997e87 < B < 5.49999999999999967e209 or 1.4000000000000001e295 < B Initial program 0.8%
Simplified5.5%
Taylor expanded in A around 0 15.6%
mul-1-neg15.6%
distribute-rgt-neg-in15.6%
*-commutative15.6%
unpow215.6%
unpow215.6%
hypot-def69.0%
Simplified69.0%
if 5.49999999999999967e209 < B < 1.4000000000000001e295Initial program 0.0%
Simplified0.0%
Taylor expanded in B around inf 0.0%
Taylor expanded in C around 0 55.0%
mul-1-neg55.0%
Simplified55.0%
Final simplification42.0%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))))
(if (<= B 8.5e+84)
(/
(* (sqrt (+ C (+ A (hypot B (- A C))))) (- (sqrt (* 2.0 (* F t_0)))))
t_0)
(if (or (<= B 1.5e+209) (not (<= B 2.2e+295)))
(* (/ (sqrt 2.0) B) (- (sqrt (* F (+ C (hypot B C))))))
(* (sqrt 2.0) (- (sqrt (/ F B))))))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (B <= 8.5e+84) {
tmp = (sqrt((C + (A + hypot(B, (A - C))))) * -sqrt((2.0 * (F * t_0)))) / t_0;
} else if ((B <= 1.5e+209) || !(B <= 2.2e+295)) {
tmp = (sqrt(2.0) / B) * -sqrt((F * (C + hypot(B, C))));
} else {
tmp = sqrt(2.0) * -sqrt((F / B));
}
return tmp;
}
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (B <= 8.5e+84) {
tmp = (Math.sqrt((C + (A + Math.hypot(B, (A - C))))) * -Math.sqrt((2.0 * (F * t_0)))) / t_0;
} else if ((B <= 1.5e+209) || !(B <= 2.2e+295)) {
tmp = (Math.sqrt(2.0) / B) * -Math.sqrt((F * (C + Math.hypot(B, C))));
} else {
tmp = Math.sqrt(2.0) * -Math.sqrt((F / B));
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) tmp = 0 if B <= 8.5e+84: tmp = (math.sqrt((C + (A + math.hypot(B, (A - C))))) * -math.sqrt((2.0 * (F * t_0)))) / t_0 elif (B <= 1.5e+209) or not (B <= 2.2e+295): tmp = (math.sqrt(2.0) / B) * -math.sqrt((F * (C + math.hypot(B, C)))) else: tmp = math.sqrt(2.0) * -math.sqrt((F / B)) return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) tmp = 0.0 if (B <= 8.5e+84) tmp = Float64(Float64(sqrt(Float64(C + Float64(A + hypot(B, Float64(A - C))))) * Float64(-sqrt(Float64(2.0 * Float64(F * t_0))))) / t_0); elseif ((B <= 1.5e+209) || !(B <= 2.2e+295)) tmp = Float64(Float64(sqrt(2.0) / B) * Float64(-sqrt(Float64(F * Float64(C + hypot(B, C)))))); else tmp = Float64(sqrt(2.0) * Float64(-sqrt(Float64(F / B)))); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); tmp = 0.0; if (B <= 8.5e+84) tmp = (sqrt((C + (A + hypot(B, (A - C))))) * -sqrt((2.0 * (F * t_0)))) / t_0; elseif ((B <= 1.5e+209) || ~((B <= 2.2e+295))) tmp = (sqrt(2.0) / B) * -sqrt((F * (C + hypot(B, C)))); else tmp = sqrt(2.0) * -sqrt((F / B)); end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[B, 8.5e+84], N[(N[(N[Sqrt[N[(C + N[(A + N[Sqrt[B ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * (-N[Sqrt[N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision])), $MachinePrecision] / t$95$0), $MachinePrecision], If[Or[LessEqual[B, 1.5e+209], N[Not[LessEqual[B, 2.2e+295]], $MachinePrecision]], N[(N[(N[Sqrt[2.0], $MachinePrecision] / B), $MachinePrecision] * (-N[Sqrt[N[(F * N[(C + N[Sqrt[B ^ 2 + C ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision])), $MachinePrecision], N[(N[Sqrt[2.0], $MachinePrecision] * (-N[Sqrt[N[(F / B), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
\mathbf{if}\;B \leq 8.5 \cdot 10^{+84}:\\
\;\;\;\;\frac{\sqrt{C + \left(A + \mathsf{hypot}\left(B, A - C\right)\right)} \cdot \left(-\sqrt{2 \cdot \left(F \cdot t_0\right)}\right)}{t_0}\\
\mathbf{elif}\;B \leq 1.5 \cdot 10^{+209} \lor \neg \left(B \leq 2.2 \cdot 10^{+295}\right):\\
\;\;\;\;\frac{\sqrt{2}}{B} \cdot \left(-\sqrt{F \cdot \left(C + \mathsf{hypot}\left(B, C\right)\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2} \cdot \left(-\sqrt{\frac{F}{B}}\right)\\
\end{array}
\end{array}
if B < 8.5000000000000008e84Initial program 24.8%
associate-*l*24.8%
unpow224.8%
+-commutative24.8%
unpow224.8%
associate-*l*24.8%
unpow224.8%
Simplified24.8%
sqrt-prod27.8%
*-commutative27.8%
*-commutative27.8%
associate-+l+28.4%
unpow228.4%
hypot-udef38.0%
associate-+r+37.3%
+-commutative37.3%
associate-+r+37.9%
Applied egg-rr37.9%
if 8.5000000000000008e84 < B < 1.49999999999999993e209 or 2.1999999999999998e295 < B Initial program 0.8%
Simplified5.5%
Taylor expanded in A around 0 15.6%
mul-1-neg15.6%
distribute-rgt-neg-in15.6%
*-commutative15.6%
unpow215.6%
unpow215.6%
hypot-def69.0%
Simplified69.0%
if 1.49999999999999993e209 < B < 2.1999999999999998e295Initial program 0.0%
Simplified0.0%
Taylor expanded in B around inf 0.0%
Taylor expanded in C around 0 55.0%
mul-1-neg55.0%
Simplified55.0%
Final simplification42.0%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C))))
(t_1 (* F t_0))
(t_2 (+ C (+ A (hypot B (- A C))))))
(if (<= B 9.6e-155)
(/ (* (sqrt t_2) (- (sqrt (* 2.0 t_1)))) (* -4.0 (* A C)))
(if (<= B 9.2e+35)
(* (sqrt (* 2.0 (* t_2 t_1))) (/ -1.0 t_0))
(if (or (<= B 1.75e+210) (not (<= B 2.05e+295)))
(* (/ (sqrt 2.0) B) (- (sqrt (* F (+ C (hypot B C))))))
(* (sqrt 2.0) (- (sqrt (/ F B)))))))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = F * t_0;
double t_2 = C + (A + hypot(B, (A - C)));
double tmp;
if (B <= 9.6e-155) {
tmp = (sqrt(t_2) * -sqrt((2.0 * t_1))) / (-4.0 * (A * C));
} else if (B <= 9.2e+35) {
tmp = sqrt((2.0 * (t_2 * t_1))) * (-1.0 / t_0);
} else if ((B <= 1.75e+210) || !(B <= 2.05e+295)) {
tmp = (sqrt(2.0) / B) * -sqrt((F * (C + hypot(B, C))));
} else {
tmp = sqrt(2.0) * -sqrt((F / B));
}
return tmp;
}
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = F * t_0;
double t_2 = C + (A + Math.hypot(B, (A - C)));
double tmp;
if (B <= 9.6e-155) {
tmp = (Math.sqrt(t_2) * -Math.sqrt((2.0 * t_1))) / (-4.0 * (A * C));
} else if (B <= 9.2e+35) {
tmp = Math.sqrt((2.0 * (t_2 * t_1))) * (-1.0 / t_0);
} else if ((B <= 1.75e+210) || !(B <= 2.05e+295)) {
tmp = (Math.sqrt(2.0) / B) * -Math.sqrt((F * (C + Math.hypot(B, C))));
} else {
tmp = Math.sqrt(2.0) * -Math.sqrt((F / B));
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) t_1 = F * t_0 t_2 = C + (A + math.hypot(B, (A - C))) tmp = 0 if B <= 9.6e-155: tmp = (math.sqrt(t_2) * -math.sqrt((2.0 * t_1))) / (-4.0 * (A * C)) elif B <= 9.2e+35: tmp = math.sqrt((2.0 * (t_2 * t_1))) * (-1.0 / t_0) elif (B <= 1.75e+210) or not (B <= 2.05e+295): tmp = (math.sqrt(2.0) / B) * -math.sqrt((F * (C + math.hypot(B, C)))) else: tmp = math.sqrt(2.0) * -math.sqrt((F / B)) return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) t_1 = Float64(F * t_0) t_2 = Float64(C + Float64(A + hypot(B, Float64(A - C)))) tmp = 0.0 if (B <= 9.6e-155) tmp = Float64(Float64(sqrt(t_2) * Float64(-sqrt(Float64(2.0 * t_1)))) / Float64(-4.0 * Float64(A * C))); elseif (B <= 9.2e+35) tmp = Float64(sqrt(Float64(2.0 * Float64(t_2 * t_1))) * Float64(-1.0 / t_0)); elseif ((B <= 1.75e+210) || !(B <= 2.05e+295)) tmp = Float64(Float64(sqrt(2.0) / B) * Float64(-sqrt(Float64(F * Float64(C + hypot(B, C)))))); else tmp = Float64(sqrt(2.0) * Float64(-sqrt(Float64(F / B)))); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); t_1 = F * t_0; t_2 = C + (A + hypot(B, (A - C))); tmp = 0.0; if (B <= 9.6e-155) tmp = (sqrt(t_2) * -sqrt((2.0 * t_1))) / (-4.0 * (A * C)); elseif (B <= 9.2e+35) tmp = sqrt((2.0 * (t_2 * t_1))) * (-1.0 / t_0); elseif ((B <= 1.75e+210) || ~((B <= 2.05e+295))) tmp = (sqrt(2.0) / B) * -sqrt((F * (C + hypot(B, C)))); else tmp = sqrt(2.0) * -sqrt((F / B)); end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(F * t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[(C + N[(A + N[Sqrt[B ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[B, 9.6e-155], N[(N[(N[Sqrt[t$95$2], $MachinePrecision] * (-N[Sqrt[N[(2.0 * t$95$1), $MachinePrecision]], $MachinePrecision])), $MachinePrecision] / N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[B, 9.2e+35], N[(N[Sqrt[N[(2.0 * N[(t$95$2 * t$95$1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(-1.0 / t$95$0), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[B, 1.75e+210], N[Not[LessEqual[B, 2.05e+295]], $MachinePrecision]], N[(N[(N[Sqrt[2.0], $MachinePrecision] / B), $MachinePrecision] * (-N[Sqrt[N[(F * N[(C + N[Sqrt[B ^ 2 + C ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision])), $MachinePrecision], N[(N[Sqrt[2.0], $MachinePrecision] * (-N[Sqrt[N[(F / B), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]]]]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
t_1 := F \cdot t_0\\
t_2 := C + \left(A + \mathsf{hypot}\left(B, A - C\right)\right)\\
\mathbf{if}\;B \leq 9.6 \cdot 10^{-155}:\\
\;\;\;\;\frac{\sqrt{t_2} \cdot \left(-\sqrt{2 \cdot t_1}\right)}{-4 \cdot \left(A \cdot C\right)}\\
\mathbf{elif}\;B \leq 9.2 \cdot 10^{+35}:\\
\;\;\;\;\sqrt{2 \cdot \left(t_2 \cdot t_1\right)} \cdot \frac{-1}{t_0}\\
\mathbf{elif}\;B \leq 1.75 \cdot 10^{+210} \lor \neg \left(B \leq 2.05 \cdot 10^{+295}\right):\\
\;\;\;\;\frac{\sqrt{2}}{B} \cdot \left(-\sqrt{F \cdot \left(C + \mathsf{hypot}\left(B, C\right)\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2} \cdot \left(-\sqrt{\frac{F}{B}}\right)\\
\end{array}
\end{array}
if B < 9.600000000000001e-155Initial program 22.1%
associate-*l*22.1%
unpow222.1%
+-commutative22.1%
unpow222.1%
associate-*l*22.1%
unpow222.1%
Simplified22.1%
sqrt-prod25.6%
*-commutative25.6%
*-commutative25.6%
associate-+l+26.1%
unpow226.1%
hypot-udef35.9%
associate-+r+35.2%
+-commutative35.2%
associate-+r+35.8%
Applied egg-rr35.8%
Taylor expanded in B around 0 21.4%
if 9.600000000000001e-155 < B < 9.1999999999999993e35Initial program 35.5%
associate-*l*35.5%
unpow235.5%
+-commutative35.5%
unpow235.5%
associate-*l*35.5%
unpow235.5%
Simplified35.5%
div-inv35.4%
Applied egg-rr40.9%
if 9.1999999999999993e35 < B < 1.75e210 or 2.04999999999999992e295 < B Initial program 7.4%
Simplified10.7%
Taylor expanded in A around 0 20.9%
mul-1-neg20.9%
distribute-rgt-neg-in20.9%
*-commutative20.9%
unpow220.9%
unpow220.9%
hypot-def57.9%
Simplified57.9%
if 1.75e210 < B < 2.04999999999999992e295Initial program 0.0%
Simplified0.0%
Taylor expanded in B around inf 0.0%
Taylor expanded in C around 0 55.0%
mul-1-neg55.0%
Simplified55.0%
Final simplification32.1%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))))
(if (<= B 1e+36)
(/ (- (sqrt (* 2.0 (* (+ C (+ A (hypot B (- A C)))) (* F t_0))))) t_0)
(if (or (<= B 9e+209) (not (<= B 2.2e+295)))
(* (/ (sqrt 2.0) B) (- (sqrt (* F (+ C (hypot B C))))))
(* (sqrt 2.0) (- (sqrt (/ F B))))))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (B <= 1e+36) {
tmp = -sqrt((2.0 * ((C + (A + hypot(B, (A - C)))) * (F * t_0)))) / t_0;
} else if ((B <= 9e+209) || !(B <= 2.2e+295)) {
tmp = (sqrt(2.0) / B) * -sqrt((F * (C + hypot(B, C))));
} else {
tmp = sqrt(2.0) * -sqrt((F / B));
}
return tmp;
}
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (B <= 1e+36) {
tmp = -Math.sqrt((2.0 * ((C + (A + Math.hypot(B, (A - C)))) * (F * t_0)))) / t_0;
} else if ((B <= 9e+209) || !(B <= 2.2e+295)) {
tmp = (Math.sqrt(2.0) / B) * -Math.sqrt((F * (C + Math.hypot(B, C))));
} else {
tmp = Math.sqrt(2.0) * -Math.sqrt((F / B));
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) tmp = 0 if B <= 1e+36: tmp = -math.sqrt((2.0 * ((C + (A + math.hypot(B, (A - C)))) * (F * t_0)))) / t_0 elif (B <= 9e+209) or not (B <= 2.2e+295): tmp = (math.sqrt(2.0) / B) * -math.sqrt((F * (C + math.hypot(B, C)))) else: tmp = math.sqrt(2.0) * -math.sqrt((F / B)) return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) tmp = 0.0 if (B <= 1e+36) tmp = Float64(Float64(-sqrt(Float64(2.0 * Float64(Float64(C + Float64(A + hypot(B, Float64(A - C)))) * Float64(F * t_0))))) / t_0); elseif ((B <= 9e+209) || !(B <= 2.2e+295)) tmp = Float64(Float64(sqrt(2.0) / B) * Float64(-sqrt(Float64(F * Float64(C + hypot(B, C)))))); else tmp = Float64(sqrt(2.0) * Float64(-sqrt(Float64(F / B)))); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); tmp = 0.0; if (B <= 1e+36) tmp = -sqrt((2.0 * ((C + (A + hypot(B, (A - C)))) * (F * t_0)))) / t_0; elseif ((B <= 9e+209) || ~((B <= 2.2e+295))) tmp = (sqrt(2.0) / B) * -sqrt((F * (C + hypot(B, C)))); else tmp = sqrt(2.0) * -sqrt((F / B)); end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[B, 1e+36], N[((-N[Sqrt[N[(2.0 * N[(N[(C + N[(A + N[Sqrt[B ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], If[Or[LessEqual[B, 9e+209], N[Not[LessEqual[B, 2.2e+295]], $MachinePrecision]], N[(N[(N[Sqrt[2.0], $MachinePrecision] / B), $MachinePrecision] * (-N[Sqrt[N[(F * N[(C + N[Sqrt[B ^ 2 + C ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision])), $MachinePrecision], N[(N[Sqrt[2.0], $MachinePrecision] * (-N[Sqrt[N[(F / B), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
\mathbf{if}\;B \leq 10^{+36}:\\
\;\;\;\;\frac{-\sqrt{2 \cdot \left(\left(C + \left(A + \mathsf{hypot}\left(B, A - C\right)\right)\right) \cdot \left(F \cdot t_0\right)\right)}}{t_0}\\
\mathbf{elif}\;B \leq 9 \cdot 10^{+209} \lor \neg \left(B \leq 2.2 \cdot 10^{+295}\right):\\
\;\;\;\;\frac{\sqrt{2}}{B} \cdot \left(-\sqrt{F \cdot \left(C + \mathsf{hypot}\left(B, C\right)\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2} \cdot \left(-\sqrt{\frac{F}{B}}\right)\\
\end{array}
\end{array}
if B < 1.00000000000000004e36Initial program 25.0%
associate-*l*25.0%
unpow225.0%
+-commutative25.0%
unpow225.0%
associate-*l*25.0%
unpow225.0%
Simplified25.0%
distribute-frac-neg25.0%
Applied egg-rr30.7%
if 1.00000000000000004e36 < B < 9.0000000000000007e209 or 2.1999999999999998e295 < B Initial program 7.4%
Simplified10.7%
Taylor expanded in A around 0 20.9%
mul-1-neg20.9%
distribute-rgt-neg-in20.9%
*-commutative20.9%
unpow220.9%
unpow220.9%
hypot-def57.9%
Simplified57.9%
if 9.0000000000000007e209 < B < 2.1999999999999998e295Initial program 0.0%
Simplified0.0%
Taylor expanded in B around inf 0.0%
Taylor expanded in C around 0 55.0%
mul-1-neg55.0%
Simplified55.0%
Final simplification36.1%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C))))
(t_1
(*
(sqrt (* 2.0 (* (+ C (+ A (hypot B (- A C)))) (* F t_0))))
(/ -1.0 t_0))))
(if (<= B 2.4e-174)
t_1
(if (<= B 5.3e-158)
(* (sqrt 2.0) (* (sqrt 0.5) (- (sqrt (/ F (- A C))))))
(if (<= B 1e+89) t_1 (* (sqrt 2.0) (- (sqrt (/ F B)))))))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = sqrt((2.0 * ((C + (A + hypot(B, (A - C)))) * (F * t_0)))) * (-1.0 / t_0);
double tmp;
if (B <= 2.4e-174) {
tmp = t_1;
} else if (B <= 5.3e-158) {
tmp = sqrt(2.0) * (sqrt(0.5) * -sqrt((F / (A - C))));
} else if (B <= 1e+89) {
tmp = t_1;
} else {
tmp = sqrt(2.0) * -sqrt((F / B));
}
return tmp;
}
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = Math.sqrt((2.0 * ((C + (A + Math.hypot(B, (A - C)))) * (F * t_0)))) * (-1.0 / t_0);
double tmp;
if (B <= 2.4e-174) {
tmp = t_1;
} else if (B <= 5.3e-158) {
tmp = Math.sqrt(2.0) * (Math.sqrt(0.5) * -Math.sqrt((F / (A - C))));
} else if (B <= 1e+89) {
tmp = t_1;
} else {
tmp = Math.sqrt(2.0) * -Math.sqrt((F / B));
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) t_1 = math.sqrt((2.0 * ((C + (A + math.hypot(B, (A - C)))) * (F * t_0)))) * (-1.0 / t_0) tmp = 0 if B <= 2.4e-174: tmp = t_1 elif B <= 5.3e-158: tmp = math.sqrt(2.0) * (math.sqrt(0.5) * -math.sqrt((F / (A - C)))) elif B <= 1e+89: tmp = t_1 else: tmp = math.sqrt(2.0) * -math.sqrt((F / B)) return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) t_1 = Float64(sqrt(Float64(2.0 * Float64(Float64(C + Float64(A + hypot(B, Float64(A - C)))) * Float64(F * t_0)))) * Float64(-1.0 / t_0)) tmp = 0.0 if (B <= 2.4e-174) tmp = t_1; elseif (B <= 5.3e-158) tmp = Float64(sqrt(2.0) * Float64(sqrt(0.5) * Float64(-sqrt(Float64(F / Float64(A - C)))))); elseif (B <= 1e+89) tmp = t_1; else tmp = Float64(sqrt(2.0) * Float64(-sqrt(Float64(F / B)))); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); t_1 = sqrt((2.0 * ((C + (A + hypot(B, (A - C)))) * (F * t_0)))) * (-1.0 / t_0); tmp = 0.0; if (B <= 2.4e-174) tmp = t_1; elseif (B <= 5.3e-158) tmp = sqrt(2.0) * (sqrt(0.5) * -sqrt((F / (A - C)))); elseif (B <= 1e+89) tmp = t_1; else tmp = sqrt(2.0) * -sqrt((F / B)); end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Sqrt[N[(2.0 * N[(N[(C + N[(A + N[Sqrt[B ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(-1.0 / t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[B, 2.4e-174], t$95$1, If[LessEqual[B, 5.3e-158], N[(N[Sqrt[2.0], $MachinePrecision] * N[(N[Sqrt[0.5], $MachinePrecision] * (-N[Sqrt[N[(F / N[(A - C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]), $MachinePrecision], If[LessEqual[B, 1e+89], t$95$1, N[(N[Sqrt[2.0], $MachinePrecision] * (-N[Sqrt[N[(F / B), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]]]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
t_1 := \sqrt{2 \cdot \left(\left(C + \left(A + \mathsf{hypot}\left(B, A - C\right)\right)\right) \cdot \left(F \cdot t_0\right)\right)} \cdot \frac{-1}{t_0}\\
\mathbf{if}\;B \leq 2.4 \cdot 10^{-174}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;B \leq 5.3 \cdot 10^{-158}:\\
\;\;\;\;\sqrt{2} \cdot \left(\sqrt{0.5} \cdot \left(-\sqrt{\frac{F}{A - C}}\right)\right)\\
\mathbf{elif}\;B \leq 10^{+89}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2} \cdot \left(-\sqrt{\frac{F}{B}}\right)\\
\end{array}
\end{array}
if B < 2.4e-174 or 5.2999999999999999e-158 < B < 9.99999999999999995e88Initial program 25.1%
associate-*l*25.1%
unpow225.1%
+-commutative25.1%
unpow225.1%
associate-*l*25.1%
unpow225.1%
Simplified25.1%
div-inv25.1%
Applied egg-rr30.7%
if 2.4e-174 < B < 5.2999999999999999e-158Initial program 4.6%
associate-*l*4.6%
unpow24.6%
+-commutative4.6%
unpow24.6%
associate-*l*4.6%
unpow24.6%
Simplified4.6%
Taylor expanded in B around 0 2.1%
associate--l+2.1%
associate-*r/2.1%
*-commutative2.1%
unpow22.1%
Simplified2.1%
Taylor expanded in B around inf 65.2%
mul-1-neg65.2%
associate-*l*65.5%
Simplified65.5%
if 9.99999999999999995e88 < B Initial program 0.4%
Simplified2.8%
Taylor expanded in B around inf 0.6%
Taylor expanded in C around 0 43.5%
mul-1-neg43.5%
Simplified43.5%
Final simplification33.2%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))))
(if (<= B 1.2e+89)
(/ (- (sqrt (* 2.0 (* (+ C (+ A (hypot B (- A C)))) (* F t_0))))) t_0)
(* (sqrt 2.0) (- (sqrt (/ F B)))))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (B <= 1.2e+89) {
tmp = -sqrt((2.0 * ((C + (A + hypot(B, (A - C)))) * (F * t_0)))) / t_0;
} else {
tmp = sqrt(2.0) * -sqrt((F / B));
}
return tmp;
}
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (B <= 1.2e+89) {
tmp = -Math.sqrt((2.0 * ((C + (A + Math.hypot(B, (A - C)))) * (F * t_0)))) / t_0;
} else {
tmp = Math.sqrt(2.0) * -Math.sqrt((F / B));
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) tmp = 0 if B <= 1.2e+89: tmp = -math.sqrt((2.0 * ((C + (A + math.hypot(B, (A - C)))) * (F * t_0)))) / t_0 else: tmp = math.sqrt(2.0) * -math.sqrt((F / B)) return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) tmp = 0.0 if (B <= 1.2e+89) tmp = Float64(Float64(-sqrt(Float64(2.0 * Float64(Float64(C + Float64(A + hypot(B, Float64(A - C)))) * Float64(F * t_0))))) / t_0); else tmp = Float64(sqrt(2.0) * Float64(-sqrt(Float64(F / B)))); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); tmp = 0.0; if (B <= 1.2e+89) tmp = -sqrt((2.0 * ((C + (A + hypot(B, (A - C)))) * (F * t_0)))) / t_0; else tmp = sqrt(2.0) * -sqrt((F / B)); end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[B, 1.2e+89], N[((-N[Sqrt[N[(2.0 * N[(N[(C + N[(A + N[Sqrt[B ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], N[(N[Sqrt[2.0], $MachinePrecision] * (-N[Sqrt[N[(F / B), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
\mathbf{if}\;B \leq 1.2 \cdot 10^{+89}:\\
\;\;\;\;\frac{-\sqrt{2 \cdot \left(\left(C + \left(A + \mathsf{hypot}\left(B, A - C\right)\right)\right) \cdot \left(F \cdot t_0\right)\right)}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2} \cdot \left(-\sqrt{\frac{F}{B}}\right)\\
\end{array}
\end{array}
if B < 1.20000000000000002e89Initial program 24.8%
associate-*l*24.8%
unpow224.8%
+-commutative24.8%
unpow224.8%
associate-*l*24.8%
unpow224.8%
Simplified24.8%
distribute-frac-neg24.8%
Applied egg-rr30.3%
if 1.20000000000000002e89 < B Initial program 0.4%
Simplified2.8%
Taylor expanded in B around inf 0.6%
Taylor expanded in C around 0 43.5%
mul-1-neg43.5%
Simplified43.5%
Final simplification32.5%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))))
(if (<= B 2.5e+35)
(/ (- (sqrt (* (* 2.0 (* F t_0)) (+ (+ A C) (hypot B A))))) t_0)
(* (sqrt 2.0) (- (sqrt (/ F B)))))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (B <= 2.5e+35) {
tmp = -sqrt(((2.0 * (F * t_0)) * ((A + C) + hypot(B, A)))) / t_0;
} else {
tmp = sqrt(2.0) * -sqrt((F / B));
}
return tmp;
}
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (B <= 2.5e+35) {
tmp = -Math.sqrt(((2.0 * (F * t_0)) * ((A + C) + Math.hypot(B, A)))) / t_0;
} else {
tmp = Math.sqrt(2.0) * -Math.sqrt((F / B));
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) tmp = 0 if B <= 2.5e+35: tmp = -math.sqrt(((2.0 * (F * t_0)) * ((A + C) + math.hypot(B, A)))) / t_0 else: tmp = math.sqrt(2.0) * -math.sqrt((F / B)) return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) tmp = 0.0 if (B <= 2.5e+35) tmp = Float64(Float64(-sqrt(Float64(Float64(2.0 * Float64(F * t_0)) * Float64(Float64(A + C) + hypot(B, A))))) / t_0); else tmp = Float64(sqrt(2.0) * Float64(-sqrt(Float64(F / B)))); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); tmp = 0.0; if (B <= 2.5e+35) tmp = -sqrt(((2.0 * (F * t_0)) * ((A + C) + hypot(B, A)))) / t_0; else tmp = sqrt(2.0) * -sqrt((F / B)); end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[B, 2.5e+35], N[((-N[Sqrt[N[(N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision] * N[(N[(A + C), $MachinePrecision] + N[Sqrt[B ^ 2 + A ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], N[(N[Sqrt[2.0], $MachinePrecision] * (-N[Sqrt[N[(F / B), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
\mathbf{if}\;B \leq 2.5 \cdot 10^{+35}:\\
\;\;\;\;\frac{-\sqrt{\left(2 \cdot \left(F \cdot t_0\right)\right) \cdot \left(\left(A + C\right) + \mathsf{hypot}\left(B, A\right)\right)}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2} \cdot \left(-\sqrt{\frac{F}{B}}\right)\\
\end{array}
\end{array}
if B < 2.50000000000000011e35Initial program 25.0%
associate-*l*25.0%
unpow225.0%
+-commutative25.0%
unpow225.0%
associate-*l*25.0%
unpow225.0%
Simplified25.0%
Taylor expanded in C around 0 20.4%
unpow220.4%
unpow220.4%
hypot-def24.3%
Simplified24.3%
if 2.50000000000000011e35 < B Initial program 4.4%
Simplified6.4%
Taylor expanded in B around inf 2.9%
Taylor expanded in C around 0 40.2%
mul-1-neg40.2%
Simplified40.2%
Final simplification27.6%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(if (<= A 5.5e+34)
(* (sqrt 2.0) (- (sqrt (/ F B))))
(/
(*
(sqrt (* 2.0 (* F (+ (* B B) (* -4.0 (* A C))))))
(- (sqrt (+ A (+ A C)))))
(- (* B B) (* 4.0 (* A C))))))B = abs(B);
double code(double A, double B, double C, double F) {
double tmp;
if (A <= 5.5e+34) {
tmp = sqrt(2.0) * -sqrt((F / B));
} else {
tmp = (sqrt((2.0 * (F * ((B * B) + (-4.0 * (A * C)))))) * -sqrt((A + (A + C)))) / ((B * B) - (4.0 * (A * C)));
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (a <= 5.5d+34) then
tmp = sqrt(2.0d0) * -sqrt((f / b))
else
tmp = (sqrt((2.0d0 * (f * ((b * b) + ((-4.0d0) * (a * c)))))) * -sqrt((a + (a + c)))) / ((b * b) - (4.0d0 * (a * c)))
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double tmp;
if (A <= 5.5e+34) {
tmp = Math.sqrt(2.0) * -Math.sqrt((F / B));
} else {
tmp = (Math.sqrt((2.0 * (F * ((B * B) + (-4.0 * (A * C)))))) * -Math.sqrt((A + (A + C)))) / ((B * B) - (4.0 * (A * C)));
}
return tmp;
}
B = abs(B) def code(A, B, C, F): tmp = 0 if A <= 5.5e+34: tmp = math.sqrt(2.0) * -math.sqrt((F / B)) else: tmp = (math.sqrt((2.0 * (F * ((B * B) + (-4.0 * (A * C)))))) * -math.sqrt((A + (A + C)))) / ((B * B) - (4.0 * (A * C))) return tmp
B = abs(B) function code(A, B, C, F) tmp = 0.0 if (A <= 5.5e+34) tmp = Float64(sqrt(2.0) * Float64(-sqrt(Float64(F / B)))); else tmp = Float64(Float64(sqrt(Float64(2.0 * Float64(F * Float64(Float64(B * B) + Float64(-4.0 * Float64(A * C)))))) * Float64(-sqrt(Float64(A + Float64(A + C))))) / Float64(Float64(B * B) - Float64(4.0 * Float64(A * C)))); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) tmp = 0.0; if (A <= 5.5e+34) tmp = sqrt(2.0) * -sqrt((F / B)); else tmp = (sqrt((2.0 * (F * ((B * B) + (-4.0 * (A * C)))))) * -sqrt((A + (A + C)))) / ((B * B) - (4.0 * (A * C))); end tmp_2 = tmp; end
NOTE: B should be positive before calling this function code[A_, B_, C_, F_] := If[LessEqual[A, 5.5e+34], N[(N[Sqrt[2.0], $MachinePrecision] * (-N[Sqrt[N[(F / B), $MachinePrecision]], $MachinePrecision])), $MachinePrecision], N[(N[(N[Sqrt[N[(2.0 * N[(F * N[(N[(B * B), $MachinePrecision] + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * (-N[Sqrt[N[(A + N[(A + C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision])), $MachinePrecision] / N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
\mathbf{if}\;A \leq 5.5 \cdot 10^{+34}:\\
\;\;\;\;\sqrt{2} \cdot \left(-\sqrt{\frac{F}{B}}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2 \cdot \left(F \cdot \left(B \cdot B + -4 \cdot \left(A \cdot C\right)\right)\right)} \cdot \left(-\sqrt{A + \left(A + C\right)}\right)}{B \cdot B - 4 \cdot \left(A \cdot C\right)}\\
\end{array}
\end{array}
if A < 5.4999999999999996e34Initial program 20.4%
Simplified23.1%
Taylor expanded in B around inf 5.3%
Taylor expanded in C around 0 15.0%
mul-1-neg15.0%
Simplified15.0%
if 5.4999999999999996e34 < A Initial program 21.5%
associate-*l*21.5%
unpow221.5%
+-commutative21.5%
unpow221.5%
associate-*l*21.5%
unpow221.5%
Simplified21.5%
Taylor expanded in A around inf 30.5%
sqrt-prod36.3%
*-commutative36.3%
cancel-sign-sub-inv36.3%
*-commutative36.3%
metadata-eval36.3%
associate-+l+36.3%
Applied egg-rr36.3%
Final simplification20.7%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))))
(if (<= A 7.5e+34)
(* (sqrt 2.0) (- (sqrt (/ F B))))
(/
(-
(sqrt
(* (* 2.0 (* F t_0)) (+ (+ A C) (+ A (- (* 0.5 (/ (* B B) A)) C))))))
t_0))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (A <= 7.5e+34) {
tmp = sqrt(2.0) * -sqrt((F / B));
} else {
tmp = -sqrt(((2.0 * (F * t_0)) * ((A + C) + (A + ((0.5 * ((B * B) / A)) - C))))) / t_0;
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: tmp
t_0 = (b * b) - (4.0d0 * (a * c))
if (a <= 7.5d+34) then
tmp = sqrt(2.0d0) * -sqrt((f / b))
else
tmp = -sqrt(((2.0d0 * (f * t_0)) * ((a + c) + (a + ((0.5d0 * ((b * b) / a)) - c))))) / t_0
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (A <= 7.5e+34) {
tmp = Math.sqrt(2.0) * -Math.sqrt((F / B));
} else {
tmp = -Math.sqrt(((2.0 * (F * t_0)) * ((A + C) + (A + ((0.5 * ((B * B) / A)) - C))))) / t_0;
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) tmp = 0 if A <= 7.5e+34: tmp = math.sqrt(2.0) * -math.sqrt((F / B)) else: tmp = -math.sqrt(((2.0 * (F * t_0)) * ((A + C) + (A + ((0.5 * ((B * B) / A)) - C))))) / t_0 return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) tmp = 0.0 if (A <= 7.5e+34) tmp = Float64(sqrt(2.0) * Float64(-sqrt(Float64(F / B)))); else tmp = Float64(Float64(-sqrt(Float64(Float64(2.0 * Float64(F * t_0)) * Float64(Float64(A + C) + Float64(A + Float64(Float64(0.5 * Float64(Float64(B * B) / A)) - C)))))) / t_0); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); tmp = 0.0; if (A <= 7.5e+34) tmp = sqrt(2.0) * -sqrt((F / B)); else tmp = -sqrt(((2.0 * (F * t_0)) * ((A + C) + (A + ((0.5 * ((B * B) / A)) - C))))) / t_0; end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[A, 7.5e+34], N[(N[Sqrt[2.0], $MachinePrecision] * (-N[Sqrt[N[(F / B), $MachinePrecision]], $MachinePrecision])), $MachinePrecision], N[((-N[Sqrt[N[(N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision] * N[(N[(A + C), $MachinePrecision] + N[(A + N[(N[(0.5 * N[(N[(B * B), $MachinePrecision] / A), $MachinePrecision]), $MachinePrecision] - C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
\mathbf{if}\;A \leq 7.5 \cdot 10^{+34}:\\
\;\;\;\;\sqrt{2} \cdot \left(-\sqrt{\frac{F}{B}}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{-\sqrt{\left(2 \cdot \left(F \cdot t_0\right)\right) \cdot \left(\left(A + C\right) + \left(A + \left(0.5 \cdot \frac{B \cdot B}{A} - C\right)\right)\right)}}{t_0}\\
\end{array}
\end{array}
if A < 7.49999999999999976e34Initial program 20.4%
Simplified23.1%
Taylor expanded in B around inf 5.3%
Taylor expanded in C around 0 15.0%
mul-1-neg15.0%
Simplified15.0%
if 7.49999999999999976e34 < A Initial program 21.5%
associate-*l*21.5%
unpow221.5%
+-commutative21.5%
unpow221.5%
associate-*l*21.5%
unpow221.5%
Simplified21.5%
Taylor expanded in B around 0 30.8%
associate--l+30.8%
associate-*r/30.8%
*-commutative30.8%
unpow230.8%
Simplified30.8%
Taylor expanded in A around inf 30.8%
unpow230.8%
Simplified30.8%
Final simplification19.2%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))) (t_1 (* 2.0 (* F t_0))))
(if (<= A -11000.0)
(/ (- (sqrt (* t_1 (* 2.0 C)))) t_0)
(if (<= A 1.7e-61)
(/ (- (sqrt (* t_1 (+ B (+ A C))))) t_0)
(/
(- (sqrt (* t_1 (+ (+ A C) (+ A (- (* 0.5 (/ (* B B) A)) C))))))
t_0)))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = 2.0 * (F * t_0);
double tmp;
if (A <= -11000.0) {
tmp = -sqrt((t_1 * (2.0 * C))) / t_0;
} else if (A <= 1.7e-61) {
tmp = -sqrt((t_1 * (B + (A + C)))) / t_0;
} else {
tmp = -sqrt((t_1 * ((A + C) + (A + ((0.5 * ((B * B) / A)) - C))))) / t_0;
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (b * b) - (4.0d0 * (a * c))
t_1 = 2.0d0 * (f * t_0)
if (a <= (-11000.0d0)) then
tmp = -sqrt((t_1 * (2.0d0 * c))) / t_0
else if (a <= 1.7d-61) then
tmp = -sqrt((t_1 * (b + (a + c)))) / t_0
else
tmp = -sqrt((t_1 * ((a + c) + (a + ((0.5d0 * ((b * b) / a)) - c))))) / t_0
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = 2.0 * (F * t_0);
double tmp;
if (A <= -11000.0) {
tmp = -Math.sqrt((t_1 * (2.0 * C))) / t_0;
} else if (A <= 1.7e-61) {
tmp = -Math.sqrt((t_1 * (B + (A + C)))) / t_0;
} else {
tmp = -Math.sqrt((t_1 * ((A + C) + (A + ((0.5 * ((B * B) / A)) - C))))) / t_0;
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) t_1 = 2.0 * (F * t_0) tmp = 0 if A <= -11000.0: tmp = -math.sqrt((t_1 * (2.0 * C))) / t_0 elif A <= 1.7e-61: tmp = -math.sqrt((t_1 * (B + (A + C)))) / t_0 else: tmp = -math.sqrt((t_1 * ((A + C) + (A + ((0.5 * ((B * B) / A)) - C))))) / t_0 return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) t_1 = Float64(2.0 * Float64(F * t_0)) tmp = 0.0 if (A <= -11000.0) tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(2.0 * C)))) / t_0); elseif (A <= 1.7e-61) tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(B + Float64(A + C))))) / t_0); else tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(Float64(A + C) + Float64(A + Float64(Float64(0.5 * Float64(Float64(B * B) / A)) - C)))))) / t_0); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); t_1 = 2.0 * (F * t_0); tmp = 0.0; if (A <= -11000.0) tmp = -sqrt((t_1 * (2.0 * C))) / t_0; elseif (A <= 1.7e-61) tmp = -sqrt((t_1 * (B + (A + C)))) / t_0; else tmp = -sqrt((t_1 * ((A + C) + (A + ((0.5 * ((B * B) / A)) - C))))) / t_0; end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[A, -11000.0], N[((-N[Sqrt[N[(t$95$1 * N[(2.0 * C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], If[LessEqual[A, 1.7e-61], N[((-N[Sqrt[N[(t$95$1 * N[(B + N[(A + C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], N[((-N[Sqrt[N[(t$95$1 * N[(N[(A + C), $MachinePrecision] + N[(A + N[(N[(0.5 * N[(N[(B * B), $MachinePrecision] / A), $MachinePrecision]), $MachinePrecision] - C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
t_1 := 2 \cdot \left(F \cdot t_0\right)\\
\mathbf{if}\;A \leq -11000:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(2 \cdot C\right)}}{t_0}\\
\mathbf{elif}\;A \leq 1.7 \cdot 10^{-61}:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(B + \left(A + C\right)\right)}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(\left(A + C\right) + \left(A + \left(0.5 \cdot \frac{B \cdot B}{A} - C\right)\right)\right)}}{t_0}\\
\end{array}
\end{array}
if A < -11000Initial program 5.6%
associate-*l*5.6%
unpow25.6%
+-commutative5.6%
unpow25.6%
associate-*l*5.6%
unpow25.6%
Simplified5.6%
Taylor expanded in A around -inf 4.9%
associate-+r+4.9%
mul-1-neg4.9%
unsub-neg4.9%
distribute-lft-out4.9%
associate-/l*4.9%
unpow24.9%
unpow24.9%
times-frac4.9%
unpow24.9%
Simplified4.9%
Taylor expanded in A around inf 17.0%
if -11000 < A < 1.6999999999999999e-61Initial program 28.7%
associate-*l*28.7%
unpow228.7%
+-commutative28.7%
unpow228.7%
associate-*l*28.7%
unpow228.7%
Simplified28.7%
unpow228.7%
hypot-udef30.8%
add-exp-log29.2%
Applied egg-rr29.2%
Taylor expanded in B around inf 9.1%
if 1.6999999999999999e-61 < A Initial program 22.7%
associate-*l*22.7%
unpow222.7%
+-commutative22.7%
unpow222.7%
associate-*l*22.7%
unpow222.7%
Simplified22.7%
Taylor expanded in B around 0 26.7%
associate--l+26.7%
associate-*r/26.7%
*-commutative26.7%
unpow226.7%
Simplified26.7%
Taylor expanded in A around inf 26.8%
unpow226.8%
Simplified26.8%
Final simplification17.2%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))) (t_1 (* 2.0 (* F t_0))))
(if (<= A -1900.0)
(/ (- (sqrt (* t_1 (* 2.0 C)))) t_0)
(if (<= A 7.5e-62)
(/ (- (sqrt (* t_1 (+ B (+ A C))))) t_0)
(/ (- (sqrt (* t_1 (+ (+ A C) (- A C))))) t_0)))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = 2.0 * (F * t_0);
double tmp;
if (A <= -1900.0) {
tmp = -sqrt((t_1 * (2.0 * C))) / t_0;
} else if (A <= 7.5e-62) {
tmp = -sqrt((t_1 * (B + (A + C)))) / t_0;
} else {
tmp = -sqrt((t_1 * ((A + C) + (A - C)))) / t_0;
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (b * b) - (4.0d0 * (a * c))
t_1 = 2.0d0 * (f * t_0)
if (a <= (-1900.0d0)) then
tmp = -sqrt((t_1 * (2.0d0 * c))) / t_0
else if (a <= 7.5d-62) then
tmp = -sqrt((t_1 * (b + (a + c)))) / t_0
else
tmp = -sqrt((t_1 * ((a + c) + (a - c)))) / t_0
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = 2.0 * (F * t_0);
double tmp;
if (A <= -1900.0) {
tmp = -Math.sqrt((t_1 * (2.0 * C))) / t_0;
} else if (A <= 7.5e-62) {
tmp = -Math.sqrt((t_1 * (B + (A + C)))) / t_0;
} else {
tmp = -Math.sqrt((t_1 * ((A + C) + (A - C)))) / t_0;
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) t_1 = 2.0 * (F * t_0) tmp = 0 if A <= -1900.0: tmp = -math.sqrt((t_1 * (2.0 * C))) / t_0 elif A <= 7.5e-62: tmp = -math.sqrt((t_1 * (B + (A + C)))) / t_0 else: tmp = -math.sqrt((t_1 * ((A + C) + (A - C)))) / t_0 return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) t_1 = Float64(2.0 * Float64(F * t_0)) tmp = 0.0 if (A <= -1900.0) tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(2.0 * C)))) / t_0); elseif (A <= 7.5e-62) tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(B + Float64(A + C))))) / t_0); else tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(Float64(A + C) + Float64(A - C))))) / t_0); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); t_1 = 2.0 * (F * t_0); tmp = 0.0; if (A <= -1900.0) tmp = -sqrt((t_1 * (2.0 * C))) / t_0; elseif (A <= 7.5e-62) tmp = -sqrt((t_1 * (B + (A + C)))) / t_0; else tmp = -sqrt((t_1 * ((A + C) + (A - C)))) / t_0; end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[A, -1900.0], N[((-N[Sqrt[N[(t$95$1 * N[(2.0 * C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], If[LessEqual[A, 7.5e-62], N[((-N[Sqrt[N[(t$95$1 * N[(B + N[(A + C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], N[((-N[Sqrt[N[(t$95$1 * N[(N[(A + C), $MachinePrecision] + N[(A - C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
t_1 := 2 \cdot \left(F \cdot t_0\right)\\
\mathbf{if}\;A \leq -1900:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(2 \cdot C\right)}}{t_0}\\
\mathbf{elif}\;A \leq 7.5 \cdot 10^{-62}:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(B + \left(A + C\right)\right)}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(\left(A + C\right) + \left(A - C\right)\right)}}{t_0}\\
\end{array}
\end{array}
if A < -1900Initial program 5.6%
associate-*l*5.6%
unpow25.6%
+-commutative5.6%
unpow25.6%
associate-*l*5.6%
unpow25.6%
Simplified5.6%
Taylor expanded in A around -inf 4.9%
associate-+r+4.9%
mul-1-neg4.9%
unsub-neg4.9%
distribute-lft-out4.9%
associate-/l*4.9%
unpow24.9%
unpow24.9%
times-frac4.9%
unpow24.9%
Simplified4.9%
Taylor expanded in A around inf 17.0%
if -1900 < A < 7.5000000000000003e-62Initial program 28.7%
associate-*l*28.7%
unpow228.7%
+-commutative28.7%
unpow228.7%
associate-*l*28.7%
unpow228.7%
Simplified28.7%
unpow228.7%
hypot-udef30.8%
add-exp-log29.2%
Applied egg-rr29.2%
Taylor expanded in B around inf 9.1%
if 7.5000000000000003e-62 < A Initial program 22.7%
associate-*l*22.7%
unpow222.7%
+-commutative22.7%
unpow222.7%
associate-*l*22.7%
unpow222.7%
Simplified22.7%
Taylor expanded in B around 0 26.5%
Final simplification17.1%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))) (t_1 (* 2.0 (* F t_0))))
(if (<= A -7000.0)
(/ (- (sqrt (* t_1 (* 2.0 C)))) t_0)
(if (<= A 8e-62)
(/ (- (sqrt (* t_1 (+ B (+ A C))))) t_0)
(/ (- (sqrt (* t_1 (- (+ C (* 2.0 A)) C)))) t_0)))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = 2.0 * (F * t_0);
double tmp;
if (A <= -7000.0) {
tmp = -sqrt((t_1 * (2.0 * C))) / t_0;
} else if (A <= 8e-62) {
tmp = -sqrt((t_1 * (B + (A + C)))) / t_0;
} else {
tmp = -sqrt((t_1 * ((C + (2.0 * A)) - C))) / t_0;
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (b * b) - (4.0d0 * (a * c))
t_1 = 2.0d0 * (f * t_0)
if (a <= (-7000.0d0)) then
tmp = -sqrt((t_1 * (2.0d0 * c))) / t_0
else if (a <= 8d-62) then
tmp = -sqrt((t_1 * (b + (a + c)))) / t_0
else
tmp = -sqrt((t_1 * ((c + (2.0d0 * a)) - c))) / t_0
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = 2.0 * (F * t_0);
double tmp;
if (A <= -7000.0) {
tmp = -Math.sqrt((t_1 * (2.0 * C))) / t_0;
} else if (A <= 8e-62) {
tmp = -Math.sqrt((t_1 * (B + (A + C)))) / t_0;
} else {
tmp = -Math.sqrt((t_1 * ((C + (2.0 * A)) - C))) / t_0;
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) t_1 = 2.0 * (F * t_0) tmp = 0 if A <= -7000.0: tmp = -math.sqrt((t_1 * (2.0 * C))) / t_0 elif A <= 8e-62: tmp = -math.sqrt((t_1 * (B + (A + C)))) / t_0 else: tmp = -math.sqrt((t_1 * ((C + (2.0 * A)) - C))) / t_0 return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) t_1 = Float64(2.0 * Float64(F * t_0)) tmp = 0.0 if (A <= -7000.0) tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(2.0 * C)))) / t_0); elseif (A <= 8e-62) tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(B + Float64(A + C))))) / t_0); else tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(Float64(C + Float64(2.0 * A)) - C)))) / t_0); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); t_1 = 2.0 * (F * t_0); tmp = 0.0; if (A <= -7000.0) tmp = -sqrt((t_1 * (2.0 * C))) / t_0; elseif (A <= 8e-62) tmp = -sqrt((t_1 * (B + (A + C)))) / t_0; else tmp = -sqrt((t_1 * ((C + (2.0 * A)) - C))) / t_0; end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[A, -7000.0], N[((-N[Sqrt[N[(t$95$1 * N[(2.0 * C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], If[LessEqual[A, 8e-62], N[((-N[Sqrt[N[(t$95$1 * N[(B + N[(A + C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], N[((-N[Sqrt[N[(t$95$1 * N[(N[(C + N[(2.0 * A), $MachinePrecision]), $MachinePrecision] - C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
t_1 := 2 \cdot \left(F \cdot t_0\right)\\
\mathbf{if}\;A \leq -7000:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(2 \cdot C\right)}}{t_0}\\
\mathbf{elif}\;A \leq 8 \cdot 10^{-62}:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(B + \left(A + C\right)\right)}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(\left(C + 2 \cdot A\right) - C\right)}}{t_0}\\
\end{array}
\end{array}
if A < -7e3Initial program 5.6%
associate-*l*5.6%
unpow25.6%
+-commutative5.6%
unpow25.6%
associate-*l*5.6%
unpow25.6%
Simplified5.6%
Taylor expanded in A around -inf 4.9%
associate-+r+4.9%
mul-1-neg4.9%
unsub-neg4.9%
distribute-lft-out4.9%
associate-/l*4.9%
unpow24.9%
unpow24.9%
times-frac4.9%
unpow24.9%
Simplified4.9%
Taylor expanded in A around inf 17.0%
if -7e3 < A < 8.0000000000000003e-62Initial program 28.7%
associate-*l*28.7%
unpow228.7%
+-commutative28.7%
unpow228.7%
associate-*l*28.7%
unpow228.7%
Simplified28.7%
unpow228.7%
hypot-udef30.8%
add-exp-log29.2%
Applied egg-rr29.2%
Taylor expanded in B around inf 9.1%
if 8.0000000000000003e-62 < A Initial program 22.7%
associate-*l*22.7%
unpow222.7%
+-commutative22.7%
unpow222.7%
associate-*l*22.7%
unpow222.7%
Simplified22.7%
Taylor expanded in A around inf 26.5%
Final simplification17.2%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C))))
(t_1 (+ (* B B) (* -4.0 (* A C))))
(t_2 (* 2.0 (* F t_0))))
(if (<= A -9000.0)
(/ (- (sqrt (* t_2 (* 2.0 C)))) t_0)
(if (<= A 8.2e+30)
(/ (- (sqrt (* t_2 (+ B (+ A C))))) t_0)
(* (sqrt (* 2.0 (* (* F t_1) (+ A (+ A C))))) (/ -1.0 t_1))))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = (B * B) + (-4.0 * (A * C));
double t_2 = 2.0 * (F * t_0);
double tmp;
if (A <= -9000.0) {
tmp = -sqrt((t_2 * (2.0 * C))) / t_0;
} else if (A <= 8.2e+30) {
tmp = -sqrt((t_2 * (B + (A + C)))) / t_0;
} else {
tmp = sqrt((2.0 * ((F * t_1) * (A + (A + C))))) * (-1.0 / t_1);
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = (b * b) - (4.0d0 * (a * c))
t_1 = (b * b) + ((-4.0d0) * (a * c))
t_2 = 2.0d0 * (f * t_0)
if (a <= (-9000.0d0)) then
tmp = -sqrt((t_2 * (2.0d0 * c))) / t_0
else if (a <= 8.2d+30) then
tmp = -sqrt((t_2 * (b + (a + c)))) / t_0
else
tmp = sqrt((2.0d0 * ((f * t_1) * (a + (a + c))))) * ((-1.0d0) / t_1)
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = (B * B) + (-4.0 * (A * C));
double t_2 = 2.0 * (F * t_0);
double tmp;
if (A <= -9000.0) {
tmp = -Math.sqrt((t_2 * (2.0 * C))) / t_0;
} else if (A <= 8.2e+30) {
tmp = -Math.sqrt((t_2 * (B + (A + C)))) / t_0;
} else {
tmp = Math.sqrt((2.0 * ((F * t_1) * (A + (A + C))))) * (-1.0 / t_1);
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) t_1 = (B * B) + (-4.0 * (A * C)) t_2 = 2.0 * (F * t_0) tmp = 0 if A <= -9000.0: tmp = -math.sqrt((t_2 * (2.0 * C))) / t_0 elif A <= 8.2e+30: tmp = -math.sqrt((t_2 * (B + (A + C)))) / t_0 else: tmp = math.sqrt((2.0 * ((F * t_1) * (A + (A + C))))) * (-1.0 / t_1) return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) t_1 = Float64(Float64(B * B) + Float64(-4.0 * Float64(A * C))) t_2 = Float64(2.0 * Float64(F * t_0)) tmp = 0.0 if (A <= -9000.0) tmp = Float64(Float64(-sqrt(Float64(t_2 * Float64(2.0 * C)))) / t_0); elseif (A <= 8.2e+30) tmp = Float64(Float64(-sqrt(Float64(t_2 * Float64(B + Float64(A + C))))) / t_0); else tmp = Float64(sqrt(Float64(2.0 * Float64(Float64(F * t_1) * Float64(A + Float64(A + C))))) * Float64(-1.0 / t_1)); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); t_1 = (B * B) + (-4.0 * (A * C)); t_2 = 2.0 * (F * t_0); tmp = 0.0; if (A <= -9000.0) tmp = -sqrt((t_2 * (2.0 * C))) / t_0; elseif (A <= 8.2e+30) tmp = -sqrt((t_2 * (B + (A + C)))) / t_0; else tmp = sqrt((2.0 * ((F * t_1) * (A + (A + C))))) * (-1.0 / t_1); end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(B * B), $MachinePrecision] + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[A, -9000.0], N[((-N[Sqrt[N[(t$95$2 * N[(2.0 * C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], If[LessEqual[A, 8.2e+30], N[((-N[Sqrt[N[(t$95$2 * N[(B + N[(A + C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], N[(N[Sqrt[N[(2.0 * N[(N[(F * t$95$1), $MachinePrecision] * N[(A + N[(A + C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(-1.0 / t$95$1), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
t_1 := B \cdot B + -4 \cdot \left(A \cdot C\right)\\
t_2 := 2 \cdot \left(F \cdot t_0\right)\\
\mathbf{if}\;A \leq -9000:\\
\;\;\;\;\frac{-\sqrt{t_2 \cdot \left(2 \cdot C\right)}}{t_0}\\
\mathbf{elif}\;A \leq 8.2 \cdot 10^{+30}:\\
\;\;\;\;\frac{-\sqrt{t_2 \cdot \left(B + \left(A + C\right)\right)}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2 \cdot \left(\left(F \cdot t_1\right) \cdot \left(A + \left(A + C\right)\right)\right)} \cdot \frac{-1}{t_1}\\
\end{array}
\end{array}
if A < -9e3Initial program 5.6%
associate-*l*5.6%
unpow25.6%
+-commutative5.6%
unpow25.6%
associate-*l*5.6%
unpow25.6%
Simplified5.6%
Taylor expanded in A around -inf 4.9%
associate-+r+4.9%
mul-1-neg4.9%
unsub-neg4.9%
distribute-lft-out4.9%
associate-/l*4.9%
unpow24.9%
unpow24.9%
times-frac4.9%
unpow24.9%
Simplified4.9%
Taylor expanded in A around inf 17.0%
if -9e3 < A < 8.20000000000000011e30Initial program 28.8%
associate-*l*28.8%
unpow228.8%
+-commutative28.8%
unpow228.8%
associate-*l*28.8%
unpow228.8%
Simplified28.8%
unpow228.8%
hypot-udef30.4%
add-exp-log29.0%
Applied egg-rr29.0%
Taylor expanded in B around inf 9.2%
if 8.20000000000000011e30 < A Initial program 21.0%
associate-*l*21.0%
unpow221.0%
+-commutative21.0%
unpow221.0%
associate-*l*21.0%
unpow221.0%
Simplified21.0%
Taylor expanded in A around inf 29.6%
div-inv29.7%
associate-*l*29.7%
*-commutative29.7%
cancel-sign-sub-inv29.7%
*-commutative29.7%
metadata-eval29.7%
associate-+l+29.7%
cancel-sign-sub-inv29.7%
Applied egg-rr29.7%
Final simplification16.8%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))) (t_1 (* 2.0 (* F t_0))))
(if (<= A -220.0)
(/ (- (sqrt (* t_1 (* 2.0 C)))) t_0)
(if (<= A 9.8e+30)
(/ (- (sqrt (* t_1 (+ B (+ A C))))) t_0)
(/ (- (sqrt (* t_1 (+ A (+ A C))))) t_0)))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = 2.0 * (F * t_0);
double tmp;
if (A <= -220.0) {
tmp = -sqrt((t_1 * (2.0 * C))) / t_0;
} else if (A <= 9.8e+30) {
tmp = -sqrt((t_1 * (B + (A + C)))) / t_0;
} else {
tmp = -sqrt((t_1 * (A + (A + C)))) / t_0;
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (b * b) - (4.0d0 * (a * c))
t_1 = 2.0d0 * (f * t_0)
if (a <= (-220.0d0)) then
tmp = -sqrt((t_1 * (2.0d0 * c))) / t_0
else if (a <= 9.8d+30) then
tmp = -sqrt((t_1 * (b + (a + c)))) / t_0
else
tmp = -sqrt((t_1 * (a + (a + c)))) / t_0
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = 2.0 * (F * t_0);
double tmp;
if (A <= -220.0) {
tmp = -Math.sqrt((t_1 * (2.0 * C))) / t_0;
} else if (A <= 9.8e+30) {
tmp = -Math.sqrt((t_1 * (B + (A + C)))) / t_0;
} else {
tmp = -Math.sqrt((t_1 * (A + (A + C)))) / t_0;
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) t_1 = 2.0 * (F * t_0) tmp = 0 if A <= -220.0: tmp = -math.sqrt((t_1 * (2.0 * C))) / t_0 elif A <= 9.8e+30: tmp = -math.sqrt((t_1 * (B + (A + C)))) / t_0 else: tmp = -math.sqrt((t_1 * (A + (A + C)))) / t_0 return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) t_1 = Float64(2.0 * Float64(F * t_0)) tmp = 0.0 if (A <= -220.0) tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(2.0 * C)))) / t_0); elseif (A <= 9.8e+30) tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(B + Float64(A + C))))) / t_0); else tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(A + Float64(A + C))))) / t_0); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); t_1 = 2.0 * (F * t_0); tmp = 0.0; if (A <= -220.0) tmp = -sqrt((t_1 * (2.0 * C))) / t_0; elseif (A <= 9.8e+30) tmp = -sqrt((t_1 * (B + (A + C)))) / t_0; else tmp = -sqrt((t_1 * (A + (A + C)))) / t_0; end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[A, -220.0], N[((-N[Sqrt[N[(t$95$1 * N[(2.0 * C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], If[LessEqual[A, 9.8e+30], N[((-N[Sqrt[N[(t$95$1 * N[(B + N[(A + C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], N[((-N[Sqrt[N[(t$95$1 * N[(A + N[(A + C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
t_1 := 2 \cdot \left(F \cdot t_0\right)\\
\mathbf{if}\;A \leq -220:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(2 \cdot C\right)}}{t_0}\\
\mathbf{elif}\;A \leq 9.8 \cdot 10^{+30}:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(B + \left(A + C\right)\right)}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(A + \left(A + C\right)\right)}}{t_0}\\
\end{array}
\end{array}
if A < -220Initial program 5.6%
associate-*l*5.6%
unpow25.6%
+-commutative5.6%
unpow25.6%
associate-*l*5.6%
unpow25.6%
Simplified5.6%
Taylor expanded in A around -inf 4.9%
associate-+r+4.9%
mul-1-neg4.9%
unsub-neg4.9%
distribute-lft-out4.9%
associate-/l*4.9%
unpow24.9%
unpow24.9%
times-frac4.9%
unpow24.9%
Simplified4.9%
Taylor expanded in A around inf 17.0%
if -220 < A < 9.79999999999999969e30Initial program 28.8%
associate-*l*28.8%
unpow228.8%
+-commutative28.8%
unpow228.8%
associate-*l*28.8%
unpow228.8%
Simplified28.8%
unpow228.8%
hypot-udef30.4%
add-exp-log29.0%
Applied egg-rr29.0%
Taylor expanded in B around inf 9.2%
if 9.79999999999999969e30 < A Initial program 21.0%
associate-*l*21.0%
unpow221.0%
+-commutative21.0%
unpow221.0%
associate-*l*21.0%
unpow221.0%
Simplified21.0%
Taylor expanded in A around inf 29.6%
Final simplification16.8%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))) (t_1 (* 2.0 (* F t_0))))
(if (<= A 4.8e-99)
(/ (- (sqrt (* t_1 (* 2.0 C)))) t_0)
(/ (- (sqrt (* t_1 (+ A (+ A C))))) t_0))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = 2.0 * (F * t_0);
double tmp;
if (A <= 4.8e-99) {
tmp = -sqrt((t_1 * (2.0 * C))) / t_0;
} else {
tmp = -sqrt((t_1 * (A + (A + C)))) / t_0;
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (b * b) - (4.0d0 * (a * c))
t_1 = 2.0d0 * (f * t_0)
if (a <= 4.8d-99) then
tmp = -sqrt((t_1 * (2.0d0 * c))) / t_0
else
tmp = -sqrt((t_1 * (a + (a + c)))) / t_0
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double t_1 = 2.0 * (F * t_0);
double tmp;
if (A <= 4.8e-99) {
tmp = -Math.sqrt((t_1 * (2.0 * C))) / t_0;
} else {
tmp = -Math.sqrt((t_1 * (A + (A + C)))) / t_0;
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) t_1 = 2.0 * (F * t_0) tmp = 0 if A <= 4.8e-99: tmp = -math.sqrt((t_1 * (2.0 * C))) / t_0 else: tmp = -math.sqrt((t_1 * (A + (A + C)))) / t_0 return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) t_1 = Float64(2.0 * Float64(F * t_0)) tmp = 0.0 if (A <= 4.8e-99) tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(2.0 * C)))) / t_0); else tmp = Float64(Float64(-sqrt(Float64(t_1 * Float64(A + Float64(A + C))))) / t_0); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); t_1 = 2.0 * (F * t_0); tmp = 0.0; if (A <= 4.8e-99) tmp = -sqrt((t_1 * (2.0 * C))) / t_0; else tmp = -sqrt((t_1 * (A + (A + C)))) / t_0; end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[A, 4.8e-99], N[((-N[Sqrt[N[(t$95$1 * N[(2.0 * C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], N[((-N[Sqrt[N[(t$95$1 * N[(A + N[(A + C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
t_1 := 2 \cdot \left(F \cdot t_0\right)\\
\mathbf{if}\;A \leq 4.8 \cdot 10^{-99}:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(2 \cdot C\right)}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{-\sqrt{t_1 \cdot \left(A + \left(A + C\right)\right)}}{t_0}\\
\end{array}
\end{array}
if A < 4.8000000000000001e-99Initial program 18.6%
associate-*l*18.6%
unpow218.6%
+-commutative18.6%
unpow218.6%
associate-*l*18.6%
unpow218.6%
Simplified18.6%
Taylor expanded in A around -inf 3.9%
associate-+r+3.9%
mul-1-neg3.9%
unsub-neg3.9%
distribute-lft-out3.9%
associate-/l*4.0%
unpow24.0%
unpow24.0%
times-frac5.9%
unpow25.9%
Simplified5.9%
Taylor expanded in A around inf 15.8%
if 4.8000000000000001e-99 < A Initial program 24.2%
associate-*l*24.2%
unpow224.2%
+-commutative24.2%
unpow224.2%
associate-*l*24.2%
unpow224.2%
Simplified24.2%
Taylor expanded in A around inf 24.0%
Final simplification18.9%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))))
(if (<= A 2.8e-67)
(/ (- (sqrt (* (* 2.0 (* F t_0)) (* 2.0 C)))) t_0)
(/ (- (sqrt (* (+ A (+ A C)) (* 2.0 (* F (* -4.0 (* A C))))))) t_0))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (A <= 2.8e-67) {
tmp = -sqrt(((2.0 * (F * t_0)) * (2.0 * C))) / t_0;
} else {
tmp = -sqrt(((A + (A + C)) * (2.0 * (F * (-4.0 * (A * C)))))) / t_0;
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: tmp
t_0 = (b * b) - (4.0d0 * (a * c))
if (a <= 2.8d-67) then
tmp = -sqrt(((2.0d0 * (f * t_0)) * (2.0d0 * c))) / t_0
else
tmp = -sqrt(((a + (a + c)) * (2.0d0 * (f * ((-4.0d0) * (a * c)))))) / t_0
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (A <= 2.8e-67) {
tmp = -Math.sqrt(((2.0 * (F * t_0)) * (2.0 * C))) / t_0;
} else {
tmp = -Math.sqrt(((A + (A + C)) * (2.0 * (F * (-4.0 * (A * C)))))) / t_0;
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) tmp = 0 if A <= 2.8e-67: tmp = -math.sqrt(((2.0 * (F * t_0)) * (2.0 * C))) / t_0 else: tmp = -math.sqrt(((A + (A + C)) * (2.0 * (F * (-4.0 * (A * C)))))) / t_0 return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) tmp = 0.0 if (A <= 2.8e-67) tmp = Float64(Float64(-sqrt(Float64(Float64(2.0 * Float64(F * t_0)) * Float64(2.0 * C)))) / t_0); else tmp = Float64(Float64(-sqrt(Float64(Float64(A + Float64(A + C)) * Float64(2.0 * Float64(F * Float64(-4.0 * Float64(A * C))))))) / t_0); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); tmp = 0.0; if (A <= 2.8e-67) tmp = -sqrt(((2.0 * (F * t_0)) * (2.0 * C))) / t_0; else tmp = -sqrt(((A + (A + C)) * (2.0 * (F * (-4.0 * (A * C)))))) / t_0; end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[A, 2.8e-67], N[((-N[Sqrt[N[(N[(2.0 * N[(F * t$95$0), $MachinePrecision]), $MachinePrecision] * N[(2.0 * C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], N[((-N[Sqrt[N[(N[(A + N[(A + C), $MachinePrecision]), $MachinePrecision] * N[(2.0 * N[(F * N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
\mathbf{if}\;A \leq 2.8 \cdot 10^{-67}:\\
\;\;\;\;\frac{-\sqrt{\left(2 \cdot \left(F \cdot t_0\right)\right) \cdot \left(2 \cdot C\right)}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{-\sqrt{\left(A + \left(A + C\right)\right) \cdot \left(2 \cdot \left(F \cdot \left(-4 \cdot \left(A \cdot C\right)\right)\right)\right)}}{t_0}\\
\end{array}
\end{array}
if A < 2.8000000000000001e-67Initial program 18.7%
associate-*l*18.7%
unpow218.7%
+-commutative18.7%
unpow218.7%
associate-*l*18.7%
unpow218.7%
Simplified18.7%
Taylor expanded in A around -inf 3.9%
associate-+r+3.9%
mul-1-neg3.9%
unsub-neg3.9%
distribute-lft-out3.9%
associate-/l*3.9%
unpow23.9%
unpow23.9%
times-frac5.8%
unpow25.8%
Simplified5.8%
Taylor expanded in A around inf 15.4%
if 2.8000000000000001e-67 < A Initial program 24.4%
associate-*l*24.4%
unpow224.4%
+-commutative24.4%
unpow224.4%
associate-*l*24.4%
unpow224.4%
Simplified24.4%
Taylor expanded in A around inf 25.2%
Taylor expanded in B around 0 18.1%
Final simplification16.4%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))))
(if (<= A 1.2e-67)
(/ (- (sqrt (* 2.0 (* -8.0 (* A (* F (* C C))))))) t_0)
(/ (- (sqrt (* (+ A (+ A C)) (* 2.0 (* F (* -4.0 (* A C))))))) t_0))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (A <= 1.2e-67) {
tmp = -sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / t_0;
} else {
tmp = -sqrt(((A + (A + C)) * (2.0 * (F * (-4.0 * (A * C)))))) / t_0;
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: tmp
t_0 = (b * b) - (4.0d0 * (a * c))
if (a <= 1.2d-67) then
tmp = -sqrt((2.0d0 * ((-8.0d0) * (a * (f * (c * c)))))) / t_0
else
tmp = -sqrt(((a + (a + c)) * (2.0d0 * (f * ((-4.0d0) * (a * c)))))) / t_0
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (A <= 1.2e-67) {
tmp = -Math.sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / t_0;
} else {
tmp = -Math.sqrt(((A + (A + C)) * (2.0 * (F * (-4.0 * (A * C)))))) / t_0;
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) tmp = 0 if A <= 1.2e-67: tmp = -math.sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / t_0 else: tmp = -math.sqrt(((A + (A + C)) * (2.0 * (F * (-4.0 * (A * C)))))) / t_0 return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) tmp = 0.0 if (A <= 1.2e-67) tmp = Float64(Float64(-sqrt(Float64(2.0 * Float64(-8.0 * Float64(A * Float64(F * Float64(C * C))))))) / t_0); else tmp = Float64(Float64(-sqrt(Float64(Float64(A + Float64(A + C)) * Float64(2.0 * Float64(F * Float64(-4.0 * Float64(A * C))))))) / t_0); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); tmp = 0.0; if (A <= 1.2e-67) tmp = -sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / t_0; else tmp = -sqrt(((A + (A + C)) * (2.0 * (F * (-4.0 * (A * C)))))) / t_0; end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[A, 1.2e-67], N[((-N[Sqrt[N[(2.0 * N[(-8.0 * N[(A * N[(F * N[(C * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], N[((-N[Sqrt[N[(N[(A + N[(A + C), $MachinePrecision]), $MachinePrecision] * N[(2.0 * N[(F * N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
\mathbf{if}\;A \leq 1.2 \cdot 10^{-67}:\\
\;\;\;\;\frac{-\sqrt{2 \cdot \left(-8 \cdot \left(A \cdot \left(F \cdot \left(C \cdot C\right)\right)\right)\right)}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{-\sqrt{\left(A + \left(A + C\right)\right) \cdot \left(2 \cdot \left(F \cdot \left(-4 \cdot \left(A \cdot C\right)\right)\right)\right)}}{t_0}\\
\end{array}
\end{array}
if A < 1.2e-67Initial program 18.7%
associate-*l*18.7%
unpow218.7%
+-commutative18.7%
unpow218.7%
associate-*l*18.7%
unpow218.7%
Simplified18.7%
*-un-lft-identity18.7%
associate-*l*18.7%
*-commutative18.7%
*-commutative18.7%
unpow218.7%
hypot-udef20.8%
+-commutative20.8%
associate-+r+21.9%
Applied egg-rr21.9%
*-lft-identity21.9%
associate-*r*21.9%
fma-neg21.9%
*-commutative21.9%
*-commutative21.9%
distribute-rgt-neg-in21.9%
*-commutative21.9%
metadata-eval21.9%
associate-*r*21.9%
associate-*r*21.9%
associate-*l*18.0%
Simplified18.0%
Taylor expanded in C around inf 9.9%
unpow29.9%
Simplified9.9%
if 1.2e-67 < A Initial program 24.4%
associate-*l*24.4%
unpow224.4%
+-commutative24.4%
unpow224.4%
associate-*l*24.4%
unpow224.4%
Simplified24.4%
Taylor expanded in A around inf 25.2%
Taylor expanded in B around 0 18.1%
Final simplification12.8%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(let* ((t_0 (- (* B B) (* 4.0 (* A C)))))
(if (<= C -1.2e-173)
(/ (- (sqrt (* 2.0 (* -8.0 (* (* A A) (* C F)))))) t_0)
(if (<= C 1.3e-151)
(* (pow (* A F) 0.5) (/ (- 2.0) B))
(/ (- (sqrt (* 2.0 (* -8.0 (* A (* F (* C C))))))) t_0)))))B = abs(B);
double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (C <= -1.2e-173) {
tmp = -sqrt((2.0 * (-8.0 * ((A * A) * (C * F))))) / t_0;
} else if (C <= 1.3e-151) {
tmp = pow((A * F), 0.5) * (-2.0 / B);
} else {
tmp = -sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / t_0;
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: tmp
t_0 = (b * b) - (4.0d0 * (a * c))
if (c <= (-1.2d-173)) then
tmp = -sqrt((2.0d0 * ((-8.0d0) * ((a * a) * (c * f))))) / t_0
else if (c <= 1.3d-151) then
tmp = ((a * f) ** 0.5d0) * (-2.0d0 / b)
else
tmp = -sqrt((2.0d0 * ((-8.0d0) * (a * (f * (c * c)))))) / t_0
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double t_0 = (B * B) - (4.0 * (A * C));
double tmp;
if (C <= -1.2e-173) {
tmp = -Math.sqrt((2.0 * (-8.0 * ((A * A) * (C * F))))) / t_0;
} else if (C <= 1.3e-151) {
tmp = Math.pow((A * F), 0.5) * (-2.0 / B);
} else {
tmp = -Math.sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / t_0;
}
return tmp;
}
B = abs(B) def code(A, B, C, F): t_0 = (B * B) - (4.0 * (A * C)) tmp = 0 if C <= -1.2e-173: tmp = -math.sqrt((2.0 * (-8.0 * ((A * A) * (C * F))))) / t_0 elif C <= 1.3e-151: tmp = math.pow((A * F), 0.5) * (-2.0 / B) else: tmp = -math.sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / t_0 return tmp
B = abs(B) function code(A, B, C, F) t_0 = Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))) tmp = 0.0 if (C <= -1.2e-173) tmp = Float64(Float64(-sqrt(Float64(2.0 * Float64(-8.0 * Float64(Float64(A * A) * Float64(C * F)))))) / t_0); elseif (C <= 1.3e-151) tmp = Float64((Float64(A * F) ^ 0.5) * Float64(Float64(-2.0) / B)); else tmp = Float64(Float64(-sqrt(Float64(2.0 * Float64(-8.0 * Float64(A * Float64(F * Float64(C * C))))))) / t_0); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) t_0 = (B * B) - (4.0 * (A * C)); tmp = 0.0; if (C <= -1.2e-173) tmp = -sqrt((2.0 * (-8.0 * ((A * A) * (C * F))))) / t_0; elseif (C <= 1.3e-151) tmp = ((A * F) ^ 0.5) * (-2.0 / B); else tmp = -sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / t_0; end tmp_2 = tmp; end
NOTE: B should be positive before calling this function
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[C, -1.2e-173], N[((-N[Sqrt[N[(2.0 * N[(-8.0 * N[(N[(A * A), $MachinePrecision] * N[(C * F), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision], If[LessEqual[C, 1.3e-151], N[(N[Power[N[(A * F), $MachinePrecision], 0.5], $MachinePrecision] * N[((-2.0) / B), $MachinePrecision]), $MachinePrecision], N[((-N[Sqrt[N[(2.0 * N[(-8.0 * N[(A * N[(F * N[(C * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$0), $MachinePrecision]]]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
t_0 := B \cdot B - 4 \cdot \left(A \cdot C\right)\\
\mathbf{if}\;C \leq -1.2 \cdot 10^{-173}:\\
\;\;\;\;\frac{-\sqrt{2 \cdot \left(-8 \cdot \left(\left(A \cdot A\right) \cdot \left(C \cdot F\right)\right)\right)}}{t_0}\\
\mathbf{elif}\;C \leq 1.3 \cdot 10^{-151}:\\
\;\;\;\;{\left(A \cdot F\right)}^{0.5} \cdot \frac{-2}{B}\\
\mathbf{else}:\\
\;\;\;\;\frac{-\sqrt{2 \cdot \left(-8 \cdot \left(A \cdot \left(F \cdot \left(C \cdot C\right)\right)\right)\right)}}{t_0}\\
\end{array}
\end{array}
if C < -1.20000000000000008e-173Initial program 14.0%
associate-*l*14.0%
unpow214.0%
+-commutative14.0%
unpow214.0%
associate-*l*14.0%
unpow214.0%
Simplified14.0%
*-un-lft-identity14.0%
associate-*l*14.0%
*-commutative14.0%
*-commutative14.0%
unpow214.0%
hypot-udef17.2%
+-commutative17.2%
associate-+r+18.0%
Applied egg-rr18.0%
*-lft-identity18.0%
associate-*r*18.0%
fma-neg18.0%
*-commutative18.0%
*-commutative18.0%
distribute-rgt-neg-in18.0%
*-commutative18.0%
metadata-eval18.0%
associate-*r*18.0%
associate-*r*18.0%
associate-*l*14.9%
Simplified14.9%
Taylor expanded in B around 0 15.6%
unpow215.6%
Simplified15.6%
if -1.20000000000000008e-173 < C < 1.3e-151Initial program 29.4%
associate-*l*29.4%
unpow229.4%
+-commutative29.4%
unpow229.4%
associate-*l*29.4%
unpow229.4%
Simplified29.4%
Taylor expanded in A around inf 15.5%
Taylor expanded in C around 0 6.9%
mul-1-neg6.9%
unpow26.9%
rem-square-sqrt6.9%
Simplified6.9%
pow1/27.0%
*-commutative7.0%
Applied egg-rr7.0%
if 1.3e-151 < C Initial program 21.4%
associate-*l*21.4%
unpow221.4%
+-commutative21.4%
unpow221.4%
associate-*l*21.4%
unpow221.4%
Simplified21.4%
*-un-lft-identity21.4%
associate-*l*21.4%
*-commutative21.4%
*-commutative21.4%
unpow221.4%
hypot-udef25.5%
+-commutative25.5%
associate-+r+25.7%
Applied egg-rr25.7%
*-lft-identity25.7%
associate-*r*25.7%
fma-neg25.7%
*-commutative25.7%
*-commutative25.7%
distribute-rgt-neg-in25.7%
*-commutative25.7%
metadata-eval25.7%
associate-*r*25.7%
associate-*r*25.7%
associate-*l*20.6%
Simplified20.6%
Taylor expanded in C around inf 10.7%
unpow210.7%
Simplified10.7%
Final simplification11.6%
NOTE: B should be positive before calling this function
(FPCore (A B C F)
:precision binary64
(if (<= A 5.3e-102)
(/
(- (sqrt (* 2.0 (* -8.0 (* A (* F (* C C)))))))
(- (* B B) (* 4.0 (* A C))))
(* (pow (* A F) 0.5) (/ (- 2.0) B))))B = abs(B);
double code(double A, double B, double C, double F) {
double tmp;
if (A <= 5.3e-102) {
tmp = -sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / ((B * B) - (4.0 * (A * C)));
} else {
tmp = pow((A * F), 0.5) * (-2.0 / B);
}
return tmp;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (a <= 5.3d-102) then
tmp = -sqrt((2.0d0 * ((-8.0d0) * (a * (f * (c * c)))))) / ((b * b) - (4.0d0 * (a * c)))
else
tmp = ((a * f) ** 0.5d0) * (-2.0d0 / b)
end if
code = tmp
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
double tmp;
if (A <= 5.3e-102) {
tmp = -Math.sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / ((B * B) - (4.0 * (A * C)));
} else {
tmp = Math.pow((A * F), 0.5) * (-2.0 / B);
}
return tmp;
}
B = abs(B) def code(A, B, C, F): tmp = 0 if A <= 5.3e-102: tmp = -math.sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / ((B * B) - (4.0 * (A * C))) else: tmp = math.pow((A * F), 0.5) * (-2.0 / B) return tmp
B = abs(B) function code(A, B, C, F) tmp = 0.0 if (A <= 5.3e-102) tmp = Float64(Float64(-sqrt(Float64(2.0 * Float64(-8.0 * Float64(A * Float64(F * Float64(C * C))))))) / Float64(Float64(B * B) - Float64(4.0 * Float64(A * C)))); else tmp = Float64((Float64(A * F) ^ 0.5) * Float64(Float64(-2.0) / B)); end return tmp end
B = abs(B) function tmp_2 = code(A, B, C, F) tmp = 0.0; if (A <= 5.3e-102) tmp = -sqrt((2.0 * (-8.0 * (A * (F * (C * C)))))) / ((B * B) - (4.0 * (A * C))); else tmp = ((A * F) ^ 0.5) * (-2.0 / B); end tmp_2 = tmp; end
NOTE: B should be positive before calling this function code[A_, B_, C_, F_] := If[LessEqual[A, 5.3e-102], N[((-N[Sqrt[N[(2.0 * N[(-8.0 * N[(A * N[(F * N[(C * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[N[(A * F), $MachinePrecision], 0.5], $MachinePrecision] * N[((-2.0) / B), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
B = |B|\\
\\
\begin{array}{l}
\mathbf{if}\;A \leq 5.3 \cdot 10^{-102}:\\
\;\;\;\;\frac{-\sqrt{2 \cdot \left(-8 \cdot \left(A \cdot \left(F \cdot \left(C \cdot C\right)\right)\right)\right)}}{B \cdot B - 4 \cdot \left(A \cdot C\right)}\\
\mathbf{else}:\\
\;\;\;\;{\left(A \cdot F\right)}^{0.5} \cdot \frac{-2}{B}\\
\end{array}
\end{array}
if A < 5.3000000000000003e-102Initial program 18.6%
associate-*l*18.6%
unpow218.6%
+-commutative18.6%
unpow218.6%
associate-*l*18.6%
unpow218.6%
Simplified18.6%
*-un-lft-identity18.6%
associate-*l*18.6%
*-commutative18.6%
*-commutative18.6%
unpow218.6%
hypot-udef20.8%
+-commutative20.8%
associate-+r+21.9%
Applied egg-rr21.9%
*-lft-identity21.9%
associate-*r*21.9%
fma-neg21.9%
*-commutative21.9%
*-commutative21.9%
distribute-rgt-neg-in21.9%
*-commutative21.9%
metadata-eval21.9%
associate-*r*21.9%
associate-*r*21.9%
associate-*l*17.8%
Simplified17.8%
Taylor expanded in C around inf 10.1%
unpow210.1%
Simplified10.1%
if 5.3000000000000003e-102 < A Initial program 24.2%
associate-*l*24.2%
unpow224.2%
+-commutative24.2%
unpow224.2%
associate-*l*24.2%
unpow224.2%
Simplified24.2%
Taylor expanded in A around inf 24.0%
Taylor expanded in C around 0 6.8%
mul-1-neg6.8%
unpow26.8%
rem-square-sqrt6.9%
Simplified6.9%
pow1/27.0%
*-commutative7.0%
Applied egg-rr7.0%
Final simplification8.9%
NOTE: B should be positive before calling this function (FPCore (A B C F) :precision binary64 (* (pow (* A F) 0.5) (/ (- 2.0) B)))
B = abs(B);
double code(double A, double B, double C, double F) {
return pow((A * F), 0.5) * (-2.0 / B);
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
code = ((a * f) ** 0.5d0) * (-2.0d0 / b)
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
return Math.pow((A * F), 0.5) * (-2.0 / B);
}
B = abs(B) def code(A, B, C, F): return math.pow((A * F), 0.5) * (-2.0 / B)
B = abs(B) function code(A, B, C, F) return Float64((Float64(A * F) ^ 0.5) * Float64(Float64(-2.0) / B)) end
B = abs(B) function tmp = code(A, B, C, F) tmp = ((A * F) ^ 0.5) * (-2.0 / B); end
NOTE: B should be positive before calling this function code[A_, B_, C_, F_] := N[(N[Power[N[(A * F), $MachinePrecision], 0.5], $MachinePrecision] * N[((-2.0) / B), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
B = |B|\\
\\
{\left(A \cdot F\right)}^{0.5} \cdot \frac{-2}{B}
\end{array}
Initial program 20.7%
associate-*l*20.7%
unpow220.7%
+-commutative20.7%
unpow220.7%
associate-*l*20.7%
unpow220.7%
Simplified20.7%
Taylor expanded in A around inf 10.9%
Taylor expanded in C around 0 3.3%
mul-1-neg3.3%
unpow23.3%
rem-square-sqrt3.4%
Simplified3.4%
pow1/23.5%
*-commutative3.5%
Applied egg-rr3.5%
Final simplification3.5%
NOTE: B should be positive before calling this function (FPCore (A B C F) :precision binary64 (* (/ 2.0 B) (- (sqrt (* A F)))))
B = abs(B);
double code(double A, double B, double C, double F) {
return (2.0 / B) * -sqrt((A * F));
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
code = (2.0d0 / b) * -sqrt((a * f))
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
return (2.0 / B) * -Math.sqrt((A * F));
}
B = abs(B) def code(A, B, C, F): return (2.0 / B) * -math.sqrt((A * F))
B = abs(B) function code(A, B, C, F) return Float64(Float64(2.0 / B) * Float64(-sqrt(Float64(A * F)))) end
B = abs(B) function tmp = code(A, B, C, F) tmp = (2.0 / B) * -sqrt((A * F)); end
NOTE: B should be positive before calling this function code[A_, B_, C_, F_] := N[(N[(2.0 / B), $MachinePrecision] * (-N[Sqrt[N[(A * F), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]
\begin{array}{l}
B = |B|\\
\\
\frac{2}{B} \cdot \left(-\sqrt{A \cdot F}\right)
\end{array}
Initial program 20.7%
associate-*l*20.7%
unpow220.7%
+-commutative20.7%
unpow220.7%
associate-*l*20.7%
unpow220.7%
Simplified20.7%
Taylor expanded in A around inf 10.9%
Taylor expanded in C around 0 3.3%
mul-1-neg3.3%
unpow23.3%
rem-square-sqrt3.4%
Simplified3.4%
Final simplification3.4%
NOTE: B should be positive before calling this function (FPCore (A B C F) :precision binary64 (/ (* 2.0 (- (sqrt (* A F)))) B))
B = abs(B);
double code(double A, double B, double C, double F) {
return (2.0 * -sqrt((A * F))) / B;
}
NOTE: B should be positive before calling this function
real(8) function code(a, b, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: f
code = (2.0d0 * -sqrt((a * f))) / b
end function
B = Math.abs(B);
public static double code(double A, double B, double C, double F) {
return (2.0 * -Math.sqrt((A * F))) / B;
}
B = abs(B) def code(A, B, C, F): return (2.0 * -math.sqrt((A * F))) / B
B = abs(B) function code(A, B, C, F) return Float64(Float64(2.0 * Float64(-sqrt(Float64(A * F)))) / B) end
B = abs(B) function tmp = code(A, B, C, F) tmp = (2.0 * -sqrt((A * F))) / B; end
NOTE: B should be positive before calling this function code[A_, B_, C_, F_] := N[(N[(2.0 * (-N[Sqrt[N[(A * F), $MachinePrecision]], $MachinePrecision])), $MachinePrecision] / B), $MachinePrecision]
\begin{array}{l}
B = |B|\\
\\
\frac{2 \cdot \left(-\sqrt{A \cdot F}\right)}{B}
\end{array}
Initial program 20.7%
associate-*l*20.7%
unpow220.7%
+-commutative20.7%
unpow220.7%
associate-*l*20.7%
unpow220.7%
Simplified20.7%
Taylor expanded in A around inf 10.9%
Taylor expanded in C around 0 3.3%
mul-1-neg3.3%
unpow23.3%
rem-square-sqrt3.4%
Simplified3.4%
associate-*l/3.4%
*-commutative3.4%
Applied egg-rr3.4%
Final simplification3.4%
herbie shell --seed 2023224
(FPCore (A B C F)
:name "ABCF->ab-angle a"
:precision binary64
(/ (- (sqrt (* (* 2.0 (* (- (pow B 2.0) (* (* 4.0 A) C)) F)) (+ (+ A C) (sqrt (+ (pow (- A C) 2.0) (pow B 2.0))))))) (- (pow B 2.0) (* (* 4.0 A) C))))