
(FPCore (c0 w h D d M) :precision binary64 (let* ((t_0 (/ (* c0 (* d d)) (* (* w h) (* D D))))) (* (/ c0 (* 2.0 w)) (+ t_0 (sqrt (- (* t_0 t_0) (* M M)))))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 * (d * d)) / ((w * h) * (D * D));
return (c0 / (2.0 * w)) * (t_0 + sqrt(((t_0 * t_0) - (M * M))));
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
real(8) :: t_0
t_0 = (c0 * (d_1 * d_1)) / ((w * h) * (d * d))
code = (c0 / (2.0d0 * w)) * (t_0 + sqrt(((t_0 * t_0) - (m * m))))
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 * (d * d)) / ((w * h) * (D * D));
return (c0 / (2.0 * w)) * (t_0 + Math.sqrt(((t_0 * t_0) - (M * M))));
}
def code(c0, w, h, D, d, M): t_0 = (c0 * (d * d)) / ((w * h) * (D * D)) return (c0 / (2.0 * w)) * (t_0 + math.sqrt(((t_0 * t_0) - (M * M))))
function code(c0, w, h, D, d, M) t_0 = Float64(Float64(c0 * Float64(d * d)) / Float64(Float64(w * h) * Float64(D * D))) return Float64(Float64(c0 / Float64(2.0 * w)) * Float64(t_0 + sqrt(Float64(Float64(t_0 * t_0) - Float64(M * M))))) end
function tmp = code(c0, w, h, D, d, M) t_0 = (c0 * (d * d)) / ((w * h) * (D * D)); tmp = (c0 / (2.0 * w)) * (t_0 + sqrt(((t_0 * t_0) - (M * M)))); end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision] / N[(N[(w * h), $MachinePrecision] * N[(D * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(N[(c0 / N[(2.0 * w), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 + N[Sqrt[N[(N[(t$95$0 * t$95$0), $MachinePrecision] - N[(M * M), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)}\\
\frac{c0}{2 \cdot w} \cdot \left(t\_0 + \sqrt{t\_0 \cdot t\_0 - M \cdot M}\right)
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 11 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (c0 w h D d M) :precision binary64 (let* ((t_0 (/ (* c0 (* d d)) (* (* w h) (* D D))))) (* (/ c0 (* 2.0 w)) (+ t_0 (sqrt (- (* t_0 t_0) (* M M)))))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 * (d * d)) / ((w * h) * (D * D));
return (c0 / (2.0 * w)) * (t_0 + sqrt(((t_0 * t_0) - (M * M))));
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
real(8) :: t_0
t_0 = (c0 * (d_1 * d_1)) / ((w * h) * (d * d))
code = (c0 / (2.0d0 * w)) * (t_0 + sqrt(((t_0 * t_0) - (m * m))))
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 * (d * d)) / ((w * h) * (D * D));
return (c0 / (2.0 * w)) * (t_0 + Math.sqrt(((t_0 * t_0) - (M * M))));
}
def code(c0, w, h, D, d, M): t_0 = (c0 * (d * d)) / ((w * h) * (D * D)) return (c0 / (2.0 * w)) * (t_0 + math.sqrt(((t_0 * t_0) - (M * M))))
function code(c0, w, h, D, d, M) t_0 = Float64(Float64(c0 * Float64(d * d)) / Float64(Float64(w * h) * Float64(D * D))) return Float64(Float64(c0 / Float64(2.0 * w)) * Float64(t_0 + sqrt(Float64(Float64(t_0 * t_0) - Float64(M * M))))) end
function tmp = code(c0, w, h, D, d, M) t_0 = (c0 * (d * d)) / ((w * h) * (D * D)); tmp = (c0 / (2.0 * w)) * (t_0 + sqrt(((t_0 * t_0) - (M * M)))); end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision] / N[(N[(w * h), $MachinePrecision] * N[(D * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(N[(c0 / N[(2.0 * w), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 + N[Sqrt[N[(N[(t$95$0 * t$95$0), $MachinePrecision] - N[(M * M), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)}\\
\frac{c0}{2 \cdot w} \cdot \left(t\_0 + \sqrt{t\_0 \cdot t\_0 - M \cdot M}\right)
\end{array}
\end{array}
(FPCore (c0 w h D d M)
:precision binary64
(let* ((t_0 (/ (* c0 (* d d)) (* (* w h) (* D D)))))
(if (<=
(* (/ c0 (* 2.0 w)) (+ t_0 (sqrt (- (* t_0 t_0) (* M M)))))
INFINITY)
(/ (/ c0 2.0) (/ (/ w (/ (* c0 d) (* D (* w (* h D))))) (* 2.0 d)))
(* (/ D d) (/ (* D (* h (* (* M M) 0.25))) d)))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 * (d * d)) / ((w * h) * (D * D));
double tmp;
if (((c0 / (2.0 * w)) * (t_0 + sqrt(((t_0 * t_0) - (M * M))))) <= ((double) INFINITY)) {
tmp = (c0 / 2.0) / ((w / ((c0 * d) / (D * (w * (h * D))))) / (2.0 * d));
} else {
tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d);
}
return tmp;
}
public static double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 * (d * d)) / ((w * h) * (D * D));
double tmp;
if (((c0 / (2.0 * w)) * (t_0 + Math.sqrt(((t_0 * t_0) - (M * M))))) <= Double.POSITIVE_INFINITY) {
tmp = (c0 / 2.0) / ((w / ((c0 * d) / (D * (w * (h * D))))) / (2.0 * d));
} else {
tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d);
}
return tmp;
}
def code(c0, w, h, D, d, M): t_0 = (c0 * (d * d)) / ((w * h) * (D * D)) tmp = 0 if ((c0 / (2.0 * w)) * (t_0 + math.sqrt(((t_0 * t_0) - (M * M))))) <= math.inf: tmp = (c0 / 2.0) / ((w / ((c0 * d) / (D * (w * (h * D))))) / (2.0 * d)) else: tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d) return tmp
function code(c0, w, h, D, d, M) t_0 = Float64(Float64(c0 * Float64(d * d)) / Float64(Float64(w * h) * Float64(D * D))) tmp = 0.0 if (Float64(Float64(c0 / Float64(2.0 * w)) * Float64(t_0 + sqrt(Float64(Float64(t_0 * t_0) - Float64(M * M))))) <= Inf) tmp = Float64(Float64(c0 / 2.0) / Float64(Float64(w / Float64(Float64(c0 * d) / Float64(D * Float64(w * Float64(h * D))))) / Float64(2.0 * d))); else tmp = Float64(Float64(D / d) * Float64(Float64(D * Float64(h * Float64(Float64(M * M) * 0.25))) / d)); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) t_0 = (c0 * (d * d)) / ((w * h) * (D * D)); tmp = 0.0; if (((c0 / (2.0 * w)) * (t_0 + sqrt(((t_0 * t_0) - (M * M))))) <= Inf) tmp = (c0 / 2.0) / ((w / ((c0 * d) / (D * (w * (h * D))))) / (2.0 * d)); else tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision] / N[(N[(w * h), $MachinePrecision] * N[(D * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(c0 / N[(2.0 * w), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 + N[Sqrt[N[(N[(t$95$0 * t$95$0), $MachinePrecision] - N[(M * M), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(c0 / 2.0), $MachinePrecision] / N[(N[(w / N[(N[(c0 * d), $MachinePrecision] / N[(D * N[(w * N[(h * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(D / d), $MachinePrecision] * N[(N[(D * N[(h * N[(N[(M * M), $MachinePrecision] * 0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)}\\
\mathbf{if}\;\frac{c0}{2 \cdot w} \cdot \left(t\_0 + \sqrt{t\_0 \cdot t\_0 - M \cdot M}\right) \leq \infty:\\
\;\;\;\;\frac{\frac{c0}{2}}{\frac{\frac{w}{\frac{c0 \cdot d}{D \cdot \left(w \cdot \left(h \cdot D\right)\right)}}}{2 \cdot d}}\\
\mathbf{else}:\\
\;\;\;\;\frac{D}{d} \cdot \frac{D \cdot \left(h \cdot \left(\left(M \cdot M\right) \cdot 0.25\right)\right)}{d}\\
\end{array}
\end{array}
if (*.f64 (/.f64 c0 (*.f64 #s(literal 2 binary64) w)) (+.f64 (/.f64 (*.f64 c0 (*.f64 d d)) (*.f64 (*.f64 w h) (*.f64 D D))) (sqrt.f64 (-.f64 (*.f64 (/.f64 (*.f64 c0 (*.f64 d d)) (*.f64 (*.f64 w h) (*.f64 D D))) (/.f64 (*.f64 c0 (*.f64 d d)) (*.f64 (*.f64 w h) (*.f64 D D)))) (*.f64 M M))))) < +inf.0Initial program 71.2%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified61.9%
Taylor expanded in c0 around inf
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6472.8%
Simplified72.8%
associate-/l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6481.5%
Applied egg-rr81.5%
associate-*r*N/A
associate-*r*N/A
associate-/r*N/A
associate-*r/N/A
associate-/l*N/A
*-commutativeN/A
associate-/l/N/A
associate-/r*N/A
clear-numN/A
frac-timesN/A
*-rgt-identityN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr74.6%
associate-/r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6483.4%
Applied egg-rr83.4%
if +inf.0 < (*.f64 (/.f64 c0 (*.f64 #s(literal 2 binary64) w)) (+.f64 (/.f64 (*.f64 c0 (*.f64 d d)) (*.f64 (*.f64 w h) (*.f64 D D))) (sqrt.f64 (-.f64 (*.f64 (/.f64 (*.f64 c0 (*.f64 d d)) (*.f64 (*.f64 w h) (*.f64 D D))) (/.f64 (*.f64 c0 (*.f64 d d)) (*.f64 (*.f64 w h) (*.f64 D D)))) (*.f64 M M))))) Initial program 0.0%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr0.7%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified17.9%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6444.3%
Simplified44.3%
associate-*l*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.2%
Applied egg-rr63.2%
Final simplification70.9%
(FPCore (c0 w h D d M)
:precision binary64
(let* ((t_0 (/ (/ c0 2.0) w)))
(if (<= d 2.5e-254)
(* (/ D d) (/ (* D (* h (* (* M M) 0.25))) d))
(if (<= d 4.2e-60)
(* t_0 (* 2.0 (/ (/ d (/ D (/ d (* w h)))) (/ D c0))))
(if (<= d 6.5e+51)
(/ (* h (* 0.25 (* M (* D (* D M))))) (* d d))
(if (<= d 1.13e+225)
(* (* 2.0 d) (* (/ (* c0 d) (* D (* w (* h D)))) t_0))
(/
(/ c0 2.0)
(/ (* 0.5 (* D (* D (/ (/ (* h (* w w)) d) d)))) c0))))))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 / 2.0) / w;
double tmp;
if (d <= 2.5e-254) {
tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d);
} else if (d <= 4.2e-60) {
tmp = t_0 * (2.0 * ((d / (D / (d / (w * h)))) / (D / c0)));
} else if (d <= 6.5e+51) {
tmp = (h * (0.25 * (M * (D * (D * M))))) / (d * d);
} else if (d <= 1.13e+225) {
tmp = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * t_0);
} else {
tmp = (c0 / 2.0) / ((0.5 * (D * (D * (((h * (w * w)) / d) / d)))) / c0);
}
return tmp;
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: tmp
t_0 = (c0 / 2.0d0) / w
if (d_1 <= 2.5d-254) then
tmp = (d / d_1) * ((d * (h * ((m * m) * 0.25d0))) / d_1)
else if (d_1 <= 4.2d-60) then
tmp = t_0 * (2.0d0 * ((d_1 / (d / (d_1 / (w * h)))) / (d / c0)))
else if (d_1 <= 6.5d+51) then
tmp = (h * (0.25d0 * (m * (d * (d * m))))) / (d_1 * d_1)
else if (d_1 <= 1.13d+225) then
tmp = (2.0d0 * d_1) * (((c0 * d_1) / (d * (w * (h * d)))) * t_0)
else
tmp = (c0 / 2.0d0) / ((0.5d0 * (d * (d * (((h * (w * w)) / d_1) / d_1)))) / c0)
end if
code = tmp
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 / 2.0) / w;
double tmp;
if (d <= 2.5e-254) {
tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d);
} else if (d <= 4.2e-60) {
tmp = t_0 * (2.0 * ((d / (D / (d / (w * h)))) / (D / c0)));
} else if (d <= 6.5e+51) {
tmp = (h * (0.25 * (M * (D * (D * M))))) / (d * d);
} else if (d <= 1.13e+225) {
tmp = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * t_0);
} else {
tmp = (c0 / 2.0) / ((0.5 * (D * (D * (((h * (w * w)) / d) / d)))) / c0);
}
return tmp;
}
def code(c0, w, h, D, d, M): t_0 = (c0 / 2.0) / w tmp = 0 if d <= 2.5e-254: tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d) elif d <= 4.2e-60: tmp = t_0 * (2.0 * ((d / (D / (d / (w * h)))) / (D / c0))) elif d <= 6.5e+51: tmp = (h * (0.25 * (M * (D * (D * M))))) / (d * d) elif d <= 1.13e+225: tmp = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * t_0) else: tmp = (c0 / 2.0) / ((0.5 * (D * (D * (((h * (w * w)) / d) / d)))) / c0) return tmp
function code(c0, w, h, D, d, M) t_0 = Float64(Float64(c0 / 2.0) / w) tmp = 0.0 if (d <= 2.5e-254) tmp = Float64(Float64(D / d) * Float64(Float64(D * Float64(h * Float64(Float64(M * M) * 0.25))) / d)); elseif (d <= 4.2e-60) tmp = Float64(t_0 * Float64(2.0 * Float64(Float64(d / Float64(D / Float64(d / Float64(w * h)))) / Float64(D / c0)))); elseif (d <= 6.5e+51) tmp = Float64(Float64(h * Float64(0.25 * Float64(M * Float64(D * Float64(D * M))))) / Float64(d * d)); elseif (d <= 1.13e+225) tmp = Float64(Float64(2.0 * d) * Float64(Float64(Float64(c0 * d) / Float64(D * Float64(w * Float64(h * D)))) * t_0)); else tmp = Float64(Float64(c0 / 2.0) / Float64(Float64(0.5 * Float64(D * Float64(D * Float64(Float64(Float64(h * Float64(w * w)) / d) / d)))) / c0)); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) t_0 = (c0 / 2.0) / w; tmp = 0.0; if (d <= 2.5e-254) tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d); elseif (d <= 4.2e-60) tmp = t_0 * (2.0 * ((d / (D / (d / (w * h)))) / (D / c0))); elseif (d <= 6.5e+51) tmp = (h * (0.25 * (M * (D * (D * M))))) / (d * d); elseif (d <= 1.13e+225) tmp = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * t_0); else tmp = (c0 / 2.0) / ((0.5 * (D * (D * (((h * (w * w)) / d) / d)))) / c0); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(N[(c0 / 2.0), $MachinePrecision] / w), $MachinePrecision]}, If[LessEqual[d, 2.5e-254], N[(N[(D / d), $MachinePrecision] * N[(N[(D * N[(h * N[(N[(M * M), $MachinePrecision] * 0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 4.2e-60], N[(t$95$0 * N[(2.0 * N[(N[(d / N[(D / N[(d / N[(w * h), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(D / c0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 6.5e+51], N[(N[(h * N[(0.25 * N[(M * N[(D * N[(D * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(d * d), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 1.13e+225], N[(N[(2.0 * d), $MachinePrecision] * N[(N[(N[(c0 * d), $MachinePrecision] / N[(D * N[(w * N[(h * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], N[(N[(c0 / 2.0), $MachinePrecision] / N[(N[(0.5 * N[(D * N[(D * N[(N[(N[(h * N[(w * w), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c0), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\frac{c0}{2}}{w}\\
\mathbf{if}\;d \leq 2.5 \cdot 10^{-254}:\\
\;\;\;\;\frac{D}{d} \cdot \frac{D \cdot \left(h \cdot \left(\left(M \cdot M\right) \cdot 0.25\right)\right)}{d}\\
\mathbf{elif}\;d \leq 4.2 \cdot 10^{-60}:\\
\;\;\;\;t\_0 \cdot \left(2 \cdot \frac{\frac{d}{\frac{D}{\frac{d}{w \cdot h}}}}{\frac{D}{c0}}\right)\\
\mathbf{elif}\;d \leq 6.5 \cdot 10^{+51}:\\
\;\;\;\;\frac{h \cdot \left(0.25 \cdot \left(M \cdot \left(D \cdot \left(D \cdot M\right)\right)\right)\right)}{d \cdot d}\\
\mathbf{elif}\;d \leq 1.13 \cdot 10^{+225}:\\
\;\;\;\;\left(2 \cdot d\right) \cdot \left(\frac{c0 \cdot d}{D \cdot \left(w \cdot \left(h \cdot D\right)\right)} \cdot t\_0\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{c0}{2}}{\frac{0.5 \cdot \left(D \cdot \left(D \cdot \frac{\frac{h \cdot \left(w \cdot w\right)}{d}}{d}\right)\right)}{c0}}\\
\end{array}
\end{array}
if d < 2.5000000000000002e-254Initial program 26.7%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr25.7%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified12.2%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6431.4%
Simplified31.4%
associate-*l*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6444.6%
Applied egg-rr44.6%
if 2.5000000000000002e-254 < d < 4.19999999999999982e-60Initial program 41.3%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified29.5%
Taylor expanded in c0 around inf
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6441.6%
Simplified41.6%
associate-*l*N/A
associate-*l*N/A
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6441.8%
Applied egg-rr41.8%
*-commutativeN/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
associate-/l/N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6472.1%
Applied egg-rr72.1%
if 4.19999999999999982e-60 < d < 6.5e51Initial program 27.4%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr27.4%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified30.3%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6457.6%
Simplified57.6%
Taylor expanded in D around 0
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6472.6%
Simplified72.6%
if 6.5e51 < d < 1.13e225Initial program 28.0%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified25.1%
Taylor expanded in c0 around inf
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6435.6%
Simplified35.6%
associate-/l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6450.0%
Applied egg-rr50.0%
Taylor expanded in d around 0
associate-*r/N/A
unpow2N/A
associate-*l*N/A
associate-/l/N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6447.7%
Simplified47.7%
associate-*r*N/A
associate-/l/N/A
associate-*r/N/A
associate-/l/N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
clear-numN/A
div-invN/A
associate-*r/N/A
associate-/r/N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr55.3%
if 1.13e225 < d Initial program 18.5%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified18.5%
Taylor expanded in c0 around inf
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6437.4%
Simplified37.4%
associate-/l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6437.9%
Applied egg-rr37.9%
associate-*r*N/A
associate-*r*N/A
associate-/r*N/A
associate-*r/N/A
associate-/l*N/A
*-commutativeN/A
associate-/l/N/A
associate-/r*N/A
clear-numN/A
frac-timesN/A
*-rgt-identityN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr41.6%
Taylor expanded in w around 0
associate-*r/N/A
associate-/l/N/A
/-lowering-/.f64N/A
Simplified63.8%
Final simplification54.6%
(FPCore (c0 w h D d M)
:precision binary64
(let* ((t_0 (/ (/ c0 2.0) w)) (t_1 (* h (* (* M M) 0.25))))
(if (<= d 8.2e-257)
(* (/ D d) (/ (* D t_1) d))
(if (<= d 2.2e-63)
(* t_0 (* 2.0 (/ (/ d (/ D (/ d (* w h)))) (/ D c0))))
(if (<= d 5.2e+51)
(/ (* h (* 0.25 (* M (* D (* D M))))) (* d d))
(if (<= d 2.3e+279)
(* (* 2.0 d) (* (/ (* c0 d) (* D (* w (* h D)))) t_0))
(/ (/ (* D D) (/ d t_1)) d)))))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 / 2.0) / w;
double t_1 = h * ((M * M) * 0.25);
double tmp;
if (d <= 8.2e-257) {
tmp = (D / d) * ((D * t_1) / d);
} else if (d <= 2.2e-63) {
tmp = t_0 * (2.0 * ((d / (D / (d / (w * h)))) / (D / c0)));
} else if (d <= 5.2e+51) {
tmp = (h * (0.25 * (M * (D * (D * M))))) / (d * d);
} else if (d <= 2.3e+279) {
tmp = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * t_0);
} else {
tmp = ((D * D) / (d / t_1)) / d;
}
return tmp;
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (c0 / 2.0d0) / w
t_1 = h * ((m * m) * 0.25d0)
if (d_1 <= 8.2d-257) then
tmp = (d / d_1) * ((d * t_1) / d_1)
else if (d_1 <= 2.2d-63) then
tmp = t_0 * (2.0d0 * ((d_1 / (d / (d_1 / (w * h)))) / (d / c0)))
else if (d_1 <= 5.2d+51) then
tmp = (h * (0.25d0 * (m * (d * (d * m))))) / (d_1 * d_1)
else if (d_1 <= 2.3d+279) then
tmp = (2.0d0 * d_1) * (((c0 * d_1) / (d * (w * (h * d)))) * t_0)
else
tmp = ((d * d) / (d_1 / t_1)) / d_1
end if
code = tmp
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 / 2.0) / w;
double t_1 = h * ((M * M) * 0.25);
double tmp;
if (d <= 8.2e-257) {
tmp = (D / d) * ((D * t_1) / d);
} else if (d <= 2.2e-63) {
tmp = t_0 * (2.0 * ((d / (D / (d / (w * h)))) / (D / c0)));
} else if (d <= 5.2e+51) {
tmp = (h * (0.25 * (M * (D * (D * M))))) / (d * d);
} else if (d <= 2.3e+279) {
tmp = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * t_0);
} else {
tmp = ((D * D) / (d / t_1)) / d;
}
return tmp;
}
def code(c0, w, h, D, d, M): t_0 = (c0 / 2.0) / w t_1 = h * ((M * M) * 0.25) tmp = 0 if d <= 8.2e-257: tmp = (D / d) * ((D * t_1) / d) elif d <= 2.2e-63: tmp = t_0 * (2.0 * ((d / (D / (d / (w * h)))) / (D / c0))) elif d <= 5.2e+51: tmp = (h * (0.25 * (M * (D * (D * M))))) / (d * d) elif d <= 2.3e+279: tmp = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * t_0) else: tmp = ((D * D) / (d / t_1)) / d return tmp
function code(c0, w, h, D, d, M) t_0 = Float64(Float64(c0 / 2.0) / w) t_1 = Float64(h * Float64(Float64(M * M) * 0.25)) tmp = 0.0 if (d <= 8.2e-257) tmp = Float64(Float64(D / d) * Float64(Float64(D * t_1) / d)); elseif (d <= 2.2e-63) tmp = Float64(t_0 * Float64(2.0 * Float64(Float64(d / Float64(D / Float64(d / Float64(w * h)))) / Float64(D / c0)))); elseif (d <= 5.2e+51) tmp = Float64(Float64(h * Float64(0.25 * Float64(M * Float64(D * Float64(D * M))))) / Float64(d * d)); elseif (d <= 2.3e+279) tmp = Float64(Float64(2.0 * d) * Float64(Float64(Float64(c0 * d) / Float64(D * Float64(w * Float64(h * D)))) * t_0)); else tmp = Float64(Float64(Float64(D * D) / Float64(d / t_1)) / d); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) t_0 = (c0 / 2.0) / w; t_1 = h * ((M * M) * 0.25); tmp = 0.0; if (d <= 8.2e-257) tmp = (D / d) * ((D * t_1) / d); elseif (d <= 2.2e-63) tmp = t_0 * (2.0 * ((d / (D / (d / (w * h)))) / (D / c0))); elseif (d <= 5.2e+51) tmp = (h * (0.25 * (M * (D * (D * M))))) / (d * d); elseif (d <= 2.3e+279) tmp = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * t_0); else tmp = ((D * D) / (d / t_1)) / d; end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(N[(c0 / 2.0), $MachinePrecision] / w), $MachinePrecision]}, Block[{t$95$1 = N[(h * N[(N[(M * M), $MachinePrecision] * 0.25), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[d, 8.2e-257], N[(N[(D / d), $MachinePrecision] * N[(N[(D * t$95$1), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 2.2e-63], N[(t$95$0 * N[(2.0 * N[(N[(d / N[(D / N[(d / N[(w * h), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(D / c0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 5.2e+51], N[(N[(h * N[(0.25 * N[(M * N[(D * N[(D * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(d * d), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 2.3e+279], N[(N[(2.0 * d), $MachinePrecision] * N[(N[(N[(c0 * d), $MachinePrecision] / N[(D * N[(w * N[(h * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(D * D), $MachinePrecision] / N[(d / t$95$1), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\frac{c0}{2}}{w}\\
t_1 := h \cdot \left(\left(M \cdot M\right) \cdot 0.25\right)\\
\mathbf{if}\;d \leq 8.2 \cdot 10^{-257}:\\
\;\;\;\;\frac{D}{d} \cdot \frac{D \cdot t\_1}{d}\\
\mathbf{elif}\;d \leq 2.2 \cdot 10^{-63}:\\
\;\;\;\;t\_0 \cdot \left(2 \cdot \frac{\frac{d}{\frac{D}{\frac{d}{w \cdot h}}}}{\frac{D}{c0}}\right)\\
\mathbf{elif}\;d \leq 5.2 \cdot 10^{+51}:\\
\;\;\;\;\frac{h \cdot \left(0.25 \cdot \left(M \cdot \left(D \cdot \left(D \cdot M\right)\right)\right)\right)}{d \cdot d}\\
\mathbf{elif}\;d \leq 2.3 \cdot 10^{+279}:\\
\;\;\;\;\left(2 \cdot d\right) \cdot \left(\frac{c0 \cdot d}{D \cdot \left(w \cdot \left(h \cdot D\right)\right)} \cdot t\_0\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{D \cdot D}{\frac{d}{t\_1}}}{d}\\
\end{array}
\end{array}
if d < 8.1999999999999994e-257Initial program 26.7%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr25.7%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified12.2%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6431.4%
Simplified31.4%
associate-*l*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6444.6%
Applied egg-rr44.6%
if 8.1999999999999994e-257 < d < 2.2e-63Initial program 41.3%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified29.5%
Taylor expanded in c0 around inf
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6441.6%
Simplified41.6%
associate-*l*N/A
associate-*l*N/A
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6441.8%
Applied egg-rr41.8%
*-commutativeN/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
associate-/l/N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6472.1%
Applied egg-rr72.1%
if 2.2e-63 < d < 5.2000000000000002e51Initial program 27.4%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr27.4%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified30.3%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6457.6%
Simplified57.6%
Taylor expanded in D around 0
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6472.6%
Simplified72.6%
if 5.2000000000000002e51 < d < 2.3e279Initial program 27.5%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified25.1%
Taylor expanded in c0 around inf
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6437.2%
Simplified37.2%
associate-/l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6450.1%
Applied egg-rr50.1%
Taylor expanded in d around 0
associate-*r/N/A
unpow2N/A
associate-*l*N/A
associate-/l/N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6449.4%
Simplified49.4%
associate-*r*N/A
associate-/l/N/A
associate-*r/N/A
associate-/l/N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
clear-numN/A
div-invN/A
associate-*r/N/A
associate-/r/N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr54.4%
if 2.3e279 < d Initial program 9.1%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr9.1%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified45.6%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6445.6%
Simplified45.6%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.8%
Applied egg-rr63.8%
Final simplification53.8%
(FPCore (c0 w h D d M)
:precision binary64
(let* ((t_0 (* h (* (* M M) 0.25)))
(t_1 (* (/ D d) (/ (* D t_0) d)))
(t_2
(* (* 2.0 d) (* (/ (* c0 d) (* D (* w (* h D)))) (/ (/ c0 2.0) w)))))
(if (<= w -1.6e-50)
(/ (/ (* D D) (/ d t_0)) d)
(if (<= w -3.3e-288)
t_2
(if (<= w 4.9e-244) t_1 (if (<= w 4.5e+65) t_2 t_1))))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = h * ((M * M) * 0.25);
double t_1 = (D / d) * ((D * t_0) / d);
double t_2 = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * ((c0 / 2.0) / w));
double tmp;
if (w <= -1.6e-50) {
tmp = ((D * D) / (d / t_0)) / d;
} else if (w <= -3.3e-288) {
tmp = t_2;
} else if (w <= 4.9e-244) {
tmp = t_1;
} else if (w <= 4.5e+65) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = h * ((m * m) * 0.25d0)
t_1 = (d / d_1) * ((d * t_0) / d_1)
t_2 = (2.0d0 * d_1) * (((c0 * d_1) / (d * (w * (h * d)))) * ((c0 / 2.0d0) / w))
if (w <= (-1.6d-50)) then
tmp = ((d * d) / (d_1 / t_0)) / d_1
else if (w <= (-3.3d-288)) then
tmp = t_2
else if (w <= 4.9d-244) then
tmp = t_1
else if (w <= 4.5d+65) then
tmp = t_2
else
tmp = t_1
end if
code = tmp
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = h * ((M * M) * 0.25);
double t_1 = (D / d) * ((D * t_0) / d);
double t_2 = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * ((c0 / 2.0) / w));
double tmp;
if (w <= -1.6e-50) {
tmp = ((D * D) / (d / t_0)) / d;
} else if (w <= -3.3e-288) {
tmp = t_2;
} else if (w <= 4.9e-244) {
tmp = t_1;
} else if (w <= 4.5e+65) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
def code(c0, w, h, D, d, M): t_0 = h * ((M * M) * 0.25) t_1 = (D / d) * ((D * t_0) / d) t_2 = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * ((c0 / 2.0) / w)) tmp = 0 if w <= -1.6e-50: tmp = ((D * D) / (d / t_0)) / d elif w <= -3.3e-288: tmp = t_2 elif w <= 4.9e-244: tmp = t_1 elif w <= 4.5e+65: tmp = t_2 else: tmp = t_1 return tmp
function code(c0, w, h, D, d, M) t_0 = Float64(h * Float64(Float64(M * M) * 0.25)) t_1 = Float64(Float64(D / d) * Float64(Float64(D * t_0) / d)) t_2 = Float64(Float64(2.0 * d) * Float64(Float64(Float64(c0 * d) / Float64(D * Float64(w * Float64(h * D)))) * Float64(Float64(c0 / 2.0) / w))) tmp = 0.0 if (w <= -1.6e-50) tmp = Float64(Float64(Float64(D * D) / Float64(d / t_0)) / d); elseif (w <= -3.3e-288) tmp = t_2; elseif (w <= 4.9e-244) tmp = t_1; elseif (w <= 4.5e+65) tmp = t_2; else tmp = t_1; end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) t_0 = h * ((M * M) * 0.25); t_1 = (D / d) * ((D * t_0) / d); t_2 = (2.0 * d) * (((c0 * d) / (D * (w * (h * D)))) * ((c0 / 2.0) / w)); tmp = 0.0; if (w <= -1.6e-50) tmp = ((D * D) / (d / t_0)) / d; elseif (w <= -3.3e-288) tmp = t_2; elseif (w <= 4.9e-244) tmp = t_1; elseif (w <= 4.5e+65) tmp = t_2; else tmp = t_1; end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(h * N[(N[(M * M), $MachinePrecision] * 0.25), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(D / d), $MachinePrecision] * N[(N[(D * t$95$0), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(2.0 * d), $MachinePrecision] * N[(N[(N[(c0 * d), $MachinePrecision] / N[(D * N[(w * N[(h * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(c0 / 2.0), $MachinePrecision] / w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[w, -1.6e-50], N[(N[(N[(D * D), $MachinePrecision] / N[(d / t$95$0), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision], If[LessEqual[w, -3.3e-288], t$95$2, If[LessEqual[w, 4.9e-244], t$95$1, If[LessEqual[w, 4.5e+65], t$95$2, t$95$1]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := h \cdot \left(\left(M \cdot M\right) \cdot 0.25\right)\\
t_1 := \frac{D}{d} \cdot \frac{D \cdot t\_0}{d}\\
t_2 := \left(2 \cdot d\right) \cdot \left(\frac{c0 \cdot d}{D \cdot \left(w \cdot \left(h \cdot D\right)\right)} \cdot \frac{\frac{c0}{2}}{w}\right)\\
\mathbf{if}\;w \leq -1.6 \cdot 10^{-50}:\\
\;\;\;\;\frac{\frac{D \cdot D}{\frac{d}{t\_0}}}{d}\\
\mathbf{elif}\;w \leq -3.3 \cdot 10^{-288}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;w \leq 4.9 \cdot 10^{-244}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;w \leq 4.5 \cdot 10^{+65}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if w < -1.6e-50Initial program 18.4%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr20.1%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified20.2%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.4%
Simplified42.4%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6452.5%
Applied egg-rr52.5%
if -1.6e-50 < w < -3.29999999999999988e-288 or 4.90000000000000015e-244 < w < 4.5e65Initial program 31.4%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified29.3%
Taylor expanded in c0 around inf
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6443.8%
Simplified43.8%
associate-/l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6460.3%
Applied egg-rr60.3%
Taylor expanded in d around 0
associate-*r/N/A
unpow2N/A
associate-*l*N/A
associate-/l/N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6458.4%
Simplified58.4%
associate-*r*N/A
associate-/l/N/A
associate-*r/N/A
associate-/l/N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
clear-numN/A
div-invN/A
associate-*r/N/A
associate-/r/N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr60.8%
if -3.29999999999999988e-288 < w < 4.90000000000000015e-244 or 4.5e65 < w Initial program 24.5%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr21.3%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified22.9%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6445.7%
Simplified45.7%
associate-*l*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6468.2%
Applied egg-rr68.2%
Final simplification60.8%
(FPCore (c0 w h D d M)
:precision binary64
(let* ((t_0 (* h (* (* M M) 0.25)))
(t_1 (* (/ D d) (/ (* D t_0) d)))
(t_2
(* (* c0 d) (* (/ c0 (* 2.0 w)) (/ (* 2.0 d) (* w (* D (* h D))))))))
(if (<= w -5.5e-48)
(/ (/ (* D D) (/ d t_0)) d)
(if (<= w -3.3e-288)
t_2
(if (<= w 9.5e-245) t_1 (if (<= w 3e+67) t_2 t_1))))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = h * ((M * M) * 0.25);
double t_1 = (D / d) * ((D * t_0) / d);
double t_2 = (c0 * d) * ((c0 / (2.0 * w)) * ((2.0 * d) / (w * (D * (h * D)))));
double tmp;
if (w <= -5.5e-48) {
tmp = ((D * D) / (d / t_0)) / d;
} else if (w <= -3.3e-288) {
tmp = t_2;
} else if (w <= 9.5e-245) {
tmp = t_1;
} else if (w <= 3e+67) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = h * ((m * m) * 0.25d0)
t_1 = (d / d_1) * ((d * t_0) / d_1)
t_2 = (c0 * d_1) * ((c0 / (2.0d0 * w)) * ((2.0d0 * d_1) / (w * (d * (h * d)))))
if (w <= (-5.5d-48)) then
tmp = ((d * d) / (d_1 / t_0)) / d_1
else if (w <= (-3.3d-288)) then
tmp = t_2
else if (w <= 9.5d-245) then
tmp = t_1
else if (w <= 3d+67) then
tmp = t_2
else
tmp = t_1
end if
code = tmp
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = h * ((M * M) * 0.25);
double t_1 = (D / d) * ((D * t_0) / d);
double t_2 = (c0 * d) * ((c0 / (2.0 * w)) * ((2.0 * d) / (w * (D * (h * D)))));
double tmp;
if (w <= -5.5e-48) {
tmp = ((D * D) / (d / t_0)) / d;
} else if (w <= -3.3e-288) {
tmp = t_2;
} else if (w <= 9.5e-245) {
tmp = t_1;
} else if (w <= 3e+67) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
def code(c0, w, h, D, d, M): t_0 = h * ((M * M) * 0.25) t_1 = (D / d) * ((D * t_0) / d) t_2 = (c0 * d) * ((c0 / (2.0 * w)) * ((2.0 * d) / (w * (D * (h * D))))) tmp = 0 if w <= -5.5e-48: tmp = ((D * D) / (d / t_0)) / d elif w <= -3.3e-288: tmp = t_2 elif w <= 9.5e-245: tmp = t_1 elif w <= 3e+67: tmp = t_2 else: tmp = t_1 return tmp
function code(c0, w, h, D, d, M) t_0 = Float64(h * Float64(Float64(M * M) * 0.25)) t_1 = Float64(Float64(D / d) * Float64(Float64(D * t_0) / d)) t_2 = Float64(Float64(c0 * d) * Float64(Float64(c0 / Float64(2.0 * w)) * Float64(Float64(2.0 * d) / Float64(w * Float64(D * Float64(h * D)))))) tmp = 0.0 if (w <= -5.5e-48) tmp = Float64(Float64(Float64(D * D) / Float64(d / t_0)) / d); elseif (w <= -3.3e-288) tmp = t_2; elseif (w <= 9.5e-245) tmp = t_1; elseif (w <= 3e+67) tmp = t_2; else tmp = t_1; end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) t_0 = h * ((M * M) * 0.25); t_1 = (D / d) * ((D * t_0) / d); t_2 = (c0 * d) * ((c0 / (2.0 * w)) * ((2.0 * d) / (w * (D * (h * D))))); tmp = 0.0; if (w <= -5.5e-48) tmp = ((D * D) / (d / t_0)) / d; elseif (w <= -3.3e-288) tmp = t_2; elseif (w <= 9.5e-245) tmp = t_1; elseif (w <= 3e+67) tmp = t_2; else tmp = t_1; end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(h * N[(N[(M * M), $MachinePrecision] * 0.25), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(D / d), $MachinePrecision] * N[(N[(D * t$95$0), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(c0 * d), $MachinePrecision] * N[(N[(c0 / N[(2.0 * w), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * d), $MachinePrecision] / N[(w * N[(D * N[(h * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[w, -5.5e-48], N[(N[(N[(D * D), $MachinePrecision] / N[(d / t$95$0), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision], If[LessEqual[w, -3.3e-288], t$95$2, If[LessEqual[w, 9.5e-245], t$95$1, If[LessEqual[w, 3e+67], t$95$2, t$95$1]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := h \cdot \left(\left(M \cdot M\right) \cdot 0.25\right)\\
t_1 := \frac{D}{d} \cdot \frac{D \cdot t\_0}{d}\\
t_2 := \left(c0 \cdot d\right) \cdot \left(\frac{c0}{2 \cdot w} \cdot \frac{2 \cdot d}{w \cdot \left(D \cdot \left(h \cdot D\right)\right)}\right)\\
\mathbf{if}\;w \leq -5.5 \cdot 10^{-48}:\\
\;\;\;\;\frac{\frac{D \cdot D}{\frac{d}{t\_0}}}{d}\\
\mathbf{elif}\;w \leq -3.3 \cdot 10^{-288}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;w \leq 9.5 \cdot 10^{-245}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;w \leq 3 \cdot 10^{+67}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if w < -5.50000000000000047e-48Initial program 18.4%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr20.1%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified20.2%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.4%
Simplified42.4%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6452.5%
Applied egg-rr52.5%
if -5.50000000000000047e-48 < w < -3.29999999999999988e-288 or 9.5000000000000002e-245 < w < 3.0000000000000001e67Initial program 31.4%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified29.3%
Taylor expanded in c0 around inf
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6443.8%
Simplified43.8%
associate-/l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6460.3%
Applied egg-rr60.3%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr58.2%
if -3.29999999999999988e-288 < w < 9.5000000000000002e-245 or 3.0000000000000001e67 < w Initial program 24.5%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr21.3%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified22.9%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6445.7%
Simplified45.7%
associate-*l*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6468.2%
Applied egg-rr68.2%
Final simplification59.3%
(FPCore (c0 w h D d M)
:precision binary64
(let* ((t_0 (* h (* (* M M) 0.25)))
(t_1 (/ (* c0 (/ (* c0 (/ (/ (/ (* d d) D) D) w)) w)) h)))
(if (<= h -3.4e-239)
(* (/ D d) (/ (* D t_0) d))
(if (<= h 3.2e-104)
t_1
(if (<= h 2.25e+43) (/ (/ (* D D) (/ d t_0)) d) t_1)))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = h * ((M * M) * 0.25);
double t_1 = (c0 * ((c0 * ((((d * d) / D) / D) / w)) / w)) / h;
double tmp;
if (h <= -3.4e-239) {
tmp = (D / d) * ((D * t_0) / d);
} else if (h <= 3.2e-104) {
tmp = t_1;
} else if (h <= 2.25e+43) {
tmp = ((D * D) / (d / t_0)) / d;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = h * ((m * m) * 0.25d0)
t_1 = (c0 * ((c0 * ((((d_1 * d_1) / d) / d) / w)) / w)) / h
if (h <= (-3.4d-239)) then
tmp = (d / d_1) * ((d * t_0) / d_1)
else if (h <= 3.2d-104) then
tmp = t_1
else if (h <= 2.25d+43) then
tmp = ((d * d) / (d_1 / t_0)) / d_1
else
tmp = t_1
end if
code = tmp
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = h * ((M * M) * 0.25);
double t_1 = (c0 * ((c0 * ((((d * d) / D) / D) / w)) / w)) / h;
double tmp;
if (h <= -3.4e-239) {
tmp = (D / d) * ((D * t_0) / d);
} else if (h <= 3.2e-104) {
tmp = t_1;
} else if (h <= 2.25e+43) {
tmp = ((D * D) / (d / t_0)) / d;
} else {
tmp = t_1;
}
return tmp;
}
def code(c0, w, h, D, d, M): t_0 = h * ((M * M) * 0.25) t_1 = (c0 * ((c0 * ((((d * d) / D) / D) / w)) / w)) / h tmp = 0 if h <= -3.4e-239: tmp = (D / d) * ((D * t_0) / d) elif h <= 3.2e-104: tmp = t_1 elif h <= 2.25e+43: tmp = ((D * D) / (d / t_0)) / d else: tmp = t_1 return tmp
function code(c0, w, h, D, d, M) t_0 = Float64(h * Float64(Float64(M * M) * 0.25)) t_1 = Float64(Float64(c0 * Float64(Float64(c0 * Float64(Float64(Float64(Float64(d * d) / D) / D) / w)) / w)) / h) tmp = 0.0 if (h <= -3.4e-239) tmp = Float64(Float64(D / d) * Float64(Float64(D * t_0) / d)); elseif (h <= 3.2e-104) tmp = t_1; elseif (h <= 2.25e+43) tmp = Float64(Float64(Float64(D * D) / Float64(d / t_0)) / d); else tmp = t_1; end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) t_0 = h * ((M * M) * 0.25); t_1 = (c0 * ((c0 * ((((d * d) / D) / D) / w)) / w)) / h; tmp = 0.0; if (h <= -3.4e-239) tmp = (D / d) * ((D * t_0) / d); elseif (h <= 3.2e-104) tmp = t_1; elseif (h <= 2.25e+43) tmp = ((D * D) / (d / t_0)) / d; else tmp = t_1; end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(h * N[(N[(M * M), $MachinePrecision] * 0.25), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(c0 * N[(N[(c0 * N[(N[(N[(N[(d * d), $MachinePrecision] / D), $MachinePrecision] / D), $MachinePrecision] / w), $MachinePrecision]), $MachinePrecision] / w), $MachinePrecision]), $MachinePrecision] / h), $MachinePrecision]}, If[LessEqual[h, -3.4e-239], N[(N[(D / d), $MachinePrecision] * N[(N[(D * t$95$0), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision], If[LessEqual[h, 3.2e-104], t$95$1, If[LessEqual[h, 2.25e+43], N[(N[(N[(D * D), $MachinePrecision] / N[(d / t$95$0), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := h \cdot \left(\left(M \cdot M\right) \cdot 0.25\right)\\
t_1 := \frac{c0 \cdot \frac{c0 \cdot \frac{\frac{\frac{d \cdot d}{D}}{D}}{w}}{w}}{h}\\
\mathbf{if}\;h \leq -3.4 \cdot 10^{-239}:\\
\;\;\;\;\frac{D}{d} \cdot \frac{D \cdot t\_0}{d}\\
\mathbf{elif}\;h \leq 3.2 \cdot 10^{-104}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;h \leq 2.25 \cdot 10^{+43}:\\
\;\;\;\;\frac{\frac{D \cdot D}{\frac{d}{t\_0}}}{d}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if h < -3.4e-239Initial program 21.4%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr22.4%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified11.4%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.2%
Simplified36.2%
associate-*l*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6454.2%
Applied egg-rr54.2%
if -3.4e-239 < h < 3.19999999999999989e-104 or 2.25e43 < h Initial program 34.2%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified30.5%
Taylor expanded in c0 around inf
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6446.2%
Simplified46.2%
associate-/l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6452.2%
Applied egg-rr52.2%
Taylor expanded in d around 0
associate-*r/N/A
unpow2N/A
associate-*l*N/A
associate-/l/N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6453.1%
Simplified53.1%
Taylor expanded in c0 around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
times-fracN/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified56.0%
if 3.19999999999999989e-104 < h < 2.25e43Initial program 12.1%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr12.1%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified21.9%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6437.4%
Simplified37.4%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6451.3%
Applied egg-rr51.3%
Final simplification54.8%
(FPCore (c0 w h D d M) :precision binary64 (if (<= h -2.2e-237) (* (/ D d) (/ (* D (* h (* (* M M) 0.25))) d)) (/ (/ (/ (/ (/ (* c0 (* c0 (* d d))) w) w) h) D) D)))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (h <= -2.2e-237) {
tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d);
} else {
tmp = (((((c0 * (c0 * (d * d))) / w) / w) / h) / D) / D;
}
return tmp;
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
real(8) :: tmp
if (h <= (-2.2d-237)) then
tmp = (d / d_1) * ((d * (h * ((m * m) * 0.25d0))) / d_1)
else
tmp = (((((c0 * (c0 * (d_1 * d_1))) / w) / w) / h) / d) / d
end if
code = tmp
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (h <= -2.2e-237) {
tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d);
} else {
tmp = (((((c0 * (c0 * (d * d))) / w) / w) / h) / D) / D;
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if h <= -2.2e-237: tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d) else: tmp = (((((c0 * (c0 * (d * d))) / w) / w) / h) / D) / D return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (h <= -2.2e-237) tmp = Float64(Float64(D / d) * Float64(Float64(D * Float64(h * Float64(Float64(M * M) * 0.25))) / d)); else tmp = Float64(Float64(Float64(Float64(Float64(Float64(c0 * Float64(c0 * Float64(d * d))) / w) / w) / h) / D) / D); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (h <= -2.2e-237) tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d); else tmp = (((((c0 * (c0 * (d * d))) / w) / w) / h) / D) / D; end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[h, -2.2e-237], N[(N[(D / d), $MachinePrecision] * N[(N[(D * N[(h * N[(N[(M * M), $MachinePrecision] * 0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(N[(c0 * N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / w), $MachinePrecision] / w), $MachinePrecision] / h), $MachinePrecision] / D), $MachinePrecision] / D), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;h \leq -2.2 \cdot 10^{-237}:\\
\;\;\;\;\frac{D}{d} \cdot \frac{D \cdot \left(h \cdot \left(\left(M \cdot M\right) \cdot 0.25\right)\right)}{d}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{\frac{\frac{c0 \cdot \left(c0 \cdot \left(d \cdot d\right)\right)}{w}}{w}}{h}}{D}}{D}\\
\end{array}
\end{array}
if h < -2.19999999999999998e-237Initial program 21.4%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr22.4%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified11.4%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.2%
Simplified36.2%
associate-*l*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6454.2%
Applied egg-rr54.2%
if -2.19999999999999998e-237 < h Initial program 30.5%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified26.8%
Taylor expanded in c0 around inf
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6440.5%
Simplified40.5%
associate-/l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6451.5%
Applied egg-rr51.5%
Taylor expanded in c0 around 0
*-commutativeN/A
times-fracN/A
unpow2N/A
associate-/r*N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified52.0%
Final simplification52.8%
(FPCore (c0 w h D d M) :precision binary64 (if (<= (* d d) 1.55e+308) (* D (/ D (/ (* d d) (* h (* (* M M) 0.25))))) 0.0))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if ((d * d) <= 1.55e+308) {
tmp = D * (D / ((d * d) / (h * ((M * M) * 0.25))));
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
real(8) :: tmp
if ((d_1 * d_1) <= 1.55d+308) then
tmp = d * (d / ((d_1 * d_1) / (h * ((m * m) * 0.25d0))))
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if ((d * d) <= 1.55e+308) {
tmp = D * (D / ((d * d) / (h * ((M * M) * 0.25))));
} else {
tmp = 0.0;
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if (d * d) <= 1.55e+308: tmp = D * (D / ((d * d) / (h * ((M * M) * 0.25)))) else: tmp = 0.0 return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (Float64(d * d) <= 1.55e+308) tmp = Float64(D * Float64(D / Float64(Float64(d * d) / Float64(h * Float64(Float64(M * M) * 0.25))))); else tmp = 0.0; end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if ((d * d) <= 1.55e+308) tmp = D * (D / ((d * d) / (h * ((M * M) * 0.25)))); else tmp = 0.0; end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[N[(d * d), $MachinePrecision], 1.55e+308], N[(D * N[(D / N[(N[(d * d), $MachinePrecision] / N[(h * N[(N[(M * M), $MachinePrecision] * 0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;d \cdot d \leq 1.55 \cdot 10^{+308}:\\
\;\;\;\;D \cdot \frac{D}{\frac{d \cdot d}{h \cdot \left(\left(M \cdot M\right) \cdot 0.25\right)}}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if (*.f64 d d) < 1.5500000000000001e308Initial program 27.9%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr27.0%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified13.2%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6437.6%
Simplified37.6%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6443.2%
Applied egg-rr43.2%
if 1.5500000000000001e308 < (*.f64 d d) Initial program 26.4%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified22.4%
Taylor expanded in c0 around -inf
associate-*r*N/A
mul-1-negN/A
distribute-lft1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval29.3%
Simplified29.3%
associate-*r*N/A
mul0-rgt31.4%
Applied egg-rr31.4%
Final simplification38.4%
(FPCore (c0 w h D d M) :precision binary64 (if (<= (* D D) 4e-300) 0.0 (* 0.25 (/ (* (* D D) (* h (* M M))) (* d d)))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if ((D * D) <= 4e-300) {
tmp = 0.0;
} else {
tmp = 0.25 * (((D * D) * (h * (M * M))) / (d * d));
}
return tmp;
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
real(8) :: tmp
if ((d * d) <= 4d-300) then
tmp = 0.0d0
else
tmp = 0.25d0 * (((d * d) * (h * (m * m))) / (d_1 * d_1))
end if
code = tmp
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if ((D * D) <= 4e-300) {
tmp = 0.0;
} else {
tmp = 0.25 * (((D * D) * (h * (M * M))) / (d * d));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if (D * D) <= 4e-300: tmp = 0.0 else: tmp = 0.25 * (((D * D) * (h * (M * M))) / (d * d)) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (Float64(D * D) <= 4e-300) tmp = 0.0; else tmp = Float64(0.25 * Float64(Float64(Float64(D * D) * Float64(h * Float64(M * M))) / Float64(d * d))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if ((D * D) <= 4e-300) tmp = 0.0; else tmp = 0.25 * (((D * D) * (h * (M * M))) / (d * d)); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[N[(D * D), $MachinePrecision], 4e-300], 0.0, N[(0.25 * N[(N[(N[(D * D), $MachinePrecision] * N[(h * N[(M * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(d * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;D \cdot D \leq 4 \cdot 10^{-300}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \frac{\left(D \cdot D\right) \cdot \left(h \cdot \left(M \cdot M\right)\right)}{d \cdot d}\\
\end{array}
\end{array}
if (*.f64 D D) < 4.0000000000000001e-300Initial program 27.1%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified25.0%
Taylor expanded in c0 around -inf
associate-*r*N/A
mul-1-negN/A
distribute-lft1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval36.1%
Simplified36.1%
associate-*r*N/A
mul0-rgt43.7%
Applied egg-rr43.7%
if 4.0000000000000001e-300 < (*.f64 D D) Initial program 27.4%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified22.9%
Taylor expanded in c0 around -inf
mul-1-negN/A
neg-sub0N/A
metadata-evalN/A
--lowering--.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
Simplified17.1%
Taylor expanded in c0 around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6431.3%
Simplified31.3%
(FPCore (c0 w h D d M) :precision binary64 (* (/ D d) (/ (* D (* h (* (* M M) 0.25))) d)))
double code(double c0, double w, double h, double D, double d, double M) {
return (D / d) * ((D * (h * ((M * M) * 0.25))) / d);
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
code = (d / d_1) * ((d * (h * ((m * m) * 0.25d0))) / d_1)
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
return (D / d) * ((D * (h * ((M * M) * 0.25))) / d);
}
def code(c0, w, h, D, d, M): return (D / d) * ((D * (h * ((M * M) * 0.25))) / d)
function code(c0, w, h, D, d, M) return Float64(Float64(D / d) * Float64(Float64(D * Float64(h * Float64(Float64(M * M) * 0.25))) / d)) end
function tmp = code(c0, w, h, D, d, M) tmp = (D / d) * ((D * (h * ((M * M) * 0.25))) / d); end
code[c0_, w_, h_, D_, d_, M_] := N[(N[(D / d), $MachinePrecision] * N[(N[(D * N[(h * N[(N[(M * M), $MachinePrecision] * 0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{D}{d} \cdot \frac{D \cdot \left(h \cdot \left(\left(M \cdot M\right) \cdot 0.25\right)\right)}{d}
\end{array}
Initial program 27.3%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr26.7%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
Simplified13.8%
Taylor expanded in c0 around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.8%
Simplified32.8%
associate-*l*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6446.4%
Applied egg-rr46.4%
Final simplification46.4%
(FPCore (c0 w h D d M) :precision binary64 0.0)
double code(double c0, double w, double h, double D, double d, double M) {
return 0.0;
}
real(8) function code(c0, w, h, d, d_1, m)
real(8), intent (in) :: c0
real(8), intent (in) :: w
real(8), intent (in) :: h
real(8), intent (in) :: d
real(8), intent (in) :: d_1
real(8), intent (in) :: m
code = 0.0d0
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
return 0.0;
}
def code(c0, w, h, D, d, M): return 0.0
function code(c0, w, h, D, d, M) return 0.0 end
function tmp = code(c0, w, h, D, d, M) tmp = 0.0; end
code[c0_, w_, h_, D_, d_, M_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 27.3%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified23.7%
Taylor expanded in c0 around -inf
associate-*r*N/A
mul-1-negN/A
distribute-lft1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval25.9%
Simplified25.9%
associate-*r*N/A
mul0-rgt29.0%
Applied egg-rr29.0%
herbie shell --seed 2024163
(FPCore (c0 w h D d M)
:name "Henrywood and Agarwal, Equation (13)"
:precision binary64
(* (/ c0 (* 2.0 w)) (+ (/ (* c0 (* d d)) (* (* w h) (* D D))) (sqrt (- (* (/ (* c0 (* d d)) (* (* w h) (* D D))) (/ (* c0 (* d d)) (* (* w h) (* D D)))) (* M M))))))