
(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 8 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))))
(t_1 (/ (* c0 (* d d)) (* (* D D) (* w h)))))
(if (<=
(* (/ c0 (* 2.0 w)) (+ t_1 (sqrt (- (* t_1 t_1) (* M M)))))
INFINITY)
(* h (fma t_0 t_0 (/ (* D (* (* M (* D M)) -0.25)) (* d d))))
(* 0.25 (* (* (* D M) (/ (* D M) d)) (/ h d))))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = (c0 * d) / (D * (w * h));
double t_1 = (c0 * (d * d)) / ((D * D) * (w * h));
double tmp;
if (((c0 / (2.0 * w)) * (t_1 + sqrt(((t_1 * t_1) - (M * M))))) <= ((double) INFINITY)) {
tmp = h * fma(t_0, t_0, ((D * ((M * (D * M)) * -0.25)) / (d * d)));
} else {
tmp = 0.25 * (((D * M) * ((D * M) / d)) * (h / d));
}
return tmp;
}
function code(c0, w, h, D, d, M) t_0 = Float64(Float64(c0 * d) / Float64(D * Float64(w * h))) t_1 = Float64(Float64(c0 * Float64(d * d)) / Float64(Float64(D * D) * Float64(w * h))) tmp = 0.0 if (Float64(Float64(c0 / Float64(2.0 * w)) * Float64(t_1 + sqrt(Float64(Float64(t_1 * t_1) - Float64(M * M))))) <= Inf) tmp = Float64(h * fma(t_0, t_0, Float64(Float64(D * Float64(Float64(M * Float64(D * M)) * -0.25)) / Float64(d * d)))); else tmp = Float64(0.25 * Float64(Float64(Float64(D * M) * Float64(Float64(D * M) / d)) * Float64(h / d))); end return tmp end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(N[(c0 * d), $MachinePrecision] / N[(D * N[(w * h), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(c0 * N[(d * d), $MachinePrecision]), $MachinePrecision] / N[(N[(D * D), $MachinePrecision] * N[(w * h), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(c0 / N[(2.0 * w), $MachinePrecision]), $MachinePrecision] * N[(t$95$1 + N[Sqrt[N[(N[(t$95$1 * t$95$1), $MachinePrecision] - N[(M * M), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(h * N[(t$95$0 * t$95$0 + N[(N[(D * N[(N[(M * N[(D * M), $MachinePrecision]), $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision] / N[(d * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(N[(N[(D * M), $MachinePrecision] * N[(N[(D * M), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision] * N[(h / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{c0 \cdot d}{D \cdot \left(w \cdot h\right)}\\
t_1 := \frac{c0 \cdot \left(d \cdot d\right)}{\left(D \cdot D\right) \cdot \left(w \cdot h\right)}\\
\mathbf{if}\;\frac{c0}{2 \cdot w} \cdot \left(t\_1 + \sqrt{t\_1 \cdot t\_1 - M \cdot M}\right) \leq \infty:\\
\;\;\;\;h \cdot \mathsf{fma}\left(t\_0, t\_0, \frac{D \cdot \left(\left(M \cdot \left(D \cdot M\right)\right) \cdot -0.25\right)}{d \cdot d}\right)\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(\left(\left(D \cdot M\right) \cdot \frac{D \cdot M}{d}\right) \cdot \frac{h}{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 79.1%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified70.0%
Taylor expanded in h around 0
/-lowering-/.f64N/A
Simplified63.2%
Taylor expanded in h around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified68.5%
+-commutativeN/A
unswap-sqrN/A
unswap-sqrN/A
times-fracN/A
fma-defineN/A
fma-lowering-fma.f64N/A
Applied egg-rr84.2%
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
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified0.0%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified19.9%
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-*.f6447.7%
Simplified47.7%
associate-*r*N/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6461.4%
Applied egg-rr61.4%
associate-*r*N/A
unswap-sqrN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6469.4%
Applied egg-rr69.4%
Final simplification73.9%
(FPCore (c0 w h D d M)
:precision binary64
(let* ((t_0 (* D (* w h))))
(if (<= c0 6.2e-127)
(* 0.25 (* D (* (/ M d) (* D (/ M (/ d h))))))
(if (<= c0 6.6e+25)
(*
h
(+
(/ (* D (* (* M (* D M)) -0.25)) (* d d))
(/ (* d (* d (* c0 c0))) (* t_0 t_0))))
(* 0.25 (* (* (* D M) (/ (* D M) d)) (/ h d)))))))
double code(double c0, double w, double h, double D, double d, double M) {
double t_0 = D * (w * h);
double tmp;
if (c0 <= 6.2e-127) {
tmp = 0.25 * (D * ((M / d) * (D * (M / (d / h)))));
} else if (c0 <= 6.6e+25) {
tmp = h * (((D * ((M * (D * M)) * -0.25)) / (d * d)) + ((d * (d * (c0 * c0))) / (t_0 * t_0)));
} else {
tmp = 0.25 * (((D * M) * ((D * M) / 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) :: t_0
real(8) :: tmp
t_0 = d * (w * h)
if (c0 <= 6.2d-127) then
tmp = 0.25d0 * (d * ((m / d_1) * (d * (m / (d_1 / h)))))
else if (c0 <= 6.6d+25) then
tmp = h * (((d * ((m * (d * m)) * (-0.25d0))) / (d_1 * d_1)) + ((d_1 * (d_1 * (c0 * c0))) / (t_0 * t_0)))
else
tmp = 0.25d0 * (((d * m) * ((d * m) / d_1)) * (h / 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 * (w * h);
double tmp;
if (c0 <= 6.2e-127) {
tmp = 0.25 * (D * ((M / d) * (D * (M / (d / h)))));
} else if (c0 <= 6.6e+25) {
tmp = h * (((D * ((M * (D * M)) * -0.25)) / (d * d)) + ((d * (d * (c0 * c0))) / (t_0 * t_0)));
} else {
tmp = 0.25 * (((D * M) * ((D * M) / d)) * (h / d));
}
return tmp;
}
def code(c0, w, h, D, d, M): t_0 = D * (w * h) tmp = 0 if c0 <= 6.2e-127: tmp = 0.25 * (D * ((M / d) * (D * (M / (d / h))))) elif c0 <= 6.6e+25: tmp = h * (((D * ((M * (D * M)) * -0.25)) / (d * d)) + ((d * (d * (c0 * c0))) / (t_0 * t_0))) else: tmp = 0.25 * (((D * M) * ((D * M) / d)) * (h / d)) return tmp
function code(c0, w, h, D, d, M) t_0 = Float64(D * Float64(w * h)) tmp = 0.0 if (c0 <= 6.2e-127) tmp = Float64(0.25 * Float64(D * Float64(Float64(M / d) * Float64(D * Float64(M / Float64(d / h)))))); elseif (c0 <= 6.6e+25) tmp = Float64(h * Float64(Float64(Float64(D * Float64(Float64(M * Float64(D * M)) * -0.25)) / Float64(d * d)) + Float64(Float64(d * Float64(d * Float64(c0 * c0))) / Float64(t_0 * t_0)))); else tmp = Float64(0.25 * Float64(Float64(Float64(D * M) * Float64(Float64(D * M) / d)) * Float64(h / d))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) t_0 = D * (w * h); tmp = 0.0; if (c0 <= 6.2e-127) tmp = 0.25 * (D * ((M / d) * (D * (M / (d / h))))); elseif (c0 <= 6.6e+25) tmp = h * (((D * ((M * (D * M)) * -0.25)) / (d * d)) + ((d * (d * (c0 * c0))) / (t_0 * t_0))); else tmp = 0.25 * (((D * M) * ((D * M) / d)) * (h / d)); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := Block[{t$95$0 = N[(D * N[(w * h), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c0, 6.2e-127], N[(0.25 * N[(D * N[(N[(M / d), $MachinePrecision] * N[(D * N[(M / N[(d / h), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c0, 6.6e+25], N[(h * N[(N[(N[(D * N[(N[(M * N[(D * M), $MachinePrecision]), $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision] / N[(d * d), $MachinePrecision]), $MachinePrecision] + N[(N[(d * N[(d * N[(c0 * c0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(N[(N[(D * M), $MachinePrecision] * N[(N[(D * M), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision] * N[(h / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := D \cdot \left(w \cdot h\right)\\
\mathbf{if}\;c0 \leq 6.2 \cdot 10^{-127}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(\frac{M}{d} \cdot \left(D \cdot \frac{M}{\frac{d}{h}}\right)\right)\right)\\
\mathbf{elif}\;c0 \leq 6.6 \cdot 10^{+25}:\\
\;\;\;\;h \cdot \left(\frac{D \cdot \left(\left(M \cdot \left(D \cdot M\right)\right) \cdot -0.25\right)}{d \cdot d} + \frac{d \cdot \left(d \cdot \left(c0 \cdot c0\right)\right)}{t\_0 \cdot t\_0}\right)\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(\left(\left(D \cdot M\right) \cdot \frac{D \cdot M}{d}\right) \cdot \frac{h}{d}\right)\\
\end{array}
\end{array}
if c0 < 6.2e-127Initial program 22.1%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified20.3%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified16.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-*.f6442.4%
Simplified42.4%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6450.3%
Applied egg-rr50.3%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6459.0%
Applied egg-rr59.0%
if 6.2e-127 < c0 < 6.6000000000000002e25Initial program 38.8%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified35.1%
Taylor expanded in h around 0
/-lowering-/.f64N/A
Simplified46.4%
Taylor expanded in h around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified46.3%
div-invN/A
fma-defineN/A
frac-2negN/A
distribute-frac-negN/A
fmm-undefN/A
div-invN/A
--lowering--.f64N/A
Applied egg-rr65.7%
if 6.6000000000000002e25 < c0 Initial program 23.1%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified18.7%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/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.8%
Simplified33.8%
associate-*r*N/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6449.1%
Applied egg-rr49.1%
associate-*r*N/A
unswap-sqrN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6456.7%
Applied egg-rr56.7%
Final simplification59.0%
(FPCore (c0 w h D d M)
:precision binary64
(if (<= (* M M) 1e-267)
0.0
(if (<= (* M M) 5e+116)
(* 0.25 (* D (* D (* (* M M) (/ h (* 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 ((M * M) <= 1e-267) {
tmp = 0.0;
} else if ((M * M) <= 5e+116) {
tmp = 0.25 * (D * (D * ((M * M) * (h / (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 ((m * m) <= 1d-267) then
tmp = 0.0d0
else if ((m * m) <= 5d+116) then
tmp = 0.25d0 * (d * (d * ((m * m) * (h / (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 ((M * M) <= 1e-267) {
tmp = 0.0;
} else if ((M * M) <= 5e+116) {
tmp = 0.25 * (D * (D * ((M * M) * (h / (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 (M * M) <= 1e-267: tmp = 0.0 elif (M * M) <= 5e+116: tmp = 0.25 * (D * (D * ((M * M) * (h / (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 (Float64(M * M) <= 1e-267) tmp = 0.0; elseif (Float64(M * M) <= 5e+116) tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(M * M) * Float64(h / Float64(d * d)))))); else tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(M * Float64(h * M)) / Float64(d * d))))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if ((M * M) <= 1e-267) tmp = 0.0; elseif ((M * M) <= 5e+116) tmp = 0.25 * (D * (D * ((M * M) * (h / (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[N[(M * M), $MachinePrecision], 1e-267], 0.0, If[LessEqual[N[(M * M), $MachinePrecision], 5e+116], N[(0.25 * N[(D * N[(D * N[(N[(M * M), $MachinePrecision] * N[(h / N[(d * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(D * N[(D * N[(N[(M * N[(h * M), $MachinePrecision]), $MachinePrecision] / N[(d * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \cdot M \leq 10^{-267}:\\
\;\;\;\;0\\
\mathbf{elif}\;M \cdot M \leq 5 \cdot 10^{+116}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \left(\left(M \cdot M\right) \cdot \frac{h}{d \cdot d}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \frac{M \cdot \left(h \cdot M\right)}{d \cdot d}\right)\right)\\
\end{array}
\end{array}
if (*.f64 M M) < 9.9999999999999998e-268Initial program 31.7%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified25.0%
Taylor expanded in c0 around -inf
associate-*r*N/A
mul-1-negN/A
distribute-lft1-inN/A
metadata-evalN/A
mul0-lftN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f6444.1%
Simplified44.1%
associate-*r*N/A
mul0-rgt54.5%
Applied egg-rr54.5%
if 9.9999999999999998e-268 < (*.f64 M M) < 5.00000000000000025e116Initial program 26.5%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified25.1%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified20.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-*.f6447.0%
Simplified47.0%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6448.6%
Applied egg-rr48.6%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6455.2%
Applied egg-rr55.2%
if 5.00000000000000025e116 < (*.f64 M M) Initial program 13.4%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified13.4%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified7.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-*.f6421.6%
Simplified21.6%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6442.1%
Applied egg-rr42.1%
Final simplification50.8%
(FPCore (c0 w h D d M) :precision binary64 (if (<= (* M M) 1.2e-266) 0.0 (* 0.25 (* D (* D (* (* M M) (/ h (* d d))))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if ((M * M) <= 1.2e-266) {
tmp = 0.0;
} else {
tmp = 0.25 * (D * (D * ((M * M) * (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 * m) <= 1.2d-266) then
tmp = 0.0d0
else
tmp = 0.25d0 * (d * (d * ((m * m) * (h / (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 * M) <= 1.2e-266) {
tmp = 0.0;
} else {
tmp = 0.25 * (D * (D * ((M * M) * (h / (d * d)))));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if (M * M) <= 1.2e-266: tmp = 0.0 else: tmp = 0.25 * (D * (D * ((M * M) * (h / (d * d))))) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (Float64(M * M) <= 1.2e-266) tmp = 0.0; else tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(M * M) * Float64(h / Float64(d * d)))))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if ((M * M) <= 1.2e-266) tmp = 0.0; else tmp = 0.25 * (D * (D * ((M * M) * (h / (d * d))))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[N[(M * M), $MachinePrecision], 1.2e-266], 0.0, N[(0.25 * N[(D * N[(D * N[(N[(M * M), $MachinePrecision] * N[(h / N[(d * d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;M \cdot M \leq 1.2 \cdot 10^{-266}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \left(\left(M \cdot M\right) \cdot \frac{h}{d \cdot d}\right)\right)\right)\\
\end{array}
\end{array}
if (*.f64 M M) < 1.2e-266Initial program 31.7%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified25.0%
Taylor expanded in c0 around -inf
associate-*r*N/A
mul-1-negN/A
distribute-lft1-inN/A
metadata-evalN/A
mul0-lftN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f6444.1%
Simplified44.1%
associate-*r*N/A
mul0-rgt54.5%
Applied egg-rr54.5%
if 1.2e-266 < (*.f64 M M) Initial program 20.1%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified19.4%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified14.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-*.f6434.5%
Simplified34.5%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6445.4%
Applied egg-rr45.4%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6442.8%
Applied egg-rr42.8%
(FPCore (c0 w h D d M) :precision binary64 (if (<= h 3.3e-163) (* 0.25 (* (* (* D M) (/ (* D M) d)) (/ h d))) (* 0.25 (* D (* D (/ (* (/ M d) (* h M)) d))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (h <= 3.3e-163) {
tmp = 0.25 * (((D * M) * ((D * M) / d)) * (h / d));
} else {
tmp = 0.25 * (D * (D * (((M / 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 (h <= 3.3d-163) then
tmp = 0.25d0 * (((d * m) * ((d * m) / d_1)) * (h / d_1))
else
tmp = 0.25d0 * (d * (d * (((m / d_1) * (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 (h <= 3.3e-163) {
tmp = 0.25 * (((D * M) * ((D * M) / d)) * (h / d));
} else {
tmp = 0.25 * (D * (D * (((M / d) * (h * M)) / d)));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if h <= 3.3e-163: tmp = 0.25 * (((D * M) * ((D * M) / d)) * (h / d)) else: tmp = 0.25 * (D * (D * (((M / d) * (h * M)) / d))) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (h <= 3.3e-163) tmp = Float64(0.25 * Float64(Float64(Float64(D * M) * Float64(Float64(D * M) / d)) * Float64(h / d))); else tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(Float64(M / d) * Float64(h * M)) / d)))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (h <= 3.3e-163) tmp = 0.25 * (((D * M) * ((D * M) / d)) * (h / d)); else tmp = 0.25 * (D * (D * (((M / d) * (h * M)) / d))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[h, 3.3e-163], N[(0.25 * N[(N[(N[(D * M), $MachinePrecision] * N[(N[(D * M), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision] * N[(h / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(D * N[(D * N[(N[(N[(M / d), $MachinePrecision] * N[(h * M), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;h \leq 3.3 \cdot 10^{-163}:\\
\;\;\;\;0.25 \cdot \left(\left(\left(D \cdot M\right) \cdot \frac{D \cdot M}{d}\right) \cdot \frac{h}{d}\right)\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \frac{\frac{M}{d} \cdot \left(h \cdot M\right)}{d}\right)\right)\\
\end{array}
\end{array}
if h < 3.30000000000000001e-163Initial program 23.3%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified20.9%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/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-*.f6437.9%
Simplified37.9%
associate-*r*N/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6449.5%
Applied egg-rr49.5%
associate-*r*N/A
unswap-sqrN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6458.1%
Applied egg-rr58.1%
if 3.30000000000000001e-163 < h Initial program 25.6%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified22.2%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified21.2%
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.9%
Simplified39.9%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6448.4%
Applied egg-rr48.4%
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6455.6%
Applied egg-rr55.6%
Final simplification57.2%
(FPCore (c0 w h D d M) :precision binary64 (if (<= c0 -4e-104) (* 0.25 (* D (* (/ M d) (* D (/ M (/ d h)))))) (* 0.25 (* D (* D (/ (* (/ M d) (* h M)) d))))))
double code(double c0, double w, double h, double D, double d, double M) {
double tmp;
if (c0 <= -4e-104) {
tmp = 0.25 * (D * ((M / d) * (D * (M / (d / h)))));
} else {
tmp = 0.25 * (D * (D * (((M / 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 (c0 <= (-4d-104)) then
tmp = 0.25d0 * (d * ((m / d_1) * (d * (m / (d_1 / h)))))
else
tmp = 0.25d0 * (d * (d * (((m / d_1) * (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 (c0 <= -4e-104) {
tmp = 0.25 * (D * ((M / d) * (D * (M / (d / h)))));
} else {
tmp = 0.25 * (D * (D * (((M / d) * (h * M)) / d)));
}
return tmp;
}
def code(c0, w, h, D, d, M): tmp = 0 if c0 <= -4e-104: tmp = 0.25 * (D * ((M / d) * (D * (M / (d / h))))) else: tmp = 0.25 * (D * (D * (((M / d) * (h * M)) / d))) return tmp
function code(c0, w, h, D, d, M) tmp = 0.0 if (c0 <= -4e-104) tmp = Float64(0.25 * Float64(D * Float64(Float64(M / d) * Float64(D * Float64(M / Float64(d / h)))))); else tmp = Float64(0.25 * Float64(D * Float64(D * Float64(Float64(Float64(M / d) * Float64(h * M)) / d)))); end return tmp end
function tmp_2 = code(c0, w, h, D, d, M) tmp = 0.0; if (c0 <= -4e-104) tmp = 0.25 * (D * ((M / d) * (D * (M / (d / h))))); else tmp = 0.25 * (D * (D * (((M / d) * (h * M)) / d))); end tmp_2 = tmp; end
code[c0_, w_, h_, D_, d_, M_] := If[LessEqual[c0, -4e-104], N[(0.25 * N[(D * N[(N[(M / d), $MachinePrecision] * N[(D * N[(M / N[(d / h), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.25 * N[(D * N[(D * N[(N[(N[(M / d), $MachinePrecision] * N[(h * M), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c0 \leq -4 \cdot 10^{-104}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(\frac{M}{d} \cdot \left(D \cdot \frac{M}{\frac{d}{h}}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;0.25 \cdot \left(D \cdot \left(D \cdot \frac{\frac{M}{d} \cdot \left(h \cdot M\right)}{d}\right)\right)\\
\end{array}
\end{array}
if c0 < -3.99999999999999971e-104Initial program 25.5%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified23.1%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified17.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-*.f6432.0%
Simplified32.0%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6442.3%
Applied egg-rr42.3%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6454.1%
Applied egg-rr54.1%
if -3.99999999999999971e-104 < c0 Initial program 23.4%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified20.5%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified16.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-*.f6441.8%
Simplified41.8%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6448.8%
Applied egg-rr48.8%
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6456.9%
Applied egg-rr56.9%
Final simplification56.0%
(FPCore (c0 w h D d M) :precision binary64 (* 0.25 (* D (* D (/ (* (/ M d) (* h M)) d)))))
double code(double c0, double w, double h, double D, double d, double M) {
return 0.25 * (D * (D * (((M / d) * (h * M)) / 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 / d_1) * (h * m)) / 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 / d) * (h * M)) / d)));
}
def code(c0, w, h, D, d, M): return 0.25 * (D * (D * (((M / d) * (h * M)) / d)))
function code(c0, w, h, D, d, M) return Float64(0.25 * Float64(D * Float64(D * Float64(Float64(Float64(M / d) * Float64(h * M)) / d)))) end
function tmp = code(c0, w, h, D, d, M) tmp = 0.25 * (D * (D * (((M / d) * (h * M)) / d))); end
code[c0_, w_, h_, D_, d_, M_] := N[(0.25 * N[(D * N[(D * N[(N[(N[(M / d), $MachinePrecision] * N[(h * M), $MachinePrecision]), $MachinePrecision] / d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.25 \cdot \left(D \cdot \left(D \cdot \frac{\frac{M}{d} \cdot \left(h \cdot M\right)}{d}\right)\right)
\end{array}
Initial program 24.1%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified21.3%
Taylor expanded in c0 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified16.9%
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.6%
Simplified38.6%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6446.7%
Applied egg-rr46.7%
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6454.1%
Applied egg-rr54.1%
Final simplification54.1%
(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.1%
*-lowering-*.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified21.3%
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.9%
Simplified30.9%
associate-*r*N/A
mul0-rgt36.7%
Applied egg-rr36.7%
herbie shell --seed 2024162
(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))))))