
(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 10 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 d) -2.0) (/ h d)) (* w D)) (/ D (/ c0 (* w -2.0))))
(* 0.25 (* (/ (* h (* D (* M M))) d) (/ D 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 * d) * -2.0) / (h / d)) / (w * D)) / (D / (c0 / (w * -2.0)));
} else {
tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / 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 * d) * -2.0) / (h / d)) / (w * D)) / (D / (c0 / (w * -2.0)));
} else {
tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / 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 * d) * -2.0) / (h / d)) / (w * D)) / (D / (c0 / (w * -2.0))) else: tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / 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(Float64(Float64(Float64(c0 * d) * -2.0) / Float64(h / d)) / Float64(w * D)) / Float64(D / Float64(c0 / Float64(w * -2.0)))); else tmp = Float64(0.25 * Float64(Float64(Float64(h * Float64(D * Float64(M * M))) / d) * Float64(D / 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 * d) * -2.0) / (h / d)) / (w * D)) / (D / (c0 / (w * -2.0))); else tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / 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[(N[(N[(N[(c0 * d), $MachinePrecision] * -2.0), $MachinePrecision] / N[(h / d), $MachinePrecision]), $MachinePrecision] / N[(w * D), $MachinePrecision]), $MachinePrecision] / N[(D / N[(c0 / N[(w * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(N[(N[(h * N[(D * N[(M * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] * N[(D / d), $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)}\\
\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{\frac{\left(c0 \cdot d\right) \cdot -2}{\frac{h}{d}}}{w \cdot D}}{\frac{D}{\frac{c0}{w \cdot -2}}}\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(\frac{h \cdot \left(D \cdot \left(M \cdot M\right)\right)}{d} \cdot \frac{D}{d}\right)\\
\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 77.7%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified77.5%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6479.6%
Simplified79.6%
associate-*r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr81.8%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6483.9%
Applied egg-rr83.9%
associate-/l/N/A
div-invN/A
clear-numN/A
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
clear-numN/A
frac-timesN/A
frac-2negN/A
*-rgt-identityN/A
clear-numN/A
div-invN/A
distribute-frac-neg2N/A
Applied egg-rr86.9%
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%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified2.0%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified20.8%
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-*.f6451.8%
Simplified51.8%
associate-*l*N/A
times-fracN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6466.7%
Applied egg-rr66.7%
(FPCore (c0 w h D d M)
:precision binary64
(let* ((t_0 (* D (* D 0.25))))
(if (<= d 1.3e-107)
(/ (/ (* (/ d h) (* (* 2.0 (* c0 d)) (/ (/ c0 (* 2.0 w)) w))) D) D)
(if (<= d 1.45e-11)
(* (/ t_0 (/ (* d d) M)) (* h M))
(if (<= d 2.8e+116)
(/ (* (/ (* 2.0 (/ (* c0 (* d d)) h)) (* w D)) (/ c0 2.0)) (* w D))
(if (<= d 1.05e+180)
(* 0.25 (* (/ (* h (* D (* M M))) d) (/ D d)))
(if (<= d 3.2e+279)
(/ (/ (/ (/ (/ (* c0 (* d (* c0 d))) D) D) w) h) w)
(* (/ (* M M) (/ d h)) (/ t_0 d)))))))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = D * (D * 0.25);
double tmp;
if (d <= 1.3e-107) {
tmp = (((d / h) * ((2.0 * (c0 * d)) * ((c0 / (2.0 * w)) / w))) / D) / D;
} else if (d <= 1.45e-11) {
tmp = (t_0 / ((d * d) / M)) * (h * M);
} else if (d <= 2.8e+116) {
tmp = (((2.0 * ((c0 * (d * d)) / h)) / (w * D)) * (c0 / 2.0)) / (w * D);
} else if (d <= 1.05e+180) {
tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d));
} else if (d <= 3.2e+279) {
tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w;
} else {
tmp = ((M * M) / (d / h)) * (t_0 / 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) :: tmp
t_0 = d * (d * 0.25d0)
if (d_1 <= 1.3d-107) then
tmp = (((d_1 / h) * ((2.0d0 * (c0 * d_1)) * ((c0 / (2.0d0 * w)) / w))) / d) / d
else if (d_1 <= 1.45d-11) then
tmp = (t_0 / ((d_1 * d_1) / m)) * (h * m)
else if (d_1 <= 2.8d+116) then
tmp = (((2.0d0 * ((c0 * (d_1 * d_1)) / h)) / (w * d)) * (c0 / 2.0d0)) / (w * d)
else if (d_1 <= 1.05d+180) then
tmp = 0.25d0 * (((h * (d * (m * m))) / d_1) * (d / d_1))
else if (d_1 <= 3.2d+279) then
tmp = (((((c0 * (d_1 * (c0 * d_1))) / d) / d) / w) / h) / w
else
tmp = ((m * m) / (d_1 / h)) * (t_0 / 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 = D * (D * 0.25);
double tmp;
if (d <= 1.3e-107) {
tmp = (((d / h) * ((2.0 * (c0 * d)) * ((c0 / (2.0 * w)) / w))) / D) / D;
} else if (d <= 1.45e-11) {
tmp = (t_0 / ((d * d) / M)) * (h * M);
} else if (d <= 2.8e+116) {
tmp = (((2.0 * ((c0 * (d * d)) / h)) / (w * D)) * (c0 / 2.0)) / (w * D);
} else if (d <= 1.05e+180) {
tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d));
} else if (d <= 3.2e+279) {
tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w;
} else {
tmp = ((M * M) / (d / h)) * (t_0 / d);
}
return tmp;
}
def code(c0, w, h, D, d, M): t_0 = D * (D * 0.25) tmp = 0 if d <= 1.3e-107: tmp = (((d / h) * ((2.0 * (c0 * d)) * ((c0 / (2.0 * w)) / w))) / D) / D elif d <= 1.45e-11: tmp = (t_0 / ((d * d) / M)) * (h * M) elif d <= 2.8e+116: tmp = (((2.0 * ((c0 * (d * d)) / h)) / (w * D)) * (c0 / 2.0)) / (w * D) elif d <= 1.05e+180: tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d)) elif d <= 3.2e+279: tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w else: tmp = ((M * M) / (d / h)) * (t_0 / d) return tmp
function code(c0, w, h, D, d, M) t_0 = Float64(D * Float64(D * 0.25)) tmp = 0.0 if (d <= 1.3e-107) tmp = Float64(Float64(Float64(Float64(d / h) * Float64(Float64(2.0 * Float64(c0 * d)) * Float64(Float64(c0 / Float64(2.0 * w)) / w))) / D) / D); elseif (d <= 1.45e-11) tmp = Float64(Float64(t_0 / Float64(Float64(d * d) / M)) * Float64(h * M)); elseif (d <= 2.8e+116) tmp = Float64(Float64(Float64(Float64(2.0 * Float64(Float64(c0 * Float64(d * d)) / h)) / Float64(w * D)) * Float64(c0 / 2.0)) / Float64(w * D)); elseif (d <= 1.05e+180) tmp = Float64(0.25 * Float64(Float64(Float64(h * Float64(D * Float64(M * M))) / d) * Float64(D / d))); elseif (d <= 3.2e+279) tmp = Float64(Float64(Float64(Float64(Float64(Float64(c0 * Float64(d * Float64(c0 * d))) / D) / D) / w) / h) / w); else tmp = Float64(Float64(Float64(M * M) / Float64(d / h)) * Float64(t_0 / d)); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) t_0 = D * (D * 0.25); tmp = 0.0; if (d <= 1.3e-107) tmp = (((d / h) * ((2.0 * (c0 * d)) * ((c0 / (2.0 * w)) / w))) / D) / D; elseif (d <= 1.45e-11) tmp = (t_0 / ((d * d) / M)) * (h * M); elseif (d <= 2.8e+116) tmp = (((2.0 * ((c0 * (d * d)) / h)) / (w * D)) * (c0 / 2.0)) / (w * D); elseif (d <= 1.05e+180) tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d)); elseif (d <= 3.2e+279) tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w; else tmp = ((M * M) / (d / h)) * (t_0 / d); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(D * N[(D * 0.25), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[d, 1.3e-107], N[(N[(N[(N[(d / h), $MachinePrecision] * N[(N[(2.0 * N[(c0 * d), $MachinePrecision]), $MachinePrecision] * N[(N[(c0 / N[(2.0 * w), $MachinePrecision]), $MachinePrecision] / w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / D), $MachinePrecision] / D), $MachinePrecision], If[LessEqual[d, 1.45e-11], N[(N[(t$95$0 / N[(N[(d * d), $MachinePrecision] / M), $MachinePrecision]), $MachinePrecision] * N[(h * M), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 2.8e+116], N[(N[(N[(N[(2.0 * N[(N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision] / h), $MachinePrecision]), $MachinePrecision] / N[(w * D), $MachinePrecision]), $MachinePrecision] * N[(c0 / 2.0), $MachinePrecision]), $MachinePrecision] / N[(w * D), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 1.05e+180], N[(0.25 * N[(N[(N[(h * N[(D * N[(M * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] * N[(D / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 3.2e+279], N[(N[(N[(N[(N[(N[(c0 * N[(d * N[(c0 * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / D), $MachinePrecision] / D), $MachinePrecision] / w), $MachinePrecision] / h), $MachinePrecision] / w), $MachinePrecision], N[(N[(N[(M * M), $MachinePrecision] / N[(d / h), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 / d), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := D \cdot \left(D \cdot 0.25\right)\\
\mathbf{if}\;d \leq 1.3 \cdot 10^{-107}:\\
\;\;\;\;\frac{\frac{\frac{d}{h} \cdot \left(\left(2 \cdot \left(c0 \cdot d\right)\right) \cdot \frac{\frac{c0}{2 \cdot w}}{w}\right)}{D}}{D}\\
\mathbf{elif}\;d \leq 1.45 \cdot 10^{-11}:\\
\;\;\;\;\frac{t\_0}{\frac{d \cdot d}{M}} \cdot \left(h \cdot M\right)\\
\mathbf{elif}\;d \leq 2.8 \cdot 10^{+116}:\\
\;\;\;\;\frac{\frac{2 \cdot \frac{c0 \cdot \left(d \cdot d\right)}{h}}{w \cdot D} \cdot \frac{c0}{2}}{w \cdot D}\\
\mathbf{elif}\;d \leq 1.05 \cdot 10^{+180}:\\
\;\;\;\;0.25 \cdot \left(\frac{h \cdot \left(D \cdot \left(M \cdot M\right)\right)}{d} \cdot \frac{D}{d}\right)\\
\mathbf{elif}\;d \leq 3.2 \cdot 10^{+279}:\\
\;\;\;\;\frac{\frac{\frac{\frac{\frac{c0 \cdot \left(d \cdot \left(c0 \cdot d\right)\right)}{D}}{D}}{w}}{h}}{w}\\
\mathbf{else}:\\
\;\;\;\;\frac{M \cdot M}{\frac{d}{h}} \cdot \frac{t\_0}{d}\\
\end{array}
\end{array}
if d < 1.3e-107Initial program 26.5%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified26.5%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.9%
Simplified32.9%
associate-*r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr42.1%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6447.9%
Applied egg-rr47.9%
associate-/r/N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6452.0%
Applied egg-rr52.0%
if 1.3e-107 < d < 1.45e-11Initial program 9.5%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified9.5%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified23.2%
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-*.f6464.3%
Simplified64.3%
associate-/l*N/A
associate-*r*N/A
associate-*r/N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.3%
Applied egg-rr73.3%
if 1.45e-11 < d < 2.80000000000000004e116Initial program 42.9%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified45.9%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6463.4%
Simplified63.4%
*-commutativeN/A
associate-/r*N/A
associate-/r*N/A
frac-timesN/A
*-commutativeN/A
/-lowering-/.f64N/A
Applied egg-rr71.8%
if 2.80000000000000004e116 < d < 1.05e180Initial program 15.4%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified15.4%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified23.5%
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-*.f6470.2%
Simplified70.2%
associate-*l*N/A
times-fracN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6462.7%
Applied egg-rr62.7%
if 1.05e180 < d < 3.19999999999999988e279Initial program 33.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified40.5%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.6%
Simplified47.6%
associate-*r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr51.1%
Taylor expanded in c0 around 0
associate-*r*N/A
associate-/r*N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
Simplified64.2%
if 3.19999999999999988e279 < d Initial program 20.0%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified20.0%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified20.2%
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-*.f6470.3%
Simplified70.3%
associate-*r/N/A
associate-*r*N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6471.1%
Applied egg-rr71.1%
(FPCore (c0 w h D d M)
:precision binary64
(if (<= d 3.7e-17)
(* 0.25 (* D (* D (/ (/ (* h (* M M)) d) d))))
(if (<= d 8.5e+114)
(/ (* (/ (* 2.0 (/ (* c0 (* d d)) h)) (* w D)) (/ c0 2.0)) (* w D))
(if (<= d 2.2e+178)
(* 0.25 (* (/ (* h (* D (* M M))) d) (/ D d)))
(if (<= d 5.6e+278)
(/ (/ (/ (/ (/ (* c0 (* d (* c0 d))) D) D) w) h) w)
(* (/ (* M M) (/ d h)) (/ (* D (* D 0.25)) d)))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (d <= 3.7e-17) {
tmp = 0.25 * (D * (D * (((h * (M * M)) / d) / d)));
} else if (d <= 8.5e+114) {
tmp = (((2.0 * ((c0 * (d * d)) / h)) / (w * D)) * (c0 / 2.0)) / (w * D);
} else if (d <= 2.2e+178) {
tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d));
} else if (d <= 5.6e+278) {
tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w;
} else {
tmp = ((M * M) / (d / h)) * ((D * (D * 0.25)) / 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_1 <= 3.7d-17) then
tmp = 0.25d0 * (d * (d * (((h * (m * m)) / d_1) / d_1)))
else if (d_1 <= 8.5d+114) then
tmp = (((2.0d0 * ((c0 * (d_1 * d_1)) / h)) / (w * d)) * (c0 / 2.0d0)) / (w * d)
else if (d_1 <= 2.2d+178) then
tmp = 0.25d0 * (((h * (d * (m * m))) / d_1) * (d / d_1))
else if (d_1 <= 5.6d+278) then
tmp = (((((c0 * (d_1 * (c0 * d_1))) / d) / d) / w) / h) / w
else
tmp = ((m * m) / (d_1 / h)) * ((d * (d * 0.25d0)) / 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 <= 3.7e-17) {
tmp = 0.25 * (D * (D * (((h * (M * M)) / d) / d)));
} else if (d <= 8.5e+114) {
tmp = (((2.0 * ((c0 * (d * d)) / h)) / (w * D)) * (c0 / 2.0)) / (w * D);
} else if (d <= 2.2e+178) {
tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d));
} else if (d <= 5.6e+278) {
tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w;
} else {
tmp = ((M * M) / (d / h)) * ((D * (D * 0.25)) / d);
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if d <= 3.7e-17: tmp = 0.25 * (D * (D * (((h * (M * M)) / d) / d))) elif d <= 8.5e+114: tmp = (((2.0 * ((c0 * (d * d)) / h)) / (w * D)) * (c0 / 2.0)) / (w * D) elif d <= 2.2e+178: tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d)) elif d <= 5.6e+278: tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w else: tmp = ((M * M) / (d / h)) * ((D * (D * 0.25)) / d) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (d <= 3.7e-17) tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(Float64(h * Float64(M * M)) / d) / d)))); elseif (d <= 8.5e+114) tmp = Float64(Float64(Float64(Float64(2.0 * Float64(Float64(c0 * Float64(d * d)) / h)) / Float64(w * D)) * Float64(c0 / 2.0)) / Float64(w * D)); elseif (d <= 2.2e+178) tmp = Float64(0.25 * Float64(Float64(Float64(h * Float64(D * Float64(M * M))) / d) * Float64(D / d))); elseif (d <= 5.6e+278) tmp = Float64(Float64(Float64(Float64(Float64(Float64(c0 * Float64(d * Float64(c0 * d))) / D) / D) / w) / h) / w); else tmp = Float64(Float64(Float64(M * M) / Float64(d / h)) * Float64(Float64(D * Float64(D * 0.25)) / d)); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (d <= 3.7e-17) tmp = 0.25 * (D * (D * (((h * (M * M)) / d) / d))); elseif (d <= 8.5e+114) tmp = (((2.0 * ((c0 * (d * d)) / h)) / (w * D)) * (c0 / 2.0)) / (w * D); elseif (d <= 2.2e+178) tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d)); elseif (d <= 5.6e+278) tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w; else tmp = ((M * M) / (d / h)) * ((D * (D * 0.25)) / d); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[d, 3.7e-17], N[(0.25 * N[(D * N[(D * N[(N[(N[(h * N[(M * M), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 8.5e+114], N[(N[(N[(N[(2.0 * N[(N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision] / h), $MachinePrecision]), $MachinePrecision] / N[(w * D), $MachinePrecision]), $MachinePrecision] * N[(c0 / 2.0), $MachinePrecision]), $MachinePrecision] / N[(w * D), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 2.2e+178], N[(0.25 * N[(N[(N[(h * N[(D * N[(M * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] * N[(D / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 5.6e+278], N[(N[(N[(N[(N[(N[(c0 * N[(d * N[(c0 * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / D), $MachinePrecision] / D), $MachinePrecision] / w), $MachinePrecision] / h), $MachinePrecision] / w), $MachinePrecision], N[(N[(N[(M * M), $MachinePrecision] / N[(d / h), $MachinePrecision]), $MachinePrecision] * N[(N[(D * N[(D * 0.25), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;d \leq 3.7 \cdot 10^{-17}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \frac{\frac{h \cdot \left(M \cdot M\right)}{d}}{d}\right)\right)\\
\mathbf{elif}\;d \leq 8.5 \cdot 10^{+114}:\\
\;\;\;\;\frac{\frac{2 \cdot \frac{c0 \cdot \left(d \cdot d\right)}{h}}{w \cdot D} \cdot \frac{c0}{2}}{w \cdot D}\\
\mathbf{elif}\;d \leq 2.2 \cdot 10^{+178}:\\
\;\;\;\;0.25 \cdot \left(\frac{h \cdot \left(D \cdot \left(M \cdot M\right)\right)}{d} \cdot \frac{D}{d}\right)\\
\mathbf{elif}\;d \leq 5.6 \cdot 10^{+278}:\\
\;\;\;\;\frac{\frac{\frac{\frac{\frac{c0 \cdot \left(d \cdot \left(c0 \cdot d\right)\right)}{D}}{D}}{w}}{h}}{w}\\
\mathbf{else}:\\
\;\;\;\;\frac{M \cdot M}{\frac{d}{h}} \cdot \frac{D \cdot \left(D \cdot 0.25\right)}{d}\\
\end{array}
\end{array}
if d < 3.6999999999999997e-17Initial program 24.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified24.6%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified14.9%
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-*.f6437.6%
Simplified37.6%
Taylor expanded in D around 0
associate-/l*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6449.2%
Simplified49.2%
if 3.6999999999999997e-17 < d < 8.5000000000000001e114Initial program 40.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified43.4%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.0%
Simplified60.0%
*-commutativeN/A
associate-/r*N/A
associate-/r*N/A
frac-timesN/A
*-commutativeN/A
/-lowering-/.f64N/A
Applied egg-rr70.7%
if 8.5000000000000001e114 < d < 2.19999999999999997e178Initial program 15.4%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified15.4%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified23.5%
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-*.f6470.2%
Simplified70.2%
associate-*l*N/A
times-fracN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6462.7%
Applied egg-rr62.7%
if 2.19999999999999997e178 < d < 5.6000000000000002e278Initial program 33.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified40.5%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.6%
Simplified47.6%
associate-*r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr51.1%
Taylor expanded in c0 around 0
associate-*r*N/A
associate-/r*N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
Simplified64.2%
if 5.6000000000000002e278 < d Initial program 20.0%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified20.0%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified20.2%
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-*.f6470.3%
Simplified70.3%
associate-*r/N/A
associate-*r*N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6471.1%
Applied egg-rr71.1%
(FPCore (c0 w h D d M)
:precision binary64
(if (<= d 1.04e-16)
(* 0.25 (* D (* D (/ (/ (* h (* M M)) d) d))))
(if (<= d 6.6e+115)
(* (/ c0 (* 2.0 w)) (/ (/ (/ (* 2.0 (* c0 d)) (/ h d)) (* w D)) D))
(if (<= d 2.4e+178)
(* 0.25 (* (/ (* h (* D (* M M))) d) (/ D d)))
(if (<= d 3.5e+279)
(/ (/ (/ (/ (/ (* c0 (* d (* c0 d))) D) D) w) h) w)
(* (/ (* M M) (/ d h)) (/ (* D (* D 0.25)) d)))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (d <= 1.04e-16) {
tmp = 0.25 * (D * (D * (((h * (M * M)) / d) / d)));
} else if (d <= 6.6e+115) {
tmp = (c0 / (2.0 * w)) * ((((2.0 * (c0 * d)) / (h / d)) / (w * D)) / D);
} else if (d <= 2.4e+178) {
tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d));
} else if (d <= 3.5e+279) {
tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w;
} else {
tmp = ((M * M) / (d / h)) * ((D * (D * 0.25)) / 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_1 <= 1.04d-16) then
tmp = 0.25d0 * (d * (d * (((h * (m * m)) / d_1) / d_1)))
else if (d_1 <= 6.6d+115) then
tmp = (c0 / (2.0d0 * w)) * ((((2.0d0 * (c0 * d_1)) / (h / d_1)) / (w * d)) / d)
else if (d_1 <= 2.4d+178) then
tmp = 0.25d0 * (((h * (d * (m * m))) / d_1) * (d / d_1))
else if (d_1 <= 3.5d+279) then
tmp = (((((c0 * (d_1 * (c0 * d_1))) / d) / d) / w) / h) / w
else
tmp = ((m * m) / (d_1 / h)) * ((d * (d * 0.25d0)) / 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 <= 1.04e-16) {
tmp = 0.25 * (D * (D * (((h * (M * M)) / d) / d)));
} else if (d <= 6.6e+115) {
tmp = (c0 / (2.0 * w)) * ((((2.0 * (c0 * d)) / (h / d)) / (w * D)) / D);
} else if (d <= 2.4e+178) {
tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d));
} else if (d <= 3.5e+279) {
tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w;
} else {
tmp = ((M * M) / (d / h)) * ((D * (D * 0.25)) / d);
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if d <= 1.04e-16: tmp = 0.25 * (D * (D * (((h * (M * M)) / d) / d))) elif d <= 6.6e+115: tmp = (c0 / (2.0 * w)) * ((((2.0 * (c0 * d)) / (h / d)) / (w * D)) / D) elif d <= 2.4e+178: tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d)) elif d <= 3.5e+279: tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w else: tmp = ((M * M) / (d / h)) * ((D * (D * 0.25)) / d) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (d <= 1.04e-16) tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(Float64(h * Float64(M * M)) / d) / d)))); elseif (d <= 6.6e+115) tmp = Float64(Float64(c0 / Float64(2.0 * w)) * Float64(Float64(Float64(Float64(2.0 * Float64(c0 * d)) / Float64(h / d)) / Float64(w * D)) / D)); elseif (d <= 2.4e+178) tmp = Float64(0.25 * Float64(Float64(Float64(h * Float64(D * Float64(M * M))) / d) * Float64(D / d))); elseif (d <= 3.5e+279) tmp = Float64(Float64(Float64(Float64(Float64(Float64(c0 * Float64(d * Float64(c0 * d))) / D) / D) / w) / h) / w); else tmp = Float64(Float64(Float64(M * M) / Float64(d / h)) * Float64(Float64(D * Float64(D * 0.25)) / d)); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (d <= 1.04e-16) tmp = 0.25 * (D * (D * (((h * (M * M)) / d) / d))); elseif (d <= 6.6e+115) tmp = (c0 / (2.0 * w)) * ((((2.0 * (c0 * d)) / (h / d)) / (w * D)) / D); elseif (d <= 2.4e+178) tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d)); elseif (d <= 3.5e+279) tmp = (((((c0 * (d * (c0 * d))) / D) / D) / w) / h) / w; else tmp = ((M * M) / (d / h)) * ((D * (D * 0.25)) / d); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[d, 1.04e-16], N[(0.25 * N[(D * N[(D * N[(N[(N[(h * N[(M * M), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 6.6e+115], N[(N[(c0 / N[(2.0 * w), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(2.0 * N[(c0 * d), $MachinePrecision]), $MachinePrecision] / N[(h / d), $MachinePrecision]), $MachinePrecision] / N[(w * D), $MachinePrecision]), $MachinePrecision] / D), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 2.4e+178], N[(0.25 * N[(N[(N[(h * N[(D * N[(M * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] * N[(D / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[d, 3.5e+279], N[(N[(N[(N[(N[(N[(c0 * N[(d * N[(c0 * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / D), $MachinePrecision] / D), $MachinePrecision] / w), $MachinePrecision] / h), $MachinePrecision] / w), $MachinePrecision], N[(N[(N[(M * M), $MachinePrecision] / N[(d / h), $MachinePrecision]), $MachinePrecision] * N[(N[(D * N[(D * 0.25), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;d \leq 1.04 \cdot 10^{-16}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \frac{\frac{h \cdot \left(M \cdot M\right)}{d}}{d}\right)\right)\\
\mathbf{elif}\;d \leq 6.6 \cdot 10^{+115}:\\
\;\;\;\;\frac{c0}{2 \cdot w} \cdot \frac{\frac{\frac{2 \cdot \left(c0 \cdot d\right)}{\frac{h}{d}}}{w \cdot D}}{D}\\
\mathbf{elif}\;d \leq 2.4 \cdot 10^{+178}:\\
\;\;\;\;0.25 \cdot \left(\frac{h \cdot \left(D \cdot \left(M \cdot M\right)\right)}{d} \cdot \frac{D}{d}\right)\\
\mathbf{elif}\;d \leq 3.5 \cdot 10^{+279}:\\
\;\;\;\;\frac{\frac{\frac{\frac{\frac{c0 \cdot \left(d \cdot \left(c0 \cdot d\right)\right)}{D}}{D}}{w}}{h}}{w}\\
\mathbf{else}:\\
\;\;\;\;\frac{M \cdot M}{\frac{d}{h}} \cdot \frac{D \cdot \left(D \cdot 0.25\right)}{d}\\
\end{array}
\end{array}
if d < 1.04000000000000001e-16Initial program 24.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified24.6%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified14.9%
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-*.f6437.6%
Simplified37.6%
Taylor expanded in D around 0
associate-/l*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6449.2%
Simplified49.2%
if 1.04000000000000001e-16 < d < 6.6000000000000001e115Initial program 40.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified43.4%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.0%
Simplified60.0%
associate-*r*N/A
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6460.1%
Applied egg-rr60.1%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
associate-*r*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6470.7%
Applied egg-rr70.7%
if 6.6000000000000001e115 < d < 2.4e178Initial program 15.4%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified15.4%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified23.5%
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-*.f6470.2%
Simplified70.2%
associate-*l*N/A
times-fracN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6462.7%
Applied egg-rr62.7%
if 2.4e178 < d < 3.50000000000000001e279Initial program 33.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified40.5%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.6%
Simplified47.6%
associate-*r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr51.1%
Taylor expanded in c0 around 0
associate-*r*N/A
associate-/r*N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
Simplified64.2%
if 3.50000000000000001e279 < d Initial program 20.0%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified20.0%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified20.2%
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-*.f6470.3%
Simplified70.3%
associate-*r/N/A
associate-*r*N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6471.1%
Applied egg-rr71.1%
(FPCore (c0 w h D d M)
:precision binary64
(let* ((t_0 (* h (* M M)))
(t_1 (/ (* d (* d (/ (/ (/ (/ (* c0 c0) w) w) h) D))) D)))
(if (<= M 2.95e-227)
0.0
(if (<= M 1e-185)
t_1
(if (<= M 2.3e-9)
(/ (/ 0.25 d) (/ d (* (* D D) t_0)))
(if (<= M 4.2e+127) t_1 (* 0.25 (* D (* D (/ (/ t_0 d) d))))))))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = h * (M * M);
double t_1 = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D;
double tmp;
if (M <= 2.95e-227) {
tmp = 0.0;
} else if (M <= 1e-185) {
tmp = t_1;
} else if (M <= 2.3e-9) {
tmp = (0.25 / d) / (d / ((D * D) * t_0));
} else if (M <= 4.2e+127) {
tmp = t_1;
} else {
tmp = 0.25 * (D * (D * ((t_0 / 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) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = h * (m * m)
t_1 = (d_1 * (d_1 * (((((c0 * c0) / w) / w) / h) / d))) / d
if (m <= 2.95d-227) then
tmp = 0.0d0
else if (m <= 1d-185) then
tmp = t_1
else if (m <= 2.3d-9) then
tmp = (0.25d0 / d_1) / (d_1 / ((d * d) * t_0))
else if (m <= 4.2d+127) then
tmp = t_1
else
tmp = 0.25d0 * (d * (d * ((t_0 / 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 t_0 = h * (M * M);
double t_1 = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D;
double tmp;
if (M <= 2.95e-227) {
tmp = 0.0;
} else if (M <= 1e-185) {
tmp = t_1;
} else if (M <= 2.3e-9) {
tmp = (0.25 / d) / (d / ((D * D) * t_0));
} else if (M <= 4.2e+127) {
tmp = t_1;
} else {
tmp = 0.25 * (D * (D * ((t_0 / d) / d)));
}
return tmp;
}
def code(c0, w, h, D, d, M): t_0 = h * (M * M) t_1 = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D tmp = 0 if M <= 2.95e-227: tmp = 0.0 elif M <= 1e-185: tmp = t_1 elif M <= 2.3e-9: tmp = (0.25 / d) / (d / ((D * D) * t_0)) elif M <= 4.2e+127: tmp = t_1 else: tmp = 0.25 * (D * (D * ((t_0 / d) / d))) return tmp
function code(c0, w, h, D, d, M) t_0 = Float64(h * Float64(M * M)) t_1 = Float64(Float64(d * Float64(d * Float64(Float64(Float64(Float64(Float64(c0 * c0) / w) / w) / h) / D))) / D) tmp = 0.0 if (M <= 2.95e-227) tmp = 0.0; elseif (M <= 1e-185) tmp = t_1; elseif (M <= 2.3e-9) tmp = Float64(Float64(0.25 / d) / Float64(d / Float64(Float64(D * D) * t_0))); elseif (M <= 4.2e+127) tmp = t_1; else tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(t_0 / d) / d)))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) t_0 = h * (M * M); t_1 = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D; tmp = 0.0; if (M <= 2.95e-227) tmp = 0.0; elseif (M <= 1e-185) tmp = t_1; elseif (M <= 2.3e-9) tmp = (0.25 / d) / (d / ((D * D) * t_0)); elseif (M <= 4.2e+127) tmp = t_1; else tmp = 0.25 * (D * (D * ((t_0 / d) / d))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(h * N[(M * M), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(d * N[(d * N[(N[(N[(N[(N[(c0 * c0), $MachinePrecision] / w), $MachinePrecision] / w), $MachinePrecision] / h), $MachinePrecision] / D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / D), $MachinePrecision]}, If[LessEqual[M, 2.95e-227], 0.0, If[LessEqual[M, 1e-185], t$95$1, If[LessEqual[M, 2.3e-9], N[(N[(0.25 / d), $MachinePrecision] / N[(d / N[(N[(D * D), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[M, 4.2e+127], t$95$1, N[(0.25 * N[(D * N[(D * N[(N[(t$95$0 / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := h \cdot \left(M \cdot M\right)\\
t_1 := \frac{d \cdot \left(d \cdot \frac{\frac{\frac{\frac{c0 \cdot c0}{w}}{w}}{h}}{D}\right)}{D}\\
\mathbf{if}\;M \leq 2.95 \cdot 10^{-227}:\\
\;\;\;\;0\\
\mathbf{elif}\;M \leq 10^{-185}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;M \leq 2.3 \cdot 10^{-9}:\\
\;\;\;\;\frac{\frac{0.25}{d}}{\frac{d}{\left(D \cdot D\right) \cdot t\_0}}\\
\mathbf{elif}\;M \leq 4.2 \cdot 10^{+127}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \frac{\frac{t\_0}{d}}{d}\right)\right)\\
\end{array}
\end{array}
if M < 2.9499999999999998e-227Initial program 28.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified29.9%
Taylor expanded in c0 around -inf
*-commutativeN/A
distribute-lft1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l/N/A
Simplified30.1%
mul0-rgtN/A
div0N/A
metadata-eval38.8%
Applied egg-rr38.8%
if 2.9499999999999998e-227 < M < 9.9999999999999999e-186 or 2.2999999999999999e-9 < M < 4.19999999999999983e127Initial program 40.3%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified43.2%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6443.7%
Simplified43.7%
associate-*r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr47.0%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6452.6%
Applied egg-rr52.6%
Taylor expanded in c0 around 0
unpow2N/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
associate-*r/N/A
/-lowering-/.f64N/A
Simplified53.3%
if 9.9999999999999999e-186 < M < 2.2999999999999999e-9Initial program 23.9%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified23.8%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified26.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-*.f6456.0%
Simplified56.0%
associate-*r/N/A
times-fracN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6457.1%
Applied egg-rr57.1%
if 4.19999999999999983e127 < M Initial program 0.0%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified0.0%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified0.9%
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-*.f6434.3%
Simplified34.3%
Taylor expanded in D around 0
associate-/l*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6462.8%
Simplified62.8%
Final simplification45.9%
(FPCore (c0 w h D d M)
:precision binary64
(if (<= M 1.18e-228)
0.0
(if (<= M 8e-185)
(/ (* d (* d (/ (/ (/ (/ (* c0 c0) w) w) h) D))) D)
(if (<= M 2.3e-10)
(/ (/ 0.25 d) (/ d (* (* D D) (* h (* M M)))))
(/ (* (/ (/ 0.5 w) (/ (/ w (* 2.0 (* c0 d))) (/ d h))) (/ c0 D)) D)))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 1.18e-228) {
tmp = 0.0;
} else if (M <= 8e-185) {
tmp = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D;
} else if (M <= 2.3e-10) {
tmp = (0.25 / d) / (d / ((D * D) * (h * (M * M))));
} else {
tmp = (((0.5 / w) / ((w / (2.0 * (c0 * d))) / (d / h))) * (c0 / 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 (m <= 1.18d-228) then
tmp = 0.0d0
else if (m <= 8d-185) then
tmp = (d_1 * (d_1 * (((((c0 * c0) / w) / w) / h) / d))) / d
else if (m <= 2.3d-10) then
tmp = (0.25d0 / d_1) / (d_1 / ((d * d) * (h * (m * m))))
else
tmp = (((0.5d0 / w) / ((w / (2.0d0 * (c0 * d_1))) / (d_1 / h))) * (c0 / 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 (M <= 1.18e-228) {
tmp = 0.0;
} else if (M <= 8e-185) {
tmp = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D;
} else if (M <= 2.3e-10) {
tmp = (0.25 / d) / (d / ((D * D) * (h * (M * M))));
} else {
tmp = (((0.5 / w) / ((w / (2.0 * (c0 * d))) / (d / h))) * (c0 / D)) / D;
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 1.18e-228: tmp = 0.0 elif M <= 8e-185: tmp = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D elif M <= 2.3e-10: tmp = (0.25 / d) / (d / ((D * D) * (h * (M * M)))) else: tmp = (((0.5 / w) / ((w / (2.0 * (c0 * d))) / (d / h))) * (c0 / D)) / D return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 1.18e-228) tmp = 0.0; elseif (M <= 8e-185) tmp = Float64(Float64(d * Float64(d * Float64(Float64(Float64(Float64(Float64(c0 * c0) / w) / w) / h) / D))) / D); elseif (M <= 2.3e-10) tmp = Float64(Float64(0.25 / d) / Float64(d / Float64(Float64(D * D) * Float64(h * Float64(M * M))))); else tmp = Float64(Float64(Float64(Float64(0.5 / w) / Float64(Float64(w / Float64(2.0 * Float64(c0 * d))) / Float64(d / h))) * Float64(c0 / D)) / D); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 1.18e-228) tmp = 0.0; elseif (M <= 8e-185) tmp = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D; elseif (M <= 2.3e-10) tmp = (0.25 / d) / (d / ((D * D) * (h * (M * M)))); else tmp = (((0.5 / w) / ((w / (2.0 * (c0 * d))) / (d / h))) * (c0 / D)) / D; end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 1.18e-228], 0.0, If[LessEqual[M, 8e-185], N[(N[(d * N[(d * N[(N[(N[(N[(N[(c0 * c0), $MachinePrecision] / w), $MachinePrecision] / w), $MachinePrecision] / h), $MachinePrecision] / D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / D), $MachinePrecision], If[LessEqual[M, 2.3e-10], N[(N[(0.25 / d), $MachinePrecision] / N[(d / N[(N[(D * D), $MachinePrecision] * N[(h * N[(M * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.5 / w), $MachinePrecision] / N[(N[(w / N[(2.0 * N[(c0 * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(d / h), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(c0 / D), $MachinePrecision]), $MachinePrecision] / D), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 1.18 \cdot 10^{-228}:\\
\;\;\;\;0\\
\mathbf{elif}\;M \leq 8 \cdot 10^{-185}:\\
\;\;\;\;\frac{d \cdot \left(d \cdot \frac{\frac{\frac{\frac{c0 \cdot c0}{w}}{w}}{h}}{D}\right)}{D}\\
\mathbf{elif}\;M \leq 2.3 \cdot 10^{-10}:\\
\;\;\;\;\frac{\frac{0.25}{d}}{\frac{d}{\left(D \cdot D\right) \cdot \left(h \cdot \left(M \cdot M\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{0.5}{w}}{\frac{\frac{w}{2 \cdot \left(c0 \cdot d\right)}}{\frac{d}{h}}} \cdot \frac{c0}{D}}{D}\\
\end{array}
\end{array}
if M < 1.1800000000000001e-228Initial program 28.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified29.9%
Taylor expanded in c0 around -inf
*-commutativeN/A
distribute-lft1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l/N/A
Simplified30.1%
mul0-rgtN/A
div0N/A
metadata-eval38.8%
Applied egg-rr38.8%
if 1.1800000000000001e-228 < M < 7.9999999999999999e-185Initial program 41.4%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified47.5%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6448.1%
Simplified48.1%
associate-*r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr53.8%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6459.3%
Applied egg-rr59.3%
Taylor expanded in c0 around 0
unpow2N/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
associate-*r/N/A
/-lowering-/.f64N/A
Simplified71.1%
if 7.9999999999999999e-185 < M < 2.30000000000000007e-10Initial program 23.9%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified23.8%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified26.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-*.f6456.0%
Simplified56.0%
associate-*r/N/A
times-fracN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6457.1%
Applied egg-rr57.1%
if 2.30000000000000007e-10 < M Initial program 19.7%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified19.6%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6440.1%
Simplified40.1%
associate-*r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr46.1%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6451.4%
Applied egg-rr51.4%
div-invN/A
associate-/l/N/A
associate-*l/N/A
div-invN/A
div-invN/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr51.4%
Final simplification46.1%
(FPCore (c0 w h D d M)
:precision binary64
(if (<= M 2.6e-227)
0.0
(if (<= M 5.5e-186)
(/ (* d (* d (/ (/ (/ (/ (* c0 c0) w) w) h) D))) D)
(if (<= M 1.9e-29)
(/ (/ 0.25 d) (/ d (* (* D D) (* h (* M M)))))
(/ (/ (/ (* c0 (/ (/ (* d (* c0 d)) w) w)) h) D) D)))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 2.6e-227) {
tmp = 0.0;
} else if (M <= 5.5e-186) {
tmp = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D;
} else if (M <= 1.9e-29) {
tmp = (0.25 / d) / (d / ((D * D) * (h * (M * M))));
} else {
tmp = (((c0 * (((d * (c0 * 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 (m <= 2.6d-227) then
tmp = 0.0d0
else if (m <= 5.5d-186) then
tmp = (d_1 * (d_1 * (((((c0 * c0) / w) / w) / h) / d))) / d
else if (m <= 1.9d-29) then
tmp = (0.25d0 / d_1) / (d_1 / ((d * d) * (h * (m * m))))
else
tmp = (((c0 * (((d_1 * (c0 * 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 (M <= 2.6e-227) {
tmp = 0.0;
} else if (M <= 5.5e-186) {
tmp = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D;
} else if (M <= 1.9e-29) {
tmp = (0.25 / d) / (d / ((D * D) * (h * (M * M))));
} else {
tmp = (((c0 * (((d * (c0 * d)) / w) / w)) / h) / D) / D;
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 2.6e-227: tmp = 0.0 elif M <= 5.5e-186: tmp = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D elif M <= 1.9e-29: tmp = (0.25 / d) / (d / ((D * D) * (h * (M * M)))) else: tmp = (((c0 * (((d * (c0 * d)) / w) / w)) / h) / D) / D return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 2.6e-227) tmp = 0.0; elseif (M <= 5.5e-186) tmp = Float64(Float64(d * Float64(d * Float64(Float64(Float64(Float64(Float64(c0 * c0) / w) / w) / h) / D))) / D); elseif (M <= 1.9e-29) tmp = Float64(Float64(0.25 / d) / Float64(d / Float64(Float64(D * D) * Float64(h * Float64(M * M))))); else tmp = Float64(Float64(Float64(Float64(c0 * Float64(Float64(Float64(d * Float64(c0 * d)) / w) / w)) / h) / D) / D); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 2.6e-227) tmp = 0.0; elseif (M <= 5.5e-186) tmp = (d * (d * (((((c0 * c0) / w) / w) / h) / D))) / D; elseif (M <= 1.9e-29) tmp = (0.25 / d) / (d / ((D * D) * (h * (M * M)))); else tmp = (((c0 * (((d * (c0 * d)) / w) / w)) / h) / D) / D; end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 2.6e-227], 0.0, If[LessEqual[M, 5.5e-186], N[(N[(d * N[(d * N[(N[(N[(N[(N[(c0 * c0), $MachinePrecision] / w), $MachinePrecision] / w), $MachinePrecision] / h), $MachinePrecision] / D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / D), $MachinePrecision], If[LessEqual[M, 1.9e-29], N[(N[(0.25 / d), $MachinePrecision] / N[(d / N[(N[(D * D), $MachinePrecision] * N[(h * N[(M * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(c0 * N[(N[(N[(d * N[(c0 * d), $MachinePrecision]), $MachinePrecision] / w), $MachinePrecision] / w), $MachinePrecision]), $MachinePrecision] / h), $MachinePrecision] / D), $MachinePrecision] / D), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 2.6 \cdot 10^{-227}:\\
\;\;\;\;0\\
\mathbf{elif}\;M \leq 5.5 \cdot 10^{-186}:\\
\;\;\;\;\frac{d \cdot \left(d \cdot \frac{\frac{\frac{\frac{c0 \cdot c0}{w}}{w}}{h}}{D}\right)}{D}\\
\mathbf{elif}\;M \leq 1.9 \cdot 10^{-29}:\\
\;\;\;\;\frac{\frac{0.25}{d}}{\frac{d}{\left(D \cdot D\right) \cdot \left(h \cdot \left(M \cdot M\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{c0 \cdot \frac{\frac{d \cdot \left(c0 \cdot d\right)}{w}}{w}}{h}}{D}}{D}\\
\end{array}
\end{array}
if M < 2.60000000000000011e-227Initial program 28.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified29.9%
Taylor expanded in c0 around -inf
*-commutativeN/A
distribute-lft1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l/N/A
Simplified30.1%
mul0-rgtN/A
div0N/A
metadata-eval38.8%
Applied egg-rr38.8%
if 2.60000000000000011e-227 < M < 5.5000000000000001e-186Initial program 41.4%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified47.5%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6448.1%
Simplified48.1%
associate-*r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr53.8%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6459.3%
Applied egg-rr59.3%
Taylor expanded in c0 around 0
unpow2N/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
associate-*r/N/A
/-lowering-/.f64N/A
Simplified71.1%
if 5.5000000000000001e-186 < M < 1.89999999999999988e-29Initial program 23.3%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified23.1%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified27.6%
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-*.f6455.2%
Simplified55.2%
associate-*r/N/A
times-fracN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6456.5%
Applied egg-rr56.5%
if 1.89999999999999988e-29 < M Initial program 20.7%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified20.6%
Taylor expanded in c0 around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6439.7%
Simplified39.7%
associate-*r/N/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr45.2%
Taylor expanded in c0 around 0
associate-/l/N/A
/-lowering-/.f64N/A
Simplified46.3%
Final simplification45.1%
(FPCore (c0 w h D d M) :precision binary64 (if (<= M 3.5e+123) (* 0.25 (* (/ (* h (* D (* M M))) d) (/ D d))) (* 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 (M <= 3.5e+123) {
tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d));
} 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 (m <= 3.5d+123) then
tmp = 0.25d0 * (((h * (d * (m * m))) / d_1) * (d / d_1))
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 (M <= 3.5e+123) {
tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d));
} else {
tmp = 0.25 * (D * (D * (((h * (M * M)) / d) / d)));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 3.5e+123: tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d)) 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 (M <= 3.5e+123) tmp = Float64(0.25 * Float64(Float64(Float64(h * Float64(D * Float64(M * M))) / d) * Float64(D / d))); else tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(Float64(h * Float64(M * M)) / d) / d)))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 3.5e+123) tmp = 0.25 * (((h * (D * (M * M))) / d) * (D / d)); 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[M, 3.5e+123], N[(0.25 * N[(N[(N[(h * N[(D * N[(M * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] * N[(D / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(D * N[(D * N[(N[(N[(h * N[(M * M), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 3.5 \cdot 10^{+123}:\\
\;\;\;\;0.25 \cdot \left(\frac{h \cdot \left(D \cdot \left(M \cdot M\right)\right)}{d} \cdot \frac{D}{d}\right)\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \frac{\frac{h \cdot \left(M \cdot M\right)}{d}}{d}\right)\right)\\
\end{array}
\end{array}
if M < 3.5e123Initial program 29.1%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified30.4%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified17.8%
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-*.f6440.6%
Simplified40.6%
associate-*l*N/A
times-fracN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6449.4%
Applied egg-rr49.4%
if 3.5e123 < M Initial program 5.3%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified5.3%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified0.9%
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-*.f6432.5%
Simplified32.5%
Taylor expanded in D around 0
associate-/l*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6464.8%
Simplified64.8%
(FPCore (c0 w h D d M) :precision binary64 (* 0.25 (* D (* D (/ (/ (* h (* M M)) d) d)))))
double code(double c0, double w, double h, double D, double d, double M) {
return 0.25 * (D * (D * (((h * (M * M)) / d) / 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 = 0.25d0 * (d * (d * (((h * (m * m)) / d_1) / d_1)))
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
return 0.25 * (D * (D * (((h * (M * M)) / d) / d)));
}
def code(c0, w, h, D, d, M): return 0.25 * (D * (D * (((h * (M * M)) / d) / d)))
function code(c0, w, h, D, d, M) return Float64(0.25 * Float64(D * Float64(D * Float64(Float64(Float64(h * Float64(M * M)) / d) / d)))) end
function tmp = code(c0, w, h, D, d, M) tmp = 0.25 * (D * (D * (((h * (M * M)) / d) / d))); end
code[c0_, w_, h_, D_, d_, M_] := N[(0.25 * N[(D * N[(D * N[(N[(N[(h * N[(M * M), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.25 \cdot \left(D \cdot \left(D \cdot \frac{\frac{h \cdot \left(M \cdot M\right)}{d}}{d}\right)\right)
\end{array}
Initial program 27.3%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified28.5%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified16.6%
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-*.f6440.0%
Simplified40.0%
Taylor expanded in D around 0
associate-/l*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.3%
Simplified50.3%
(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
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified28.5%
Taylor expanded in c0 around -inf
*-commutativeN/A
distribute-lft1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l/N/A
Simplified27.1%
mul0-rgtN/A
div0N/A
metadata-eval36.6%
Applied egg-rr36.6%
herbie shell --seed 2024138
(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))))))