
(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 18 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}
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(let* ((t_0 (+ (+ A C) (hypot B_m (- A C))))
(t_1 (+ (* B_m B_m) (* -4.0 (* A C))))
(t_2 (* (* 4.0 A) C))
(t_3
(/
(sqrt
(*
(* 2.0 (* (- (pow B_m 2.0) t_2) F))
(+ (+ A C) (sqrt (+ (pow B_m 2.0) (pow (- A C) 2.0))))))
(- t_2 (pow B_m 2.0)))))
(if (<= t_3 -2e-208)
(*
(sqrt (/ t_0 (+ (* B_m B_m) (* A (* C -4.0)))))
(- 0.0 (sqrt (* 2.0 F))))
(if (<= t_3 4e-86)
(/
(sqrt (* t_1 (* (* 2.0 F) (+ (* -0.5 (/ (* B_m B_m) A)) (* 2.0 C)))))
(- (* 4.0 (* A C)) (* B_m B_m)))
(if (<= t_3 INFINITY)
(* (/ (sqrt t_0) (- t_2 (* B_m B_m))) (sqrt (* (* 2.0 F) t_1)))
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m))))))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double t_0 = (A + C) + hypot(B_m, (A - C));
double t_1 = (B_m * B_m) + (-4.0 * (A * C));
double t_2 = (4.0 * A) * C;
double t_3 = sqrt(((2.0 * ((pow(B_m, 2.0) - t_2) * F)) * ((A + C) + sqrt((pow(B_m, 2.0) + pow((A - C), 2.0)))))) / (t_2 - pow(B_m, 2.0));
double tmp;
if (t_3 <= -2e-208) {
tmp = sqrt((t_0 / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - sqrt((2.0 * F)));
} else if (t_3 <= 4e-86) {
tmp = sqrt((t_1 * ((2.0 * F) * ((-0.5 * ((B_m * B_m) / A)) + (2.0 * C))))) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (t_3 <= ((double) INFINITY)) {
tmp = (sqrt(t_0) / (t_2 - (B_m * B_m))) * sqrt(((2.0 * F) * t_1));
} else {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
}
return tmp;
}
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double t_0 = (A + C) + Math.hypot(B_m, (A - C));
double t_1 = (B_m * B_m) + (-4.0 * (A * C));
double t_2 = (4.0 * A) * C;
double t_3 = Math.sqrt(((2.0 * ((Math.pow(B_m, 2.0) - t_2) * F)) * ((A + C) + Math.sqrt((Math.pow(B_m, 2.0) + Math.pow((A - C), 2.0)))))) / (t_2 - Math.pow(B_m, 2.0));
double tmp;
if (t_3 <= -2e-208) {
tmp = Math.sqrt((t_0 / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - Math.sqrt((2.0 * F)));
} else if (t_3 <= 4e-86) {
tmp = Math.sqrt((t_1 * ((2.0 * F) * ((-0.5 * ((B_m * B_m) / A)) + (2.0 * C))))) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (t_3 <= Double.POSITIVE_INFINITY) {
tmp = (Math.sqrt(t_0) / (t_2 - (B_m * B_m))) * Math.sqrt(((2.0 * F) * t_1));
} else {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): t_0 = (A + C) + math.hypot(B_m, (A - C)) t_1 = (B_m * B_m) + (-4.0 * (A * C)) t_2 = (4.0 * A) * C t_3 = math.sqrt(((2.0 * ((math.pow(B_m, 2.0) - t_2) * F)) * ((A + C) + math.sqrt((math.pow(B_m, 2.0) + math.pow((A - C), 2.0)))))) / (t_2 - math.pow(B_m, 2.0)) tmp = 0 if t_3 <= -2e-208: tmp = math.sqrt((t_0 / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - math.sqrt((2.0 * F))) elif t_3 <= 4e-86: tmp = math.sqrt((t_1 * ((2.0 * F) * ((-0.5 * ((B_m * B_m) / A)) + (2.0 * C))))) / ((4.0 * (A * C)) - (B_m * B_m)) elif t_3 <= math.inf: tmp = (math.sqrt(t_0) / (t_2 - (B_m * B_m))) * math.sqrt(((2.0 * F) * t_1)) else: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) return tmp
B_m = abs(B) function code(A, B_m, C, F) t_0 = Float64(Float64(A + C) + hypot(B_m, Float64(A - C))) t_1 = Float64(Float64(B_m * B_m) + Float64(-4.0 * Float64(A * C))) t_2 = Float64(Float64(4.0 * A) * C) t_3 = Float64(sqrt(Float64(Float64(2.0 * Float64(Float64((B_m ^ 2.0) - t_2) * F)) * Float64(Float64(A + C) + sqrt(Float64((B_m ^ 2.0) + (Float64(A - C) ^ 2.0)))))) / Float64(t_2 - (B_m ^ 2.0))) tmp = 0.0 if (t_3 <= -2e-208) tmp = Float64(sqrt(Float64(t_0 / Float64(Float64(B_m * B_m) + Float64(A * Float64(C * -4.0))))) * Float64(0.0 - sqrt(Float64(2.0 * F)))); elseif (t_3 <= 4e-86) tmp = Float64(sqrt(Float64(t_1 * Float64(Float64(2.0 * F) * Float64(Float64(-0.5 * Float64(Float64(B_m * B_m) / A)) + Float64(2.0 * C))))) / Float64(Float64(4.0 * Float64(A * C)) - Float64(B_m * B_m))); elseif (t_3 <= Inf) tmp = Float64(Float64(sqrt(t_0) / Float64(t_2 - Float64(B_m * B_m))) * sqrt(Float64(Float64(2.0 * F) * t_1))); else tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) t_0 = (A + C) + hypot(B_m, (A - C)); t_1 = (B_m * B_m) + (-4.0 * (A * C)); t_2 = (4.0 * A) * C; t_3 = sqrt(((2.0 * (((B_m ^ 2.0) - t_2) * F)) * ((A + C) + sqrt(((B_m ^ 2.0) + ((A - C) ^ 2.0)))))) / (t_2 - (B_m ^ 2.0)); tmp = 0.0; if (t_3 <= -2e-208) tmp = sqrt((t_0 / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - sqrt((2.0 * F))); elseif (t_3 <= 4e-86) tmp = sqrt((t_1 * ((2.0 * F) * ((-0.5 * ((B_m * B_m) / A)) + (2.0 * C))))) / ((4.0 * (A * C)) - (B_m * B_m)); elseif (t_3 <= Inf) tmp = (sqrt(t_0) / (t_2 - (B_m * B_m))) * sqrt(((2.0 * F) * t_1)); else tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision]
code[A_, B$95$m_, C_, F_] := Block[{t$95$0 = N[(N[(A + C), $MachinePrecision] + N[Sqrt[B$95$m ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(4.0 * A), $MachinePrecision] * C), $MachinePrecision]}, Block[{t$95$3 = N[(N[Sqrt[N[(N[(2.0 * N[(N[(N[Power[B$95$m, 2.0], $MachinePrecision] - t$95$2), $MachinePrecision] * F), $MachinePrecision]), $MachinePrecision] * N[(N[(A + C), $MachinePrecision] + N[Sqrt[N[(N[Power[B$95$m, 2.0], $MachinePrecision] + N[Power[N[(A - C), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(t$95$2 - N[Power[B$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, -2e-208], N[(N[Sqrt[N[(t$95$0 / N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(A * N[(C * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(0.0 - N[Sqrt[N[(2.0 * F), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 4e-86], N[(N[Sqrt[N[(t$95$1 * N[(N[(2.0 * F), $MachinePrecision] * N[(N[(-0.5 * N[(N[(B$95$m * B$95$m), $MachinePrecision] / A), $MachinePrecision]), $MachinePrecision] + N[(2.0 * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, Infinity], N[(N[(N[Sqrt[t$95$0], $MachinePrecision] / N[(t$95$2 - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[(2.0 * F), $MachinePrecision] * t$95$1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
t_0 := \left(A + C\right) + \mathsf{hypot}\left(B\_m, A - C\right)\\
t_1 := B\_m \cdot B\_m + -4 \cdot \left(A \cdot C\right)\\
t_2 := \left(4 \cdot A\right) \cdot C\\
t_3 := \frac{\sqrt{\left(2 \cdot \left(\left({B\_m}^{2} - t\_2\right) \cdot F\right)\right) \cdot \left(\left(A + C\right) + \sqrt{{B\_m}^{2} + {\left(A - C\right)}^{2}}\right)}}{t\_2 - {B\_m}^{2}}\\
\mathbf{if}\;t\_3 \leq -2 \cdot 10^{-208}:\\
\;\;\;\;\sqrt{\frac{t\_0}{B\_m \cdot B\_m + A \cdot \left(C \cdot -4\right)}} \cdot \left(0 - \sqrt{2 \cdot F}\right)\\
\mathbf{elif}\;t\_3 \leq 4 \cdot 10^{-86}:\\
\;\;\;\;\frac{\sqrt{t\_1 \cdot \left(\left(2 \cdot F\right) \cdot \left(-0.5 \cdot \frac{B\_m \cdot B\_m}{A} + 2 \cdot C\right)\right)}}{4 \cdot \left(A \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{elif}\;t\_3 \leq \infty:\\
\;\;\;\;\frac{\sqrt{t\_0}}{t\_2 - B\_m \cdot B\_m} \cdot \sqrt{\left(2 \cdot F\right) \cdot t\_1}\\
\mathbf{else}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\end{array}
\end{array}
if (/.f64 (neg.f64 (sqrt.f64 (*.f64 (*.f64 #s(literal 2 binary64) (*.f64 (-.f64 (pow.f64 B #s(literal 2 binary64)) (*.f64 (*.f64 #s(literal 4 binary64) A) C)) F)) (+.f64 (+.f64 A C) (sqrt.f64 (+.f64 (pow.f64 (-.f64 A C) #s(literal 2 binary64)) (pow.f64 B #s(literal 2 binary64)))))))) (-.f64 (pow.f64 B #s(literal 2 binary64)) (*.f64 (*.f64 #s(literal 4 binary64) A) C))) < -2.0000000000000002e-208Initial program 35.6%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified63.0%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr84.6%
if -2.0000000000000002e-208 < (/.f64 (neg.f64 (sqrt.f64 (*.f64 (*.f64 #s(literal 2 binary64) (*.f64 (-.f64 (pow.f64 B #s(literal 2 binary64)) (*.f64 (*.f64 #s(literal 4 binary64) A) C)) F)) (+.f64 (+.f64 A C) (sqrt.f64 (+.f64 (pow.f64 (-.f64 A C) #s(literal 2 binary64)) (pow.f64 B #s(literal 2 binary64)))))))) (-.f64 (pow.f64 B #s(literal 2 binary64)) (*.f64 (*.f64 #s(literal 4 binary64) A) C))) < 4.00000000000000034e-86Initial program 10.8%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified10.8%
associate-*l*N/A
associate-*r*N/A
pow2N/A
*-lowering-*.f64N/A
pow2N/A
associate-*r*N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr10.8%
Taylor expanded in A around -inf
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6443.9%
Simplified43.9%
if 4.00000000000000034e-86 < (/.f64 (neg.f64 (sqrt.f64 (*.f64 (*.f64 #s(literal 2 binary64) (*.f64 (-.f64 (pow.f64 B #s(literal 2 binary64)) (*.f64 (*.f64 #s(literal 4 binary64) A) C)) F)) (+.f64 (+.f64 A C) (sqrt.f64 (+.f64 (pow.f64 (-.f64 A C) #s(literal 2 binary64)) (pow.f64 B #s(literal 2 binary64)))))))) (-.f64 (pow.f64 B #s(literal 2 binary64)) (*.f64 (*.f64 #s(literal 4 binary64) A) C))) < +inf.0Initial program 36.1%
Applied egg-rr74.5%
if +inf.0 < (/.f64 (neg.f64 (sqrt.f64 (*.f64 (*.f64 #s(literal 2 binary64) (*.f64 (-.f64 (pow.f64 B #s(literal 2 binary64)) (*.f64 (*.f64 #s(literal 4 binary64) A) C)) F)) (+.f64 (+.f64 A C) (sqrt.f64 (+.f64 (pow.f64 (-.f64 A C) #s(literal 2 binary64)) (pow.f64 B #s(literal 2 binary64)))))))) (-.f64 (pow.f64 B #s(literal 2 binary64)) (*.f64 (*.f64 #s(literal 4 binary64) A) C))) Initial program 0.0%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified3.0%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr3.9%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6426.6%
Simplified26.6%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6426.6%
Applied egg-rr26.6%
Final simplification54.7%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(let* ((t_0 (+ (+ A C) (hypot B_m (- A C)))))
(if (<= B_m 1.4e-105)
(*
(sqrt t_0)
(/
(sqrt (* (* 2.0 F) (+ (* B_m B_m) (* -4.0 (* A C)))))
(- (* (* 4.0 A) C) (* B_m B_m))))
(if (<= B_m 2.8e+146)
(*
(sqrt (/ t_0 (+ (* B_m B_m) (* A (* C -4.0)))))
(- 0.0 (sqrt (* 2.0 F))))
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m)))))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double t_0 = (A + C) + hypot(B_m, (A - C));
double tmp;
if (B_m <= 1.4e-105) {
tmp = sqrt(t_0) * (sqrt(((2.0 * F) * ((B_m * B_m) + (-4.0 * (A * C))))) / (((4.0 * A) * C) - (B_m * B_m)));
} else if (B_m <= 2.8e+146) {
tmp = sqrt((t_0 / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - sqrt((2.0 * F)));
} else {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
}
return tmp;
}
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double t_0 = (A + C) + Math.hypot(B_m, (A - C));
double tmp;
if (B_m <= 1.4e-105) {
tmp = Math.sqrt(t_0) * (Math.sqrt(((2.0 * F) * ((B_m * B_m) + (-4.0 * (A * C))))) / (((4.0 * A) * C) - (B_m * B_m)));
} else if (B_m <= 2.8e+146) {
tmp = Math.sqrt((t_0 / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - Math.sqrt((2.0 * F)));
} else {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): t_0 = (A + C) + math.hypot(B_m, (A - C)) tmp = 0 if B_m <= 1.4e-105: tmp = math.sqrt(t_0) * (math.sqrt(((2.0 * F) * ((B_m * B_m) + (-4.0 * (A * C))))) / (((4.0 * A) * C) - (B_m * B_m))) elif B_m <= 2.8e+146: tmp = math.sqrt((t_0 / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - math.sqrt((2.0 * F))) else: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) return tmp
B_m = abs(B) function code(A, B_m, C, F) t_0 = Float64(Float64(A + C) + hypot(B_m, Float64(A - C))) tmp = 0.0 if (B_m <= 1.4e-105) tmp = Float64(sqrt(t_0) * Float64(sqrt(Float64(Float64(2.0 * F) * Float64(Float64(B_m * B_m) + Float64(-4.0 * Float64(A * C))))) / Float64(Float64(Float64(4.0 * A) * C) - Float64(B_m * B_m)))); elseif (B_m <= 2.8e+146) tmp = Float64(sqrt(Float64(t_0 / Float64(Float64(B_m * B_m) + Float64(A * Float64(C * -4.0))))) * Float64(0.0 - sqrt(Float64(2.0 * F)))); else tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) t_0 = (A + C) + hypot(B_m, (A - C)); tmp = 0.0; if (B_m <= 1.4e-105) tmp = sqrt(t_0) * (sqrt(((2.0 * F) * ((B_m * B_m) + (-4.0 * (A * C))))) / (((4.0 * A) * C) - (B_m * B_m))); elseif (B_m <= 2.8e+146) tmp = sqrt((t_0 / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - sqrt((2.0 * F))); else tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision]
code[A_, B$95$m_, C_, F_] := Block[{t$95$0 = N[(N[(A + C), $MachinePrecision] + N[Sqrt[B$95$m ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[B$95$m, 1.4e-105], N[(N[Sqrt[t$95$0], $MachinePrecision] * N[(N[Sqrt[N[(N[(2.0 * F), $MachinePrecision] * N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(N[(N[(4.0 * A), $MachinePrecision] * C), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[B$95$m, 2.8e+146], N[(N[Sqrt[N[(t$95$0 / N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(A * N[(C * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(0.0 - N[Sqrt[N[(2.0 * F), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
t_0 := \left(A + C\right) + \mathsf{hypot}\left(B\_m, A - C\right)\\
\mathbf{if}\;B\_m \leq 1.4 \cdot 10^{-105}:\\
\;\;\;\;\sqrt{t\_0} \cdot \frac{\sqrt{\left(2 \cdot F\right) \cdot \left(B\_m \cdot B\_m + -4 \cdot \left(A \cdot C\right)\right)}}{\left(4 \cdot A\right) \cdot C - B\_m \cdot B\_m}\\
\mathbf{elif}\;B\_m \leq 2.8 \cdot 10^{+146}:\\
\;\;\;\;\sqrt{\frac{t\_0}{B\_m \cdot B\_m + A \cdot \left(C \cdot -4\right)}} \cdot \left(0 - \sqrt{2 \cdot F}\right)\\
\mathbf{else}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\end{array}
\end{array}
if B < 1.4e-105Initial program 22.0%
Applied egg-rr41.1%
if 1.4e-105 < B < 2.8000000000000001e146Initial program 17.6%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified39.8%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr50.5%
if 2.8000000000000001e146 < B Initial program 0.2%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified5.2%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr8.9%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6474.6%
Simplified74.6%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6474.6%
Applied egg-rr74.6%
Final simplification48.0%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(if (<= B_m 1.4e-139)
(/ (* (sqrt (* F (* C -16.0))) (fabs A)) (- (* 4.0 (* A C)) (* B_m B_m)))
(if (<= B_m 2.8e+146)
(*
(sqrt
(/ (+ (+ A C) (hypot B_m (- A C))) (+ (* B_m B_m) (* A (* C -4.0)))))
(- 0.0 (sqrt (* 2.0 F))))
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m))))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 1.4e-139) {
tmp = (sqrt((F * (C * -16.0))) * fabs(A)) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 2.8e+146) {
tmp = sqrt((((A + C) + hypot(B_m, (A - C))) / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - sqrt((2.0 * F)));
} else {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
}
return tmp;
}
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 1.4e-139) {
tmp = (Math.sqrt((F * (C * -16.0))) * Math.abs(A)) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 2.8e+146) {
tmp = Math.sqrt((((A + C) + Math.hypot(B_m, (A - C))) / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - Math.sqrt((2.0 * F)));
} else {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if B_m <= 1.4e-139: tmp = (math.sqrt((F * (C * -16.0))) * math.fabs(A)) / ((4.0 * (A * C)) - (B_m * B_m)) elif B_m <= 2.8e+146: tmp = math.sqrt((((A + C) + math.hypot(B_m, (A - C))) / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - math.sqrt((2.0 * F))) else: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (B_m <= 1.4e-139) tmp = Float64(Float64(sqrt(Float64(F * Float64(C * -16.0))) * abs(A)) / Float64(Float64(4.0 * Float64(A * C)) - Float64(B_m * B_m))); elseif (B_m <= 2.8e+146) tmp = Float64(sqrt(Float64(Float64(Float64(A + C) + hypot(B_m, Float64(A - C))) / Float64(Float64(B_m * B_m) + Float64(A * Float64(C * -4.0))))) * Float64(0.0 - sqrt(Float64(2.0 * F)))); else tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (B_m <= 1.4e-139) tmp = (sqrt((F * (C * -16.0))) * abs(A)) / ((4.0 * (A * C)) - (B_m * B_m)); elseif (B_m <= 2.8e+146) tmp = sqrt((((A + C) + hypot(B_m, (A - C))) / ((B_m * B_m) + (A * (C * -4.0))))) * (0.0 - sqrt((2.0 * F))); else tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[B$95$m, 1.4e-139], N[(N[(N[Sqrt[N[(F * N[(C * -16.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Abs[A], $MachinePrecision]), $MachinePrecision] / N[(N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[B$95$m, 2.8e+146], N[(N[Sqrt[N[(N[(N[(A + C), $MachinePrecision] + N[Sqrt[B$95$m ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision] / N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(A * N[(C * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(0.0 - N[Sqrt[N[(2.0 * F), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;B\_m \leq 1.4 \cdot 10^{-139}:\\
\;\;\;\;\frac{\sqrt{F \cdot \left(C \cdot -16\right)} \cdot \left|A\right|}{4 \cdot \left(A \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{elif}\;B\_m \leq 2.8 \cdot 10^{+146}:\\
\;\;\;\;\sqrt{\frac{\left(A + C\right) + \mathsf{hypot}\left(B\_m, A - C\right)}{B\_m \cdot B\_m + A \cdot \left(C \cdot -4\right)}} \cdot \left(0 - \sqrt{2 \cdot F}\right)\\
\mathbf{else}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\end{array}
\end{array}
if B < 1.3999999999999999e-139Initial program 21.9%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified28.0%
Taylor expanded in B around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6410.1%
Simplified10.1%
*-commutativeN/A
associate-*r*N/A
sqrt-prodN/A
pow1/2N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
pow1/2N/A
rem-sqrt-squareN/A
fabs-lowering-fabs.f6417.9%
Applied egg-rr17.9%
if 1.3999999999999999e-139 < B < 2.8000000000000001e146Initial program 18.8%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified42.5%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr51.4%
if 2.8000000000000001e146 < B Initial program 0.2%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified5.2%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr8.9%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6474.6%
Simplified74.6%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6474.6%
Applied egg-rr74.6%
Final simplification33.9%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(if (<= B_m 2.8e-136)
(/ (* (sqrt (* F (* C -16.0))) (fabs A)) (- (* 4.0 (* A C)) (* B_m B_m)))
(if (<= B_m 9e+58)
(-
0.0
(sqrt
(/
(* (* 2.0 F) (+ (+ A C) (hypot B_m (- A C))))
(+ (* B_m B_m) (* A (* C -4.0))))))
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m))))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 2.8e-136) {
tmp = (sqrt((F * (C * -16.0))) * fabs(A)) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 9e+58) {
tmp = 0.0 - sqrt((((2.0 * F) * ((A + C) + hypot(B_m, (A - C)))) / ((B_m * B_m) + (A * (C * -4.0)))));
} else {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
}
return tmp;
}
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 2.8e-136) {
tmp = (Math.sqrt((F * (C * -16.0))) * Math.abs(A)) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 9e+58) {
tmp = 0.0 - Math.sqrt((((2.0 * F) * ((A + C) + Math.hypot(B_m, (A - C)))) / ((B_m * B_m) + (A * (C * -4.0)))));
} else {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if B_m <= 2.8e-136: tmp = (math.sqrt((F * (C * -16.0))) * math.fabs(A)) / ((4.0 * (A * C)) - (B_m * B_m)) elif B_m <= 9e+58: tmp = 0.0 - math.sqrt((((2.0 * F) * ((A + C) + math.hypot(B_m, (A - C)))) / ((B_m * B_m) + (A * (C * -4.0))))) else: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (B_m <= 2.8e-136) tmp = Float64(Float64(sqrt(Float64(F * Float64(C * -16.0))) * abs(A)) / Float64(Float64(4.0 * Float64(A * C)) - Float64(B_m * B_m))); elseif (B_m <= 9e+58) tmp = Float64(0.0 - sqrt(Float64(Float64(Float64(2.0 * F) * Float64(Float64(A + C) + hypot(B_m, Float64(A - C)))) / Float64(Float64(B_m * B_m) + Float64(A * Float64(C * -4.0)))))); else tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (B_m <= 2.8e-136) tmp = (sqrt((F * (C * -16.0))) * abs(A)) / ((4.0 * (A * C)) - (B_m * B_m)); elseif (B_m <= 9e+58) tmp = 0.0 - sqrt((((2.0 * F) * ((A + C) + hypot(B_m, (A - C)))) / ((B_m * B_m) + (A * (C * -4.0))))); else tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[B$95$m, 2.8e-136], N[(N[(N[Sqrt[N[(F * N[(C * -16.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Abs[A], $MachinePrecision]), $MachinePrecision] / N[(N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[B$95$m, 9e+58], N[(0.0 - N[Sqrt[N[(N[(N[(2.0 * F), $MachinePrecision] * N[(N[(A + C), $MachinePrecision] + N[Sqrt[B$95$m ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(A * N[(C * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;B\_m \leq 2.8 \cdot 10^{-136}:\\
\;\;\;\;\frac{\sqrt{F \cdot \left(C \cdot -16\right)} \cdot \left|A\right|}{4 \cdot \left(A \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{elif}\;B\_m \leq 9 \cdot 10^{+58}:\\
\;\;\;\;0 - \sqrt{\frac{\left(2 \cdot F\right) \cdot \left(\left(A + C\right) + \mathsf{hypot}\left(B\_m, A - C\right)\right)}{B\_m \cdot B\_m + A \cdot \left(C \cdot -4\right)}}\\
\mathbf{else}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\end{array}
\end{array}
if B < 2.8000000000000001e-136Initial program 21.9%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified28.0%
Taylor expanded in B around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6410.1%
Simplified10.1%
*-commutativeN/A
associate-*r*N/A
sqrt-prodN/A
pow1/2N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
pow1/2N/A
rem-sqrt-squareN/A
fabs-lowering-fabs.f6417.9%
Applied egg-rr17.9%
if 2.8000000000000001e-136 < B < 8.9999999999999996e58Initial program 17.3%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified40.3%
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
*-commutativeN/A
associate-*r/N/A
associate-+r+N/A
associate-*l*N/A
/-lowering-/.f64N/A
Applied egg-rr37.2%
if 8.9999999999999996e58 < B Initial program 6.1%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified16.5%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr25.2%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6472.3%
Simplified72.3%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6472.2%
Applied egg-rr72.2%
Final simplification32.4%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(let* ((t_0 (- 0.0 (sqrt (* 2.0 F)))))
(if (<= C -4.2e-60)
(* (sqrt (/ -0.5 C)) t_0)
(if (<= C 7.5e+99)
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m)))
(* (sqrt (/ (* 2.0 C) (+ (* B_m B_m) (* A (* C -4.0))))) t_0)))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double t_0 = 0.0 - sqrt((2.0 * F));
double tmp;
if (C <= -4.2e-60) {
tmp = sqrt((-0.5 / C)) * t_0;
} else if (C <= 7.5e+99) {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
} else {
tmp = sqrt(((2.0 * C) / ((B_m * B_m) + (A * (C * -4.0))))) * t_0;
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: tmp
t_0 = 0.0d0 - sqrt((2.0d0 * f))
if (c <= (-4.2d-60)) then
tmp = sqrt(((-0.5d0) / c)) * t_0
else if (c <= 7.5d+99) then
tmp = 0.0d0 - (((2.0d0 * f) ** 0.5d0) / sqrt(b_m))
else
tmp = sqrt(((2.0d0 * c) / ((b_m * b_m) + (a * (c * (-4.0d0)))))) * t_0
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double t_0 = 0.0 - Math.sqrt((2.0 * F));
double tmp;
if (C <= -4.2e-60) {
tmp = Math.sqrt((-0.5 / C)) * t_0;
} else if (C <= 7.5e+99) {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
} else {
tmp = Math.sqrt(((2.0 * C) / ((B_m * B_m) + (A * (C * -4.0))))) * t_0;
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): t_0 = 0.0 - math.sqrt((2.0 * F)) tmp = 0 if C <= -4.2e-60: tmp = math.sqrt((-0.5 / C)) * t_0 elif C <= 7.5e+99: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) else: tmp = math.sqrt(((2.0 * C) / ((B_m * B_m) + (A * (C * -4.0))))) * t_0 return tmp
B_m = abs(B) function code(A, B_m, C, F) t_0 = Float64(0.0 - sqrt(Float64(2.0 * F))) tmp = 0.0 if (C <= -4.2e-60) tmp = Float64(sqrt(Float64(-0.5 / C)) * t_0); elseif (C <= 7.5e+99) tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); else tmp = Float64(sqrt(Float64(Float64(2.0 * C) / Float64(Float64(B_m * B_m) + Float64(A * Float64(C * -4.0))))) * t_0); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) t_0 = 0.0 - sqrt((2.0 * F)); tmp = 0.0; if (C <= -4.2e-60) tmp = sqrt((-0.5 / C)) * t_0; elseif (C <= 7.5e+99) tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); else tmp = sqrt(((2.0 * C) / ((B_m * B_m) + (A * (C * -4.0))))) * t_0; end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision]
code[A_, B$95$m_, C_, F_] := Block[{t$95$0 = N[(0.0 - N[Sqrt[N[(2.0 * F), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[C, -4.2e-60], N[(N[Sqrt[N[(-0.5 / C), $MachinePrecision]], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[C, 7.5e+99], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Sqrt[N[(N[(2.0 * C), $MachinePrecision] / N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(A * N[(C * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * t$95$0), $MachinePrecision]]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
t_0 := 0 - \sqrt{2 \cdot F}\\
\mathbf{if}\;C \leq -4.2 \cdot 10^{-60}:\\
\;\;\;\;\sqrt{\frac{-0.5}{C}} \cdot t\_0\\
\mathbf{elif}\;C \leq 7.5 \cdot 10^{+99}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{2 \cdot C}{B\_m \cdot B\_m + A \cdot \left(C \cdot -4\right)}} \cdot t\_0\\
\end{array}
\end{array}
if C < -4.19999999999999982e-60Initial program 6.5%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified19.1%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr20.7%
Taylor expanded in A around inf
/-lowering-/.f6446.7%
Simplified46.7%
if -4.19999999999999982e-60 < C < 7.49999999999999963e99Initial program 28.1%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified31.0%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr41.4%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6428.3%
Simplified28.3%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6428.3%
Applied egg-rr28.3%
if 7.49999999999999963e99 < C Initial program 7.8%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified26.8%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr40.2%
Taylor expanded in A around -inf
*-lowering-*.f6442.1%
Simplified42.1%
Final simplification36.4%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(let* ((t_0 (- 0.0 (sqrt (* 2.0 F)))))
(if (<= C -4.2e-60)
(* (sqrt (/ -0.5 C)) t_0)
(if (<= C 1.3e+104)
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m)))
(if (<= C 1.45e+228)
(/
-1.0
(/ (- (* B_m B_m) (* A (* 4.0 C))) (* C (sqrt (* F (* A -16.0))))))
(* (sqrt (/ -0.5 A)) t_0))))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double t_0 = 0.0 - sqrt((2.0 * F));
double tmp;
if (C <= -4.2e-60) {
tmp = sqrt((-0.5 / C)) * t_0;
} else if (C <= 1.3e+104) {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
} else if (C <= 1.45e+228) {
tmp = -1.0 / (((B_m * B_m) - (A * (4.0 * C))) / (C * sqrt((F * (A * -16.0)))));
} else {
tmp = sqrt((-0.5 / A)) * t_0;
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: t_0
real(8) :: tmp
t_0 = 0.0d0 - sqrt((2.0d0 * f))
if (c <= (-4.2d-60)) then
tmp = sqrt(((-0.5d0) / c)) * t_0
else if (c <= 1.3d+104) then
tmp = 0.0d0 - (((2.0d0 * f) ** 0.5d0) / sqrt(b_m))
else if (c <= 1.45d+228) then
tmp = (-1.0d0) / (((b_m * b_m) - (a * (4.0d0 * c))) / (c * sqrt((f * (a * (-16.0d0))))))
else
tmp = sqrt(((-0.5d0) / a)) * t_0
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double t_0 = 0.0 - Math.sqrt((2.0 * F));
double tmp;
if (C <= -4.2e-60) {
tmp = Math.sqrt((-0.5 / C)) * t_0;
} else if (C <= 1.3e+104) {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
} else if (C <= 1.45e+228) {
tmp = -1.0 / (((B_m * B_m) - (A * (4.0 * C))) / (C * Math.sqrt((F * (A * -16.0)))));
} else {
tmp = Math.sqrt((-0.5 / A)) * t_0;
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): t_0 = 0.0 - math.sqrt((2.0 * F)) tmp = 0 if C <= -4.2e-60: tmp = math.sqrt((-0.5 / C)) * t_0 elif C <= 1.3e+104: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) elif C <= 1.45e+228: tmp = -1.0 / (((B_m * B_m) - (A * (4.0 * C))) / (C * math.sqrt((F * (A * -16.0))))) else: tmp = math.sqrt((-0.5 / A)) * t_0 return tmp
B_m = abs(B) function code(A, B_m, C, F) t_0 = Float64(0.0 - sqrt(Float64(2.0 * F))) tmp = 0.0 if (C <= -4.2e-60) tmp = Float64(sqrt(Float64(-0.5 / C)) * t_0); elseif (C <= 1.3e+104) tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); elseif (C <= 1.45e+228) tmp = Float64(-1.0 / Float64(Float64(Float64(B_m * B_m) - Float64(A * Float64(4.0 * C))) / Float64(C * sqrt(Float64(F * Float64(A * -16.0)))))); else tmp = Float64(sqrt(Float64(-0.5 / A)) * t_0); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) t_0 = 0.0 - sqrt((2.0 * F)); tmp = 0.0; if (C <= -4.2e-60) tmp = sqrt((-0.5 / C)) * t_0; elseif (C <= 1.3e+104) tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); elseif (C <= 1.45e+228) tmp = -1.0 / (((B_m * B_m) - (A * (4.0 * C))) / (C * sqrt((F * (A * -16.0))))); else tmp = sqrt((-0.5 / A)) * t_0; end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision]
code[A_, B$95$m_, C_, F_] := Block[{t$95$0 = N[(0.0 - N[Sqrt[N[(2.0 * F), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[C, -4.2e-60], N[(N[Sqrt[N[(-0.5 / C), $MachinePrecision]], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[C, 1.3e+104], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[C, 1.45e+228], N[(-1.0 / N[(N[(N[(B$95$m * B$95$m), $MachinePrecision] - N[(A * N[(4.0 * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(C * N[Sqrt[N[(F * N[(A * -16.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Sqrt[N[(-0.5 / A), $MachinePrecision]], $MachinePrecision] * t$95$0), $MachinePrecision]]]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
t_0 := 0 - \sqrt{2 \cdot F}\\
\mathbf{if}\;C \leq -4.2 \cdot 10^{-60}:\\
\;\;\;\;\sqrt{\frac{-0.5}{C}} \cdot t\_0\\
\mathbf{elif}\;C \leq 1.3 \cdot 10^{+104}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\mathbf{elif}\;C \leq 1.45 \cdot 10^{+228}:\\
\;\;\;\;\frac{-1}{\frac{B\_m \cdot B\_m - A \cdot \left(4 \cdot C\right)}{C \cdot \sqrt{F \cdot \left(A \cdot -16\right)}}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{-0.5}{A}} \cdot t\_0\\
\end{array}
\end{array}
if C < -4.19999999999999982e-60Initial program 6.5%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified19.1%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr20.7%
Taylor expanded in A around inf
/-lowering-/.f6446.7%
Simplified46.7%
if -4.19999999999999982e-60 < C < 1.3e104Initial program 28.6%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified31.5%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr41.9%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6428.2%
Simplified28.2%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6428.2%
Applied egg-rr28.2%
if 1.3e104 < C < 1.45000000000000001e228Initial program 9.9%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified32.9%
Taylor expanded in A around -inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f646.3%
Simplified6.3%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
Applied egg-rr2.3%
associate-*r*N/A
associate-*r*N/A
sqrt-prodN/A
pow2N/A
sqrt-pow1N/A
metadata-evalN/A
metadata-evalN/A
pow-flipN/A
inv-powN/A
remove-double-divN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6449.1%
Applied egg-rr49.1%
if 1.45000000000000001e228 < C Initial program 2.1%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified25.7%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr45.0%
Taylor expanded in A around -inf
/-lowering-/.f6448.0%
Simplified48.0%
Final simplification37.5%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(if (<= C -6.5e+63)
(* (sqrt (* -0.5 (/ F C))) (- 0.0 (sqrt 2.0)))
(if (<= C 4.5e+104)
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m)))
(if (<= C 3.5e+228)
(/
-1.0
(/ (- (* B_m B_m) (* A (* 4.0 C))) (* C (sqrt (* F (* A -16.0))))))
(- 0.0 (* (sqrt 2.0) (sqrt (* -0.5 (/ F A)))))))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (C <= -6.5e+63) {
tmp = sqrt((-0.5 * (F / C))) * (0.0 - sqrt(2.0));
} else if (C <= 4.5e+104) {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
} else if (C <= 3.5e+228) {
tmp = -1.0 / (((B_m * B_m) - (A * (4.0 * C))) / (C * sqrt((F * (A * -16.0)))));
} else {
tmp = 0.0 - (sqrt(2.0) * sqrt((-0.5 * (F / A))));
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (c <= (-6.5d+63)) then
tmp = sqrt(((-0.5d0) * (f / c))) * (0.0d0 - sqrt(2.0d0))
else if (c <= 4.5d+104) then
tmp = 0.0d0 - (((2.0d0 * f) ** 0.5d0) / sqrt(b_m))
else if (c <= 3.5d+228) then
tmp = (-1.0d0) / (((b_m * b_m) - (a * (4.0d0 * c))) / (c * sqrt((f * (a * (-16.0d0))))))
else
tmp = 0.0d0 - (sqrt(2.0d0) * sqrt(((-0.5d0) * (f / a))))
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (C <= -6.5e+63) {
tmp = Math.sqrt((-0.5 * (F / C))) * (0.0 - Math.sqrt(2.0));
} else if (C <= 4.5e+104) {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
} else if (C <= 3.5e+228) {
tmp = -1.0 / (((B_m * B_m) - (A * (4.0 * C))) / (C * Math.sqrt((F * (A * -16.0)))));
} else {
tmp = 0.0 - (Math.sqrt(2.0) * Math.sqrt((-0.5 * (F / A))));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if C <= -6.5e+63: tmp = math.sqrt((-0.5 * (F / C))) * (0.0 - math.sqrt(2.0)) elif C <= 4.5e+104: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) elif C <= 3.5e+228: tmp = -1.0 / (((B_m * B_m) - (A * (4.0 * C))) / (C * math.sqrt((F * (A * -16.0))))) else: tmp = 0.0 - (math.sqrt(2.0) * math.sqrt((-0.5 * (F / A)))) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (C <= -6.5e+63) tmp = Float64(sqrt(Float64(-0.5 * Float64(F / C))) * Float64(0.0 - sqrt(2.0))); elseif (C <= 4.5e+104) tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); elseif (C <= 3.5e+228) tmp = Float64(-1.0 / Float64(Float64(Float64(B_m * B_m) - Float64(A * Float64(4.0 * C))) / Float64(C * sqrt(Float64(F * Float64(A * -16.0)))))); else tmp = Float64(0.0 - Float64(sqrt(2.0) * sqrt(Float64(-0.5 * Float64(F / A))))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (C <= -6.5e+63) tmp = sqrt((-0.5 * (F / C))) * (0.0 - sqrt(2.0)); elseif (C <= 4.5e+104) tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); elseif (C <= 3.5e+228) tmp = -1.0 / (((B_m * B_m) - (A * (4.0 * C))) / (C * sqrt((F * (A * -16.0))))); else tmp = 0.0 - (sqrt(2.0) * sqrt((-0.5 * (F / A)))); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[C, -6.5e+63], N[(N[Sqrt[N[(-0.5 * N[(F / C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(0.0 - N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[C, 4.5e+104], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[C, 3.5e+228], N[(-1.0 / N[(N[(N[(B$95$m * B$95$m), $MachinePrecision] - N[(A * N[(4.0 * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(C * N[Sqrt[N[(F * N[(A * -16.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(N[Sqrt[2.0], $MachinePrecision] * N[Sqrt[N[(-0.5 * N[(F / A), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;C \leq -6.5 \cdot 10^{+63}:\\
\;\;\;\;\sqrt{-0.5 \cdot \frac{F}{C}} \cdot \left(0 - \sqrt{2}\right)\\
\mathbf{elif}\;C \leq 4.5 \cdot 10^{+104}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\mathbf{elif}\;C \leq 3.5 \cdot 10^{+228}:\\
\;\;\;\;\frac{-1}{\frac{B\_m \cdot B\_m - A \cdot \left(4 \cdot C\right)}{C \cdot \sqrt{F \cdot \left(A \cdot -16\right)}}}\\
\mathbf{else}:\\
\;\;\;\;0 - \sqrt{2} \cdot \sqrt{-0.5 \cdot \frac{F}{A}}\\
\end{array}
\end{array}
if C < -6.49999999999999992e63Initial program 0.9%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified8.9%
Taylor expanded in A around inf
*-lowering-*.f64N/A
/-lowering-/.f6434.8%
Simplified34.8%
if -6.49999999999999992e63 < C < 4.4999999999999998e104Initial program 27.2%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified33.0%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr43.0%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6427.0%
Simplified27.0%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6427.0%
Applied egg-rr27.0%
if 4.4999999999999998e104 < C < 3.5000000000000002e228Initial program 9.9%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified32.9%
Taylor expanded in A around -inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f646.3%
Simplified6.3%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
Applied egg-rr2.3%
associate-*r*N/A
associate-*r*N/A
sqrt-prodN/A
pow2N/A
sqrt-pow1N/A
metadata-evalN/A
metadata-evalN/A
pow-flipN/A
inv-powN/A
remove-double-divN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6449.1%
Applied egg-rr49.1%
if 3.5000000000000002e228 < C Initial program 2.1%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified25.7%
Taylor expanded in A around -inf
*-lowering-*.f64N/A
/-lowering-/.f6444.9%
Simplified44.9%
Final simplification32.5%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(if (<= B_m 3.4e-137)
(/ (* (sqrt (* F (* C -16.0))) (fabs A)) (- (* 4.0 (* A C)) (* B_m B_m)))
(if (<= B_m 95000000.0)
(* (sqrt (/ -0.5 A)) (- 0.0 (sqrt (* 2.0 F))))
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m))))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 3.4e-137) {
tmp = (sqrt((F * (C * -16.0))) * fabs(A)) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 95000000.0) {
tmp = sqrt((-0.5 / A)) * (0.0 - sqrt((2.0 * F)));
} else {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (b_m <= 3.4d-137) then
tmp = (sqrt((f * (c * (-16.0d0)))) * abs(a)) / ((4.0d0 * (a * c)) - (b_m * b_m))
else if (b_m <= 95000000.0d0) then
tmp = sqrt(((-0.5d0) / a)) * (0.0d0 - sqrt((2.0d0 * f)))
else
tmp = 0.0d0 - (((2.0d0 * f) ** 0.5d0) / sqrt(b_m))
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 3.4e-137) {
tmp = (Math.sqrt((F * (C * -16.0))) * Math.abs(A)) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 95000000.0) {
tmp = Math.sqrt((-0.5 / A)) * (0.0 - Math.sqrt((2.0 * F)));
} else {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if B_m <= 3.4e-137: tmp = (math.sqrt((F * (C * -16.0))) * math.fabs(A)) / ((4.0 * (A * C)) - (B_m * B_m)) elif B_m <= 95000000.0: tmp = math.sqrt((-0.5 / A)) * (0.0 - math.sqrt((2.0 * F))) else: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (B_m <= 3.4e-137) tmp = Float64(Float64(sqrt(Float64(F * Float64(C * -16.0))) * abs(A)) / Float64(Float64(4.0 * Float64(A * C)) - Float64(B_m * B_m))); elseif (B_m <= 95000000.0) tmp = Float64(sqrt(Float64(-0.5 / A)) * Float64(0.0 - sqrt(Float64(2.0 * F)))); else tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (B_m <= 3.4e-137) tmp = (sqrt((F * (C * -16.0))) * abs(A)) / ((4.0 * (A * C)) - (B_m * B_m)); elseif (B_m <= 95000000.0) tmp = sqrt((-0.5 / A)) * (0.0 - sqrt((2.0 * F))); else tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[B$95$m, 3.4e-137], N[(N[(N[Sqrt[N[(F * N[(C * -16.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Abs[A], $MachinePrecision]), $MachinePrecision] / N[(N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[B$95$m, 95000000.0], N[(N[Sqrt[N[(-0.5 / A), $MachinePrecision]], $MachinePrecision] * N[(0.0 - N[Sqrt[N[(2.0 * F), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;B\_m \leq 3.4 \cdot 10^{-137}:\\
\;\;\;\;\frac{\sqrt{F \cdot \left(C \cdot -16\right)} \cdot \left|A\right|}{4 \cdot \left(A \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{elif}\;B\_m \leq 95000000:\\
\;\;\;\;\sqrt{\frac{-0.5}{A}} \cdot \left(0 - \sqrt{2 \cdot F}\right)\\
\mathbf{else}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\end{array}
\end{array}
if B < 3.40000000000000014e-137Initial program 21.9%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified28.0%
Taylor expanded in B around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6410.1%
Simplified10.1%
*-commutativeN/A
associate-*r*N/A
sqrt-prodN/A
pow1/2N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
pow1/2N/A
rem-sqrt-squareN/A
fabs-lowering-fabs.f6417.9%
Applied egg-rr17.9%
if 3.40000000000000014e-137 < B < 9.5e7Initial program 17.8%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified36.4%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr42.4%
Taylor expanded in A around -inf
/-lowering-/.f6438.6%
Simplified38.6%
if 9.5e7 < B Initial program 7.3%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified21.5%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr28.6%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6463.8%
Simplified63.8%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6463.8%
Applied egg-rr63.8%
Final simplification31.6%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(if (<= B_m 1.32e-276)
(/
(sqrt (* (+ (* B_m B_m) (* -4.0 (* A C))) (* (* 2.0 F) (* 2.0 A))))
(- (* 4.0 (* A C)) (* B_m B_m)))
(if (<= B_m 64000000.0)
(* (sqrt (/ -0.5 A)) (- 0.0 (sqrt (* 2.0 F))))
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m))))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 1.32e-276) {
tmp = sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * (2.0 * A)))) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 64000000.0) {
tmp = sqrt((-0.5 / A)) * (0.0 - sqrt((2.0 * F)));
} else {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (b_m <= 1.32d-276) then
tmp = sqrt((((b_m * b_m) + ((-4.0d0) * (a * c))) * ((2.0d0 * f) * (2.0d0 * a)))) / ((4.0d0 * (a * c)) - (b_m * b_m))
else if (b_m <= 64000000.0d0) then
tmp = sqrt(((-0.5d0) / a)) * (0.0d0 - sqrt((2.0d0 * f)))
else
tmp = 0.0d0 - (((2.0d0 * f) ** 0.5d0) / sqrt(b_m))
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 1.32e-276) {
tmp = Math.sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * (2.0 * A)))) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 64000000.0) {
tmp = Math.sqrt((-0.5 / A)) * (0.0 - Math.sqrt((2.0 * F)));
} else {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if B_m <= 1.32e-276: tmp = math.sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * (2.0 * A)))) / ((4.0 * (A * C)) - (B_m * B_m)) elif B_m <= 64000000.0: tmp = math.sqrt((-0.5 / A)) * (0.0 - math.sqrt((2.0 * F))) else: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (B_m <= 1.32e-276) tmp = Float64(sqrt(Float64(Float64(Float64(B_m * B_m) + Float64(-4.0 * Float64(A * C))) * Float64(Float64(2.0 * F) * Float64(2.0 * A)))) / Float64(Float64(4.0 * Float64(A * C)) - Float64(B_m * B_m))); elseif (B_m <= 64000000.0) tmp = Float64(sqrt(Float64(-0.5 / A)) * Float64(0.0 - sqrt(Float64(2.0 * F)))); else tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (B_m <= 1.32e-276) tmp = sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * (2.0 * A)))) / ((4.0 * (A * C)) - (B_m * B_m)); elseif (B_m <= 64000000.0) tmp = sqrt((-0.5 / A)) * (0.0 - sqrt((2.0 * F))); else tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[B$95$m, 1.32e-276], N[(N[Sqrt[N[(N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * F), $MachinePrecision] * N[(2.0 * A), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[B$95$m, 64000000.0], N[(N[Sqrt[N[(-0.5 / A), $MachinePrecision]], $MachinePrecision] * N[(0.0 - N[Sqrt[N[(2.0 * F), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;B\_m \leq 1.32 \cdot 10^{-276}:\\
\;\;\;\;\frac{\sqrt{\left(B\_m \cdot B\_m + -4 \cdot \left(A \cdot C\right)\right) \cdot \left(\left(2 \cdot F\right) \cdot \left(2 \cdot A\right)\right)}}{4 \cdot \left(A \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{elif}\;B\_m \leq 64000000:\\
\;\;\;\;\sqrt{\frac{-0.5}{A}} \cdot \left(0 - \sqrt{2 \cdot F}\right)\\
\mathbf{else}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\end{array}
\end{array}
if B < 1.31999999999999992e-276Initial program 21.2%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified27.8%
associate-*l*N/A
associate-*r*N/A
pow2N/A
*-lowering-*.f64N/A
pow2N/A
associate-*r*N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr27.2%
Taylor expanded in A around inf
*-lowering-*.f6412.2%
Simplified12.2%
if 1.31999999999999992e-276 < B < 6.4e7Initial program 21.1%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified27.5%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr32.7%
Taylor expanded in A around -inf
/-lowering-/.f6432.5%
Simplified32.5%
if 6.4e7 < B Initial program 7.3%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified21.5%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr28.6%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6463.8%
Simplified63.8%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6463.8%
Applied egg-rr63.8%
Final simplification29.5%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(if (<= B_m 7.6e-137)
(/
(sqrt
(*
(+ (* B_m B_m) (* -4.0 (* A C)))
(* (* 2.0 F) (+ (* 2.0 A) (* -0.5 (/ (* B_m B_m) C))))))
(- (* 4.0 (* A C)) (* B_m B_m)))
(if (<= B_m 6.5e-91)
(- 0.0 (* (sqrt 2.0) (sqrt (* -0.5 (/ F A)))))
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m))))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 7.6e-137) {
tmp = sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 6.5e-91) {
tmp = 0.0 - (sqrt(2.0) * sqrt((-0.5 * (F / A))));
} else {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (b_m <= 7.6d-137) then
tmp = sqrt((((b_m * b_m) + ((-4.0d0) * (a * c))) * ((2.0d0 * f) * ((2.0d0 * a) + ((-0.5d0) * ((b_m * b_m) / c)))))) / ((4.0d0 * (a * c)) - (b_m * b_m))
else if (b_m <= 6.5d-91) then
tmp = 0.0d0 - (sqrt(2.0d0) * sqrt(((-0.5d0) * (f / a))))
else
tmp = 0.0d0 - (((2.0d0 * f) ** 0.5d0) / sqrt(b_m))
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 7.6e-137) {
tmp = Math.sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 6.5e-91) {
tmp = 0.0 - (Math.sqrt(2.0) * Math.sqrt((-0.5 * (F / A))));
} else {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if B_m <= 7.6e-137: tmp = math.sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m)) elif B_m <= 6.5e-91: tmp = 0.0 - (math.sqrt(2.0) * math.sqrt((-0.5 * (F / A)))) else: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (B_m <= 7.6e-137) tmp = Float64(sqrt(Float64(Float64(Float64(B_m * B_m) + Float64(-4.0 * Float64(A * C))) * Float64(Float64(2.0 * F) * Float64(Float64(2.0 * A) + Float64(-0.5 * Float64(Float64(B_m * B_m) / C)))))) / Float64(Float64(4.0 * Float64(A * C)) - Float64(B_m * B_m))); elseif (B_m <= 6.5e-91) tmp = Float64(0.0 - Float64(sqrt(2.0) * sqrt(Float64(-0.5 * Float64(F / A))))); else tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (B_m <= 7.6e-137) tmp = sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m)); elseif (B_m <= 6.5e-91) tmp = 0.0 - (sqrt(2.0) * sqrt((-0.5 * (F / A)))); else tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[B$95$m, 7.6e-137], N[(N[Sqrt[N[(N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * F), $MachinePrecision] * N[(N[(2.0 * A), $MachinePrecision] + N[(-0.5 * N[(N[(B$95$m * B$95$m), $MachinePrecision] / C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[B$95$m, 6.5e-91], N[(0.0 - N[(N[Sqrt[2.0], $MachinePrecision] * N[Sqrt[N[(-0.5 * N[(F / A), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;B\_m \leq 7.6 \cdot 10^{-137}:\\
\;\;\;\;\frac{\sqrt{\left(B\_m \cdot B\_m + -4 \cdot \left(A \cdot C\right)\right) \cdot \left(\left(2 \cdot F\right) \cdot \left(2 \cdot A + -0.5 \cdot \frac{B\_m \cdot B\_m}{C}\right)\right)}}{4 \cdot \left(A \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{elif}\;B\_m \leq 6.5 \cdot 10^{-91}:\\
\;\;\;\;0 - \sqrt{2} \cdot \sqrt{-0.5 \cdot \frac{F}{A}}\\
\mathbf{else}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\end{array}
\end{array}
if B < 7.59999999999999997e-137Initial program 21.9%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified28.0%
associate-*l*N/A
associate-*r*N/A
pow2N/A
*-lowering-*.f64N/A
pow2N/A
associate-*r*N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr27.5%
Taylor expanded in C around -inf
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6413.9%
Simplified13.9%
if 7.59999999999999997e-137 < B < 6.5000000000000001e-91Initial program 19.2%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified42.9%
Taylor expanded in A around -inf
*-lowering-*.f64N/A
/-lowering-/.f6435.1%
Simplified35.1%
if 6.5000000000000001e-91 < B Initial program 9.6%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified24.2%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr30.7%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6453.0%
Simplified53.0%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6453.0%
Applied egg-rr53.0%
Final simplification27.4%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(if (<= B_m 3.8e-133)
(/
(sqrt
(*
(+ (* B_m B_m) (* -4.0 (* A C)))
(* (* 2.0 F) (+ (* 2.0 A) (* -0.5 (/ (* B_m B_m) C))))))
(- (* 4.0 (* A C)) (* B_m B_m)))
(- 0.0 (/ (pow (* 2.0 F) 0.5) (sqrt B_m)))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 3.8e-133) {
tmp = sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m));
} else {
tmp = 0.0 - (pow((2.0 * F), 0.5) / sqrt(B_m));
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (b_m <= 3.8d-133) then
tmp = sqrt((((b_m * b_m) + ((-4.0d0) * (a * c))) * ((2.0d0 * f) * ((2.0d0 * a) + ((-0.5d0) * ((b_m * b_m) / c)))))) / ((4.0d0 * (a * c)) - (b_m * b_m))
else
tmp = 0.0d0 - (((2.0d0 * f) ** 0.5d0) / sqrt(b_m))
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 3.8e-133) {
tmp = Math.sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m));
} else {
tmp = 0.0 - (Math.pow((2.0 * F), 0.5) / Math.sqrt(B_m));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if B_m <= 3.8e-133: tmp = math.sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m)) else: tmp = 0.0 - (math.pow((2.0 * F), 0.5) / math.sqrt(B_m)) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (B_m <= 3.8e-133) tmp = Float64(sqrt(Float64(Float64(Float64(B_m * B_m) + Float64(-4.0 * Float64(A * C))) * Float64(Float64(2.0 * F) * Float64(Float64(2.0 * A) + Float64(-0.5 * Float64(Float64(B_m * B_m) / C)))))) / Float64(Float64(4.0 * Float64(A * C)) - Float64(B_m * B_m))); else tmp = Float64(0.0 - Float64((Float64(2.0 * F) ^ 0.5) / sqrt(B_m))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (B_m <= 3.8e-133) tmp = sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m)); else tmp = 0.0 - (((2.0 * F) ^ 0.5) / sqrt(B_m)); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[B$95$m, 3.8e-133], N[(N[Sqrt[N[(N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * F), $MachinePrecision] * N[(N[(2.0 * A), $MachinePrecision] + N[(-0.5 * N[(N[(B$95$m * B$95$m), $MachinePrecision] / C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(N[Power[N[(2.0 * F), $MachinePrecision], 0.5], $MachinePrecision] / N[Sqrt[B$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;B\_m \leq 3.8 \cdot 10^{-133}:\\
\;\;\;\;\frac{\sqrt{\left(B\_m \cdot B\_m + -4 \cdot \left(A \cdot C\right)\right) \cdot \left(\left(2 \cdot F\right) \cdot \left(2 \cdot A + -0.5 \cdot \frac{B\_m \cdot B\_m}{C}\right)\right)}}{4 \cdot \left(A \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{else}:\\
\;\;\;\;0 - \frac{{\left(2 \cdot F\right)}^{0.5}}{\sqrt{B\_m}}\\
\end{array}
\end{array}
if B < 3.8000000000000003e-133Initial program 21.7%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified27.8%
associate-*l*N/A
associate-*r*N/A
pow2N/A
*-lowering-*.f64N/A
pow2N/A
associate-*r*N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr27.3%
Taylor expanded in C around -inf
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6414.5%
Simplified14.5%
if 3.8000000000000003e-133 < B Initial program 11.0%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified26.9%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr33.7%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6449.3%
Simplified49.3%
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6449.3%
Applied egg-rr49.3%
Final simplification27.1%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(if (<= B_m 4.5e-133)
(/
(sqrt
(*
(+ (* B_m B_m) (* -4.0 (* A C)))
(* (* 2.0 F) (+ (* 2.0 A) (* -0.5 (/ (* B_m B_m) C))))))
(- (* 4.0 (* A C)) (* B_m B_m)))
(* (sqrt (* 2.0 F)) (- 0.0 (pow B_m -0.5)))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 4.5e-133) {
tmp = sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m));
} else {
tmp = sqrt((2.0 * F)) * (0.0 - pow(B_m, -0.5));
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (b_m <= 4.5d-133) then
tmp = sqrt((((b_m * b_m) + ((-4.0d0) * (a * c))) * ((2.0d0 * f) * ((2.0d0 * a) + ((-0.5d0) * ((b_m * b_m) / c)))))) / ((4.0d0 * (a * c)) - (b_m * b_m))
else
tmp = sqrt((2.0d0 * f)) * (0.0d0 - (b_m ** (-0.5d0)))
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 4.5e-133) {
tmp = Math.sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m));
} else {
tmp = Math.sqrt((2.0 * F)) * (0.0 - Math.pow(B_m, -0.5));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if B_m <= 4.5e-133: tmp = math.sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m)) else: tmp = math.sqrt((2.0 * F)) * (0.0 - math.pow(B_m, -0.5)) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (B_m <= 4.5e-133) tmp = Float64(sqrt(Float64(Float64(Float64(B_m * B_m) + Float64(-4.0 * Float64(A * C))) * Float64(Float64(2.0 * F) * Float64(Float64(2.0 * A) + Float64(-0.5 * Float64(Float64(B_m * B_m) / C)))))) / Float64(Float64(4.0 * Float64(A * C)) - Float64(B_m * B_m))); else tmp = Float64(sqrt(Float64(2.0 * F)) * Float64(0.0 - (B_m ^ -0.5))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (B_m <= 4.5e-133) tmp = sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * ((2.0 * A) + (-0.5 * ((B_m * B_m) / C)))))) / ((4.0 * (A * C)) - (B_m * B_m)); else tmp = sqrt((2.0 * F)) * (0.0 - (B_m ^ -0.5)); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[B$95$m, 4.5e-133], N[(N[Sqrt[N[(N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * F), $MachinePrecision] * N[(N[(2.0 * A), $MachinePrecision] + N[(-0.5 * N[(N[(B$95$m * B$95$m), $MachinePrecision] / C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Sqrt[N[(2.0 * F), $MachinePrecision]], $MachinePrecision] * N[(0.0 - N[Power[B$95$m, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;B\_m \leq 4.5 \cdot 10^{-133}:\\
\;\;\;\;\frac{\sqrt{\left(B\_m \cdot B\_m + -4 \cdot \left(A \cdot C\right)\right) \cdot \left(\left(2 \cdot F\right) \cdot \left(2 \cdot A + -0.5 \cdot \frac{B\_m \cdot B\_m}{C}\right)\right)}}{4 \cdot \left(A \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2 \cdot F} \cdot \left(0 - {B\_m}^{-0.5}\right)\\
\end{array}
\end{array}
if B < 4.50000000000000009e-133Initial program 21.7%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified27.8%
associate-*l*N/A
associate-*r*N/A
pow2N/A
*-lowering-*.f64N/A
pow2N/A
associate-*r*N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr27.3%
Taylor expanded in C around -inf
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6414.5%
Simplified14.5%
if 4.50000000000000009e-133 < B Initial program 11.0%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified26.9%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr33.7%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6449.3%
Simplified49.3%
inv-powN/A
sqrt-pow1N/A
metadata-evalN/A
pow-lowering-pow.f6449.3%
Applied egg-rr49.3%
Final simplification27.1%
B_m = (fabs.f64 B)
(FPCore (A B_m C F)
:precision binary64
(if (<= B_m 6.2e-282)
(/
(sqrt (* (+ (* B_m B_m) (* -4.0 (* A C))) (* (* 2.0 F) (* 2.0 A))))
(- (* 4.0 (* A C)) (* B_m B_m)))
(if (<= B_m 6e+19)
(/ (* C (pow (* F (* A -16.0)) 0.5)) (- (* A (* 4.0 C)) (* B_m B_m)))
(- 0.0 (pow (/ (* 2.0 F) B_m) 0.5)))))B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 6.2e-282) {
tmp = sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * (2.0 * A)))) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 6e+19) {
tmp = (C * pow((F * (A * -16.0)), 0.5)) / ((A * (4.0 * C)) - (B_m * B_m));
} else {
tmp = 0.0 - pow(((2.0 * F) / B_m), 0.5);
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (b_m <= 6.2d-282) then
tmp = sqrt((((b_m * b_m) + ((-4.0d0) * (a * c))) * ((2.0d0 * f) * (2.0d0 * a)))) / ((4.0d0 * (a * c)) - (b_m * b_m))
else if (b_m <= 6d+19) then
tmp = (c * ((f * (a * (-16.0d0))) ** 0.5d0)) / ((a * (4.0d0 * c)) - (b_m * b_m))
else
tmp = 0.0d0 - (((2.0d0 * f) / b_m) ** 0.5d0)
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 6.2e-282) {
tmp = Math.sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * (2.0 * A)))) / ((4.0 * (A * C)) - (B_m * B_m));
} else if (B_m <= 6e+19) {
tmp = (C * Math.pow((F * (A * -16.0)), 0.5)) / ((A * (4.0 * C)) - (B_m * B_m));
} else {
tmp = 0.0 - Math.pow(((2.0 * F) / B_m), 0.5);
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if B_m <= 6.2e-282: tmp = math.sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * (2.0 * A)))) / ((4.0 * (A * C)) - (B_m * B_m)) elif B_m <= 6e+19: tmp = (C * math.pow((F * (A * -16.0)), 0.5)) / ((A * (4.0 * C)) - (B_m * B_m)) else: tmp = 0.0 - math.pow(((2.0 * F) / B_m), 0.5) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (B_m <= 6.2e-282) tmp = Float64(sqrt(Float64(Float64(Float64(B_m * B_m) + Float64(-4.0 * Float64(A * C))) * Float64(Float64(2.0 * F) * Float64(2.0 * A)))) / Float64(Float64(4.0 * Float64(A * C)) - Float64(B_m * B_m))); elseif (B_m <= 6e+19) tmp = Float64(Float64(C * (Float64(F * Float64(A * -16.0)) ^ 0.5)) / Float64(Float64(A * Float64(4.0 * C)) - Float64(B_m * B_m))); else tmp = Float64(0.0 - (Float64(Float64(2.0 * F) / B_m) ^ 0.5)); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (B_m <= 6.2e-282) tmp = sqrt((((B_m * B_m) + (-4.0 * (A * C))) * ((2.0 * F) * (2.0 * A)))) / ((4.0 * (A * C)) - (B_m * B_m)); elseif (B_m <= 6e+19) tmp = (C * ((F * (A * -16.0)) ^ 0.5)) / ((A * (4.0 * C)) - (B_m * B_m)); else tmp = 0.0 - (((2.0 * F) / B_m) ^ 0.5); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[B$95$m, 6.2e-282], N[(N[Sqrt[N[(N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * F), $MachinePrecision] * N[(2.0 * A), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[B$95$m, 6e+19], N[(N[(C * N[Power[N[(F * N[(A * -16.0), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision]), $MachinePrecision] / N[(N[(A * N[(4.0 * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[Power[N[(N[(2.0 * F), $MachinePrecision] / B$95$m), $MachinePrecision], 0.5], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;B\_m \leq 6.2 \cdot 10^{-282}:\\
\;\;\;\;\frac{\sqrt{\left(B\_m \cdot B\_m + -4 \cdot \left(A \cdot C\right)\right) \cdot \left(\left(2 \cdot F\right) \cdot \left(2 \cdot A\right)\right)}}{4 \cdot \left(A \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{elif}\;B\_m \leq 6 \cdot 10^{+19}:\\
\;\;\;\;\frac{C \cdot {\left(F \cdot \left(A \cdot -16\right)\right)}^{0.5}}{A \cdot \left(4 \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{else}:\\
\;\;\;\;0 - {\left(\frac{2 \cdot F}{B\_m}\right)}^{0.5}\\
\end{array}
\end{array}
if B < 6.20000000000000027e-282Initial program 21.3%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified28.0%
associate-*l*N/A
associate-*r*N/A
pow2N/A
*-lowering-*.f64N/A
pow2N/A
associate-*r*N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr27.3%
Taylor expanded in A around inf
*-lowering-*.f6412.3%
Simplified12.3%
if 6.20000000000000027e-282 < B < 6e19Initial program 19.7%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified25.6%
Taylor expanded in A around -inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6417.8%
Simplified17.8%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
Applied egg-rr17.3%
clear-numN/A
/-lowering-/.f64N/A
pow1/2N/A
associate-*r*N/A
associate-*r*N/A
unpow-prod-downN/A
pow1/2N/A
pow2N/A
sqrt-pow1N/A
metadata-evalN/A
metadata-evalN/A
pow-flipN/A
inv-powN/A
remove-double-divN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
Applied egg-rr28.9%
if 6e19 < B Initial program 7.6%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified20.5%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr28.7%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6467.7%
Simplified67.7%
sqrt-unprodN/A
pow1/2N/A
pow-lowering-pow.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6456.3%
Applied egg-rr56.3%
Final simplification26.5%
B_m = (fabs.f64 B) (FPCore (A B_m C F) :precision binary64 (if (<= B_m 2.3e+19) (/ (* C (pow (* F (* A -16.0)) 0.5)) (- (* A (* 4.0 C)) (* B_m B_m))) (- 0.0 (pow (/ (* 2.0 F) B_m) 0.5))))
B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 2.3e+19) {
tmp = (C * pow((F * (A * -16.0)), 0.5)) / ((A * (4.0 * C)) - (B_m * B_m));
} else {
tmp = 0.0 - pow(((2.0 * F) / B_m), 0.5);
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (b_m <= 2.3d+19) then
tmp = (c * ((f * (a * (-16.0d0))) ** 0.5d0)) / ((a * (4.0d0 * c)) - (b_m * b_m))
else
tmp = 0.0d0 - (((2.0d0 * f) / b_m) ** 0.5d0)
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 2.3e+19) {
tmp = (C * Math.pow((F * (A * -16.0)), 0.5)) / ((A * (4.0 * C)) - (B_m * B_m));
} else {
tmp = 0.0 - Math.pow(((2.0 * F) / B_m), 0.5);
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if B_m <= 2.3e+19: tmp = (C * math.pow((F * (A * -16.0)), 0.5)) / ((A * (4.0 * C)) - (B_m * B_m)) else: tmp = 0.0 - math.pow(((2.0 * F) / B_m), 0.5) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (B_m <= 2.3e+19) tmp = Float64(Float64(C * (Float64(F * Float64(A * -16.0)) ^ 0.5)) / Float64(Float64(A * Float64(4.0 * C)) - Float64(B_m * B_m))); else tmp = Float64(0.0 - (Float64(Float64(2.0 * F) / B_m) ^ 0.5)); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (B_m <= 2.3e+19) tmp = (C * ((F * (A * -16.0)) ^ 0.5)) / ((A * (4.0 * C)) - (B_m * B_m)); else tmp = 0.0 - (((2.0 * F) / B_m) ^ 0.5); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[B$95$m, 2.3e+19], N[(N[(C * N[Power[N[(F * N[(A * -16.0), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision]), $MachinePrecision] / N[(N[(A * N[(4.0 * C), $MachinePrecision]), $MachinePrecision] - N[(B$95$m * B$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[Power[N[(N[(2.0 * F), $MachinePrecision] / B$95$m), $MachinePrecision], 0.5], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;B\_m \leq 2.3 \cdot 10^{+19}:\\
\;\;\;\;\frac{C \cdot {\left(F \cdot \left(A \cdot -16\right)\right)}^{0.5}}{A \cdot \left(4 \cdot C\right) - B\_m \cdot B\_m}\\
\mathbf{else}:\\
\;\;\;\;0 - {\left(\frac{2 \cdot F}{B\_m}\right)}^{0.5}\\
\end{array}
\end{array}
if B < 2.3e19Initial program 20.8%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified27.2%
Taylor expanded in A around -inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f649.7%
Simplified9.7%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
Applied egg-rr9.6%
clear-numN/A
/-lowering-/.f64N/A
pow1/2N/A
associate-*r*N/A
associate-*r*N/A
unpow-prod-downN/A
pow1/2N/A
pow2N/A
sqrt-pow1N/A
metadata-evalN/A
metadata-evalN/A
pow-flipN/A
inv-powN/A
remove-double-divN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
Applied egg-rr17.1%
if 2.3e19 < B Initial program 7.6%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified20.5%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr28.7%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6467.7%
Simplified67.7%
sqrt-unprodN/A
pow1/2N/A
pow-lowering-pow.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6456.3%
Applied egg-rr56.3%
Final simplification26.0%
B_m = (fabs.f64 B) (FPCore (A B_m C F) :precision binary64 (if (<= A 1.5e+58) (- 0.0 (pow (/ (* 2.0 F) B_m) 0.5)) (- 0.0 (sqrt (/ (* 4.0 F) (/ (+ (* B_m B_m) (/ A (/ -0.25 C))) A))))))
B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (A <= 1.5e+58) {
tmp = 0.0 - pow(((2.0 * F) / B_m), 0.5);
} else {
tmp = 0.0 - sqrt(((4.0 * F) / (((B_m * B_m) + (A / (-0.25 / C))) / A)));
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (a <= 1.5d+58) then
tmp = 0.0d0 - (((2.0d0 * f) / b_m) ** 0.5d0)
else
tmp = 0.0d0 - sqrt(((4.0d0 * f) / (((b_m * b_m) + (a / ((-0.25d0) / c))) / a)))
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (A <= 1.5e+58) {
tmp = 0.0 - Math.pow(((2.0 * F) / B_m), 0.5);
} else {
tmp = 0.0 - Math.sqrt(((4.0 * F) / (((B_m * B_m) + (A / (-0.25 / C))) / A)));
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if A <= 1.5e+58: tmp = 0.0 - math.pow(((2.0 * F) / B_m), 0.5) else: tmp = 0.0 - math.sqrt(((4.0 * F) / (((B_m * B_m) + (A / (-0.25 / C))) / A))) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (A <= 1.5e+58) tmp = Float64(0.0 - (Float64(Float64(2.0 * F) / B_m) ^ 0.5)); else tmp = Float64(0.0 - sqrt(Float64(Float64(4.0 * F) / Float64(Float64(Float64(B_m * B_m) + Float64(A / Float64(-0.25 / C))) / A)))); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (A <= 1.5e+58) tmp = 0.0 - (((2.0 * F) / B_m) ^ 0.5); else tmp = 0.0 - sqrt(((4.0 * F) / (((B_m * B_m) + (A / (-0.25 / C))) / A))); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[A, 1.5e+58], N[(0.0 - N[Power[N[(N[(2.0 * F), $MachinePrecision] / B$95$m), $MachinePrecision], 0.5], $MachinePrecision]), $MachinePrecision], N[(0.0 - N[Sqrt[N[(N[(4.0 * F), $MachinePrecision] / N[(N[(N[(B$95$m * B$95$m), $MachinePrecision] + N[(A / N[(-0.25 / C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / A), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;A \leq 1.5 \cdot 10^{+58}:\\
\;\;\;\;0 - {\left(\frac{2 \cdot F}{B\_m}\right)}^{0.5}\\
\mathbf{else}:\\
\;\;\;\;0 - \sqrt{\frac{4 \cdot F}{\frac{B\_m \cdot B\_m + \frac{A}{\frac{-0.25}{C}}}{A}}}\\
\end{array}
\end{array}
if A < 1.5000000000000001e58Initial program 17.1%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified25.7%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr33.6%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6421.1%
Simplified21.1%
sqrt-unprodN/A
pow1/2N/A
pow-lowering-pow.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6417.8%
Applied egg-rr17.8%
if 1.5000000000000001e58 < A Initial program 20.8%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified30.6%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr41.4%
Taylor expanded in A around inf
*-lowering-*.f6441.4%
Simplified41.4%
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-/l*N/A
associate-*r*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr34.9%
Final simplification21.3%
B_m = (fabs.f64 B) (FPCore (A B_m C F) :precision binary64 (if (<= B_m 1.3e-139) (/ (sqrt (* (* -16.0 (* A A)) (* C F))) (* 4.0 (* A C))) (- 0.0 (pow (/ (* 2.0 F) B_m) 0.5))))
B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 1.3e-139) {
tmp = sqrt(((-16.0 * (A * A)) * (C * F))) / (4.0 * (A * C));
} else {
tmp = 0.0 - pow(((2.0 * F) / B_m), 0.5);
}
return tmp;
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
real(8) :: tmp
if (b_m <= 1.3d-139) then
tmp = sqrt((((-16.0d0) * (a * a)) * (c * f))) / (4.0d0 * (a * c))
else
tmp = 0.0d0 - (((2.0d0 * f) / b_m) ** 0.5d0)
end if
code = tmp
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
double tmp;
if (B_m <= 1.3e-139) {
tmp = Math.sqrt(((-16.0 * (A * A)) * (C * F))) / (4.0 * (A * C));
} else {
tmp = 0.0 - Math.pow(((2.0 * F) / B_m), 0.5);
}
return tmp;
}
B_m = math.fabs(B) def code(A, B_m, C, F): tmp = 0 if B_m <= 1.3e-139: tmp = math.sqrt(((-16.0 * (A * A)) * (C * F))) / (4.0 * (A * C)) else: tmp = 0.0 - math.pow(((2.0 * F) / B_m), 0.5) return tmp
B_m = abs(B) function code(A, B_m, C, F) tmp = 0.0 if (B_m <= 1.3e-139) tmp = Float64(sqrt(Float64(Float64(-16.0 * Float64(A * A)) * Float64(C * F))) / Float64(4.0 * Float64(A * C))); else tmp = Float64(0.0 - (Float64(Float64(2.0 * F) / B_m) ^ 0.5)); end return tmp end
B_m = abs(B); function tmp_2 = code(A, B_m, C, F) tmp = 0.0; if (B_m <= 1.3e-139) tmp = sqrt(((-16.0 * (A * A)) * (C * F))) / (4.0 * (A * C)); else tmp = 0.0 - (((2.0 * F) / B_m) ^ 0.5); end tmp_2 = tmp; end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := If[LessEqual[B$95$m, 1.3e-139], N[(N[Sqrt[N[(N[(-16.0 * N[(A * A), $MachinePrecision]), $MachinePrecision] * N[(C * F), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[Power[N[(N[(2.0 * F), $MachinePrecision] / B$95$m), $MachinePrecision], 0.5], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
B_m = \left|B\right|
\\
\begin{array}{l}
\mathbf{if}\;B\_m \leq 1.3 \cdot 10^{-139}:\\
\;\;\;\;\frac{\sqrt{\left(-16 \cdot \left(A \cdot A\right)\right) \cdot \left(C \cdot F\right)}}{4 \cdot \left(A \cdot C\right)}\\
\mathbf{else}:\\
\;\;\;\;0 - {\left(\frac{2 \cdot F}{B\_m}\right)}^{0.5}\\
\end{array}
\end{array}
if B < 1.2999999999999999e-139Initial program 21.9%
distribute-frac-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified28.0%
Taylor expanded in B around 0
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6410.1%
Simplified10.1%
Taylor expanded in A around inf
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6410.4%
Simplified10.4%
if 1.2999999999999999e-139 < B Initial program 10.9%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified26.6%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr33.3%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6448.9%
Simplified48.9%
sqrt-unprodN/A
pow1/2N/A
pow-lowering-pow.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6440.8%
Applied egg-rr40.8%
Final simplification21.6%
B_m = (fabs.f64 B) (FPCore (A B_m C F) :precision binary64 (- 0.0 (pow (/ (* 2.0 F) B_m) 0.5)))
B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
return 0.0 - pow(((2.0 * F) / B_m), 0.5);
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
code = 0.0d0 - (((2.0d0 * f) / b_m) ** 0.5d0)
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
return 0.0 - Math.pow(((2.0 * F) / B_m), 0.5);
}
B_m = math.fabs(B) def code(A, B_m, C, F): return 0.0 - math.pow(((2.0 * F) / B_m), 0.5)
B_m = abs(B) function code(A, B_m, C, F) return Float64(0.0 - (Float64(Float64(2.0 * F) / B_m) ^ 0.5)) end
B_m = abs(B); function tmp = code(A, B_m, C, F) tmp = 0.0 - (((2.0 * F) / B_m) ^ 0.5); end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := N[(0.0 - N[Power[N[(N[(2.0 * F), $MachinePrecision] / B$95$m), $MachinePrecision], 0.5], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
B_m = \left|B\right|
\\
0 - {\left(\frac{2 \cdot F}{B\_m}\right)}^{0.5}
\end{array}
Initial program 17.8%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified26.7%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr35.2%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6419.1%
Simplified19.1%
sqrt-unprodN/A
pow1/2N/A
pow-lowering-pow.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6415.9%
Applied egg-rr15.9%
Final simplification15.9%
B_m = (fabs.f64 B) (FPCore (A B_m C F) :precision binary64 (- 0.0 (sqrt (/ (* 2.0 F) B_m))))
B_m = fabs(B);
double code(double A, double B_m, double C, double F) {
return 0.0 - sqrt(((2.0 * F) / B_m));
}
B_m = abs(b)
real(8) function code(a, b_m, c, f)
real(8), intent (in) :: a
real(8), intent (in) :: b_m
real(8), intent (in) :: c
real(8), intent (in) :: f
code = 0.0d0 - sqrt(((2.0d0 * f) / b_m))
end function
B_m = Math.abs(B);
public static double code(double A, double B_m, double C, double F) {
return 0.0 - Math.sqrt(((2.0 * F) / B_m));
}
B_m = math.fabs(B) def code(A, B_m, C, F): return 0.0 - math.sqrt(((2.0 * F) / B_m))
B_m = abs(B) function code(A, B_m, C, F) return Float64(0.0 - sqrt(Float64(Float64(2.0 * F) / B_m))) end
B_m = abs(B); function tmp = code(A, B_m, C, F) tmp = 0.0 - sqrt(((2.0 * F) / B_m)); end
B_m = N[Abs[B], $MachinePrecision] code[A_, B$95$m_, C_, F_] := N[(0.0 - N[Sqrt[N[(N[(2.0 * F), $MachinePrecision] / B$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
B_m = \left|B\right|
\\
0 - \sqrt{\frac{2 \cdot F}{B\_m}}
\end{array}
Initial program 17.8%
Taylor expanded in F around 0
mul-1-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
Simplified26.7%
*-commutativeN/A
pow1/2N/A
associate-+r+N/A
associate-/l*N/A
unpow-prod-downN/A
associate-*r*N/A
Applied egg-rr35.2%
Taylor expanded in B around inf
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f6419.1%
Simplified19.1%
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6415.7%
Applied egg-rr15.7%
Final simplification15.7%
herbie shell --seed 2024164
(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))))