
(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 12 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 w) D)) (/ (* w (* D (* h 0.5))) (* c0 d)))
0.0)))
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 * w) * D)) / ((w * (D * (h * 0.5))) / (c0 * d));
} else {
tmp = 0.0;
}
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 * w) * D)) / ((w * (D * (h * 0.5))) / (c0 * d));
} else {
tmp = 0.0;
}
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 * w) * D)) / ((w * (D * (h * 0.5))) / (c0 * d)) else: tmp = 0.0 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(c0 * d) / Float64(Float64(2.0 * w) * D)) / Float64(Float64(w * Float64(D * Float64(h * 0.5))) / Float64(c0 * d))); else tmp = 0.0; 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 * w) * D)) / ((w * (D * (h * 0.5))) / (c0 * d)); else tmp = 0.0; 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[(c0 * d), $MachinePrecision] / N[(N[(2.0 * w), $MachinePrecision] * D), $MachinePrecision]), $MachinePrecision] / N[(N[(w * N[(D * N[(h * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(c0 * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)}\\
\mathbf{if}\;\frac{c0}{2 \cdot w} \cdot \left(t\_0 + \sqrt{t\_0 \cdot t\_0 - M \cdot M}\right) \leq \infty:\\
\;\;\;\;\frac{\frac{c0 \cdot d}{\left(2 \cdot w\right) \cdot D}}{\frac{w \cdot \left(D \cdot \left(h \cdot 0.5\right)\right)}{c0 \cdot d}}\\
\mathbf{else}:\\
\;\;\;\;0\\
\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 75.9%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6470.8
Simplified70.8%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6469.7
Applied egg-rr69.7%
associate-/r*N/A
associate-/l/N/A
associate-*r/N/A
associate-/l/N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l/N/A
*-commutativeN/A
associate-*l/N/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr73.0%
associate-*r*N/A
associate-*l/N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
frac-timesN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
associate-*l/N/A
frac-timesN/A
*-commutativeN/A
clear-numN/A
/-lowering-/.f64N/A
Applied egg-rr78.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%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval36.1
Simplified36.1%
associate-*r*N/A
mul0-rgt41.7
Applied egg-rr41.7%
Final simplification53.7%
(FPCore (c0 w h D d M)
:precision binary64
(if (<= M 1.2e-209)
0.0
(/
(* (* c0 d) (/ (* c0 d) (* (* 2.0 w) D)))
(* (* w D) (- 0.0 (* h -0.5))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 1.2e-209) {
tmp = 0.0;
} else {
tmp = ((c0 * d) * ((c0 * d) / ((2.0 * w) * D))) / ((w * D) * (0.0 - (h * -0.5)));
}
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.2d-209) then
tmp = 0.0d0
else
tmp = ((c0 * d_1) * ((c0 * d_1) / ((2.0d0 * w) * d))) / ((w * d) * (0.0d0 - (h * (-0.5d0))))
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.2e-209) {
tmp = 0.0;
} else {
tmp = ((c0 * d) * ((c0 * d) / ((2.0 * w) * D))) / ((w * D) * (0.0 - (h * -0.5)));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 1.2e-209: tmp = 0.0 else: tmp = ((c0 * d) * ((c0 * d) / ((2.0 * w) * D))) / ((w * D) * (0.0 - (h * -0.5))) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 1.2e-209) tmp = 0.0; else tmp = Float64(Float64(Float64(c0 * d) * Float64(Float64(c0 * d) / Float64(Float64(2.0 * w) * D))) / Float64(Float64(w * D) * Float64(0.0 - Float64(h * -0.5)))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 1.2e-209) tmp = 0.0; else tmp = ((c0 * d) * ((c0 * d) / ((2.0 * w) * D))) / ((w * D) * (0.0 - (h * -0.5))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 1.2e-209], 0.0, N[(N[(N[(c0 * d), $MachinePrecision] * N[(N[(c0 * d), $MachinePrecision] / N[(N[(2.0 * w), $MachinePrecision] * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(w * D), $MachinePrecision] * N[(0.0 - N[(h * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 1.2 \cdot 10^{-209}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(c0 \cdot d\right) \cdot \frac{c0 \cdot d}{\left(2 \cdot w\right) \cdot D}}{\left(w \cdot D\right) \cdot \left(0 - h \cdot -0.5\right)}\\
\end{array}
\end{array}
if M < 1.2000000000000001e-209Initial program 29.0%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval29.2
Simplified29.2%
associate-*r*N/A
mul0-rgt33.3
Applied egg-rr33.3%
if 1.2000000000000001e-209 < M Initial program 17.6%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6428.7
Simplified28.7%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6431.0
Applied egg-rr31.0%
associate-/r*N/A
associate-/l/N/A
associate-*r/N/A
associate-/l/N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l/N/A
*-commutativeN/A
associate-*l/N/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr38.2%
associate-*r*N/A
frac-timesN/A
associate-*r/N/A
*-commutativeN/A
frac-2negN/A
/-lowering-/.f64N/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
frac-timesN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
Applied egg-rr43.3%
Final simplification37.2%
(FPCore (c0 w h D d M) :precision binary64 (if (<= M 1.8e-209) 0.0 (/ (* (* c0 d) (/ (* c0 d) (* (* 2.0 w) D))) (* w (* D (* h 0.5))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 1.8e-209) {
tmp = 0.0;
} else {
tmp = ((c0 * d) * ((c0 * d) / ((2.0 * w) * D))) / (w * (D * (h * 0.5)));
}
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.8d-209) then
tmp = 0.0d0
else
tmp = ((c0 * d_1) * ((c0 * d_1) / ((2.0d0 * w) * d))) / (w * (d * (h * 0.5d0)))
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.8e-209) {
tmp = 0.0;
} else {
tmp = ((c0 * d) * ((c0 * d) / ((2.0 * w) * D))) / (w * (D * (h * 0.5)));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 1.8e-209: tmp = 0.0 else: tmp = ((c0 * d) * ((c0 * d) / ((2.0 * w) * D))) / (w * (D * (h * 0.5))) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 1.8e-209) tmp = 0.0; else tmp = Float64(Float64(Float64(c0 * d) * Float64(Float64(c0 * d) / Float64(Float64(2.0 * w) * D))) / Float64(w * Float64(D * Float64(h * 0.5)))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 1.8e-209) tmp = 0.0; else tmp = ((c0 * d) * ((c0 * d) / ((2.0 * w) * D))) / (w * (D * (h * 0.5))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 1.8e-209], 0.0, N[(N[(N[(c0 * d), $MachinePrecision] * N[(N[(c0 * d), $MachinePrecision] / N[(N[(2.0 * w), $MachinePrecision] * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(w * N[(D * N[(h * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 1.8 \cdot 10^{-209}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(c0 \cdot d\right) \cdot \frac{c0 \cdot d}{\left(2 \cdot w\right) \cdot D}}{w \cdot \left(D \cdot \left(h \cdot 0.5\right)\right)}\\
\end{array}
\end{array}
if M < 1.80000000000000008e-209Initial program 29.0%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval29.2
Simplified29.2%
associate-*r*N/A
mul0-rgt33.3
Applied egg-rr33.3%
if 1.80000000000000008e-209 < M Initial program 17.6%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6428.7
Simplified28.7%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6431.0
Applied egg-rr31.0%
associate-/r*N/A
associate-/l/N/A
associate-*r/N/A
associate-/l/N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l/N/A
*-commutativeN/A
associate-*l/N/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr38.2%
associate-*r*N/A
frac-timesN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
frac-timesN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6442.4
Applied egg-rr42.4%
Final simplification36.8%
(FPCore (c0 w h D d M) :precision binary64 (if (<= M 5.1e-219) 0.0 (* (/ (* c0 d) (* D (* w (* D (* h 0.5))))) (* d (/ (* c0 0.5) w)))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 5.1e-219) {
tmp = 0.0;
} else {
tmp = ((c0 * d) / (D * (w * (D * (h * 0.5))))) * (d * ((c0 * 0.5) / w));
}
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 <= 5.1d-219) then
tmp = 0.0d0
else
tmp = ((c0 * d_1) / (d * (w * (d * (h * 0.5d0))))) * (d_1 * ((c0 * 0.5d0) / w))
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 <= 5.1e-219) {
tmp = 0.0;
} else {
tmp = ((c0 * d) / (D * (w * (D * (h * 0.5))))) * (d * ((c0 * 0.5) / w));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 5.1e-219: tmp = 0.0 else: tmp = ((c0 * d) / (D * (w * (D * (h * 0.5))))) * (d * ((c0 * 0.5) / w)) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 5.1e-219) tmp = 0.0; else tmp = Float64(Float64(Float64(c0 * d) / Float64(D * Float64(w * Float64(D * Float64(h * 0.5))))) * Float64(d * Float64(Float64(c0 * 0.5) / w))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 5.1e-219) tmp = 0.0; else tmp = ((c0 * d) / (D * (w * (D * (h * 0.5))))) * (d * ((c0 * 0.5) / w)); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 5.1e-219], 0.0, N[(N[(N[(c0 * d), $MachinePrecision] / N[(D * N[(w * N[(D * N[(h * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(d * N[(N[(c0 * 0.5), $MachinePrecision] / w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 5.1 \cdot 10^{-219}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\frac{c0 \cdot d}{D \cdot \left(w \cdot \left(D \cdot \left(h \cdot 0.5\right)\right)\right)} \cdot \left(d \cdot \frac{c0 \cdot 0.5}{w}\right)\\
\end{array}
\end{array}
if M < 5.0999999999999998e-219Initial program 28.3%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval29.9
Simplified29.9%
associate-*r*N/A
mul0-rgt34.1
Applied egg-rr34.1%
if 5.0999999999999998e-219 < M Initial program 19.0%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6429.6
Simplified29.6%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6431.8
Applied egg-rr31.8%
associate-/r*N/A
associate-/l/N/A
associate-*r/N/A
associate-/l/N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l/N/A
*-commutativeN/A
associate-*l/N/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr38.7%
*-commutativeN/A
associate-*l/N/A
*-commutativeN/A
*-commutativeN/A
associate-*l/N/A
associate-*l*N/A
*-lowering-*.f64N/A
Applied egg-rr41.0%
Final simplification36.8%
(FPCore (c0 w h D d M) :precision binary64 (if (<= M 7.6e-218) 0.0 (* (* d (* c0 (/ (* c0 d) (* D (* w (* D (* h 0.5))))))) (/ 0.5 w))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 7.6e-218) {
tmp = 0.0;
} else {
tmp = (d * (c0 * ((c0 * d) / (D * (w * (D * (h * 0.5))))))) * (0.5 / w);
}
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 <= 7.6d-218) then
tmp = 0.0d0
else
tmp = (d_1 * (c0 * ((c0 * d_1) / (d * (w * (d * (h * 0.5d0))))))) * (0.5d0 / w)
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 <= 7.6e-218) {
tmp = 0.0;
} else {
tmp = (d * (c0 * ((c0 * d) / (D * (w * (D * (h * 0.5))))))) * (0.5 / w);
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 7.6e-218: tmp = 0.0 else: tmp = (d * (c0 * ((c0 * d) / (D * (w * (D * (h * 0.5))))))) * (0.5 / w) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 7.6e-218) tmp = 0.0; else tmp = Float64(Float64(d * Float64(c0 * Float64(Float64(c0 * d) / Float64(D * Float64(w * Float64(D * Float64(h * 0.5))))))) * Float64(0.5 / w)); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 7.6e-218) tmp = 0.0; else tmp = (d * (c0 * ((c0 * d) / (D * (w * (D * (h * 0.5))))))) * (0.5 / w); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 7.6e-218], 0.0, N[(N[(d * N[(c0 * N[(N[(c0 * d), $MachinePrecision] / N[(D * N[(w * N[(D * N[(h * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(0.5 / w), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 7.6 \cdot 10^{-218}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\left(d \cdot \left(c0 \cdot \frac{c0 \cdot d}{D \cdot \left(w \cdot \left(D \cdot \left(h \cdot 0.5\right)\right)\right)}\right)\right) \cdot \frac{0.5}{w}\\
\end{array}
\end{array}
if M < 7.5999999999999997e-218Initial program 28.3%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval29.9
Simplified29.9%
associate-*r*N/A
mul0-rgt34.1
Applied egg-rr34.1%
if 7.5999999999999997e-218 < M Initial program 19.0%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6429.6
Simplified29.6%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6431.8
Applied egg-rr31.8%
associate-/r*N/A
associate-/l/N/A
associate-*r/N/A
associate-/l/N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l/N/A
*-commutativeN/A
associate-*l/N/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
Applied egg-rr38.7%
associate-*l/N/A
*-commutativeN/A
*-commutativeN/A
clear-numN/A
div-invN/A
div-invN/A
*-lowering-*.f64N/A
Applied egg-rr42.7%
Final simplification37.5%
(FPCore (c0 w h D d M) :precision binary64 (if (<= M 9.5e-17) 0.0 (/ (/ (* c0 (/ (* c0 (* d d)) w)) (* w D)) (* h D))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 9.5e-17) {
tmp = 0.0;
} else {
tmp = ((c0 * ((c0 * (d * d)) / w)) / (w * D)) / (h * 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 <= 9.5d-17) then
tmp = 0.0d0
else
tmp = ((c0 * ((c0 * (d_1 * d_1)) / w)) / (w * d)) / (h * 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 <= 9.5e-17) {
tmp = 0.0;
} else {
tmp = ((c0 * ((c0 * (d * d)) / w)) / (w * D)) / (h * D);
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 9.5e-17: tmp = 0.0 else: tmp = ((c0 * ((c0 * (d * d)) / w)) / (w * D)) / (h * D) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 9.5e-17) tmp = 0.0; else tmp = Float64(Float64(Float64(c0 * Float64(Float64(c0 * Float64(d * d)) / w)) / Float64(w * D)) / Float64(h * D)); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 9.5e-17) tmp = 0.0; else tmp = ((c0 * ((c0 * (d * d)) / w)) / (w * D)) / (h * D); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 9.5e-17], 0.0, N[(N[(N[(c0 * N[(N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision] / w), $MachinePrecision]), $MachinePrecision] / N[(w * D), $MachinePrecision]), $MachinePrecision] / N[(h * D), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 9.5 \cdot 10^{-17}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{c0 \cdot \frac{c0 \cdot \left(d \cdot d\right)}{w}}{w \cdot D}}{h \cdot D}\\
\end{array}
\end{array}
if M < 9.50000000000000029e-17Initial program 27.4%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval30.7
Simplified30.7%
associate-*r*N/A
mul0-rgt35.0
Applied egg-rr35.0%
if 9.50000000000000029e-17 < M Initial program 13.6%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.1
Simplified32.1%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6436.4
Applied egg-rr36.4%
Taylor expanded in c0 around 0
/-lowering-/.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6436.0
Simplified36.0%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6438.5
Applied egg-rr38.5%
Final simplification35.7%
(FPCore (c0 w h D d M) :precision binary64 (if (<= M 1.8e-15) 0.0 (/ (/ (* c0 (* c0 (* d d))) (* h (* w D))) (* w D))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 1.8e-15) {
tmp = 0.0;
} else {
tmp = ((c0 * (c0 * (d * d))) / (h * (w * D))) / (w * 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.8d-15) then
tmp = 0.0d0
else
tmp = ((c0 * (c0 * (d_1 * d_1))) / (h * (w * d))) / (w * 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.8e-15) {
tmp = 0.0;
} else {
tmp = ((c0 * (c0 * (d * d))) / (h * (w * D))) / (w * D);
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 1.8e-15: tmp = 0.0 else: tmp = ((c0 * (c0 * (d * d))) / (h * (w * D))) / (w * D) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 1.8e-15) tmp = 0.0; else tmp = Float64(Float64(Float64(c0 * Float64(c0 * Float64(d * d))) / Float64(h * Float64(w * D))) / Float64(w * D)); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 1.8e-15) tmp = 0.0; else tmp = ((c0 * (c0 * (d * d))) / (h * (w * D))) / (w * D); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 1.8e-15], 0.0, N[(N[(N[(c0 * N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(h * N[(w * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(w * D), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 1.8 \cdot 10^{-15}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{c0 \cdot \left(c0 \cdot \left(d \cdot d\right)\right)}{h \cdot \left(w \cdot D\right)}}{w \cdot D}\\
\end{array}
\end{array}
if M < 1.8000000000000001e-15Initial program 27.4%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval30.7
Simplified30.7%
associate-*r*N/A
mul0-rgt35.0
Applied egg-rr35.0%
if 1.8000000000000001e-15 < M Initial program 13.6%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.1
Simplified32.1%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6436.4
Applied egg-rr36.4%
*-commutativeN/A
associate-/r*N/A
frac-timesN/A
*-commutativeN/A
/-lowering-/.f64N/A
Applied egg-rr38.5%
Taylor expanded in c0 around 0
/-lowering-/.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6438.5
Simplified38.5%
Final simplification35.7%
(FPCore (c0 w h D d M) :precision binary64 (if (<= M 1.5e-16) 0.0 (/ (* (* c0 d) (* c0 d)) (* w (* w (* D (* h D)))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 1.5e-16) {
tmp = 0.0;
} else {
tmp = ((c0 * d) * (c0 * d)) / (w * (w * (D * (h * 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.5d-16) then
tmp = 0.0d0
else
tmp = ((c0 * d_1) * (c0 * d_1)) / (w * (w * (d * (h * 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.5e-16) {
tmp = 0.0;
} else {
tmp = ((c0 * d) * (c0 * d)) / (w * (w * (D * (h * D))));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 1.5e-16: tmp = 0.0 else: tmp = ((c0 * d) * (c0 * d)) / (w * (w * (D * (h * D)))) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 1.5e-16) tmp = 0.0; else tmp = Float64(Float64(Float64(c0 * d) * Float64(c0 * d)) / Float64(w * Float64(w * Float64(D * Float64(h * D))))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 1.5e-16) tmp = 0.0; else tmp = ((c0 * d) * (c0 * d)) / (w * (w * (D * (h * D)))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 1.5e-16], 0.0, N[(N[(N[(c0 * d), $MachinePrecision] * N[(c0 * d), $MachinePrecision]), $MachinePrecision] / N[(w * N[(w * N[(D * N[(h * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 1.5 \cdot 10^{-16}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(c0 \cdot d\right) \cdot \left(c0 \cdot d\right)}{w \cdot \left(w \cdot \left(D \cdot \left(h \cdot D\right)\right)\right)}\\
\end{array}
\end{array}
if M < 1.49999999999999997e-16Initial program 27.4%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval30.7
Simplified30.7%
associate-*r*N/A
mul0-rgt35.0
Applied egg-rr35.0%
if 1.49999999999999997e-16 < M Initial program 13.6%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.1
Simplified32.1%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6436.4
Applied egg-rr36.4%
Taylor expanded in c0 around 0
/-lowering-/.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6436.0
Simplified36.0%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6438.4
Applied egg-rr38.4%
Final simplification35.7%
(FPCore (c0 w h D d M) :precision binary64 (if (<= M 9e-16) 0.0 (/ (* c0 (* c0 (* d d))) (* w (* (* w D) (* h D))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 9e-16) {
tmp = 0.0;
} else {
tmp = (c0 * (c0 * (d * d))) / (w * ((w * D) * (h * 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 <= 9d-16) then
tmp = 0.0d0
else
tmp = (c0 * (c0 * (d_1 * d_1))) / (w * ((w * d) * (h * 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 <= 9e-16) {
tmp = 0.0;
} else {
tmp = (c0 * (c0 * (d * d))) / (w * ((w * D) * (h * D)));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 9e-16: tmp = 0.0 else: tmp = (c0 * (c0 * (d * d))) / (w * ((w * D) * (h * D))) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 9e-16) tmp = 0.0; else tmp = Float64(Float64(c0 * Float64(c0 * Float64(d * d))) / Float64(w * Float64(Float64(w * D) * Float64(h * D)))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 9e-16) tmp = 0.0; else tmp = (c0 * (c0 * (d * d))) / (w * ((w * D) * (h * D))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 9e-16], 0.0, N[(N[(c0 * N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(w * N[(N[(w * D), $MachinePrecision] * N[(h * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 9 \cdot 10^{-16}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\frac{c0 \cdot \left(c0 \cdot \left(d \cdot d\right)\right)}{w \cdot \left(\left(w \cdot D\right) \cdot \left(h \cdot D\right)\right)}\\
\end{array}
\end{array}
if M < 9.0000000000000003e-16Initial program 27.4%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval30.7
Simplified30.7%
associate-*r*N/A
mul0-rgt35.0
Applied egg-rr35.0%
if 9.0000000000000003e-16 < M Initial program 13.6%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.1
Simplified32.1%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6436.4
Applied egg-rr36.4%
Taylor expanded in c0 around 0
/-lowering-/.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6436.0
Simplified36.0%
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6437.9
Applied egg-rr37.9%
Final simplification35.6%
(FPCore (c0 w h D d M) :precision binary64 (if (<= M 2.9e-15) 0.0 (/ (* c0 (* c0 (* d d))) (* w (* w (* D (* h D)))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 2.9e-15) {
tmp = 0.0;
} else {
tmp = (c0 * (c0 * (d * d))) / (w * (w * (D * (h * 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.9d-15) then
tmp = 0.0d0
else
tmp = (c0 * (c0 * (d_1 * d_1))) / (w * (w * (d * (h * 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.9e-15) {
tmp = 0.0;
} else {
tmp = (c0 * (c0 * (d * d))) / (w * (w * (D * (h * D))));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 2.9e-15: tmp = 0.0 else: tmp = (c0 * (c0 * (d * d))) / (w * (w * (D * (h * D)))) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 2.9e-15) tmp = 0.0; else tmp = Float64(Float64(c0 * Float64(c0 * Float64(d * d))) / Float64(w * Float64(w * Float64(D * Float64(h * D))))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 2.9e-15) tmp = 0.0; else tmp = (c0 * (c0 * (d * d))) / (w * (w * (D * (h * D)))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 2.9e-15], 0.0, N[(N[(c0 * N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(w * N[(w * N[(D * N[(h * D), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 2.9 \cdot 10^{-15}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\frac{c0 \cdot \left(c0 \cdot \left(d \cdot d\right)\right)}{w \cdot \left(w \cdot \left(D \cdot \left(h \cdot D\right)\right)\right)}\\
\end{array}
\end{array}
if M < 2.90000000000000019e-15Initial program 27.4%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval30.7
Simplified30.7%
associate-*r*N/A
mul0-rgt35.0
Applied egg-rr35.0%
if 2.90000000000000019e-15 < M Initial program 13.6%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.1
Simplified32.1%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6436.4
Applied egg-rr36.4%
Taylor expanded in c0 around 0
/-lowering-/.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6436.0
Simplified36.0%
Final simplification35.2%
(FPCore (c0 w h D d M) :precision binary64 (if (<= M 1.55e-15) 0.0 (* c0 (/ (* c0 (* d d)) (* (* D (* h D)) (* w w))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (M <= 1.55e-15) {
tmp = 0.0;
} else {
tmp = c0 * ((c0 * (d * d)) / ((D * (h * D)) * (w * w)));
}
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.55d-15) then
tmp = 0.0d0
else
tmp = c0 * ((c0 * (d_1 * d_1)) / ((d * (h * d)) * (w * w)))
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.55e-15) {
tmp = 0.0;
} else {
tmp = c0 * ((c0 * (d * d)) / ((D * (h * D)) * (w * w)));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if M <= 1.55e-15: tmp = 0.0 else: tmp = c0 * ((c0 * (d * d)) / ((D * (h * D)) * (w * w))) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (M <= 1.55e-15) tmp = 0.0; else tmp = Float64(c0 * Float64(Float64(c0 * Float64(d * d)) / Float64(Float64(D * Float64(h * D)) * Float64(w * w)))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (M <= 1.55e-15) tmp = 0.0; else tmp = c0 * ((c0 * (d * d)) / ((D * (h * D)) * (w * w))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[M, 1.55e-15], 0.0, N[(c0 * N[(N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision] / N[(N[(D * N[(h * D), $MachinePrecision]), $MachinePrecision] * N[(w * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \leq 1.55 \cdot 10^{-15}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;c0 \cdot \frac{c0 \cdot \left(d \cdot d\right)}{\left(D \cdot \left(h \cdot D\right)\right) \cdot \left(w \cdot w\right)}\\
\end{array}
\end{array}
if M < 1.5499999999999999e-15Initial program 27.4%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval30.7
Simplified30.7%
associate-*r*N/A
mul0-rgt35.0
Applied egg-rr35.0%
if 1.5499999999999999e-15 < M Initial program 13.6%
Taylor expanded in c0 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.1
Simplified32.1%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6436.4
Applied egg-rr36.4%
Taylor expanded in c0 around 0
/-lowering-/.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6436.0
Simplified36.0%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6433.9
Applied egg-rr33.9%
Final simplification34.8%
(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 24.6%
Taylor expanded in c0 around -inf
associate-*r*N/A
distribute-lft1-inN/A
mul-1-negN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-eval29.2
Simplified29.2%
associate-*r*N/A
mul0-rgt33.0
Applied egg-rr33.0%
herbie shell --seed 2024197
(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))))))