
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) - 1.0));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) - 1.0d0))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) - 1.0));
}
def code(t, l, k): return 2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t), 2.0)) - 1.0))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) - 1.0))) end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) - 1.0)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 16 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) - 1.0));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) - 1.0d0))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) - 1.0));
}
def code(t, l, k): return 2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t), 2.0)) - 1.0))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) - 1.0))) end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) - 1.0)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}
\end{array}
(FPCore (t l k) :precision binary64 (/ (/ l (tan k)) (* (* k (/ (sin k) l)) (* (* k t) 0.5))))
double code(double t, double l, double k) {
return (l / tan(k)) / ((k * (sin(k) / l)) * ((k * t) * 0.5));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = (l / tan(k)) / ((k * (sin(k) / l)) * ((k * t) * 0.5d0))
end function
public static double code(double t, double l, double k) {
return (l / Math.tan(k)) / ((k * (Math.sin(k) / l)) * ((k * t) * 0.5));
}
def code(t, l, k): return (l / math.tan(k)) / ((k * (math.sin(k) / l)) * ((k * t) * 0.5))
function code(t, l, k) return Float64(Float64(l / tan(k)) / Float64(Float64(k * Float64(sin(k) / l)) * Float64(Float64(k * t) * 0.5))) end
function tmp = code(t, l, k) tmp = (l / tan(k)) / ((k * (sin(k) / l)) * ((k * t) * 0.5)); end
code[t_, l_, k_] := N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(N[(k * N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(k * t), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{\ell}{\tan k}}{\left(k \cdot \frac{\sin k}{\ell}\right) \cdot \left(\left(k \cdot t\right) \cdot 0.5\right)}
\end{array}
Initial program 40.2%
Applied egg-rr30.2%
Applied egg-rr87.6%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-/.f64N/A
lift-tan.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-invN/A
lift-/.f64N/A
clear-numN/A
*-commutativeN/A
lift-/.f64N/A
clear-numN/A
div-invN/A
lower-/.f64N/A
Applied egg-rr92.5%
lift-sin.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f6496.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6496.9
Applied egg-rr96.9%
Final simplification96.9%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (* k (* k t))))
(if (<= (* l l) 5e-167)
(/ (/ l k) (* (/ (sin k) l) (* 0.5 t_1)))
(if (<= (* l l) 8e+279)
(* (/ 2.0 (* (sin k) t_1)) (* l (/ l (tan k))))
(/ l (* (tan k) (* (sin k) (/ (* 0.5 (* t (* k k))) l))))))))
double code(double t, double l, double k) {
double t_1 = k * (k * t);
double tmp;
if ((l * l) <= 5e-167) {
tmp = (l / k) / ((sin(k) / l) * (0.5 * t_1));
} else if ((l * l) <= 8e+279) {
tmp = (2.0 / (sin(k) * t_1)) * (l * (l / tan(k)));
} else {
tmp = l / (tan(k) * (sin(k) * ((0.5 * (t * (k * k))) / l)));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = k * (k * t)
if ((l * l) <= 5d-167) then
tmp = (l / k) / ((sin(k) / l) * (0.5d0 * t_1))
else if ((l * l) <= 8d+279) then
tmp = (2.0d0 / (sin(k) * t_1)) * (l * (l / tan(k)))
else
tmp = l / (tan(k) * (sin(k) * ((0.5d0 * (t * (k * k))) / l)))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = k * (k * t);
double tmp;
if ((l * l) <= 5e-167) {
tmp = (l / k) / ((Math.sin(k) / l) * (0.5 * t_1));
} else if ((l * l) <= 8e+279) {
tmp = (2.0 / (Math.sin(k) * t_1)) * (l * (l / Math.tan(k)));
} else {
tmp = l / (Math.tan(k) * (Math.sin(k) * ((0.5 * (t * (k * k))) / l)));
}
return tmp;
}
def code(t, l, k): t_1 = k * (k * t) tmp = 0 if (l * l) <= 5e-167: tmp = (l / k) / ((math.sin(k) / l) * (0.5 * t_1)) elif (l * l) <= 8e+279: tmp = (2.0 / (math.sin(k) * t_1)) * (l * (l / math.tan(k))) else: tmp = l / (math.tan(k) * (math.sin(k) * ((0.5 * (t * (k * k))) / l))) return tmp
function code(t, l, k) t_1 = Float64(k * Float64(k * t)) tmp = 0.0 if (Float64(l * l) <= 5e-167) tmp = Float64(Float64(l / k) / Float64(Float64(sin(k) / l) * Float64(0.5 * t_1))); elseif (Float64(l * l) <= 8e+279) tmp = Float64(Float64(2.0 / Float64(sin(k) * t_1)) * Float64(l * Float64(l / tan(k)))); else tmp = Float64(l / Float64(tan(k) * Float64(sin(k) * Float64(Float64(0.5 * Float64(t * Float64(k * k))) / l)))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = k * (k * t); tmp = 0.0; if ((l * l) <= 5e-167) tmp = (l / k) / ((sin(k) / l) * (0.5 * t_1)); elseif ((l * l) <= 8e+279) tmp = (2.0 / (sin(k) * t_1)) * (l * (l / tan(k))); else tmp = l / (tan(k) * (sin(k) * ((0.5 * (t * (k * k))) / l))); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(l * l), $MachinePrecision], 5e-167], N[(N[(l / k), $MachinePrecision] / N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[(0.5 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 8e+279], N[(N[(2.0 / N[(N[Sin[k], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] * N[(l * N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l / N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(0.5 * N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := k \cdot \left(k \cdot t\right)\\
\mathbf{if}\;\ell \cdot \ell \leq 5 \cdot 10^{-167}:\\
\;\;\;\;\frac{\frac{\ell}{k}}{\frac{\sin k}{\ell} \cdot \left(0.5 \cdot t\_1\right)}\\
\mathbf{elif}\;\ell \cdot \ell \leq 8 \cdot 10^{+279}:\\
\;\;\;\;\frac{2}{\sin k \cdot t\_1} \cdot \left(\ell \cdot \frac{\ell}{\tan k}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{\tan k \cdot \left(\sin k \cdot \frac{0.5 \cdot \left(t \cdot \left(k \cdot k\right)\right)}{\ell}\right)}\\
\end{array}
\end{array}
if (*.f64 l l) < 5.0000000000000002e-167Initial program 25.7%
Applied egg-rr26.8%
Applied egg-rr90.2%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-/.f64N/A
lift-tan.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-invN/A
lift-/.f64N/A
clear-numN/A
*-commutativeN/A
lift-/.f64N/A
clear-numN/A
div-invN/A
lower-/.f64N/A
Applied egg-rr95.2%
Taylor expanded in k around 0
lower-/.f6494.3
Simplified94.3%
if 5.0000000000000002e-167 < (*.f64 l l) < 8.00000000000000046e279Initial program 54.5%
Applied egg-rr36.8%
Applied egg-rr91.4%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
associate-/r/N/A
lift-tan.f64N/A
lift-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
Applied egg-rr96.7%
if 8.00000000000000046e279 < (*.f64 l l) Initial program 40.6%
Applied egg-rr26.1%
Applied egg-rr79.0%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-/.f64N/A
lift-tan.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-invN/A
lift-/.f64N/A
clear-numN/A
*-commutativeN/A
lift-/.f64N/A
clear-numN/A
div-invN/A
lower-/.f64N/A
Applied egg-rr83.2%
lift-tan.f64N/A
lift-sin.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6483.3
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6483.3
Applied egg-rr79.1%
Final simplification91.1%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (* k (* k t))))
(if (<= k 4.5e-6)
(/ (/ l k) (* (/ (sin k) l) (* 0.5 t_1)))
(* (/ 2.0 (* (sin k) t_1)) (* l (/ l (tan k)))))))
double code(double t, double l, double k) {
double t_1 = k * (k * t);
double tmp;
if (k <= 4.5e-6) {
tmp = (l / k) / ((sin(k) / l) * (0.5 * t_1));
} else {
tmp = (2.0 / (sin(k) * t_1)) * (l * (l / tan(k)));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = k * (k * t)
if (k <= 4.5d-6) then
tmp = (l / k) / ((sin(k) / l) * (0.5d0 * t_1))
else
tmp = (2.0d0 / (sin(k) * t_1)) * (l * (l / tan(k)))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = k * (k * t);
double tmp;
if (k <= 4.5e-6) {
tmp = (l / k) / ((Math.sin(k) / l) * (0.5 * t_1));
} else {
tmp = (2.0 / (Math.sin(k) * t_1)) * (l * (l / Math.tan(k)));
}
return tmp;
}
def code(t, l, k): t_1 = k * (k * t) tmp = 0 if k <= 4.5e-6: tmp = (l / k) / ((math.sin(k) / l) * (0.5 * t_1)) else: tmp = (2.0 / (math.sin(k) * t_1)) * (l * (l / math.tan(k))) return tmp
function code(t, l, k) t_1 = Float64(k * Float64(k * t)) tmp = 0.0 if (k <= 4.5e-6) tmp = Float64(Float64(l / k) / Float64(Float64(sin(k) / l) * Float64(0.5 * t_1))); else tmp = Float64(Float64(2.0 / Float64(sin(k) * t_1)) * Float64(l * Float64(l / tan(k)))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = k * (k * t); tmp = 0.0; if (k <= 4.5e-6) tmp = (l / k) / ((sin(k) / l) * (0.5 * t_1)); else tmp = (2.0 / (sin(k) * t_1)) * (l * (l / tan(k))); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, 4.5e-6], N[(N[(l / k), $MachinePrecision] / N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[(0.5 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(N[Sin[k], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] * N[(l * N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := k \cdot \left(k \cdot t\right)\\
\mathbf{if}\;k \leq 4.5 \cdot 10^{-6}:\\
\;\;\;\;\frac{\frac{\ell}{k}}{\frac{\sin k}{\ell} \cdot \left(0.5 \cdot t\_1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\sin k \cdot t\_1} \cdot \left(\ell \cdot \frac{\ell}{\tan k}\right)\\
\end{array}
\end{array}
if k < 4.50000000000000011e-6Initial program 40.4%
Applied egg-rr30.0%
Applied egg-rr88.7%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-/.f64N/A
lift-tan.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-invN/A
lift-/.f64N/A
clear-numN/A
*-commutativeN/A
lift-/.f64N/A
clear-numN/A
div-invN/A
lower-/.f64N/A
Applied egg-rr94.3%
Taylor expanded in k around 0
lower-/.f6483.0
Simplified83.0%
if 4.50000000000000011e-6 < k Initial program 39.2%
Applied egg-rr31.2%
Applied egg-rr83.9%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
associate-/r/N/A
lift-tan.f64N/A
lift-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
Applied egg-rr77.1%
Final simplification81.7%
(FPCore (t l k) :precision binary64 (/ (/ l (tan k)) (* k (* (/ (sin k) l) (* (* k t) 0.5)))))
double code(double t, double l, double k) {
return (l / tan(k)) / (k * ((sin(k) / l) * ((k * t) * 0.5)));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = (l / tan(k)) / (k * ((sin(k) / l) * ((k * t) * 0.5d0)))
end function
public static double code(double t, double l, double k) {
return (l / Math.tan(k)) / (k * ((Math.sin(k) / l) * ((k * t) * 0.5)));
}
def code(t, l, k): return (l / math.tan(k)) / (k * ((math.sin(k) / l) * ((k * t) * 0.5)))
function code(t, l, k) return Float64(Float64(l / tan(k)) / Float64(k * Float64(Float64(sin(k) / l) * Float64(Float64(k * t) * 0.5)))) end
function tmp = code(t, l, k) tmp = (l / tan(k)) / (k * ((sin(k) / l) * ((k * t) * 0.5))); end
code[t_, l_, k_] := N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(k * N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[(N[(k * t), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{\ell}{\tan k}}{k \cdot \left(\frac{\sin k}{\ell} \cdot \left(\left(k \cdot t\right) \cdot 0.5\right)\right)}
\end{array}
Initial program 40.2%
Applied egg-rr30.2%
Applied egg-rr87.6%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-/.f64N/A
lift-tan.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-invN/A
lift-/.f64N/A
clear-numN/A
*-commutativeN/A
lift-/.f64N/A
clear-numN/A
div-invN/A
lower-/.f64N/A
Applied egg-rr92.5%
lift-sin.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f6496.8
lift-*.f64N/A
*-commutativeN/A
lower-*.f6496.8
Applied egg-rr96.8%
Final simplification96.8%
(FPCore (t l k) :precision binary64 (if (<= k 2.5e-6) (/ (/ l k) (* (/ (sin k) l) (* 0.5 (* k (* k t))))) (* (* l l) (/ 2.0 (* (sin k) (* (tan k) (* t (* k k))))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 2.5e-6) {
tmp = (l / k) / ((sin(k) / l) * (0.5 * (k * (k * t))));
} else {
tmp = (l * l) * (2.0 / (sin(k) * (tan(k) * (t * (k * k)))));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 2.5d-6) then
tmp = (l / k) / ((sin(k) / l) * (0.5d0 * (k * (k * t))))
else
tmp = (l * l) * (2.0d0 / (sin(k) * (tan(k) * (t * (k * k)))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 2.5e-6) {
tmp = (l / k) / ((Math.sin(k) / l) * (0.5 * (k * (k * t))));
} else {
tmp = (l * l) * (2.0 / (Math.sin(k) * (Math.tan(k) * (t * (k * k)))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 2.5e-6: tmp = (l / k) / ((math.sin(k) / l) * (0.5 * (k * (k * t)))) else: tmp = (l * l) * (2.0 / (math.sin(k) * (math.tan(k) * (t * (k * k))))) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 2.5e-6) tmp = Float64(Float64(l / k) / Float64(Float64(sin(k) / l) * Float64(0.5 * Float64(k * Float64(k * t))))); else tmp = Float64(Float64(l * l) * Float64(2.0 / Float64(sin(k) * Float64(tan(k) * Float64(t * Float64(k * k)))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 2.5e-6) tmp = (l / k) / ((sin(k) / l) * (0.5 * (k * (k * t)))); else tmp = (l * l) * (2.0 / (sin(k) * (tan(k) * (t * (k * k))))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 2.5e-6], N[(N[(l / k), $MachinePrecision] / N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[(0.5 * N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[Sin[k], $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 2.5 \cdot 10^{-6}:\\
\;\;\;\;\frac{\frac{\ell}{k}}{\frac{\sin k}{\ell} \cdot \left(0.5 \cdot \left(k \cdot \left(k \cdot t\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{\sin k \cdot \left(\tan k \cdot \left(t \cdot \left(k \cdot k\right)\right)\right)}\\
\end{array}
\end{array}
if k < 2.5000000000000002e-6Initial program 40.4%
Applied egg-rr30.0%
Applied egg-rr88.7%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-/.f64N/A
lift-tan.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-invN/A
lift-/.f64N/A
clear-numN/A
*-commutativeN/A
lift-/.f64N/A
clear-numN/A
div-invN/A
lower-/.f64N/A
Applied egg-rr94.3%
Taylor expanded in k around 0
lower-/.f6483.0
Simplified83.0%
if 2.5000000000000002e-6 < k Initial program 39.2%
Applied egg-rr31.2%
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-sin.f64N/A
lift-*.f64N/A
associate-*r/N/A
associate-/r/N/A
lower-*.f64N/A
Applied egg-rr76.5%
Final simplification81.6%
(FPCore (t l k) :precision binary64 (/ l (* (tan k) (* (/ (sin k) l) (* 0.5 (* k (* k t)))))))
double code(double t, double l, double k) {
return l / (tan(k) * ((sin(k) / l) * (0.5 * (k * (k * t)))));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = l / (tan(k) * ((sin(k) / l) * (0.5d0 * (k * (k * t)))))
end function
public static double code(double t, double l, double k) {
return l / (Math.tan(k) * ((Math.sin(k) / l) * (0.5 * (k * (k * t)))));
}
def code(t, l, k): return l / (math.tan(k) * ((math.sin(k) / l) * (0.5 * (k * (k * t)))))
function code(t, l, k) return Float64(l / Float64(tan(k) * Float64(Float64(sin(k) / l) * Float64(0.5 * Float64(k * Float64(k * t)))))) end
function tmp = code(t, l, k) tmp = l / (tan(k) * ((sin(k) / l) * (0.5 * (k * (k * t))))); end
code[t_, l_, k_] := N[(l / N[(N[Tan[k], $MachinePrecision] * N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[(0.5 * N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell}{\tan k \cdot \left(\frac{\sin k}{\ell} \cdot \left(0.5 \cdot \left(k \cdot \left(k \cdot t\right)\right)\right)\right)}
\end{array}
Initial program 40.2%
Applied egg-rr30.2%
Applied egg-rr87.6%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-/.f64N/A
lift-/.f64N/A
lift-tan.f64N/A
associate-/r/N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
div-invN/A
lift-/.f64N/A
Applied egg-rr91.2%
Final simplification91.2%
(FPCore (t l k)
:precision binary64
(if (<= k 17000.0)
(/
(/ (fma l (* (* k k) -0.3333333333333333) l) k)
(* (/ (sin k) l) (* 0.5 (* k (* k t)))))
(/ (* 2.0 (* l (* l (cos k)))) (* k (* k (* t (* k k)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 17000.0) {
tmp = (fma(l, ((k * k) * -0.3333333333333333), l) / k) / ((sin(k) / l) * (0.5 * (k * (k * t))));
} else {
tmp = (2.0 * (l * (l * cos(k)))) / (k * (k * (t * (k * k))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (k <= 17000.0) tmp = Float64(Float64(fma(l, Float64(Float64(k * k) * -0.3333333333333333), l) / k) / Float64(Float64(sin(k) / l) * Float64(0.5 * Float64(k * Float64(k * t))))); else tmp = Float64(Float64(2.0 * Float64(l * Float64(l * cos(k)))) / Float64(k * Float64(k * Float64(t * Float64(k * k))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[k, 17000.0], N[(N[(N[(l * N[(N[(k * k), $MachinePrecision] * -0.3333333333333333), $MachinePrecision] + l), $MachinePrecision] / k), $MachinePrecision] / N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[(0.5 * N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(l * N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * N[(k * N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 17000:\\
\;\;\;\;\frac{\frac{\mathsf{fma}\left(\ell, \left(k \cdot k\right) \cdot -0.3333333333333333, \ell\right)}{k}}{\frac{\sin k}{\ell} \cdot \left(0.5 \cdot \left(k \cdot \left(k \cdot t\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \left(\ell \cdot \left(\ell \cdot \cos k\right)\right)}{k \cdot \left(k \cdot \left(t \cdot \left(k \cdot k\right)\right)\right)}\\
\end{array}
\end{array}
if k < 17000Initial program 40.7%
Applied egg-rr30.3%
Applied egg-rr88.7%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-/.f64N/A
lift-tan.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-invN/A
lift-/.f64N/A
clear-numN/A
*-commutativeN/A
lift-/.f64N/A
clear-numN/A
div-invN/A
lower-/.f64N/A
Applied egg-rr94.3%
Taylor expanded in k around 0
lower-/.f64N/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6474.3
Simplified74.3%
if 17000 < k Initial program 38.1%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6476.6
Simplified76.6%
Taylor expanded in k around 0
lower-*.f64N/A
unpow2N/A
lower-*.f6465.6
Simplified65.6%
Final simplification72.5%
(FPCore (t l k) :precision binary64 (if (<= (* l l) 1e+192) (/ (/ l k) (* (/ (sin k) l) (* 0.5 (* k (* k t))))) (/ (* 2.0 (* l (* l (cos k)))) (* k (* k (* t (* k k)))))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e+192) {
tmp = (l / k) / ((sin(k) / l) * (0.5 * (k * (k * t))));
} else {
tmp = (2.0 * (l * (l * cos(k)))) / (k * (k * (t * (k * k))));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if ((l * l) <= 1d+192) then
tmp = (l / k) / ((sin(k) / l) * (0.5d0 * (k * (k * t))))
else
tmp = (2.0d0 * (l * (l * cos(k)))) / (k * (k * (t * (k * k))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e+192) {
tmp = (l / k) / ((Math.sin(k) / l) * (0.5 * (k * (k * t))));
} else {
tmp = (2.0 * (l * (l * Math.cos(k)))) / (k * (k * (t * (k * k))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (l * l) <= 1e+192: tmp = (l / k) / ((math.sin(k) / l) * (0.5 * (k * (k * t)))) else: tmp = (2.0 * (l * (l * math.cos(k)))) / (k * (k * (t * (k * k)))) return tmp
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e+192) tmp = Float64(Float64(l / k) / Float64(Float64(sin(k) / l) * Float64(0.5 * Float64(k * Float64(k * t))))); else tmp = Float64(Float64(2.0 * Float64(l * Float64(l * cos(k)))) / Float64(k * Float64(k * Float64(t * Float64(k * k))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((l * l) <= 1e+192) tmp = (l / k) / ((sin(k) / l) * (0.5 * (k * (k * t)))); else tmp = (2.0 * (l * (l * cos(k)))) / (k * (k * (t * (k * k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e+192], N[(N[(l / k), $MachinePrecision] / N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[(0.5 * N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(l * N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * N[(k * N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{+192}:\\
\;\;\;\;\frac{\frac{\ell}{k}}{\frac{\sin k}{\ell} \cdot \left(0.5 \cdot \left(k \cdot \left(k \cdot t\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \left(\ell \cdot \left(\ell \cdot \cos k\right)\right)}{k \cdot \left(k \cdot \left(t \cdot \left(k \cdot k\right)\right)\right)}\\
\end{array}
\end{array}
if (*.f64 l l) < 1.00000000000000004e192Initial program 36.0%
Applied egg-rr28.6%
Applied egg-rr91.0%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-/.f64N/A
lift-tan.f64N/A
lift-/.f64N/A
lift-/.f64N/A
div-invN/A
lift-/.f64N/A
clear-numN/A
*-commutativeN/A
lift-/.f64N/A
clear-numN/A
div-invN/A
lower-/.f64N/A
Applied egg-rr96.1%
Taylor expanded in k around 0
lower-/.f6485.5
Simplified85.5%
if 1.00000000000000004e192 < (*.f64 l l) Initial program 48.4%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6475.4
Simplified75.4%
Taylor expanded in k around 0
lower-*.f64N/A
unpow2N/A
lower-*.f6470.7
Simplified70.7%
Final simplification80.5%
(FPCore (t l k) :precision binary64 (if (<= (* l l) 1e-190) (* (* l 2.0) (/ 1.0 (/ (* k (* k t)) (/ (/ l k) k)))) (/ (* 2.0 (* l (* l (cos k)))) (* k (* k (* t (* k k)))))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-190) {
tmp = (l * 2.0) * (1.0 / ((k * (k * t)) / ((l / k) / k)));
} else {
tmp = (2.0 * (l * (l * cos(k)))) / (k * (k * (t * (k * k))));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if ((l * l) <= 1d-190) then
tmp = (l * 2.0d0) * (1.0d0 / ((k * (k * t)) / ((l / k) / k)))
else
tmp = (2.0d0 * (l * (l * cos(k)))) / (k * (k * (t * (k * k))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-190) {
tmp = (l * 2.0) * (1.0 / ((k * (k * t)) / ((l / k) / k)));
} else {
tmp = (2.0 * (l * (l * Math.cos(k)))) / (k * (k * (t * (k * k))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (l * l) <= 1e-190: tmp = (l * 2.0) * (1.0 / ((k * (k * t)) / ((l / k) / k))) else: tmp = (2.0 * (l * (l * math.cos(k)))) / (k * (k * (t * (k * k)))) return tmp
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-190) tmp = Float64(Float64(l * 2.0) * Float64(1.0 / Float64(Float64(k * Float64(k * t)) / Float64(Float64(l / k) / k)))); else tmp = Float64(Float64(2.0 * Float64(l * Float64(l * cos(k)))) / Float64(k * Float64(k * Float64(t * Float64(k * k))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((l * l) <= 1e-190) tmp = (l * 2.0) * (1.0 / ((k * (k * t)) / ((l / k) / k))); else tmp = (2.0 * (l * (l * cos(k)))) / (k * (k * (t * (k * k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e-190], N[(N[(l * 2.0), $MachinePrecision] * N[(1.0 / N[(N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision] / N[(N[(l / k), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(l * N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * N[(k * N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-190}:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{1}{\frac{k \cdot \left(k \cdot t\right)}{\frac{\frac{\ell}{k}}{k}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \left(\ell \cdot \left(\ell \cdot \cos k\right)\right)}{k \cdot \left(k \cdot \left(t \cdot \left(k \cdot k\right)\right)\right)}\\
\end{array}
\end{array}
if (*.f64 l l) < 1e-190Initial program 23.7%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6453.1
Simplified53.1%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6478.2
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6482.1
Applied egg-rr82.1%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
associate-/r*N/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6489.7
Applied egg-rr89.7%
if 1e-190 < (*.f64 l l) Initial program 48.9%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6485.9
Simplified85.9%
Taylor expanded in k around 0
lower-*.f64N/A
unpow2N/A
lower-*.f6473.2
Simplified73.2%
Final simplification79.0%
(FPCore (t l k) :precision binary64 (if (<= k 175000000000.0) (* (* l 2.0) (/ 1.0 (/ (* k (* k t)) (/ (/ l k) k)))) (* (/ l k) (/ (* l -0.3333333333333333) (* k t)))))
double code(double t, double l, double k) {
double tmp;
if (k <= 175000000000.0) {
tmp = (l * 2.0) * (1.0 / ((k * (k * t)) / ((l / k) / k)));
} else {
tmp = (l / k) * ((l * -0.3333333333333333) / (k * t));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 175000000000.0d0) then
tmp = (l * 2.0d0) * (1.0d0 / ((k * (k * t)) / ((l / k) / k)))
else
tmp = (l / k) * ((l * (-0.3333333333333333d0)) / (k * t))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 175000000000.0) {
tmp = (l * 2.0) * (1.0 / ((k * (k * t)) / ((l / k) / k)));
} else {
tmp = (l / k) * ((l * -0.3333333333333333) / (k * t));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 175000000000.0: tmp = (l * 2.0) * (1.0 / ((k * (k * t)) / ((l / k) / k))) else: tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 175000000000.0) tmp = Float64(Float64(l * 2.0) * Float64(1.0 / Float64(Float64(k * Float64(k * t)) / Float64(Float64(l / k) / k)))); else tmp = Float64(Float64(l / k) * Float64(Float64(l * -0.3333333333333333) / Float64(k * t))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 175000000000.0) tmp = (l * 2.0) * (1.0 / ((k * (k * t)) / ((l / k) / k))); else tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 175000000000.0], N[(N[(l * 2.0), $MachinePrecision] * N[(1.0 / N[(N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision] / N[(N[(l / k), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(N[(l * -0.3333333333333333), $MachinePrecision] / N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 175000000000:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{1}{\frac{k \cdot \left(k \cdot t\right)}{\frac{\frac{\ell}{k}}{k}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \frac{\ell \cdot -0.3333333333333333}{k \cdot t}\\
\end{array}
\end{array}
if k < 1.75e11Initial program 40.7%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6464.4
Simplified64.4%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6474.9
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6477.2
Applied egg-rr77.2%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
associate-/r*N/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6480.8
Applied egg-rr80.8%
if 1.75e11 < k Initial program 37.9%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6474.3
Simplified74.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
lower-/.f64N/A
Simplified30.9%
Taylor expanded in k around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6463.7
Simplified63.7%
associate-*l*N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f6467.4
lift-*.f64N/A
*-commutativeN/A
lower-*.f6467.4
Applied egg-rr67.4%
Final simplification78.2%
(FPCore (t l k) :precision binary64 (if (<= k 175000000000.0) (* (/ (* l 2.0) (* t (* k k))) (/ l (* k k))) (* (/ l k) (/ (* l -0.3333333333333333) (* k t)))))
double code(double t, double l, double k) {
double tmp;
if (k <= 175000000000.0) {
tmp = ((l * 2.0) / (t * (k * k))) * (l / (k * k));
} else {
tmp = (l / k) * ((l * -0.3333333333333333) / (k * t));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 175000000000.0d0) then
tmp = ((l * 2.0d0) / (t * (k * k))) * (l / (k * k))
else
tmp = (l / k) * ((l * (-0.3333333333333333d0)) / (k * t))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 175000000000.0) {
tmp = ((l * 2.0) / (t * (k * k))) * (l / (k * k));
} else {
tmp = (l / k) * ((l * -0.3333333333333333) / (k * t));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 175000000000.0: tmp = ((l * 2.0) / (t * (k * k))) * (l / (k * k)) else: tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 175000000000.0) tmp = Float64(Float64(Float64(l * 2.0) / Float64(t * Float64(k * k))) * Float64(l / Float64(k * k))); else tmp = Float64(Float64(l / k) * Float64(Float64(l * -0.3333333333333333) / Float64(k * t))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 175000000000.0) tmp = ((l * 2.0) / (t * (k * k))) * (l / (k * k)); else tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 175000000000.0], N[(N[(N[(l * 2.0), $MachinePrecision] / N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(N[(l * -0.3333333333333333), $MachinePrecision] / N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 175000000000:\\
\;\;\;\;\frac{\ell \cdot 2}{t \cdot \left(k \cdot k\right)} \cdot \frac{\ell}{k \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \frac{\ell \cdot -0.3333333333333333}{k \cdot t}\\
\end{array}
\end{array}
if k < 1.75e11Initial program 40.7%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6464.4
Simplified64.4%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6480.5
Applied egg-rr80.5%
if 1.75e11 < k Initial program 37.9%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6474.3
Simplified74.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
lower-/.f64N/A
Simplified30.9%
Taylor expanded in k around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6463.7
Simplified63.7%
associate-*l*N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f6467.4
lift-*.f64N/A
*-commutativeN/A
lower-*.f6467.4
Applied egg-rr67.4%
Final simplification78.0%
(FPCore (t l k) :precision binary64 (if (<= k 175000000000.0) (* (* l 2.0) (/ (/ l k) (* k (* k (* k t))))) (* (/ l k) (/ (* l -0.3333333333333333) (* k t)))))
double code(double t, double l, double k) {
double tmp;
if (k <= 175000000000.0) {
tmp = (l * 2.0) * ((l / k) / (k * (k * (k * t))));
} else {
tmp = (l / k) * ((l * -0.3333333333333333) / (k * t));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 175000000000.0d0) then
tmp = (l * 2.0d0) * ((l / k) / (k * (k * (k * t))))
else
tmp = (l / k) * ((l * (-0.3333333333333333d0)) / (k * t))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 175000000000.0) {
tmp = (l * 2.0) * ((l / k) / (k * (k * (k * t))));
} else {
tmp = (l / k) * ((l * -0.3333333333333333) / (k * t));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 175000000000.0: tmp = (l * 2.0) * ((l / k) / (k * (k * (k * t)))) else: tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 175000000000.0) tmp = Float64(Float64(l * 2.0) * Float64(Float64(l / k) / Float64(k * Float64(k * Float64(k * t))))); else tmp = Float64(Float64(l / k) * Float64(Float64(l * -0.3333333333333333) / Float64(k * t))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 175000000000.0) tmp = (l * 2.0) * ((l / k) / (k * (k * (k * t)))); else tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 175000000000.0], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(l / k), $MachinePrecision] / N[(k * N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(N[(l * -0.3333333333333333), $MachinePrecision] / N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 175000000000:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\ell}{k}}{k \cdot \left(k \cdot \left(k \cdot t\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \frac{\ell \cdot -0.3333333333333333}{k \cdot t}\\
\end{array}
\end{array}
if k < 1.75e11Initial program 40.7%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6464.4
Simplified64.4%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6474.9
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6477.2
Applied egg-rr77.2%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6477.5
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
lower-*.f6478.7
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6479.2
Applied egg-rr79.2%
if 1.75e11 < k Initial program 37.9%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6474.3
Simplified74.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
lower-/.f64N/A
Simplified30.9%
Taylor expanded in k around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6463.7
Simplified63.7%
associate-*l*N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f6467.4
lift-*.f64N/A
*-commutativeN/A
lower-*.f6467.4
Applied egg-rr67.4%
Final simplification76.9%
(FPCore (t l k) :precision binary64 (if (<= k 175000000000.0) (* (* l 2.0) (/ l (* k (* (* k t) (* k k))))) (* (/ l k) (/ (* l -0.3333333333333333) (* k t)))))
double code(double t, double l, double k) {
double tmp;
if (k <= 175000000000.0) {
tmp = (l * 2.0) * (l / (k * ((k * t) * (k * k))));
} else {
tmp = (l / k) * ((l * -0.3333333333333333) / (k * t));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 175000000000.0d0) then
tmp = (l * 2.0d0) * (l / (k * ((k * t) * (k * k))))
else
tmp = (l / k) * ((l * (-0.3333333333333333d0)) / (k * t))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 175000000000.0) {
tmp = (l * 2.0) * (l / (k * ((k * t) * (k * k))));
} else {
tmp = (l / k) * ((l * -0.3333333333333333) / (k * t));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 175000000000.0: tmp = (l * 2.0) * (l / (k * ((k * t) * (k * k)))) else: tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 175000000000.0) tmp = Float64(Float64(l * 2.0) * Float64(l / Float64(k * Float64(Float64(k * t) * Float64(k * k))))); else tmp = Float64(Float64(l / k) * Float64(Float64(l * -0.3333333333333333) / Float64(k * t))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 175000000000.0) tmp = (l * 2.0) * (l / (k * ((k * t) * (k * k)))); else tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 175000000000.0], N[(N[(l * 2.0), $MachinePrecision] * N[(l / N[(k * N[(N[(k * t), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(N[(l * -0.3333333333333333), $MachinePrecision] / N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 175000000000:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\ell}{k \cdot \left(\left(k \cdot t\right) \cdot \left(k \cdot k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \frac{\ell \cdot -0.3333333333333333}{k \cdot t}\\
\end{array}
\end{array}
if k < 1.75e11Initial program 40.7%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6464.4
Simplified64.4%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6474.9
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6477.2
Applied egg-rr77.2%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6478.5
Applied egg-rr78.5%
if 1.75e11 < k Initial program 37.9%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6474.3
Simplified74.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
lower-/.f64N/A
Simplified30.9%
Taylor expanded in k around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6463.7
Simplified63.7%
associate-*l*N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f6467.4
lift-*.f64N/A
*-commutativeN/A
lower-*.f6467.4
Applied egg-rr67.4%
Final simplification76.4%
(FPCore (t l k) :precision binary64 (if (<= k 175000000000.0) (* (* l 2.0) (/ l (* k (* t (* k (* k k)))))) (* (/ l k) (/ (* l -0.3333333333333333) (* k t)))))
double code(double t, double l, double k) {
double tmp;
if (k <= 175000000000.0) {
tmp = (l * 2.0) * (l / (k * (t * (k * (k * k)))));
} else {
tmp = (l / k) * ((l * -0.3333333333333333) / (k * t));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 175000000000.0d0) then
tmp = (l * 2.0d0) * (l / (k * (t * (k * (k * k)))))
else
tmp = (l / k) * ((l * (-0.3333333333333333d0)) / (k * t))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 175000000000.0) {
tmp = (l * 2.0) * (l / (k * (t * (k * (k * k)))));
} else {
tmp = (l / k) * ((l * -0.3333333333333333) / (k * t));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 175000000000.0: tmp = (l * 2.0) * (l / (k * (t * (k * (k * k))))) else: tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 175000000000.0) tmp = Float64(Float64(l * 2.0) * Float64(l / Float64(k * Float64(t * Float64(k * Float64(k * k)))))); else tmp = Float64(Float64(l / k) * Float64(Float64(l * -0.3333333333333333) / Float64(k * t))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 175000000000.0) tmp = (l * 2.0) * (l / (k * (t * (k * (k * k))))); else tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 175000000000.0], N[(N[(l * 2.0), $MachinePrecision] * N[(l / N[(k * N[(t * N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(N[(l * -0.3333333333333333), $MachinePrecision] / N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 175000000000:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\ell}{k \cdot \left(t \cdot \left(k \cdot \left(k \cdot k\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \frac{\ell \cdot -0.3333333333333333}{k \cdot t}\\
\end{array}
\end{array}
if k < 1.75e11Initial program 40.7%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6464.4
Simplified64.4%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6474.9
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6477.2
Applied egg-rr77.2%
if 1.75e11 < k Initial program 37.9%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6474.3
Simplified74.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
lower-/.f64N/A
Simplified30.9%
Taylor expanded in k around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6463.7
Simplified63.7%
associate-*l*N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f6467.4
lift-*.f64N/A
*-commutativeN/A
lower-*.f6467.4
Applied egg-rr67.4%
Final simplification75.3%
(FPCore (t l k) :precision binary64 (* (/ l k) (/ (* l -0.3333333333333333) (* k t))))
double code(double t, double l, double k) {
return (l / k) * ((l * -0.3333333333333333) / (k * t));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = (l / k) * ((l * (-0.3333333333333333d0)) / (k * t))
end function
public static double code(double t, double l, double k) {
return (l / k) * ((l * -0.3333333333333333) / (k * t));
}
def code(t, l, k): return (l / k) * ((l * -0.3333333333333333) / (k * t))
function code(t, l, k) return Float64(Float64(l / k) * Float64(Float64(l * -0.3333333333333333) / Float64(k * t))) end
function tmp = code(t, l, k) tmp = (l / k) * ((l * -0.3333333333333333) / (k * t)); end
code[t_, l_, k_] := N[(N[(l / k), $MachinePrecision] * N[(N[(l * -0.3333333333333333), $MachinePrecision] / N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell}{k} \cdot \frac{\ell \cdot -0.3333333333333333}{k \cdot t}
\end{array}
Initial program 40.2%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6475.7
Simplified75.7%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
lower-/.f64N/A
Simplified46.6%
Taylor expanded in k around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6428.3
Simplified28.3%
associate-*l*N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f6429.6
lift-*.f64N/A
*-commutativeN/A
lower-*.f6429.6
Applied egg-rr29.6%
(FPCore (t l k) :precision binary64 (* -0.3333333333333333 (* l (/ l (* t (* k k))))))
double code(double t, double l, double k) {
return -0.3333333333333333 * (l * (l / (t * (k * k))));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = (-0.3333333333333333d0) * (l * (l / (t * (k * k))))
end function
public static double code(double t, double l, double k) {
return -0.3333333333333333 * (l * (l / (t * (k * k))));
}
def code(t, l, k): return -0.3333333333333333 * (l * (l / (t * (k * k))))
function code(t, l, k) return Float64(-0.3333333333333333 * Float64(l * Float64(l / Float64(t * Float64(k * k))))) end
function tmp = code(t, l, k) tmp = -0.3333333333333333 * (l * (l / (t * (k * k)))); end
code[t_, l_, k_] := N[(-0.3333333333333333 * N[(l * N[(l / N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.3333333333333333 \cdot \left(\ell \cdot \frac{\ell}{t \cdot \left(k \cdot k\right)}\right)
\end{array}
Initial program 40.2%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6475.7
Simplified75.7%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
lower-/.f64N/A
Simplified46.6%
Taylor expanded in k around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6428.3
Simplified28.3%
Taylor expanded in l around 0
lower-*.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6429.0
Simplified29.0%
Final simplification29.0%
herbie shell --seed 2024207
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10-)"
:precision binary64
(/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))