
(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 9 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))))
(t_1 (* (/ c0 (* 2.0 w)) (+ t_0 (sqrt (- (* t_0 t_0) (* M M)))))))
(if (<= t_1 INFINITY) t_1 (* 0.25 (/ (* D (/ M d)) (/ d (* D (* h 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));
double t_1 = (c0 / (2.0 * w)) * (t_0 + sqrt(((t_0 * t_0) - (M * M))));
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = 0.25 * ((D * (M / d)) / (d / (D * (h * M))));
}
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 t_1 = (c0 / (2.0 * w)) * (t_0 + Math.sqrt(((t_0 * t_0) - (M * M))));
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = 0.25 * ((D * (M / d)) / (d / (D * (h * M))));
}
return tmp;
}
def code(c0, w, h, D, d, M): t_0 = (c0 * (d * d)) / ((w * h) * (D * D)) t_1 = (c0 / (2.0 * w)) * (t_0 + math.sqrt(((t_0 * t_0) - (M * M)))) tmp = 0 if t_1 <= math.inf: tmp = t_1 else: tmp = 0.25 * ((D * (M / d)) / (d / (D * (h * M)))) 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))) t_1 = Float64(Float64(c0 / Float64(2.0 * w)) * Float64(t_0 + sqrt(Float64(Float64(t_0 * t_0) - Float64(M * M))))) tmp = 0.0 if (t_1 <= Inf) tmp = t_1; else tmp = Float64(0.25 * Float64(Float64(D * Float64(M / d)) / Float64(d / Float64(D * Float64(h * M))))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) t_0 = (c0 * (d * d)) / ((w * h) * (D * D)); t_1 = (c0 / (2.0 * w)) * (t_0 + sqrt(((t_0 * t_0) - (M * M)))); tmp = 0.0; if (t_1 <= Inf) tmp = t_1; else tmp = 0.25 * ((D * (M / d)) / (d / (D * (h * M)))); 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]}, Block[{t$95$1 = 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]}, If[LessEqual[t$95$1, Infinity], t$95$1, N[(0.25 * N[(N[(D * N[(M / d), $MachinePrecision]), $MachinePrecision] / N[(d / N[(D * N[(h * 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)}\\
t_1 := \frac{c0}{2 \cdot w} \cdot \left(t\_0 + \sqrt{t\_0 \cdot t\_0 - M \cdot M}\right)\\
\mathbf{if}\;t\_1 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \frac{D \cdot \frac{M}{d}}{\frac{d}{D \cdot \left(h \cdot M\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 72.7%
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
Simplified1.1%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified21.0%
Taylor expanded in c0 around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6445.0%
Simplified45.0%
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6471.4%
Applied egg-rr71.4%
*-commutativeN/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.0%
Applied egg-rr73.0%
Final simplification72.9%
(FPCore (c0 w h D d M)
:precision binary64
(if (<= (* d d) 5e-182)
(* 0.25 (* D (* D (/ (* M (/ (* h M) d)) d))))
(if (<= (* d d) 5e-105)
(/
(/
(/
c0
(/
w
(+
(* (/ -0.5 (* d d)) (* (* M (* D (* D M))) (/ (* w (* w h)) c0)))
(* (/ 2.0 (* D D)) (/ (* c0 d) (/ h d))))))
w)
2.0)
(if (<= (* d d) INFINITY)
(* 0.25 (* (/ (* M (* h D)) d) (/ (* D M) d)))
(* 0.25 (/ (* M (/ (* D (* D (* h M))) d)) d))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if ((d * d) <= 5e-182) {
tmp = 0.25 * (D * (D * ((M * ((h * M) / d)) / d)));
} else if ((d * d) <= 5e-105) {
tmp = ((c0 / (w / (((-0.5 / (d * d)) * ((M * (D * (D * M))) * ((w * (w * h)) / c0))) + ((2.0 / (D * D)) * ((c0 * d) / (h / d)))))) / w) / 2.0;
} else if ((d * d) <= ((double) INFINITY)) {
tmp = 0.25 * (((M * (h * D)) / d) * ((D * M) / d));
} else {
tmp = 0.25 * ((M * ((D * (D * (h * M))) / d)) / d);
}
return tmp;
}
public static double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if ((d * d) <= 5e-182) {
tmp = 0.25 * (D * (D * ((M * ((h * M) / d)) / d)));
} else if ((d * d) <= 5e-105) {
tmp = ((c0 / (w / (((-0.5 / (d * d)) * ((M * (D * (D * M))) * ((w * (w * h)) / c0))) + ((2.0 / (D * D)) * ((c0 * d) / (h / d)))))) / w) / 2.0;
} else if ((d * d) <= Double.POSITIVE_INFINITY) {
tmp = 0.25 * (((M * (h * D)) / d) * ((D * M) / d));
} else {
tmp = 0.25 * ((M * ((D * (D * (h * M))) / d)) / d);
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if (d * d) <= 5e-182: tmp = 0.25 * (D * (D * ((M * ((h * M) / d)) / d))) elif (d * d) <= 5e-105: tmp = ((c0 / (w / (((-0.5 / (d * d)) * ((M * (D * (D * M))) * ((w * (w * h)) / c0))) + ((2.0 / (D * D)) * ((c0 * d) / (h / d)))))) / w) / 2.0 elif (d * d) <= math.inf: tmp = 0.25 * (((M * (h * D)) / d) * ((D * M) / d)) else: tmp = 0.25 * ((M * ((D * (D * (h * M))) / d)) / d) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (Float64(d * d) <= 5e-182) tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(M * Float64(Float64(h * M) / d)) / d)))); elseif (Float64(d * d) <= 5e-105) tmp = Float64(Float64(Float64(c0 / Float64(w / Float64(Float64(Float64(-0.5 / Float64(d * d)) * Float64(Float64(M * Float64(D * Float64(D * M))) * Float64(Float64(w * Float64(w * h)) / c0))) + Float64(Float64(2.0 / Float64(D * D)) * Float64(Float64(c0 * d) / Float64(h / d)))))) / w) / 2.0); elseif (Float64(d * d) <= Inf) tmp = Float64(0.25 * Float64(Float64(Float64(M * Float64(h * D)) / d) * Float64(Float64(D * M) / d))); else tmp = Float64(0.25 * Float64(Float64(M * Float64(Float64(D * Float64(D * Float64(h * M))) / d)) / d)); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if ((d * d) <= 5e-182) tmp = 0.25 * (D * (D * ((M * ((h * M) / d)) / d))); elseif ((d * d) <= 5e-105) tmp = ((c0 / (w / (((-0.5 / (d * d)) * ((M * (D * (D * M))) * ((w * (w * h)) / c0))) + ((2.0 / (D * D)) * ((c0 * d) / (h / d)))))) / w) / 2.0; elseif ((d * d) <= Inf) tmp = 0.25 * (((M * (h * D)) / d) * ((D * M) / d)); else tmp = 0.25 * ((M * ((D * (D * (h * M))) / d)) / d); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[N[(d * d), $MachinePrecision], 5e-182], N[(0.25 * N[(D * N[(D * N[(N[(M * N[(N[(h * M), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(d * d), $MachinePrecision], 5e-105], N[(N[(N[(c0 / N[(w / N[(N[(N[(-0.5 / N[(d * d), $MachinePrecision]), $MachinePrecision] * N[(N[(M * N[(D * N[(D * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(w * N[(w * h), $MachinePrecision]), $MachinePrecision] / c0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(2.0 / N[(D * D), $MachinePrecision]), $MachinePrecision] * N[(N[(c0 * d), $MachinePrecision] / N[(h / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / w), $MachinePrecision] / 2.0), $MachinePrecision], If[LessEqual[N[(d * d), $MachinePrecision], Infinity], N[(0.25 * N[(N[(N[(M * N[(h * D), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision] * N[(N[(D * M), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(N[(M * N[(N[(D * N[(D * N[(h * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;d \cdot d \leq 5 \cdot 10^{-182}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \frac{M \cdot \frac{h \cdot M}{d}}{d}\right)\right)\\
\mathbf{elif}\;d \cdot d \leq 5 \cdot 10^{-105}:\\
\;\;\;\;\frac{\frac{\frac{c0}{\frac{w}{\frac{-0.5}{d \cdot d} \cdot \left(\left(M \cdot \left(D \cdot \left(D \cdot M\right)\right)\right) \cdot \frac{w \cdot \left(w \cdot h\right)}{c0}\right) + \frac{2}{D \cdot D} \cdot \frac{c0 \cdot d}{\frac{h}{d}}}}}{w}}{2}\\
\mathbf{elif}\;d \cdot d \leq \infty:\\
\;\;\;\;0.25 \cdot \left(\frac{M \cdot \left(h \cdot D\right)}{d} \cdot \frac{D \cdot M}{d}\right)\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \frac{M \cdot \frac{D \cdot \left(D \cdot \left(h \cdot M\right)\right)}{d}}{d}\\
\end{array}
\end{array}
if (*.f64 d d) < 5.00000000000000024e-182Initial program 12.7%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified14.2%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified12.8%
Taylor expanded in c0 around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6422.9%
Simplified22.9%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6461.7%
Applied egg-rr61.7%
if 5.00000000000000024e-182 < (*.f64 d d) < 4.99999999999999963e-105Initial program 46.8%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified46.2%
Taylor expanded in w around 0
/-lowering-/.f64N/A
Simplified53.9%
Applied egg-rr70.0%
if 4.99999999999999963e-105 < (*.f64 d d) < +inf.0Initial program 22.2%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified21.7%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified18.7%
Taylor expanded in c0 around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6440.7%
Simplified40.7%
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6460.4%
Applied egg-rr60.4%
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6459.8%
Applied egg-rr59.8%
if +inf.0 < (*.f64 d d) Initial program 21.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified21.5%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.2%
Simplified36.2%
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6457.3%
Applied egg-rr57.3%
frac-timesN/A
associate-*l*N/A
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6454.6%
Applied egg-rr54.6%
Final simplification60.7%
(FPCore (c0 w h D d M) :precision binary64 (if (<= (* d d) 5e+182) (* 0.25 (* D (* M (* D (/ (/ (* h M) d) d))))) (* 0.25 (/ (* M (/ (* D (* D (* h M))) d)) d))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if ((d * d) <= 5e+182) {
tmp = 0.25 * (D * (M * (D * (((h * M) / d) / d))));
} else {
tmp = 0.25 * ((M * ((D * (D * (h * 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_1 * d_1) <= 5d+182) then
tmp = 0.25d0 * (d * (m * (d * (((h * m) / d_1) / d_1))))
else
tmp = 0.25d0 * ((m * ((d * (d * (h * 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) <= 5e+182) {
tmp = 0.25 * (D * (M * (D * (((h * M) / d) / d))));
} else {
tmp = 0.25 * ((M * ((D * (D * (h * M))) / d)) / d);
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if (d * d) <= 5e+182: tmp = 0.25 * (D * (M * (D * (((h * M) / d) / d)))) else: tmp = 0.25 * ((M * ((D * (D * (h * M))) / d)) / d) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (Float64(d * d) <= 5e+182) tmp = Float64(0.25 * Float64(D * Float64(M * Float64(D * Float64(Float64(Float64(h * M) / d) / d))))); else tmp = Float64(0.25 * Float64(Float64(M * Float64(Float64(D * Float64(D * Float64(h * M))) / d)) / d)); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if ((d * d) <= 5e+182) tmp = 0.25 * (D * (M * (D * (((h * M) / d) / d)))); else tmp = 0.25 * ((M * ((D * (D * (h * M))) / d)) / d); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[N[(d * d), $MachinePrecision], 5e+182], N[(0.25 * N[(D * N[(M * N[(D * N[(N[(N[(h * M), $MachinePrecision] / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(N[(M * N[(N[(D * N[(D * N[(h * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;d \cdot d \leq 5 \cdot 10^{+182}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(M \cdot \left(D \cdot \frac{\frac{h \cdot M}{d}}{d}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \frac{M \cdot \frac{D \cdot \left(D \cdot \left(h \cdot M\right)\right)}{d}}{d}\\
\end{array}
\end{array}
if (*.f64 d d) < 4.99999999999999973e182Initial program 24.2%
*-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
Simplified14.7%
Taylor expanded in c0 around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6433.1%
Simplified33.1%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6455.0%
Applied egg-rr55.0%
*-commutativeN/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6459.5%
Applied egg-rr59.5%
if 4.99999999999999973e182 < (*.f64 d d) Initial program 19.2%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified19.9%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified19.3%
Taylor expanded in c0 around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6439.1%
Simplified39.1%
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6458.1%
Applied egg-rr58.1%
frac-timesN/A
associate-*l*N/A
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6461.1%
Applied egg-rr61.1%
Final simplification60.3%
(FPCore (c0 w h D d M) :precision binary64 (if (<= D 5e-203) (* 0.25 (* D (* D (/ (* M (/ (* h M) d)) d)))) (* 0.25 (* (/ (* D M) d) (/ (* D (* h M)) d)))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (D <= 5e-203) {
tmp = 0.25 * (D * (D * ((M * ((h * M) / d)) / d)));
} else {
tmp = 0.25 * (((D * M) / d) * ((D * (h * M)) / 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 <= 5d-203) then
tmp = 0.25d0 * (d * (d * ((m * ((h * m) / d_1)) / d_1)))
else
tmp = 0.25d0 * (((d * m) / d_1) * ((d * (h * m)) / 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 <= 5e-203) {
tmp = 0.25 * (D * (D * ((M * ((h * M) / d)) / d)));
} else {
tmp = 0.25 * (((D * M) / d) * ((D * (h * M)) / d));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if D <= 5e-203: tmp = 0.25 * (D * (D * ((M * ((h * M) / d)) / d))) else: tmp = 0.25 * (((D * M) / d) * ((D * (h * M)) / d)) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (D <= 5e-203) tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(M * Float64(Float64(h * M) / d)) / d)))); else tmp = Float64(0.25 * Float64(Float64(Float64(D * M) / d) * Float64(Float64(D * Float64(h * M)) / d))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (D <= 5e-203) tmp = 0.25 * (D * (D * ((M * ((h * M) / d)) / d))); else tmp = 0.25 * (((D * M) / d) * ((D * (h * M)) / d)); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[D, 5e-203], N[(0.25 * N[(D * N[(D * N[(N[(M * N[(N[(h * M), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(N[(N[(D * M), $MachinePrecision] / d), $MachinePrecision] * N[(N[(D * N[(h * M), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;D \leq 5 \cdot 10^{-203}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \frac{M \cdot \frac{h \cdot M}{d}}{d}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(\frac{D \cdot M}{d} \cdot \frac{D \cdot \left(h \cdot M\right)}{d}\right)\\
\end{array}
\end{array}
if D < 5.0000000000000002e-203Initial program 19.3%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified18.1%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified19.0%
Taylor expanded in c0 around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6438.9%
Simplified38.9%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6460.5%
Applied egg-rr60.5%
if 5.0000000000000002e-203 < D Initial program 25.9%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified27.9%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified13.4%
Taylor expanded in c0 around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6431.1%
Simplified31.1%
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6455.6%
Applied egg-rr55.6%
Final simplification58.8%
(FPCore (c0 w h D d M) :precision binary64 (if (<= h -3e+65) (* 0.25 (* D (* M (* D (/ (/ (* h M) d) d))))) (* 0.25 (* D (* D (/ (* M (* h (/ M d))) d))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (h <= -3e+65) {
tmp = 0.25 * (D * (M * (D * (((h * M) / d) / d))));
} else {
tmp = 0.25 * (D * (D * ((M * (h * (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 (h <= (-3d+65)) then
tmp = 0.25d0 * (d * (m * (d * (((h * m) / d_1) / d_1))))
else
tmp = 0.25d0 * (d * (d * ((m * (h * (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 (h <= -3e+65) {
tmp = 0.25 * (D * (M * (D * (((h * M) / d) / d))));
} else {
tmp = 0.25 * (D * (D * ((M * (h * (M / d))) / d)));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if h <= -3e+65: tmp = 0.25 * (D * (M * (D * (((h * M) / d) / d)))) else: tmp = 0.25 * (D * (D * ((M * (h * (M / d))) / d))) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (h <= -3e+65) tmp = Float64(0.25 * Float64(D * Float64(M * Float64(D * Float64(Float64(Float64(h * M) / d) / d))))); else tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(M * Float64(h * Float64(M / d))) / d)))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (h <= -3e+65) tmp = 0.25 * (D * (M * (D * (((h * M) / d) / d)))); else tmp = 0.25 * (D * (D * ((M * (h * (M / d))) / d))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[h, -3e+65], N[(0.25 * N[(D * N[(M * N[(D * N[(N[(N[(h * M), $MachinePrecision] / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(D * N[(D * N[(N[(M * N[(h * N[(M / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;h \leq -3 \cdot 10^{+65}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(M \cdot \left(D \cdot \frac{\frac{h \cdot M}{d}}{d}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \frac{M \cdot \left(h \cdot \frac{M}{d}\right)}{d}\right)\right)\\
\end{array}
\end{array}
if h < -3.0000000000000002e65Initial program 18.1%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified23.7%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified18.5%
Taylor expanded in c0 around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.6%
Simplified36.6%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6448.5%
Applied egg-rr48.5%
*-commutativeN/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6462.5%
Applied egg-rr62.5%
if -3.0000000000000002e65 < h Initial program 22.1%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified21.1%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified16.8%
Taylor expanded in c0 around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.2%
Simplified36.2%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6456.6%
Applied egg-rr56.6%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6457.5%
Applied egg-rr57.5%
Final simplification58.1%
(FPCore (c0 w h D d M) :precision binary64 (* 0.25 (/ (* D (/ M d)) (/ d (* D (* h M))))))
double code(double c0, double w, double h, double D, double d, double M) {
return 0.25 * ((D * (M / d)) / (d / (D * (h * 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
code = 0.25d0 * ((d * (m / d_1)) / (d_1 / (d * (h * m))))
end function
public static double code(double c0, double w, double h, double D, double d, double M) {
return 0.25 * ((D * (M / d)) / (d / (D * (h * M))));
}
def code(c0, w, h, D, d, M): return 0.25 * ((D * (M / d)) / (d / (D * (h * M))))
function code(c0, w, h, D, d, M) return Float64(0.25 * Float64(Float64(D * Float64(M / d)) / Float64(d / Float64(D * Float64(h * M))))) end
function tmp = code(c0, w, h, D, d, M) tmp = 0.25 * ((D * (M / d)) / (d / (D * (h * M)))); end
code[c0_, w_, h_, D_, d_, M_] := N[(0.25 * N[(N[(D * N[(M / d), $MachinePrecision]), $MachinePrecision] / N[(d / N[(D * N[(h * M), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.25 \cdot \frac{D \cdot \frac{M}{d}}{\frac{d}{D \cdot \left(h \cdot M\right)}}
\end{array}
Initial program 21.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified21.5%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.2%
Simplified36.2%
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6457.3%
Applied egg-rr57.3%
*-commutativeN/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-*.f6458.6%
Applied egg-rr58.6%
Final simplification58.6%
(FPCore (c0 w h D d M) :precision binary64 (* 0.25 (* D (* D (/ (* M (* h (/ M d))) d)))))
double code(double c0, double w, double h, double D, double d, double M) {
return 0.25 * (D * (D * ((M * (h * (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 * ((m * (h * (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 * ((M * (h * (M / d))) / d)));
}
def code(c0, w, h, D, d, M): return 0.25 * (D * (D * ((M * (h * (M / d))) / d)))
function code(c0, w, h, D, d, M) return Float64(0.25 * Float64(D * Float64(D * Float64(Float64(M * Float64(h * Float64(M / d))) / d)))) end
function tmp = code(c0, w, h, D, d, M) tmp = 0.25 * (D * (D * ((M * (h * (M / d))) / d))); end
code[c0_, w_, h_, D_, d_, M_] := N[(0.25 * N[(D * N[(D * N[(N[(M * N[(h * N[(M / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.25 \cdot \left(D \cdot \left(D \cdot \frac{M \cdot \left(h \cdot \frac{M}{d}\right)}{d}\right)\right)
\end{array}
Initial program 21.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified21.5%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.2%
Simplified36.2%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6455.5%
Applied egg-rr55.5%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6455.9%
Applied egg-rr55.9%
(FPCore (c0 w h D d M) :precision binary64 (* 0.25 (* D (* D (* M (/ (/ (* h M) d) d))))))
double code(double c0, double w, double h, double D, double d, double M) {
return 0.25 * (D * (D * (M * (((h * 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 * (m * (((h * 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 * (M * (((h * M) / d) / d))));
}
def code(c0, w, h, D, d, M): return 0.25 * (D * (D * (M * (((h * M) / d) / d))))
function code(c0, w, h, D, d, M) return Float64(0.25 * Float64(D * Float64(D * Float64(M * Float64(Float64(Float64(h * M) / d) / d))))) end
function tmp = code(c0, w, h, D, d, M) tmp = 0.25 * (D * (D * (M * (((h * M) / d) / d)))); end
code[c0_, w_, h_, D_, d_, M_] := N[(0.25 * N[(D * N[(D * N[(M * N[(N[(N[(h * M), $MachinePrecision] / d), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.25 \cdot \left(D \cdot \left(D \cdot \left(M \cdot \frac{\frac{h \cdot M}{d}}{d}\right)\right)\right)
\end{array}
Initial program 21.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified21.5%
Taylor expanded in c0 around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.2%
Simplified36.2%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6455.5%
Applied egg-rr55.5%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6453.9%
Applied egg-rr53.9%
Final simplification53.9%
(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 21.6%
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified21.5%
Taylor expanded in c0 around -inf
associate-*r*N/A
mul-1-negN/A
distribute-lft1-inN/A
metadata-evalN/A
mul0-lftN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f6430.8%
Simplified30.8%
associate-*r*N/A
mul0-rgt35.4%
Applied egg-rr35.4%
herbie shell --seed 2024149
(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))))))