
(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 22 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
(let* ((t_1 (cbrt (sin k))) (t_2 (cbrt (tan k))))
(if (<= t -2.9e-29)
(/
2.0
(*
(pow (* (* (* t (pow (cbrt l) -2.0)) t_1) t_2) 3.0)
(+ 1.0 (+ 1.0 (pow (/ k t) 2.0)))))
(if (<= t 4.2e-48)
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))
(/
2.0
(*
(pow (* t_2 (* t_1 (/ t (pow (cbrt l) 2.0)))) 3.0)
(+ 1.0 (+ 1.0 (/ (/ k t) (/ t k))))))))))
double code(double t, double l, double k) {
double t_1 = cbrt(sin(k));
double t_2 = cbrt(tan(k));
double tmp;
if (t <= -2.9e-29) {
tmp = 2.0 / (pow((((t * pow(cbrt(l), -2.0)) * t_1) * t_2), 3.0) * (1.0 + (1.0 + pow((k / t), 2.0))));
} else if (t <= 4.2e-48) {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / tan(k);
} else {
tmp = 2.0 / (pow((t_2 * (t_1 * (t / pow(cbrt(l), 2.0)))), 3.0) * (1.0 + (1.0 + ((k / t) / (t / k)))));
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = Math.cbrt(Math.sin(k));
double t_2 = Math.cbrt(Math.tan(k));
double tmp;
if (t <= -2.9e-29) {
tmp = 2.0 / (Math.pow((((t * Math.pow(Math.cbrt(l), -2.0)) * t_1) * t_2), 3.0) * (1.0 + (1.0 + Math.pow((k / t), 2.0))));
} else if (t <= 4.2e-48) {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
} else {
tmp = 2.0 / (Math.pow((t_2 * (t_1 * (t / Math.pow(Math.cbrt(l), 2.0)))), 3.0) * (1.0 + (1.0 + ((k / t) / (t / k)))));
}
return tmp;
}
function code(t, l, k) t_1 = cbrt(sin(k)) t_2 = cbrt(tan(k)) tmp = 0.0 if (t <= -2.9e-29) tmp = Float64(2.0 / Float64((Float64(Float64(Float64(t * (cbrt(l) ^ -2.0)) * t_1) * t_2) ^ 3.0) * Float64(1.0 + Float64(1.0 + (Float64(k / t) ^ 2.0))))); elseif (t <= 4.2e-48) tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); else tmp = Float64(2.0 / Float64((Float64(t_2 * Float64(t_1 * Float64(t / (cbrt(l) ^ 2.0)))) ^ 3.0) * Float64(1.0 + Float64(1.0 + Float64(Float64(k / t) / Float64(t / k)))))); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]}, Block[{t$95$2 = N[Power[N[Tan[k], $MachinePrecision], 1/3], $MachinePrecision]}, If[LessEqual[t, -2.9e-29], N[(2.0 / N[(N[Power[N[(N[(N[(t * N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision] * t$95$2), $MachinePrecision], 3.0], $MachinePrecision] * N[(1.0 + N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4.2e-48], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Power[N[(t$95$2 * N[(t$95$1 * N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(1.0 + N[(1.0 + N[(N[(k / t), $MachinePrecision] / N[(t / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sqrt[3]{\sin k}\\
t_2 := \sqrt[3]{\tan k}\\
\mathbf{if}\;t \leq -2.9 \cdot 10^{-29}:\\
\;\;\;\;\frac{2}{{\left(\left(\left(t \cdot {\left(\sqrt[3]{\ell}\right)}^{-2}\right) \cdot t_1\right) \cdot t_2\right)}^{3} \cdot \left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right)}\\
\mathbf{elif}\;t \leq 4.2 \cdot 10^{-48}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(t_2 \cdot \left(t_1 \cdot \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)\right)}^{3} \cdot \left(1 + \left(1 + \frac{\frac{k}{t}}{\frac{t}{k}}\right)\right)}\\
\end{array}
\end{array}
if t < -2.90000000000000024e-29Initial program 58.1%
associate-/r*64.7%
add-cube-cbrt64.7%
*-un-lft-identity64.7%
times-frac64.7%
pow264.7%
cbrt-div64.7%
rem-cbrt-cube64.7%
cbrt-div64.7%
rem-cbrt-cube79.2%
Applied egg-rr79.2%
add-cube-cbrt79.0%
pow379.0%
cbrt-prod79.0%
frac-times73.3%
unpow273.3%
*-un-lft-identity73.3%
cbrt-div74.7%
add-cbrt-cube87.7%
Applied egg-rr87.7%
add-cube-cbrt87.6%
pow387.7%
cbrt-prod87.6%
unpow387.6%
add-cbrt-cube91.6%
associate-/l/91.5%
pow291.5%
Applied egg-rr91.5%
expm1-log1p-u61.1%
expm1-udef49.2%
div-inv49.2%
pow-flip49.2%
metadata-eval49.2%
Applied egg-rr49.2%
expm1-def61.1%
expm1-log1p91.6%
Simplified91.6%
if -2.90000000000000024e-29 < t < 4.19999999999999977e-48Initial program 36.3%
Simplified35.4%
associate-*r*41.8%
associate-*r*41.8%
times-frac43.5%
div-inv43.5%
pow-flip43.5%
metadata-eval43.5%
Applied egg-rr43.5%
times-frac43.6%
Simplified43.6%
Taylor expanded in t around 0 80.2%
associate-*r/80.2%
associate-*r*80.2%
Simplified80.2%
associate-*r/80.3%
times-frac82.5%
*-commutative82.5%
Applied egg-rr82.5%
if 4.19999999999999977e-48 < t Initial program 67.1%
associate-/r*72.7%
add-cube-cbrt72.6%
*-un-lft-identity72.6%
times-frac72.6%
pow272.6%
cbrt-div72.6%
rem-cbrt-cube72.6%
cbrt-div72.6%
rem-cbrt-cube85.0%
Applied egg-rr85.0%
add-cube-cbrt84.9%
pow384.8%
cbrt-prod84.8%
frac-times78.6%
unpow278.6%
*-un-lft-identity78.6%
cbrt-div78.6%
add-cbrt-cube89.5%
Applied egg-rr89.5%
add-cube-cbrt89.5%
pow389.5%
cbrt-prod89.2%
unpow389.3%
add-cbrt-cube95.2%
associate-/l/95.3%
pow295.3%
Applied egg-rr95.3%
unpow295.3%
clear-num95.3%
un-div-inv95.3%
Applied egg-rr95.3%
Final simplification88.6%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (pow (/ k t) 2.0)))
(if (<=
(* (+ 1.0 (+ 1.0 t_1)) (* (tan k) (* (sin k) (/ (pow t 3.0) (* l l)))))
INFINITY)
(/ (* (/ (* 2.0 (pow t -3.0)) (sin k)) (/ l (+ 2.0 t_1))) (/ (tan k) l))
(* (/ l (tan k)) (* 2.0 (/ (/ l (sin k)) (* t (pow k 2.0))))))))
double code(double t, double l, double k) {
double t_1 = pow((k / t), 2.0);
double tmp;
if (((1.0 + (1.0 + t_1)) * (tan(k) * (sin(k) * (pow(t, 3.0) / (l * l))))) <= ((double) INFINITY)) {
tmp = (((2.0 * pow(t, -3.0)) / sin(k)) * (l / (2.0 + t_1))) / (tan(k) / l);
} else {
tmp = (l / tan(k)) * (2.0 * ((l / sin(k)) / (t * pow(k, 2.0))));
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = Math.pow((k / t), 2.0);
double tmp;
if (((1.0 + (1.0 + t_1)) * (Math.tan(k) * (Math.sin(k) * (Math.pow(t, 3.0) / (l * l))))) <= Double.POSITIVE_INFINITY) {
tmp = (((2.0 * Math.pow(t, -3.0)) / Math.sin(k)) * (l / (2.0 + t_1))) / (Math.tan(k) / l);
} else {
tmp = (l / Math.tan(k)) * (2.0 * ((l / Math.sin(k)) / (t * Math.pow(k, 2.0))));
}
return tmp;
}
def code(t, l, k): t_1 = math.pow((k / t), 2.0) tmp = 0 if ((1.0 + (1.0 + t_1)) * (math.tan(k) * (math.sin(k) * (math.pow(t, 3.0) / (l * l))))) <= math.inf: tmp = (((2.0 * math.pow(t, -3.0)) / math.sin(k)) * (l / (2.0 + t_1))) / (math.tan(k) / l) else: tmp = (l / math.tan(k)) * (2.0 * ((l / math.sin(k)) / (t * math.pow(k, 2.0)))) return tmp
function code(t, l, k) t_1 = Float64(k / t) ^ 2.0 tmp = 0.0 if (Float64(Float64(1.0 + Float64(1.0 + t_1)) * Float64(tan(k) * Float64(sin(k) * Float64((t ^ 3.0) / Float64(l * l))))) <= Inf) tmp = Float64(Float64(Float64(Float64(2.0 * (t ^ -3.0)) / sin(k)) * Float64(l / Float64(2.0 + t_1))) / Float64(tan(k) / l)); else tmp = Float64(Float64(l / tan(k)) * Float64(2.0 * Float64(Float64(l / sin(k)) / Float64(t * (k ^ 2.0))))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = (k / t) ^ 2.0; tmp = 0.0; if (((1.0 + (1.0 + t_1)) * (tan(k) * (sin(k) * ((t ^ 3.0) / (l * l))))) <= Inf) tmp = (((2.0 * (t ^ -3.0)) / sin(k)) * (l / (2.0 + t_1))) / (tan(k) / l); else tmp = (l / tan(k)) * (2.0 * ((l / sin(k)) / (t * (k ^ 2.0)))); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[N[(N[(1.0 + N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(N[(2.0 * N[Power[t, -3.0], $MachinePrecision]), $MachinePrecision] / N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[(2.0 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Tan[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(2.0 * N[(N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision] / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {\left(\frac{k}{t}\right)}^{2}\\
\mathbf{if}\;\left(1 + \left(1 + t_1\right)\right) \cdot \left(\tan k \cdot \left(\sin k \cdot \frac{{t}^{3}}{\ell \cdot \ell}\right)\right) \leq \infty:\\
\;\;\;\;\frac{\frac{2 \cdot {t}^{-3}}{\sin k} \cdot \frac{\ell}{2 + t_1}}{\frac{\tan k}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{\tan k} \cdot \left(2 \cdot \frac{\frac{\ell}{\sin k}}{t \cdot {k}^{2}}\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t 3) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 1 (pow.f64 (/.f64 k t) 2)) 1)) < +inf.0Initial program 80.8%
Simplified74.4%
associate-*r*75.6%
associate-*r*75.6%
times-frac87.4%
div-inv87.4%
pow-flip87.9%
metadata-eval87.9%
Applied egg-rr87.9%
times-frac86.2%
Simplified86.2%
associate-*r/84.0%
Applied egg-rr84.0%
associate-/l*86.3%
associate-*l/89.2%
associate-*r/87.9%
associate-/r*87.9%
times-frac88.7%
Simplified88.7%
if +inf.0 < (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t 3) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 1 (pow.f64 (/.f64 k t) 2)) 1)) Initial program 0.0%
Simplified0.0%
associate-*r*14.7%
associate-*r*14.7%
times-frac14.7%
div-inv14.7%
pow-flip16.8%
metadata-eval16.8%
Applied egg-rr16.8%
times-frac16.9%
Simplified16.9%
Taylor expanded in t around 0 59.6%
associate-*r/59.6%
associate-*r*59.6%
Simplified59.6%
Taylor expanded in l around 0 59.6%
*-commutative59.6%
*-commutative59.6%
associate-*r*59.6%
associate-/r*61.3%
*-commutative61.3%
Simplified61.3%
Final simplification78.6%
(FPCore (t l k)
:precision binary64
(if (or (<= t -1.12e-30) (not (<= t 3.9e-48)))
(/
2.0
(*
(pow (* (cbrt (tan k)) (* (cbrt (sin k)) (/ t (pow (cbrt l) 2.0)))) 3.0)
(+ 1.0 (+ 1.0 (/ (/ k t) (/ t k))))))
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))))
double code(double t, double l, double k) {
double tmp;
if ((t <= -1.12e-30) || !(t <= 3.9e-48)) {
tmp = 2.0 / (pow((cbrt(tan(k)) * (cbrt(sin(k)) * (t / pow(cbrt(l), 2.0)))), 3.0) * (1.0 + (1.0 + ((k / t) / (t / k)))));
} else {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / tan(k);
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if ((t <= -1.12e-30) || !(t <= 3.9e-48)) {
tmp = 2.0 / (Math.pow((Math.cbrt(Math.tan(k)) * (Math.cbrt(Math.sin(k)) * (t / Math.pow(Math.cbrt(l), 2.0)))), 3.0) * (1.0 + (1.0 + ((k / t) / (t / k)))));
} else {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if ((t <= -1.12e-30) || !(t <= 3.9e-48)) tmp = Float64(2.0 / Float64((Float64(cbrt(tan(k)) * Float64(cbrt(sin(k)) * Float64(t / (cbrt(l) ^ 2.0)))) ^ 3.0) * Float64(1.0 + Float64(1.0 + Float64(Float64(k / t) / Float64(t / k)))))); else tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); end return tmp end
code[t_, l_, k_] := If[Or[LessEqual[t, -1.12e-30], N[Not[LessEqual[t, 3.9e-48]], $MachinePrecision]], N[(2.0 / N[(N[Power[N[(N[Power[N[Tan[k], $MachinePrecision], 1/3], $MachinePrecision] * N[(N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision] * N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(1.0 + N[(1.0 + N[(N[(k / t), $MachinePrecision] / N[(t / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.12 \cdot 10^{-30} \lor \neg \left(t \leq 3.9 \cdot 10^{-48}\right):\\
\;\;\;\;\frac{2}{{\left(\sqrt[3]{\tan k} \cdot \left(\sqrt[3]{\sin k} \cdot \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)\right)}^{3} \cdot \left(1 + \left(1 + \frac{\frac{k}{t}}{\frac{t}{k}}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\end{array}
\end{array}
if t < -1.12e-30 or 3.9e-48 < t Initial program 63.0%
associate-/r*69.1%
add-cube-cbrt69.0%
*-un-lft-identity69.0%
times-frac69.0%
pow269.0%
cbrt-div69.0%
rem-cbrt-cube69.0%
cbrt-div69.0%
rem-cbrt-cube82.3%
Applied egg-rr82.3%
add-cube-cbrt82.2%
pow382.2%
cbrt-prod82.1%
frac-times76.2%
unpow276.2%
*-un-lft-identity76.2%
cbrt-div76.8%
add-cbrt-cube88.6%
Applied egg-rr88.6%
add-cube-cbrt88.6%
pow388.7%
cbrt-prod88.5%
unpow388.5%
add-cbrt-cube93.5%
associate-/l/93.6%
pow293.6%
Applied egg-rr93.6%
unpow293.6%
clear-num93.6%
un-div-inv93.6%
Applied egg-rr93.6%
if -1.12e-30 < t < 3.9e-48Initial program 36.3%
Simplified35.4%
associate-*r*41.8%
associate-*r*41.8%
times-frac43.5%
div-inv43.5%
pow-flip43.5%
metadata-eval43.5%
Applied egg-rr43.5%
times-frac43.6%
Simplified43.6%
Taylor expanded in t around 0 80.2%
associate-*r/80.2%
associate-*r*80.2%
Simplified80.2%
associate-*r/80.3%
times-frac82.5%
*-commutative82.5%
Applied egg-rr82.5%
Final simplification88.6%
(FPCore (t l k)
:precision binary64
(if (or (<= t -5.5e+68) (not (<= t 2.4e-48)))
(/
2.0
(*
(+ 1.0 (+ 1.0 (pow (/ k t) 2.0)))
(* (tan k) (pow (* (cbrt (sin k)) (/ t (pow (cbrt l) 2.0))) 3.0))))
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))))
double code(double t, double l, double k) {
double tmp;
if ((t <= -5.5e+68) || !(t <= 2.4e-48)) {
tmp = 2.0 / ((1.0 + (1.0 + pow((k / t), 2.0))) * (tan(k) * pow((cbrt(sin(k)) * (t / pow(cbrt(l), 2.0))), 3.0)));
} else {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / tan(k);
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if ((t <= -5.5e+68) || !(t <= 2.4e-48)) {
tmp = 2.0 / ((1.0 + (1.0 + Math.pow((k / t), 2.0))) * (Math.tan(k) * Math.pow((Math.cbrt(Math.sin(k)) * (t / Math.pow(Math.cbrt(l), 2.0))), 3.0)));
} else {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if ((t <= -5.5e+68) || !(t <= 2.4e-48)) tmp = Float64(2.0 / Float64(Float64(1.0 + Float64(1.0 + (Float64(k / t) ^ 2.0))) * Float64(tan(k) * (Float64(cbrt(sin(k)) * Float64(t / (cbrt(l) ^ 2.0))) ^ 3.0)))); else tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); end return tmp end
code[t_, l_, k_] := If[Or[LessEqual[t, -5.5e+68], N[Not[LessEqual[t, 2.4e-48]], $MachinePrecision]], N[(2.0 / N[(N[(1.0 + N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[Power[N[(N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision] * N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -5.5 \cdot 10^{+68} \lor \neg \left(t \leq 2.4 \cdot 10^{-48}\right):\\
\;\;\;\;\frac{2}{\left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right) \cdot \left(\tan k \cdot {\left(\sqrt[3]{\sin k} \cdot \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\end{array}
\end{array}
if t < -5.5000000000000004e68 or 2.4e-48 < t Initial program 61.0%
associate-/r*68.4%
add-cube-cbrt68.3%
*-un-lft-identity68.3%
times-frac68.3%
pow268.3%
cbrt-div68.3%
rem-cbrt-cube68.3%
cbrt-div68.3%
rem-cbrt-cube84.4%
Applied egg-rr84.4%
add-cube-cbrt84.2%
pow384.2%
cbrt-prod84.2%
frac-times77.1%
unpow277.1%
*-un-lft-identity77.1%
cbrt-div77.0%
add-cbrt-cube91.2%
Applied egg-rr91.2%
expm1-log1p-u66.0%
expm1-udef49.9%
associate-/l/49.9%
pow249.9%
Applied egg-rr49.9%
expm1-def66.1%
expm1-log1p91.2%
Simplified91.2%
if -5.5000000000000004e68 < t < 2.4e-48Initial program 42.6%
Simplified41.1%
associate-*r*46.4%
associate-*r*46.4%
times-frac51.9%
div-inv51.9%
pow-flip51.9%
metadata-eval51.9%
Applied egg-rr51.9%
times-frac51.9%
Simplified51.9%
Taylor expanded in t around 0 80.1%
associate-*r/80.1%
associate-*r*80.2%
Simplified80.2%
associate-*r/80.2%
times-frac82.0%
*-commutative82.0%
Applied egg-rr82.0%
Final simplification86.2%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (cbrt (sin k))) (t_2 (+ 1.0 (+ 1.0 (pow (/ k t) 2.0)))))
(if (<= t -5.5e+68)
(/ 2.0 (* t_2 (* (tan k) (pow (* t_1 (/ (/ t (cbrt l)) (cbrt l))) 3.0))))
(if (<= t 3.8e-48)
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))
(/
2.0
(* t_2 (* (tan k) (pow (* t_1 (/ t (pow (cbrt l) 2.0))) 3.0))))))))
double code(double t, double l, double k) {
double t_1 = cbrt(sin(k));
double t_2 = 1.0 + (1.0 + pow((k / t), 2.0));
double tmp;
if (t <= -5.5e+68) {
tmp = 2.0 / (t_2 * (tan(k) * pow((t_1 * ((t / cbrt(l)) / cbrt(l))), 3.0)));
} else if (t <= 3.8e-48) {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / tan(k);
} else {
tmp = 2.0 / (t_2 * (tan(k) * pow((t_1 * (t / pow(cbrt(l), 2.0))), 3.0)));
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = Math.cbrt(Math.sin(k));
double t_2 = 1.0 + (1.0 + Math.pow((k / t), 2.0));
double tmp;
if (t <= -5.5e+68) {
tmp = 2.0 / (t_2 * (Math.tan(k) * Math.pow((t_1 * ((t / Math.cbrt(l)) / Math.cbrt(l))), 3.0)));
} else if (t <= 3.8e-48) {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
} else {
tmp = 2.0 / (t_2 * (Math.tan(k) * Math.pow((t_1 * (t / Math.pow(Math.cbrt(l), 2.0))), 3.0)));
}
return tmp;
}
function code(t, l, k) t_1 = cbrt(sin(k)) t_2 = Float64(1.0 + Float64(1.0 + (Float64(k / t) ^ 2.0))) tmp = 0.0 if (t <= -5.5e+68) tmp = Float64(2.0 / Float64(t_2 * Float64(tan(k) * (Float64(t_1 * Float64(Float64(t / cbrt(l)) / cbrt(l))) ^ 3.0)))); elseif (t <= 3.8e-48) tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); else tmp = Float64(2.0 / Float64(t_2 * Float64(tan(k) * (Float64(t_1 * Float64(t / (cbrt(l) ^ 2.0))) ^ 3.0)))); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]}, Block[{t$95$2 = N[(1.0 + N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -5.5e+68], N[(2.0 / N[(t$95$2 * N[(N[Tan[k], $MachinePrecision] * N[Power[N[(t$95$1 * N[(N[(t / N[Power[l, 1/3], $MachinePrecision]), $MachinePrecision] / N[Power[l, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 3.8e-48], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(t$95$2 * N[(N[Tan[k], $MachinePrecision] * N[Power[N[(t$95$1 * N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sqrt[3]{\sin k}\\
t_2 := 1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\\
\mathbf{if}\;t \leq -5.5 \cdot 10^{+68}:\\
\;\;\;\;\frac{2}{t_2 \cdot \left(\tan k \cdot {\left(t_1 \cdot \frac{\frac{t}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}\right)}^{3}\right)}\\
\mathbf{elif}\;t \leq 3.8 \cdot 10^{-48}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t_2 \cdot \left(\tan k \cdot {\left(t_1 \cdot \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)}\\
\end{array}
\end{array}
if t < -5.5000000000000004e68Initial program 49.7%
associate-/r*60.2%
add-cube-cbrt60.2%
*-un-lft-identity60.2%
times-frac60.2%
pow260.2%
cbrt-div60.2%
rem-cbrt-cube60.2%
cbrt-div60.2%
rem-cbrt-cube83.2%
Applied egg-rr83.2%
add-cube-cbrt83.1%
pow383.1%
cbrt-prod83.1%
frac-times74.1%
unpow274.1%
*-un-lft-identity74.1%
cbrt-div74.0%
add-cbrt-cube94.5%
Applied egg-rr94.5%
if -5.5000000000000004e68 < t < 3.80000000000000002e-48Initial program 42.6%
Simplified41.1%
associate-*r*46.4%
associate-*r*46.4%
times-frac51.9%
div-inv51.9%
pow-flip51.9%
metadata-eval51.9%
Applied egg-rr51.9%
times-frac51.9%
Simplified51.9%
Taylor expanded in t around 0 80.1%
associate-*r/80.1%
associate-*r*80.2%
Simplified80.2%
associate-*r/80.2%
times-frac82.0%
*-commutative82.0%
Applied egg-rr82.0%
if 3.80000000000000002e-48 < t Initial program 67.1%
associate-/r*72.7%
add-cube-cbrt72.6%
*-un-lft-identity72.6%
times-frac72.6%
pow272.6%
cbrt-div72.6%
rem-cbrt-cube72.6%
cbrt-div72.6%
rem-cbrt-cube85.0%
Applied egg-rr85.0%
add-cube-cbrt84.9%
pow384.8%
cbrt-prod84.8%
frac-times78.6%
unpow278.6%
*-un-lft-identity78.6%
cbrt-div78.6%
add-cbrt-cube89.5%
Applied egg-rr89.5%
expm1-log1p-u72.4%
expm1-udef52.4%
associate-/l/52.4%
pow252.4%
Applied egg-rr52.4%
expm1-def72.4%
expm1-log1p89.5%
Simplified89.5%
Final simplification86.2%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (pow (/ k t) 2.0))
(t_2
(/
2.0
(*
(+ 1.0 (+ 1.0 t_1))
(*
(tan k)
(pow (* (cbrt (sin k)) (* t (pow l -0.6666666666666666))) 3.0)))))
(t_3 (* (/ (* 2.0 (pow t -3.0)) (sin k)) (/ l (+ 2.0 t_1)))))
(if (<= t -6.7e+103)
t_2
(if (<= t -3e-30)
(/ t_3 (/ (tan k) l))
(if (<= t 2.1e-48)
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))
(if (<= t 2.6e+104) (* t_3 (/ l (tan k))) t_2))))))
double code(double t, double l, double k) {
double t_1 = pow((k / t), 2.0);
double t_2 = 2.0 / ((1.0 + (1.0 + t_1)) * (tan(k) * pow((cbrt(sin(k)) * (t * pow(l, -0.6666666666666666))), 3.0)));
double t_3 = ((2.0 * pow(t, -3.0)) / sin(k)) * (l / (2.0 + t_1));
double tmp;
if (t <= -6.7e+103) {
tmp = t_2;
} else if (t <= -3e-30) {
tmp = t_3 / (tan(k) / l);
} else if (t <= 2.1e-48) {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / tan(k);
} else if (t <= 2.6e+104) {
tmp = t_3 * (l / tan(k));
} else {
tmp = t_2;
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = Math.pow((k / t), 2.0);
double t_2 = 2.0 / ((1.0 + (1.0 + t_1)) * (Math.tan(k) * Math.pow((Math.cbrt(Math.sin(k)) * (t * Math.pow(l, -0.6666666666666666))), 3.0)));
double t_3 = ((2.0 * Math.pow(t, -3.0)) / Math.sin(k)) * (l / (2.0 + t_1));
double tmp;
if (t <= -6.7e+103) {
tmp = t_2;
} else if (t <= -3e-30) {
tmp = t_3 / (Math.tan(k) / l);
} else if (t <= 2.1e-48) {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
} else if (t <= 2.6e+104) {
tmp = t_3 * (l / Math.tan(k));
} else {
tmp = t_2;
}
return tmp;
}
function code(t, l, k) t_1 = Float64(k / t) ^ 2.0 t_2 = Float64(2.0 / Float64(Float64(1.0 + Float64(1.0 + t_1)) * Float64(tan(k) * (Float64(cbrt(sin(k)) * Float64(t * (l ^ -0.6666666666666666))) ^ 3.0)))) t_3 = Float64(Float64(Float64(2.0 * (t ^ -3.0)) / sin(k)) * Float64(l / Float64(2.0 + t_1))) tmp = 0.0 if (t <= -6.7e+103) tmp = t_2; elseif (t <= -3e-30) tmp = Float64(t_3 / Float64(tan(k) / l)); elseif (t <= 2.1e-48) tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); elseif (t <= 2.6e+104) tmp = Float64(t_3 * Float64(l / tan(k))); else tmp = t_2; end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[(2.0 / N[(N[(1.0 + N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[Power[N[(N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision] * N[(t * N[Power[l, -0.6666666666666666], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(2.0 * N[Power[t, -3.0], $MachinePrecision]), $MachinePrecision] / N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[(2.0 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -6.7e+103], t$95$2, If[LessEqual[t, -3e-30], N[(t$95$3 / N[(N[Tan[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.1e-48], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.6e+104], N[(t$95$3 * N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {\left(\frac{k}{t}\right)}^{2}\\
t_2 := \frac{2}{\left(1 + \left(1 + t_1\right)\right) \cdot \left(\tan k \cdot {\left(\sqrt[3]{\sin k} \cdot \left(t \cdot {\ell}^{-0.6666666666666666}\right)\right)}^{3}\right)}\\
t_3 := \frac{2 \cdot {t}^{-3}}{\sin k} \cdot \frac{\ell}{2 + t_1}\\
\mathbf{if}\;t \leq -6.7 \cdot 10^{+103}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq -3 \cdot 10^{-30}:\\
\;\;\;\;\frac{t_3}{\frac{\tan k}{\ell}}\\
\mathbf{elif}\;t \leq 2.1 \cdot 10^{-48}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\mathbf{elif}\;t \leq 2.6 \cdot 10^{+104}:\\
\;\;\;\;t_3 \cdot \frac{\ell}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;t_2\\
\end{array}
\end{array}
if t < -6.70000000000000033e103 or 2.6e104 < t Initial program 52.7%
add-cube-cbrt52.7%
pow352.7%
cbrt-prod52.7%
div-inv52.7%
cbrt-prod52.7%
rem-cbrt-cube67.2%
pow267.2%
pow-flip68.5%
metadata-eval68.5%
Applied egg-rr68.5%
pow1/367.8%
pow-pow37.8%
metadata-eval37.8%
Applied egg-rr37.8%
if -6.70000000000000033e103 < t < -2.9999999999999999e-30Initial program 67.0%
Simplified59.9%
associate-*r*63.2%
associate-*r*63.2%
times-frac88.7%
div-inv88.7%
pow-flip92.1%
metadata-eval92.1%
Applied egg-rr92.1%
times-frac89.6%
Simplified89.6%
associate-*r/86.5%
Applied egg-rr86.5%
associate-/l*89.7%
associate-*l/96.3%
associate-*r/91.9%
associate-/r*92.2%
times-frac96.3%
Simplified96.3%
if -2.9999999999999999e-30 < t < 2.09999999999999989e-48Initial program 36.3%
Simplified35.4%
associate-*r*41.8%
associate-*r*41.8%
times-frac43.5%
div-inv43.5%
pow-flip43.5%
metadata-eval43.5%
Applied egg-rr43.5%
times-frac43.6%
Simplified43.6%
Taylor expanded in t around 0 80.2%
associate-*r/80.2%
associate-*r*80.2%
Simplified80.2%
associate-*r/80.3%
times-frac82.5%
*-commutative82.5%
Applied egg-rr82.5%
if 2.09999999999999989e-48 < t < 2.6e104Initial program 81.7%
Simplified81.7%
associate-*r*84.4%
associate-*r*84.4%
times-frac91.9%
div-inv91.9%
pow-flip91.9%
metadata-eval91.9%
Applied egg-rr91.9%
times-frac86.5%
Simplified86.5%
expm1-log1p-u75.1%
expm1-udef62.3%
Applied egg-rr62.3%
expm1-def75.1%
expm1-log1p86.5%
associate-*l/94.3%
associate-*r/91.9%
associate-/r*91.9%
times-frac94.4%
Simplified94.4%
Final simplification72.3%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (+ 2.0 (pow (/ k t) 2.0)))
(t_2
(/
2.0
(* t_1 (* (sin k) (* (tan k) (pow (/ t (pow (cbrt l) 2.0)) 3.0)))))))
(if (<= t -6.8e+103)
t_2
(if (<= t -1e-30)
(/ (* (/ (* 2.0 (pow t -3.0)) (sin k)) (/ l t_1)) (/ (tan k) l))
(if (<= t 3.7e-48)
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))
t_2)))))
double code(double t, double l, double k) {
double t_1 = 2.0 + pow((k / t), 2.0);
double t_2 = 2.0 / (t_1 * (sin(k) * (tan(k) * pow((t / pow(cbrt(l), 2.0)), 3.0))));
double tmp;
if (t <= -6.8e+103) {
tmp = t_2;
} else if (t <= -1e-30) {
tmp = (((2.0 * pow(t, -3.0)) / sin(k)) * (l / t_1)) / (tan(k) / l);
} else if (t <= 3.7e-48) {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / tan(k);
} else {
tmp = t_2;
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = 2.0 + Math.pow((k / t), 2.0);
double t_2 = 2.0 / (t_1 * (Math.sin(k) * (Math.tan(k) * Math.pow((t / Math.pow(Math.cbrt(l), 2.0)), 3.0))));
double tmp;
if (t <= -6.8e+103) {
tmp = t_2;
} else if (t <= -1e-30) {
tmp = (((2.0 * Math.pow(t, -3.0)) / Math.sin(k)) * (l / t_1)) / (Math.tan(k) / l);
} else if (t <= 3.7e-48) {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
} else {
tmp = t_2;
}
return tmp;
}
function code(t, l, k) t_1 = Float64(2.0 + (Float64(k / t) ^ 2.0)) t_2 = Float64(2.0 / Float64(t_1 * Float64(sin(k) * Float64(tan(k) * (Float64(t / (cbrt(l) ^ 2.0)) ^ 3.0))))) tmp = 0.0 if (t <= -6.8e+103) tmp = t_2; elseif (t <= -1e-30) tmp = Float64(Float64(Float64(Float64(2.0 * (t ^ -3.0)) / sin(k)) * Float64(l / t_1)) / Float64(tan(k) / l)); elseif (t <= 3.7e-48) tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); else tmp = t_2; end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(2.0 / N[(t$95$1 * N[(N[Sin[k], $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[Power[N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -6.8e+103], t$95$2, If[LessEqual[t, -1e-30], N[(N[(N[(N[(2.0 * N[Power[t, -3.0], $MachinePrecision]), $MachinePrecision] / N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[(l / t$95$1), $MachinePrecision]), $MachinePrecision] / N[(N[Tan[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 3.7e-48], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 2 + {\left(\frac{k}{t}\right)}^{2}\\
t_2 := \frac{2}{t_1 \cdot \left(\sin k \cdot \left(\tan k \cdot {\left(\frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)\right)}\\
\mathbf{if}\;t \leq -6.8 \cdot 10^{+103}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq -1 \cdot 10^{-30}:\\
\;\;\;\;\frac{\frac{2 \cdot {t}^{-3}}{\sin k} \cdot \frac{\ell}{t_1}}{\frac{\tan k}{\ell}}\\
\mathbf{elif}\;t \leq 3.7 \cdot 10^{-48}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;t_2\\
\end{array}
\end{array}
if t < -6.7999999999999997e103 or 3.6999999999999998e-48 < t Initial program 61.9%
associate-/r*68.8%
add-cube-cbrt68.7%
*-un-lft-identity68.7%
times-frac68.7%
pow268.7%
cbrt-div68.7%
rem-cbrt-cube68.7%
cbrt-div68.7%
rem-cbrt-cube84.6%
Applied egg-rr84.6%
add-cube-cbrt84.5%
pow384.5%
cbrt-prod84.4%
frac-times77.0%
unpow277.0%
*-un-lft-identity77.0%
cbrt-div76.9%
add-cbrt-cube90.9%
Applied egg-rr90.9%
distribute-lft-in90.9%
Applied egg-rr87.7%
distribute-lft-out87.7%
associate-*l*87.7%
+-commutative87.7%
associate-+l+87.7%
metadata-eval87.7%
+-commutative87.7%
Simplified87.7%
if -6.7999999999999997e103 < t < -1e-30Initial program 67.0%
Simplified59.9%
associate-*r*63.2%
associate-*r*63.2%
times-frac88.7%
div-inv88.7%
pow-flip92.1%
metadata-eval92.1%
Applied egg-rr92.1%
times-frac89.6%
Simplified89.6%
associate-*r/86.5%
Applied egg-rr86.5%
associate-/l*89.7%
associate-*l/96.3%
associate-*r/91.9%
associate-/r*92.2%
times-frac96.3%
Simplified96.3%
if -1e-30 < t < 3.6999999999999998e-48Initial program 36.3%
Simplified35.4%
associate-*r*41.8%
associate-*r*41.8%
times-frac43.5%
div-inv43.5%
pow-flip43.5%
metadata-eval43.5%
Applied egg-rr43.5%
times-frac43.6%
Simplified43.6%
Taylor expanded in t around 0 80.2%
associate-*r/80.2%
associate-*r*80.2%
Simplified80.2%
associate-*r/80.3%
times-frac82.5%
*-commutative82.5%
Applied egg-rr82.5%
Final simplification86.3%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (pow (/ k t) 2.0)) (t_2 (+ 1.0 (+ 1.0 t_1))))
(if (<= t -6.4e+103)
(/ 2.0 (* t_2 (* (tan k) (* (sin k) (* (/ (pow t 2.0) l) (/ t l))))))
(if (<= t -7e-32)
(/ (* (/ (* 2.0 (pow t -3.0)) (sin k)) (/ l (+ 2.0 t_1))) (/ (tan k) l))
(if (<= t 1.9e-48)
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))
(/
2.0
(* t_2 (* (tan k) (* (sin k) (pow (/ (pow t 1.5) l) 2.0))))))))))
double code(double t, double l, double k) {
double t_1 = pow((k / t), 2.0);
double t_2 = 1.0 + (1.0 + t_1);
double tmp;
if (t <= -6.4e+103) {
tmp = 2.0 / (t_2 * (tan(k) * (sin(k) * ((pow(t, 2.0) / l) * (t / l)))));
} else if (t <= -7e-32) {
tmp = (((2.0 * pow(t, -3.0)) / sin(k)) * (l / (2.0 + t_1))) / (tan(k) / l);
} else if (t <= 1.9e-48) {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / tan(k);
} else {
tmp = 2.0 / (t_2 * (tan(k) * (sin(k) * pow((pow(t, 1.5) / l), 2.0))));
}
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) :: t_2
real(8) :: tmp
t_1 = (k / t) ** 2.0d0
t_2 = 1.0d0 + (1.0d0 + t_1)
if (t <= (-6.4d+103)) then
tmp = 2.0d0 / (t_2 * (tan(k) * (sin(k) * (((t ** 2.0d0) / l) * (t / l)))))
else if (t <= (-7d-32)) then
tmp = (((2.0d0 * (t ** (-3.0d0))) / sin(k)) * (l / (2.0d0 + t_1))) / (tan(k) / l)
else if (t <= 1.9d-48) then
tmp = (l * ((2.0d0 / (t * (k ** 2.0d0))) * (l / sin(k)))) / tan(k)
else
tmp = 2.0d0 / (t_2 * (tan(k) * (sin(k) * (((t ** 1.5d0) / l) ** 2.0d0))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = Math.pow((k / t), 2.0);
double t_2 = 1.0 + (1.0 + t_1);
double tmp;
if (t <= -6.4e+103) {
tmp = 2.0 / (t_2 * (Math.tan(k) * (Math.sin(k) * ((Math.pow(t, 2.0) / l) * (t / l)))));
} else if (t <= -7e-32) {
tmp = (((2.0 * Math.pow(t, -3.0)) / Math.sin(k)) * (l / (2.0 + t_1))) / (Math.tan(k) / l);
} else if (t <= 1.9e-48) {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
} else {
tmp = 2.0 / (t_2 * (Math.tan(k) * (Math.sin(k) * Math.pow((Math.pow(t, 1.5) / l), 2.0))));
}
return tmp;
}
def code(t, l, k): t_1 = math.pow((k / t), 2.0) t_2 = 1.0 + (1.0 + t_1) tmp = 0 if t <= -6.4e+103: tmp = 2.0 / (t_2 * (math.tan(k) * (math.sin(k) * ((math.pow(t, 2.0) / l) * (t / l))))) elif t <= -7e-32: tmp = (((2.0 * math.pow(t, -3.0)) / math.sin(k)) * (l / (2.0 + t_1))) / (math.tan(k) / l) elif t <= 1.9e-48: tmp = (l * ((2.0 / (t * math.pow(k, 2.0))) * (l / math.sin(k)))) / math.tan(k) else: tmp = 2.0 / (t_2 * (math.tan(k) * (math.sin(k) * math.pow((math.pow(t, 1.5) / l), 2.0)))) return tmp
function code(t, l, k) t_1 = Float64(k / t) ^ 2.0 t_2 = Float64(1.0 + Float64(1.0 + t_1)) tmp = 0.0 if (t <= -6.4e+103) tmp = Float64(2.0 / Float64(t_2 * Float64(tan(k) * Float64(sin(k) * Float64(Float64((t ^ 2.0) / l) * Float64(t / l)))))); elseif (t <= -7e-32) tmp = Float64(Float64(Float64(Float64(2.0 * (t ^ -3.0)) / sin(k)) * Float64(l / Float64(2.0 + t_1))) / Float64(tan(k) / l)); elseif (t <= 1.9e-48) tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); else tmp = Float64(2.0 / Float64(t_2 * Float64(tan(k) * Float64(sin(k) * (Float64((t ^ 1.5) / l) ^ 2.0))))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = (k / t) ^ 2.0; t_2 = 1.0 + (1.0 + t_1); tmp = 0.0; if (t <= -6.4e+103) tmp = 2.0 / (t_2 * (tan(k) * (sin(k) * (((t ^ 2.0) / l) * (t / l))))); elseif (t <= -7e-32) tmp = (((2.0 * (t ^ -3.0)) / sin(k)) * (l / (2.0 + t_1))) / (tan(k) / l); elseif (t <= 1.9e-48) tmp = (l * ((2.0 / (t * (k ^ 2.0))) * (l / sin(k)))) / tan(k); else tmp = 2.0 / (t_2 * (tan(k) * (sin(k) * (((t ^ 1.5) / l) ^ 2.0)))); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[(1.0 + N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -6.4e+103], N[(2.0 / N[(t$95$2 * N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, -7e-32], N[(N[(N[(N[(2.0 * N[Power[t, -3.0], $MachinePrecision]), $MachinePrecision] / N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[(2.0 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Tan[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.9e-48], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(t$95$2 * N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[Power[N[(N[Power[t, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {\left(\frac{k}{t}\right)}^{2}\\
t_2 := 1 + \left(1 + t_1\right)\\
\mathbf{if}\;t \leq -6.4 \cdot 10^{+103}:\\
\;\;\;\;\frac{2}{t_2 \cdot \left(\tan k \cdot \left(\sin k \cdot \left(\frac{{t}^{2}}{\ell} \cdot \frac{t}{\ell}\right)\right)\right)}\\
\mathbf{elif}\;t \leq -7 \cdot 10^{-32}:\\
\;\;\;\;\frac{\frac{2 \cdot {t}^{-3}}{\sin k} \cdot \frac{\ell}{2 + t_1}}{\frac{\tan k}{\ell}}\\
\mathbf{elif}\;t \leq 1.9 \cdot 10^{-48}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t_2 \cdot \left(\tan k \cdot \left(\sin k \cdot {\left(\frac{{t}^{1.5}}{\ell}\right)}^{2}\right)\right)}\\
\end{array}
\end{array}
if t < -6.39999999999999985e103Initial program 50.9%
unpow350.9%
times-frac73.5%
pow273.5%
Applied egg-rr73.5%
if -6.39999999999999985e103 < t < -6.9999999999999997e-32Initial program 67.0%
Simplified59.9%
associate-*r*63.2%
associate-*r*63.2%
times-frac88.7%
div-inv88.7%
pow-flip92.1%
metadata-eval92.1%
Applied egg-rr92.1%
times-frac89.6%
Simplified89.6%
associate-*r/86.5%
Applied egg-rr86.5%
associate-/l*89.7%
associate-*l/96.3%
associate-*r/91.9%
associate-/r*92.2%
times-frac96.3%
Simplified96.3%
if -6.9999999999999997e-32 < t < 1.90000000000000001e-48Initial program 36.3%
Simplified35.4%
associate-*r*41.8%
associate-*r*41.8%
times-frac43.5%
div-inv43.5%
pow-flip43.5%
metadata-eval43.5%
Applied egg-rr43.5%
times-frac43.6%
Simplified43.6%
Taylor expanded in t around 0 80.2%
associate-*r/80.2%
associate-*r*80.2%
Simplified80.2%
associate-*r/80.3%
times-frac82.5%
*-commutative82.5%
Applied egg-rr82.5%
if 1.90000000000000001e-48 < t Initial program 67.1%
add-sqr-sqrt67.0%
pow267.0%
sqrt-div67.1%
sqrt-pow171.0%
metadata-eval71.0%
sqrt-prod36.9%
add-sqr-sqrt82.6%
Applied egg-rr82.6%
Final simplification82.8%
(FPCore (t l k)
:precision binary64
(if (or (<= t -9e-30) (not (<= t 3.2e-48)))
(*
(* (/ (* 2.0 (pow t -3.0)) (sin k)) (/ l (+ 2.0 (pow (/ k t) 2.0))))
(/ l (tan k)))
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))))
double code(double t, double l, double k) {
double tmp;
if ((t <= -9e-30) || !(t <= 3.2e-48)) {
tmp = (((2.0 * pow(t, -3.0)) / sin(k)) * (l / (2.0 + pow((k / t), 2.0)))) * (l / tan(k));
} else {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / 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) :: tmp
if ((t <= (-9d-30)) .or. (.not. (t <= 3.2d-48))) then
tmp = (((2.0d0 * (t ** (-3.0d0))) / sin(k)) * (l / (2.0d0 + ((k / t) ** 2.0d0)))) * (l / tan(k))
else
tmp = (l * ((2.0d0 / (t * (k ** 2.0d0))) * (l / sin(k)))) / tan(k)
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((t <= -9e-30) || !(t <= 3.2e-48)) {
tmp = (((2.0 * Math.pow(t, -3.0)) / Math.sin(k)) * (l / (2.0 + Math.pow((k / t), 2.0)))) * (l / Math.tan(k));
} else {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
}
return tmp;
}
def code(t, l, k): tmp = 0 if (t <= -9e-30) or not (t <= 3.2e-48): tmp = (((2.0 * math.pow(t, -3.0)) / math.sin(k)) * (l / (2.0 + math.pow((k / t), 2.0)))) * (l / math.tan(k)) else: tmp = (l * ((2.0 / (t * math.pow(k, 2.0))) * (l / math.sin(k)))) / math.tan(k) return tmp
function code(t, l, k) tmp = 0.0 if ((t <= -9e-30) || !(t <= 3.2e-48)) tmp = Float64(Float64(Float64(Float64(2.0 * (t ^ -3.0)) / sin(k)) * Float64(l / Float64(2.0 + (Float64(k / t) ^ 2.0)))) * Float64(l / tan(k))); else tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((t <= -9e-30) || ~((t <= 3.2e-48))) tmp = (((2.0 * (t ^ -3.0)) / sin(k)) * (l / (2.0 + ((k / t) ^ 2.0)))) * (l / tan(k)); else tmp = (l * ((2.0 / (t * (k ^ 2.0))) * (l / sin(k)))) / tan(k); end tmp_2 = tmp; end
code[t_, l_, k_] := If[Or[LessEqual[t, -9e-30], N[Not[LessEqual[t, 3.2e-48]], $MachinePrecision]], N[(N[(N[(N[(2.0 * N[Power[t, -3.0], $MachinePrecision]), $MachinePrecision] / N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -9 \cdot 10^{-30} \lor \neg \left(t \leq 3.2 \cdot 10^{-48}\right):\\
\;\;\;\;\left(\frac{2 \cdot {t}^{-3}}{\sin k} \cdot \frac{\ell}{2 + {\left(\frac{k}{t}\right)}^{2}}\right) \cdot \frac{\ell}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\end{array}
\end{array}
if t < -8.99999999999999935e-30 or 3.1999999999999998e-48 < t Initial program 63.0%
Simplified56.4%
associate-*r*62.4%
associate-*r*62.4%
times-frac74.5%
div-inv74.5%
pow-flip76.5%
metadata-eval76.5%
Applied egg-rr76.5%
times-frac74.6%
Simplified74.6%
expm1-log1p-u65.3%
expm1-udef57.3%
Applied egg-rr57.3%
expm1-def65.3%
expm1-log1p74.6%
associate-*l/77.9%
associate-*r/76.5%
associate-/r*76.5%
times-frac78.0%
Simplified78.0%
if -8.99999999999999935e-30 < t < 3.1999999999999998e-48Initial program 36.3%
Simplified35.4%
associate-*r*41.8%
associate-*r*41.8%
times-frac43.5%
div-inv43.5%
pow-flip43.5%
metadata-eval43.5%
Applied egg-rr43.5%
times-frac43.6%
Simplified43.6%
Taylor expanded in t around 0 80.2%
associate-*r/80.2%
associate-*r*80.2%
Simplified80.2%
associate-*r/80.3%
times-frac82.5%
*-commutative82.5%
Applied egg-rr82.5%
Final simplification80.0%
(FPCore (t l k)
:precision binary64
(if (or (<= t -5.5e+68) (not (<= t 2e-49)))
(/
(/ 2.0 (* (tan k) (* (sin k) (/ (* t (/ (pow t 2.0) l)) l))))
(+ 2.0 (/ (/ k t) (/ t k))))
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))))
double code(double t, double l, double k) {
double tmp;
if ((t <= -5.5e+68) || !(t <= 2e-49)) {
tmp = (2.0 / (tan(k) * (sin(k) * ((t * (pow(t, 2.0) / l)) / l)))) / (2.0 + ((k / t) / (t / k)));
} else {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / 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) :: tmp
if ((t <= (-5.5d+68)) .or. (.not. (t <= 2d-49))) then
tmp = (2.0d0 / (tan(k) * (sin(k) * ((t * ((t ** 2.0d0) / l)) / l)))) / (2.0d0 + ((k / t) / (t / k)))
else
tmp = (l * ((2.0d0 / (t * (k ** 2.0d0))) * (l / sin(k)))) / tan(k)
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((t <= -5.5e+68) || !(t <= 2e-49)) {
tmp = (2.0 / (Math.tan(k) * (Math.sin(k) * ((t * (Math.pow(t, 2.0) / l)) / l)))) / (2.0 + ((k / t) / (t / k)));
} else {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
}
return tmp;
}
def code(t, l, k): tmp = 0 if (t <= -5.5e+68) or not (t <= 2e-49): tmp = (2.0 / (math.tan(k) * (math.sin(k) * ((t * (math.pow(t, 2.0) / l)) / l)))) / (2.0 + ((k / t) / (t / k))) else: tmp = (l * ((2.0 / (t * math.pow(k, 2.0))) * (l / math.sin(k)))) / math.tan(k) return tmp
function code(t, l, k) tmp = 0.0 if ((t <= -5.5e+68) || !(t <= 2e-49)) tmp = Float64(Float64(2.0 / Float64(tan(k) * Float64(sin(k) * Float64(Float64(t * Float64((t ^ 2.0) / l)) / l)))) / Float64(2.0 + Float64(Float64(k / t) / Float64(t / k)))); else tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((t <= -5.5e+68) || ~((t <= 2e-49))) tmp = (2.0 / (tan(k) * (sin(k) * ((t * ((t ^ 2.0) / l)) / l)))) / (2.0 + ((k / t) / (t / k))); else tmp = (l * ((2.0 / (t * (k ^ 2.0))) * (l / sin(k)))) / tan(k); end tmp_2 = tmp; end
code[t_, l_, k_] := If[Or[LessEqual[t, -5.5e+68], N[Not[LessEqual[t, 2e-49]], $MachinePrecision]], N[(N[(2.0 / N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(t * N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[(N[(k / t), $MachinePrecision] / N[(t / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -5.5 \cdot 10^{+68} \lor \neg \left(t \leq 2 \cdot 10^{-49}\right):\\
\;\;\;\;\frac{\frac{2}{\tan k \cdot \left(\sin k \cdot \frac{t \cdot \frac{{t}^{2}}{\ell}}{\ell}\right)}}{2 + \frac{\frac{k}{t}}{\frac{t}{k}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\end{array}
\end{array}
if t < -5.5000000000000004e68 or 1.99999999999999987e-49 < t Initial program 61.4%
Simplified68.7%
unpow296.4%
clear-num96.4%
un-div-inv96.4%
Applied egg-rr68.7%
cube-mult68.7%
*-un-lft-identity68.7%
times-frac76.0%
pow276.0%
Applied egg-rr76.0%
if -5.5000000000000004e68 < t < 1.99999999999999987e-49Initial program 42.2%
Simplified40.7%
associate-*r*46.0%
associate-*r*46.0%
times-frac51.6%
div-inv51.6%
pow-flip51.6%
metadata-eval51.6%
Applied egg-rr51.6%
times-frac51.6%
Simplified51.6%
Taylor expanded in t around 0 80.0%
associate-*r/80.0%
associate-*r*80.0%
Simplified80.0%
associate-*r/80.1%
times-frac81.8%
*-commutative81.8%
Applied egg-rr81.8%
Final simplification79.1%
(FPCore (t l k)
:precision binary64
(if (<= t -1.95e+77)
(* (/ l (tan k)) (/ l (pow (* t (cbrt k)) 3.0)))
(if (<= t 2.55e-48)
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))
(/
(/ 2.0 (* (tan k) (* (sin k) (/ (/ (pow t 3.0) l) l))))
(+ 2.0 (* (/ k t) (/ k t)))))))
double code(double t, double l, double k) {
double tmp;
if (t <= -1.95e+77) {
tmp = (l / tan(k)) * (l / pow((t * cbrt(k)), 3.0));
} else if (t <= 2.55e-48) {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / tan(k);
} else {
tmp = (2.0 / (tan(k) * (sin(k) * ((pow(t, 3.0) / l) / l)))) / (2.0 + ((k / t) * (k / t)));
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if (t <= -1.95e+77) {
tmp = (l / Math.tan(k)) * (l / Math.pow((t * Math.cbrt(k)), 3.0));
} else if (t <= 2.55e-48) {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
} else {
tmp = (2.0 / (Math.tan(k) * (Math.sin(k) * ((Math.pow(t, 3.0) / l) / l)))) / (2.0 + ((k / t) * (k / t)));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (t <= -1.95e+77) tmp = Float64(Float64(l / tan(k)) * Float64(l / (Float64(t * cbrt(k)) ^ 3.0))); elseif (t <= 2.55e-48) tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); else tmp = Float64(Float64(2.0 / Float64(tan(k) * Float64(sin(k) * Float64(Float64((t ^ 3.0) / l) / l)))) / Float64(2.0 + Float64(Float64(k / t) * Float64(k / t)))); end return tmp end
code[t_, l_, k_] := If[LessEqual[t, -1.95e+77], N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[Power[N[(t * N[Power[k, 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.55e-48], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[(N[(k / t), $MachinePrecision] * N[(k / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.95 \cdot 10^{+77}:\\
\;\;\;\;\frac{\ell}{\tan k} \cdot \frac{\ell}{{\left(t \cdot \sqrt[3]{k}\right)}^{3}}\\
\mathbf{elif}\;t \leq 2.55 \cdot 10^{-48}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\tan k \cdot \left(\sin k \cdot \frac{\frac{{t}^{3}}{\ell}}{\ell}\right)}}{2 + \frac{k}{t} \cdot \frac{k}{t}}\\
\end{array}
\end{array}
if t < -1.9499999999999999e77Initial program 50.9%
Simplified35.5%
associate-*r*43.8%
associate-*r*43.8%
times-frac60.7%
div-inv60.7%
pow-flip65.2%
metadata-eval65.2%
Applied egg-rr65.2%
times-frac65.8%
Simplified65.8%
Taylor expanded in k around 0 60.4%
add-cube-cbrt60.4%
pow360.4%
*-commutative60.4%
cbrt-prod60.4%
unpow360.4%
add-cbrt-cube70.2%
Applied egg-rr70.2%
if -1.9499999999999999e77 < t < 2.55000000000000006e-48Initial program 42.3%
Simplified40.9%
associate-*r*46.8%
associate-*r*46.8%
times-frac52.3%
div-inv52.3%
pow-flip52.2%
metadata-eval52.2%
Applied egg-rr52.2%
times-frac51.6%
Simplified51.6%
Taylor expanded in t around 0 79.6%
associate-*r/79.6%
associate-*r*79.6%
Simplified79.6%
associate-*r/79.6%
times-frac81.4%
*-commutative81.4%
Applied egg-rr81.4%
if 2.55000000000000006e-48 < t Initial program 67.1%
Simplified72.8%
unpow272.8%
Applied egg-rr72.8%
Final simplification77.1%
(FPCore (t l k)
:precision binary64
(if (<= t -1.95e+77)
(* (/ l (tan k)) (/ l (pow (* t (cbrt k)) 3.0)))
(if (<= t 1.4e-48)
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))
(/
(/ 2.0 (* (tan k) (* (sin k) (/ (/ (pow t 3.0) l) l))))
(+ 2.0 (/ (/ k t) (/ t k)))))))
double code(double t, double l, double k) {
double tmp;
if (t <= -1.95e+77) {
tmp = (l / tan(k)) * (l / pow((t * cbrt(k)), 3.0));
} else if (t <= 1.4e-48) {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / tan(k);
} else {
tmp = (2.0 / (tan(k) * (sin(k) * ((pow(t, 3.0) / l) / l)))) / (2.0 + ((k / t) / (t / k)));
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if (t <= -1.95e+77) {
tmp = (l / Math.tan(k)) * (l / Math.pow((t * Math.cbrt(k)), 3.0));
} else if (t <= 1.4e-48) {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
} else {
tmp = (2.0 / (Math.tan(k) * (Math.sin(k) * ((Math.pow(t, 3.0) / l) / l)))) / (2.0 + ((k / t) / (t / k)));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (t <= -1.95e+77) tmp = Float64(Float64(l / tan(k)) * Float64(l / (Float64(t * cbrt(k)) ^ 3.0))); elseif (t <= 1.4e-48) tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); else tmp = Float64(Float64(2.0 / Float64(tan(k) * Float64(sin(k) * Float64(Float64((t ^ 3.0) / l) / l)))) / Float64(2.0 + Float64(Float64(k / t) / Float64(t / k)))); end return tmp end
code[t_, l_, k_] := If[LessEqual[t, -1.95e+77], N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[Power[N[(t * N[Power[k, 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.4e-48], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[(N[(k / t), $MachinePrecision] / N[(t / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.95 \cdot 10^{+77}:\\
\;\;\;\;\frac{\ell}{\tan k} \cdot \frac{\ell}{{\left(t \cdot \sqrt[3]{k}\right)}^{3}}\\
\mathbf{elif}\;t \leq 1.4 \cdot 10^{-48}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\tan k \cdot \left(\sin k \cdot \frac{\frac{{t}^{3}}{\ell}}{\ell}\right)}}{2 + \frac{\frac{k}{t}}{\frac{t}{k}}}\\
\end{array}
\end{array}
if t < -1.9499999999999999e77Initial program 50.9%
Simplified35.5%
associate-*r*43.8%
associate-*r*43.8%
times-frac60.7%
div-inv60.7%
pow-flip65.2%
metadata-eval65.2%
Applied egg-rr65.2%
times-frac65.8%
Simplified65.8%
Taylor expanded in k around 0 60.4%
add-cube-cbrt60.4%
pow360.4%
*-commutative60.4%
cbrt-prod60.4%
unpow360.4%
add-cbrt-cube70.2%
Applied egg-rr70.2%
if -1.9499999999999999e77 < t < 1.40000000000000002e-48Initial program 41.9%
Simplified40.4%
associate-*r*46.4%
associate-*r*46.4%
times-frac51.9%
div-inv51.9%
pow-flip51.9%
metadata-eval51.9%
Applied egg-rr51.9%
times-frac51.3%
Simplified51.3%
Taylor expanded in t around 0 79.4%
associate-*r/79.4%
associate-*r*79.5%
Simplified79.5%
associate-*r/79.5%
times-frac81.3%
*-commutative81.3%
Applied egg-rr81.3%
if 1.40000000000000002e-48 < t Initial program 67.5%
Simplified73.1%
unpow295.3%
clear-num95.4%
un-div-inv95.4%
Applied egg-rr73.1%
Final simplification77.1%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ l (tan k))) (t_2 (/ l (pow (* t (cbrt k)) 3.0))))
(if (<= t -1.95e+77)
(* t_1 t_2)
(if (<= t 4.1e-48)
(* t_1 (* 2.0 (/ l (* (pow k 2.0) (* t (sin k))))))
(* t_2 (/ l k))))))
double code(double t, double l, double k) {
double t_1 = l / tan(k);
double t_2 = l / pow((t * cbrt(k)), 3.0);
double tmp;
if (t <= -1.95e+77) {
tmp = t_1 * t_2;
} else if (t <= 4.1e-48) {
tmp = t_1 * (2.0 * (l / (pow(k, 2.0) * (t * sin(k)))));
} else {
tmp = t_2 * (l / k);
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = l / Math.tan(k);
double t_2 = l / Math.pow((t * Math.cbrt(k)), 3.0);
double tmp;
if (t <= -1.95e+77) {
tmp = t_1 * t_2;
} else if (t <= 4.1e-48) {
tmp = t_1 * (2.0 * (l / (Math.pow(k, 2.0) * (t * Math.sin(k)))));
} else {
tmp = t_2 * (l / k);
}
return tmp;
}
function code(t, l, k) t_1 = Float64(l / tan(k)) t_2 = Float64(l / (Float64(t * cbrt(k)) ^ 3.0)) tmp = 0.0 if (t <= -1.95e+77) tmp = Float64(t_1 * t_2); elseif (t <= 4.1e-48) tmp = Float64(t_1 * Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t * sin(k)))))); else tmp = Float64(t_2 * Float64(l / k)); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(l / N[Power[N[(t * N[Power[k, 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.95e+77], N[(t$95$1 * t$95$2), $MachinePrecision], If[LessEqual[t, 4.1e-48], N[(t$95$1 * N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$2 * N[(l / k), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\ell}{\tan k}\\
t_2 := \frac{\ell}{{\left(t \cdot \sqrt[3]{k}\right)}^{3}}\\
\mathbf{if}\;t \leq -1.95 \cdot 10^{+77}:\\
\;\;\;\;t_1 \cdot t_2\\
\mathbf{elif}\;t \leq 4.1 \cdot 10^{-48}:\\
\;\;\;\;t_1 \cdot \left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t \cdot \sin k\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;t_2 \cdot \frac{\ell}{k}\\
\end{array}
\end{array}
if t < -1.9499999999999999e77Initial program 50.9%
Simplified35.5%
associate-*r*43.8%
associate-*r*43.8%
times-frac60.7%
div-inv60.7%
pow-flip65.2%
metadata-eval65.2%
Applied egg-rr65.2%
times-frac65.8%
Simplified65.8%
Taylor expanded in k around 0 60.4%
add-cube-cbrt60.4%
pow360.4%
*-commutative60.4%
cbrt-prod60.4%
unpow360.4%
add-cbrt-cube70.2%
Applied egg-rr70.2%
if -1.9499999999999999e77 < t < 4.10000000000000014e-48Initial program 42.3%
Simplified40.9%
associate-*r*46.8%
associate-*r*46.8%
times-frac52.3%
div-inv52.3%
pow-flip52.2%
metadata-eval52.2%
Applied egg-rr52.2%
times-frac51.6%
Simplified51.6%
Taylor expanded in t around 0 79.6%
if 4.10000000000000014e-48 < t Initial program 67.1%
Simplified64.3%
associate-*r*69.9%
associate-*r*69.9%
times-frac76.0%
div-inv76.0%
pow-flip77.3%
metadata-eval77.3%
Applied egg-rr77.3%
times-frac74.8%
Simplified74.8%
Taylor expanded in k around 0 65.9%
Taylor expanded in k around 0 69.9%
add-cube-cbrt69.8%
pow369.8%
*-commutative69.8%
cbrt-prod69.8%
unpow369.8%
add-cbrt-cube72.4%
Applied egg-rr72.4%
Final simplification75.9%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ l (tan k))) (t_2 (/ l (pow (* t (cbrt k)) 3.0))))
(if (<= t -1.95e+77)
(* t_1 t_2)
(if (<= t 4.4e-48)
(* t_1 (* 2.0 (/ (/ l (sin k)) (* t (pow k 2.0)))))
(* t_2 (/ l k))))))
double code(double t, double l, double k) {
double t_1 = l / tan(k);
double t_2 = l / pow((t * cbrt(k)), 3.0);
double tmp;
if (t <= -1.95e+77) {
tmp = t_1 * t_2;
} else if (t <= 4.4e-48) {
tmp = t_1 * (2.0 * ((l / sin(k)) / (t * pow(k, 2.0))));
} else {
tmp = t_2 * (l / k);
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = l / Math.tan(k);
double t_2 = l / Math.pow((t * Math.cbrt(k)), 3.0);
double tmp;
if (t <= -1.95e+77) {
tmp = t_1 * t_2;
} else if (t <= 4.4e-48) {
tmp = t_1 * (2.0 * ((l / Math.sin(k)) / (t * Math.pow(k, 2.0))));
} else {
tmp = t_2 * (l / k);
}
return tmp;
}
function code(t, l, k) t_1 = Float64(l / tan(k)) t_2 = Float64(l / (Float64(t * cbrt(k)) ^ 3.0)) tmp = 0.0 if (t <= -1.95e+77) tmp = Float64(t_1 * t_2); elseif (t <= 4.4e-48) tmp = Float64(t_1 * Float64(2.0 * Float64(Float64(l / sin(k)) / Float64(t * (k ^ 2.0))))); else tmp = Float64(t_2 * Float64(l / k)); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(l / N[Power[N[(t * N[Power[k, 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.95e+77], N[(t$95$1 * t$95$2), $MachinePrecision], If[LessEqual[t, 4.4e-48], N[(t$95$1 * N[(2.0 * N[(N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision] / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$2 * N[(l / k), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\ell}{\tan k}\\
t_2 := \frac{\ell}{{\left(t \cdot \sqrt[3]{k}\right)}^{3}}\\
\mathbf{if}\;t \leq -1.95 \cdot 10^{+77}:\\
\;\;\;\;t_1 \cdot t_2\\
\mathbf{elif}\;t \leq 4.4 \cdot 10^{-48}:\\
\;\;\;\;t_1 \cdot \left(2 \cdot \frac{\frac{\ell}{\sin k}}{t \cdot {k}^{2}}\right)\\
\mathbf{else}:\\
\;\;\;\;t_2 \cdot \frac{\ell}{k}\\
\end{array}
\end{array}
if t < -1.9499999999999999e77Initial program 50.9%
Simplified35.5%
associate-*r*43.8%
associate-*r*43.8%
times-frac60.7%
div-inv60.7%
pow-flip65.2%
metadata-eval65.2%
Applied egg-rr65.2%
times-frac65.8%
Simplified65.8%
Taylor expanded in k around 0 60.4%
add-cube-cbrt60.4%
pow360.4%
*-commutative60.4%
cbrt-prod60.4%
unpow360.4%
add-cbrt-cube70.2%
Applied egg-rr70.2%
if -1.9499999999999999e77 < t < 4.40000000000000025e-48Initial program 42.3%
Simplified40.9%
associate-*r*46.8%
associate-*r*46.8%
times-frac52.3%
div-inv52.3%
pow-flip52.2%
metadata-eval52.2%
Applied egg-rr52.2%
times-frac51.6%
Simplified51.6%
Taylor expanded in t around 0 79.6%
associate-*r/79.6%
associate-*r*79.6%
Simplified79.6%
Taylor expanded in l around 0 79.6%
*-commutative79.6%
*-commutative79.6%
associate-*r*79.6%
associate-/r*81.4%
*-commutative81.4%
Simplified81.4%
if 4.40000000000000025e-48 < t Initial program 67.1%
Simplified64.3%
associate-*r*69.9%
associate-*r*69.9%
times-frac76.0%
div-inv76.0%
pow-flip77.3%
metadata-eval77.3%
Applied egg-rr77.3%
times-frac74.8%
Simplified74.8%
Taylor expanded in k around 0 65.9%
Taylor expanded in k around 0 69.9%
add-cube-cbrt69.8%
pow369.8%
*-commutative69.8%
cbrt-prod69.8%
unpow369.8%
add-cbrt-cube72.4%
Applied egg-rr72.4%
Final simplification76.9%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ l (pow (* t (cbrt k)) 3.0))))
(if (<= t -1.95e+77)
(* (/ l (tan k)) t_1)
(if (<= t 4.4e-48)
(/ (* l (* (/ 2.0 (* t (pow k 2.0))) (/ l (sin k)))) (tan k))
(* t_1 (/ l k))))))
double code(double t, double l, double k) {
double t_1 = l / pow((t * cbrt(k)), 3.0);
double tmp;
if (t <= -1.95e+77) {
tmp = (l / tan(k)) * t_1;
} else if (t <= 4.4e-48) {
tmp = (l * ((2.0 / (t * pow(k, 2.0))) * (l / sin(k)))) / tan(k);
} else {
tmp = t_1 * (l / k);
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = l / Math.pow((t * Math.cbrt(k)), 3.0);
double tmp;
if (t <= -1.95e+77) {
tmp = (l / Math.tan(k)) * t_1;
} else if (t <= 4.4e-48) {
tmp = (l * ((2.0 / (t * Math.pow(k, 2.0))) * (l / Math.sin(k)))) / Math.tan(k);
} else {
tmp = t_1 * (l / k);
}
return tmp;
}
function code(t, l, k) t_1 = Float64(l / (Float64(t * cbrt(k)) ^ 3.0)) tmp = 0.0 if (t <= -1.95e+77) tmp = Float64(Float64(l / tan(k)) * t_1); elseif (t <= 4.4e-48) tmp = Float64(Float64(l * Float64(Float64(2.0 / Float64(t * (k ^ 2.0))) * Float64(l / sin(k)))) / tan(k)); else tmp = Float64(t_1 * Float64(l / k)); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[(l / N[Power[N[(t * N[Power[k, 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.95e+77], N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision], If[LessEqual[t, 4.4e-48], N[(N[(l * N[(N[(2.0 / N[(t * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision], N[(t$95$1 * N[(l / k), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\ell}{{\left(t \cdot \sqrt[3]{k}\right)}^{3}}\\
\mathbf{if}\;t \leq -1.95 \cdot 10^{+77}:\\
\;\;\;\;\frac{\ell}{\tan k} \cdot t_1\\
\mathbf{elif}\;t \leq 4.4 \cdot 10^{-48}:\\
\;\;\;\;\frac{\ell \cdot \left(\frac{2}{t \cdot {k}^{2}} \cdot \frac{\ell}{\sin k}\right)}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;t_1 \cdot \frac{\ell}{k}\\
\end{array}
\end{array}
if t < -1.9499999999999999e77Initial program 50.9%
Simplified35.5%
associate-*r*43.8%
associate-*r*43.8%
times-frac60.7%
div-inv60.7%
pow-flip65.2%
metadata-eval65.2%
Applied egg-rr65.2%
times-frac65.8%
Simplified65.8%
Taylor expanded in k around 0 60.4%
add-cube-cbrt60.4%
pow360.4%
*-commutative60.4%
cbrt-prod60.4%
unpow360.4%
add-cbrt-cube70.2%
Applied egg-rr70.2%
if -1.9499999999999999e77 < t < 4.40000000000000025e-48Initial program 42.3%
Simplified40.9%
associate-*r*46.8%
associate-*r*46.8%
times-frac52.3%
div-inv52.3%
pow-flip52.2%
metadata-eval52.2%
Applied egg-rr52.2%
times-frac51.6%
Simplified51.6%
Taylor expanded in t around 0 79.6%
associate-*r/79.6%
associate-*r*79.6%
Simplified79.6%
associate-*r/79.6%
times-frac81.4%
*-commutative81.4%
Applied egg-rr81.4%
if 4.40000000000000025e-48 < t Initial program 67.1%
Simplified64.3%
associate-*r*69.9%
associate-*r*69.9%
times-frac76.0%
div-inv76.0%
pow-flip77.3%
metadata-eval77.3%
Applied egg-rr77.3%
times-frac74.8%
Simplified74.8%
Taylor expanded in k around 0 65.9%
Taylor expanded in k around 0 69.9%
add-cube-cbrt69.8%
pow369.8%
*-commutative69.8%
cbrt-prod69.8%
unpow369.8%
add-cbrt-cube72.4%
Applied egg-rr72.4%
Final simplification76.9%
(FPCore (t l k) :precision binary64 (if (or (<= t -1.5e-56) (not (<= t 2.55e-48))) (* (/ l (pow (* t (cbrt k)) 3.0)) (/ l k)) (* (/ l (tan k)) (* 2.0 (/ l (* t (pow k 3.0)))))))
double code(double t, double l, double k) {
double tmp;
if ((t <= -1.5e-56) || !(t <= 2.55e-48)) {
tmp = (l / pow((t * cbrt(k)), 3.0)) * (l / k);
} else {
tmp = (l / tan(k)) * (2.0 * (l / (t * pow(k, 3.0))));
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if ((t <= -1.5e-56) || !(t <= 2.55e-48)) {
tmp = (l / Math.pow((t * Math.cbrt(k)), 3.0)) * (l / k);
} else {
tmp = (l / Math.tan(k)) * (2.0 * (l / (t * Math.pow(k, 3.0))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if ((t <= -1.5e-56) || !(t <= 2.55e-48)) tmp = Float64(Float64(l / (Float64(t * cbrt(k)) ^ 3.0)) * Float64(l / k)); else tmp = Float64(Float64(l / tan(k)) * Float64(2.0 * Float64(l / Float64(t * (k ^ 3.0))))); end return tmp end
code[t_, l_, k_] := If[Or[LessEqual[t, -1.5e-56], N[Not[LessEqual[t, 2.55e-48]], $MachinePrecision]], N[(N[(l / N[Power[N[(t * N[Power[k, 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision] * N[(l / k), $MachinePrecision]), $MachinePrecision], N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(2.0 * N[(l / N[(t * N[Power[k, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.5 \cdot 10^{-56} \lor \neg \left(t \leq 2.55 \cdot 10^{-48}\right):\\
\;\;\;\;\frac{\ell}{{\left(t \cdot \sqrt[3]{k}\right)}^{3}} \cdot \frac{\ell}{k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{\tan k} \cdot \left(2 \cdot \frac{\ell}{t \cdot {k}^{3}}\right)\\
\end{array}
\end{array}
if t < -1.49999999999999995e-56 or 2.55000000000000006e-48 < t Initial program 62.6%
Simplified56.2%
associate-*r*62.1%
associate-*r*62.1%
times-frac73.9%
div-inv73.9%
pow-flip75.8%
metadata-eval75.8%
Applied egg-rr75.8%
times-frac74.0%
Simplified74.0%
Taylor expanded in k around 0 62.5%
Taylor expanded in k around 0 65.4%
add-cube-cbrt65.3%
pow365.3%
*-commutative65.3%
cbrt-prod65.4%
unpow365.4%
add-cbrt-cube69.4%
Applied egg-rr69.4%
if -1.49999999999999995e-56 < t < 2.55000000000000006e-48Initial program 35.9%
Simplified34.9%
associate-*r*41.5%
associate-*r*41.5%
times-frac43.3%
div-inv43.3%
pow-flip43.3%
metadata-eval43.3%
Applied egg-rr43.3%
times-frac43.3%
Simplified43.3%
Taylor expanded in t around 0 81.3%
associate-*r/81.3%
associate-*r*81.3%
Simplified81.3%
Taylor expanded in k around 0 63.3%
Final simplification66.8%
(FPCore (t l k) :precision binary64 (if (<= k 1.35e+152) (* (/ l (tan k)) (/ l (* k (pow t 3.0)))) (* 2.0 (/ (pow l 2.0) (* t (pow k 4.0))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 1.35e+152) {
tmp = (l / tan(k)) * (l / (k * pow(t, 3.0)));
} else {
tmp = 2.0 * (pow(l, 2.0) / (t * pow(k, 4.0)));
}
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 <= 1.35d+152) then
tmp = (l / tan(k)) * (l / (k * (t ** 3.0d0)))
else
tmp = 2.0d0 * ((l ** 2.0d0) / (t * (k ** 4.0d0)))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 1.35e+152) {
tmp = (l / Math.tan(k)) * (l / (k * Math.pow(t, 3.0)));
} else {
tmp = 2.0 * (Math.pow(l, 2.0) / (t * Math.pow(k, 4.0)));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 1.35e+152: tmp = (l / math.tan(k)) * (l / (k * math.pow(t, 3.0))) else: tmp = 2.0 * (math.pow(l, 2.0) / (t * math.pow(k, 4.0))) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 1.35e+152) tmp = Float64(Float64(l / tan(k)) * Float64(l / Float64(k * (t ^ 3.0)))); else tmp = Float64(2.0 * Float64((l ^ 2.0) / Float64(t * (k ^ 4.0)))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 1.35e+152) tmp = (l / tan(k)) * (l / (k * (t ^ 3.0))); else tmp = 2.0 * ((l ^ 2.0) / (t * (k ^ 4.0))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 1.35e+152], N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[(k * N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] / N[(t * N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.35 \cdot 10^{+152}:\\
\;\;\;\;\frac{\ell}{\tan k} \cdot \frac{\ell}{k \cdot {t}^{3}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{{\ell}^{2}}{t \cdot {k}^{4}}\\
\end{array}
\end{array}
if k < 1.35000000000000007e152Initial program 50.8%
Simplified46.6%
associate-*r*52.8%
associate-*r*52.8%
times-frac61.4%
div-inv61.4%
pow-flip62.6%
metadata-eval62.6%
Applied egg-rr62.6%
times-frac61.4%
Simplified61.4%
Taylor expanded in k around 0 56.2%
if 1.35000000000000007e152 < k Initial program 52.9%
associate-/r*59.2%
add-cube-cbrt59.2%
*-un-lft-identity59.2%
times-frac59.2%
pow259.2%
cbrt-div59.2%
rem-cbrt-cube59.2%
cbrt-div59.2%
rem-cbrt-cube73.4%
Applied egg-rr73.4%
Taylor expanded in t around 0 68.8%
Taylor expanded in k around 0 68.8%
Final simplification57.8%
(FPCore (t l k) :precision binary64 (if (<= k 1.05e+76) (* (/ l (* k (pow t 3.0))) (+ (/ l k) (* -0.3333333333333333 (* l k)))) (* 2.0 (/ (pow l 2.0) (* t (pow k 4.0))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 1.05e+76) {
tmp = (l / (k * pow(t, 3.0))) * ((l / k) + (-0.3333333333333333 * (l * k)));
} else {
tmp = 2.0 * (pow(l, 2.0) / (t * pow(k, 4.0)));
}
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 <= 1.05d+76) then
tmp = (l / (k * (t ** 3.0d0))) * ((l / k) + ((-0.3333333333333333d0) * (l * k)))
else
tmp = 2.0d0 * ((l ** 2.0d0) / (t * (k ** 4.0d0)))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 1.05e+76) {
tmp = (l / (k * Math.pow(t, 3.0))) * ((l / k) + (-0.3333333333333333 * (l * k)));
} else {
tmp = 2.0 * (Math.pow(l, 2.0) / (t * Math.pow(k, 4.0)));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 1.05e+76: tmp = (l / (k * math.pow(t, 3.0))) * ((l / k) + (-0.3333333333333333 * (l * k))) else: tmp = 2.0 * (math.pow(l, 2.0) / (t * math.pow(k, 4.0))) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 1.05e+76) tmp = Float64(Float64(l / Float64(k * (t ^ 3.0))) * Float64(Float64(l / k) + Float64(-0.3333333333333333 * Float64(l * k)))); else tmp = Float64(2.0 * Float64((l ^ 2.0) / Float64(t * (k ^ 4.0)))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 1.05e+76) tmp = (l / (k * (t ^ 3.0))) * ((l / k) + (-0.3333333333333333 * (l * k))); else tmp = 2.0 * ((l ^ 2.0) / (t * (k ^ 4.0))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 1.05e+76], N[(N[(l / N[(k * N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(l / k), $MachinePrecision] + N[(-0.3333333333333333 * N[(l * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] / N[(t * N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.05 \cdot 10^{+76}:\\
\;\;\;\;\frac{\ell}{k \cdot {t}^{3}} \cdot \left(\frac{\ell}{k} + -0.3333333333333333 \cdot \left(\ell \cdot k\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{{\ell}^{2}}{t \cdot {k}^{4}}\\
\end{array}
\end{array}
if k < 1.05000000000000003e76Initial program 52.0%
Simplified47.5%
associate-*r*53.6%
associate-*r*53.6%
times-frac61.9%
div-inv61.9%
pow-flip63.3%
metadata-eval63.3%
Applied egg-rr63.3%
times-frac62.4%
Simplified62.4%
Taylor expanded in k around 0 57.2%
Taylor expanded in k around 0 56.3%
if 1.05000000000000003e76 < k Initial program 47.5%
associate-/r*53.9%
add-cube-cbrt53.9%
*-un-lft-identity53.9%
times-frac53.9%
pow253.9%
cbrt-div53.9%
rem-cbrt-cube53.9%
cbrt-div53.9%
rem-cbrt-cube65.6%
Applied egg-rr65.6%
Taylor expanded in t around 0 66.5%
Taylor expanded in k around 0 56.3%
Final simplification56.3%
(FPCore (t l k) :precision binary64 (* (/ l (pow (* t (cbrt k)) 3.0)) (/ l k)))
double code(double t, double l, double k) {
return (l / pow((t * cbrt(k)), 3.0)) * (l / k);
}
public static double code(double t, double l, double k) {
return (l / Math.pow((t * Math.cbrt(k)), 3.0)) * (l / k);
}
function code(t, l, k) return Float64(Float64(l / (Float64(t * cbrt(k)) ^ 3.0)) * Float64(l / k)) end
code[t_, l_, k_] := N[(N[(l / N[Power[N[(t * N[Power[k, 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision] * N[(l / k), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell}{{\left(t \cdot \sqrt[3]{k}\right)}^{3}} \cdot \frac{\ell}{k}
\end{array}
Initial program 51.1%
Simplified47.1%
associate-*r*53.2%
associate-*r*53.3%
times-frac60.7%
div-inv60.7%
pow-flip61.8%
metadata-eval61.8%
Applied egg-rr61.8%
times-frac60.8%
Simplified60.8%
Taylor expanded in k around 0 56.7%
Taylor expanded in k around 0 57.8%
add-cube-cbrt57.8%
pow357.8%
*-commutative57.8%
cbrt-prod57.8%
unpow357.8%
add-cbrt-cube60.9%
Applied egg-rr60.9%
Final simplification60.9%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ l (* k (pow t 3.0)))))
(if (<= k 1.9e+100)
(* t_1 (+ (/ l k) (* -0.3333333333333333 (* l k))))
(* (/ l k) t_1))))
double code(double t, double l, double k) {
double t_1 = l / (k * pow(t, 3.0));
double tmp;
if (k <= 1.9e+100) {
tmp = t_1 * ((l / k) + (-0.3333333333333333 * (l * k)));
} else {
tmp = (l / k) * t_1;
}
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 = l / (k * (t ** 3.0d0))
if (k <= 1.9d+100) then
tmp = t_1 * ((l / k) + ((-0.3333333333333333d0) * (l * k)))
else
tmp = (l / k) * t_1
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = l / (k * Math.pow(t, 3.0));
double tmp;
if (k <= 1.9e+100) {
tmp = t_1 * ((l / k) + (-0.3333333333333333 * (l * k)));
} else {
tmp = (l / k) * t_1;
}
return tmp;
}
def code(t, l, k): t_1 = l / (k * math.pow(t, 3.0)) tmp = 0 if k <= 1.9e+100: tmp = t_1 * ((l / k) + (-0.3333333333333333 * (l * k))) else: tmp = (l / k) * t_1 return tmp
function code(t, l, k) t_1 = Float64(l / Float64(k * (t ^ 3.0))) tmp = 0.0 if (k <= 1.9e+100) tmp = Float64(t_1 * Float64(Float64(l / k) + Float64(-0.3333333333333333 * Float64(l * k)))); else tmp = Float64(Float64(l / k) * t_1); end return tmp end
function tmp_2 = code(t, l, k) t_1 = l / (k * (t ^ 3.0)); tmp = 0.0; if (k <= 1.9e+100) tmp = t_1 * ((l / k) + (-0.3333333333333333 * (l * k))); else tmp = (l / k) * t_1; end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(l / N[(k * N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, 1.9e+100], N[(t$95$1 * N[(N[(l / k), $MachinePrecision] + N[(-0.3333333333333333 * N[(l * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * t$95$1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\ell}{k \cdot {t}^{3}}\\
\mathbf{if}\;k \leq 1.9 \cdot 10^{+100}:\\
\;\;\;\;t_1 \cdot \left(\frac{\ell}{k} + -0.3333333333333333 \cdot \left(\ell \cdot k\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot t_1\\
\end{array}
\end{array}
if k < 1.89999999999999982e100Initial program 52.0%
Simplified47.5%
associate-*r*53.5%
associate-*r*53.6%
times-frac61.8%
div-inv61.8%
pow-flip63.1%
metadata-eval63.1%
Applied egg-rr63.1%
times-frac62.3%
Simplified62.3%
Taylor expanded in k around 0 57.1%
Taylor expanded in k around 0 56.3%
if 1.89999999999999982e100 < k Initial program 47.3%
Simplified45.2%
associate-*r*51.9%
associate-*r*51.9%
times-frac55.9%
div-inv55.9%
pow-flip55.9%
metadata-eval55.9%
Applied egg-rr55.9%
times-frac54.0%
Simplified54.0%
Taylor expanded in k around 0 54.7%
Taylor expanded in k around 0 56.5%
Final simplification56.3%
(FPCore (t l k) :precision binary64 (* (/ l k) (/ l (* k (* t (pow t 2.0))))))
double code(double t, double l, double k) {
return (l / k) * (l / (k * (t * pow(t, 2.0))));
}
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 / (k * (t * (t ** 2.0d0))))
end function
public static double code(double t, double l, double k) {
return (l / k) * (l / (k * (t * Math.pow(t, 2.0))));
}
def code(t, l, k): return (l / k) * (l / (k * (t * math.pow(t, 2.0))))
function code(t, l, k) return Float64(Float64(l / k) * Float64(l / Float64(k * Float64(t * (t ^ 2.0))))) end
function tmp = code(t, l, k) tmp = (l / k) * (l / (k * (t * (t ^ 2.0)))); end
code[t_, l_, k_] := N[(N[(l / k), $MachinePrecision] * N[(l / N[(k * N[(t * N[Power[t, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell}{k} \cdot \frac{\ell}{k \cdot \left(t \cdot {t}^{2}\right)}
\end{array}
Initial program 51.1%
Simplified47.1%
associate-*r*53.2%
associate-*r*53.3%
times-frac60.7%
div-inv60.7%
pow-flip61.8%
metadata-eval61.8%
Applied egg-rr61.8%
times-frac60.8%
Simplified60.8%
Taylor expanded in k around 0 56.7%
Taylor expanded in k around 0 57.8%
unpow357.8%
pow257.8%
Applied egg-rr57.8%
Final simplification57.8%
(FPCore (t l k) :precision binary64 (* (/ l k) (/ l (* k (pow t 3.0)))))
double code(double t, double l, double k) {
return (l / k) * (l / (k * pow(t, 3.0)));
}
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 / (k * (t ** 3.0d0)))
end function
public static double code(double t, double l, double k) {
return (l / k) * (l / (k * Math.pow(t, 3.0)));
}
def code(t, l, k): return (l / k) * (l / (k * math.pow(t, 3.0)))
function code(t, l, k) return Float64(Float64(l / k) * Float64(l / Float64(k * (t ^ 3.0)))) end
function tmp = code(t, l, k) tmp = (l / k) * (l / (k * (t ^ 3.0))); end
code[t_, l_, k_] := N[(N[(l / k), $MachinePrecision] * N[(l / N[(k * N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell}{k} \cdot \frac{\ell}{k \cdot {t}^{3}}
\end{array}
Initial program 51.1%
Simplified47.1%
associate-*r*53.2%
associate-*r*53.3%
times-frac60.7%
div-inv60.7%
pow-flip61.8%
metadata-eval61.8%
Applied egg-rr61.8%
times-frac60.8%
Simplified60.8%
Taylor expanded in k around 0 56.7%
Taylor expanded in k around 0 57.8%
Final simplification57.8%
herbie shell --seed 2023320
(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))))