
(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 21 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 (pow (/ k t) 2.0)))
(if (<=
(/
2.0
(*
(* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k))
(+ 1.0 (+ 1.0 t_1))))
1e+242)
(/ (* 2.0 (* (/ l (tan k)) (/ (/ l (pow t 3.0)) (sin k)))) (+ 2.0 t_1))
(/ 2.0 (* (/ k l) (* (/ t (/ l k)) (* (sin k) (tan k))))))))
double code(double t, double l, double k) {
double t_1 = pow((k / t), 2.0);
double tmp;
if ((2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * (1.0 + (1.0 + t_1)))) <= 1e+242) {
tmp = (2.0 * ((l / tan(k)) * ((l / pow(t, 3.0)) / sin(k)))) / (2.0 + t_1);
} else {
tmp = 2.0 / ((k / l) * ((t / (l / k)) * (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) :: t_1
real(8) :: tmp
t_1 = (k / t) ** 2.0d0
if ((2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * (1.0d0 + (1.0d0 + t_1)))) <= 1d+242) then
tmp = (2.0d0 * ((l / tan(k)) * ((l / (t ** 3.0d0)) / sin(k)))) / (2.0d0 + t_1)
else
tmp = 2.0d0 / ((k / l) * ((t / (l / k)) * (sin(k) * tan(k))))
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 tmp;
if ((2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * (1.0 + (1.0 + t_1)))) <= 1e+242) {
tmp = (2.0 * ((l / Math.tan(k)) * ((l / Math.pow(t, 3.0)) / Math.sin(k)))) / (2.0 + t_1);
} else {
tmp = 2.0 / ((k / l) * ((t / (l / k)) * (Math.sin(k) * Math.tan(k))));
}
return tmp;
}
def code(t, l, k): t_1 = math.pow((k / t), 2.0) tmp = 0 if (2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * (1.0 + (1.0 + t_1)))) <= 1e+242: tmp = (2.0 * ((l / math.tan(k)) * ((l / math.pow(t, 3.0)) / math.sin(k)))) / (2.0 + t_1) else: tmp = 2.0 / ((k / l) * ((t / (l / k)) * (math.sin(k) * math.tan(k)))) return tmp
function code(t, l, k) t_1 = Float64(k / t) ^ 2.0 tmp = 0.0 if (Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(1.0 + Float64(1.0 + t_1)))) <= 1e+242) tmp = Float64(Float64(2.0 * Float64(Float64(l / tan(k)) * Float64(Float64(l / (t ^ 3.0)) / sin(k)))) / Float64(2.0 + t_1)); else tmp = Float64(2.0 / Float64(Float64(k / l) * Float64(Float64(t / Float64(l / k)) * Float64(sin(k) * tan(k))))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = (k / t) ^ 2.0; tmp = 0.0; if ((2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * (1.0 + (1.0 + t_1)))) <= 1e+242) tmp = (2.0 * ((l / tan(k)) * ((l / (t ^ 3.0)) / sin(k)))) / (2.0 + t_1); else tmp = 2.0 / ((k / l) * ((t / (l / k)) * (sin(k) * tan(k)))); 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[(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[(1.0 + N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1e+242], N[(N[(2.0 * N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(l / N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision] / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + t$95$1), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k / l), $MachinePrecision] * N[(N[(t / N[(l / k), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {\left(\frac{k}{t}\right)}^{2}\\
\mathbf{if}\;\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(1 + \left(1 + t_1\right)\right)} \leq 10^{+242}:\\
\;\;\;\;\frac{2 \cdot \left(\frac{\ell}{\tan k} \cdot \frac{\frac{\ell}{{t}^{3}}}{\sin k}\right)}{2 + t_1}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{k}{\ell} \cdot \left(\frac{t}{\frac{\ell}{k}} \cdot \left(\sin k \cdot \tan k\right)\right)}\\
\end{array}
\end{array}
if (/.f64 2 (*.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))) < 1.00000000000000005e242Initial program 81.4%
associate-*l*81.4%
associate-/l/81.4%
*-commutative81.4%
associate-*r/82.4%
associate-/l*82.4%
associate-/r/73.1%
Simplified75.0%
expm1-log1p-u58.3%
expm1-udef53.8%
associate-*l/53.8%
*-commutative53.8%
Applied egg-rr53.8%
expm1-def58.3%
expm1-log1p75.0%
times-frac75.0%
associate-*l/75.0%
associate-*r/73.1%
unpow273.1%
associate-/l/73.3%
Simplified86.6%
if 1.00000000000000005e242 < (/.f64 2 (*.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 17.5%
*-commutative17.5%
associate-*l*17.5%
associate-*r*17.5%
+-commutative17.5%
associate-+r+17.5%
metadata-eval17.5%
Simplified17.5%
Taylor expanded in k around inf 58.7%
associate-/l*55.2%
associate-/r/55.3%
unpow255.3%
unpow255.3%
Simplified55.3%
pow155.3%
associate-*l*55.2%
times-frac77.2%
Applied egg-rr77.2%
unpow177.2%
associate-*r*78.8%
*-commutative78.8%
associate-*r*89.9%
*-commutative89.9%
Simplified89.9%
pow189.9%
Applied egg-rr89.9%
unpow189.9%
associate-*l*89.9%
associate-*r/86.5%
associate-/l*89.9%
Simplified89.9%
Final simplification88.0%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ t (/ l k))) (t_2 (* (sin k) (tan k))))
(if (<= k 1.7e-131)
(/
2.0
(*
(+ 1.0 (+ 1.0 (pow (/ k t) 2.0)))
(* (tan k) (/ (* (pow t 3.0) (/ k l)) l))))
(if (<= k 3.3e+27)
(/ 2.0 (* t_2 (* (/ k l) (cbrt (* t_1 (* t (* (/ k l) t_1)))))))
(if (<= k 5.5e+170)
(/ 2.0 (* t_2 (* (/ k l) (/ (* t k) l))))
(/ 2.0 (* t_2 (* (/ k l) (* t (/ k l))))))))))
double code(double t, double l, double k) {
double t_1 = t / (l / k);
double t_2 = sin(k) * tan(k);
double tmp;
if (k <= 1.7e-131) {
tmp = 2.0 / ((1.0 + (1.0 + pow((k / t), 2.0))) * (tan(k) * ((pow(t, 3.0) * (k / l)) / l)));
} else if (k <= 3.3e+27) {
tmp = 2.0 / (t_2 * ((k / l) * cbrt((t_1 * (t * ((k / l) * t_1))))));
} else if (k <= 5.5e+170) {
tmp = 2.0 / (t_2 * ((k / l) * ((t * k) / l)));
} else {
tmp = 2.0 / (t_2 * ((k / l) * (t * (k / l))));
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = t / (l / k);
double t_2 = Math.sin(k) * Math.tan(k);
double tmp;
if (k <= 1.7e-131) {
tmp = 2.0 / ((1.0 + (1.0 + Math.pow((k / t), 2.0))) * (Math.tan(k) * ((Math.pow(t, 3.0) * (k / l)) / l)));
} else if (k <= 3.3e+27) {
tmp = 2.0 / (t_2 * ((k / l) * Math.cbrt((t_1 * (t * ((k / l) * t_1))))));
} else if (k <= 5.5e+170) {
tmp = 2.0 / (t_2 * ((k / l) * ((t * k) / l)));
} else {
tmp = 2.0 / (t_2 * ((k / l) * (t * (k / l))));
}
return tmp;
}
function code(t, l, k) t_1 = Float64(t / Float64(l / k)) t_2 = Float64(sin(k) * tan(k)) tmp = 0.0 if (k <= 1.7e-131) tmp = Float64(2.0 / Float64(Float64(1.0 + Float64(1.0 + (Float64(k / t) ^ 2.0))) * Float64(tan(k) * Float64(Float64((t ^ 3.0) * Float64(k / l)) / l)))); elseif (k <= 3.3e+27) tmp = Float64(2.0 / Float64(t_2 * Float64(Float64(k / l) * cbrt(Float64(t_1 * Float64(t * Float64(Float64(k / l) * t_1))))))); elseif (k <= 5.5e+170) tmp = Float64(2.0 / Float64(t_2 * Float64(Float64(k / l) * Float64(Float64(t * k) / l)))); else tmp = Float64(2.0 / Float64(t_2 * Float64(Float64(k / l) * Float64(t * Float64(k / l))))); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[(t / N[(l / k), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, 1.7e-131], 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[(N[(N[Power[t, 3.0], $MachinePrecision] * N[(k / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 3.3e+27], N[(2.0 / N[(t$95$2 * N[(N[(k / l), $MachinePrecision] * N[Power[N[(t$95$1 * N[(t * N[(N[(k / l), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 5.5e+170], N[(2.0 / N[(t$95$2 * N[(N[(k / l), $MachinePrecision] * N[(N[(t * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(t$95$2 * N[(N[(k / l), $MachinePrecision] * N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{\frac{\ell}{k}}\\
t_2 := \sin k \cdot \tan k\\
\mathbf{if}\;k \leq 1.7 \cdot 10^{-131}:\\
\;\;\;\;\frac{2}{\left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right) \cdot \left(\tan k \cdot \frac{{t}^{3} \cdot \frac{k}{\ell}}{\ell}\right)}\\
\mathbf{elif}\;k \leq 3.3 \cdot 10^{+27}:\\
\;\;\;\;\frac{2}{t_2 \cdot \left(\frac{k}{\ell} \cdot \sqrt[3]{t_1 \cdot \left(t \cdot \left(\frac{k}{\ell} \cdot t_1\right)\right)}\right)}\\
\mathbf{elif}\;k \leq 5.5 \cdot 10^{+170}:\\
\;\;\;\;\frac{2}{t_2 \cdot \left(\frac{k}{\ell} \cdot \frac{t \cdot k}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t_2 \cdot \left(\frac{k}{\ell} \cdot \left(t \cdot \frac{k}{\ell}\right)\right)}\\
\end{array}
\end{array}
if k < 1.69999999999999998e-131Initial program 53.1%
Taylor expanded in k around 0 52.8%
*-commutative52.8%
unpow252.8%
times-frac63.5%
Simplified63.5%
associate-*l/63.5%
Applied egg-rr63.5%
if 1.69999999999999998e-131 < k < 3.2999999999999998e27Initial program 66.4%
*-commutative66.4%
associate-*l*66.3%
associate-*r*66.2%
+-commutative66.2%
associate-+r+66.2%
metadata-eval66.2%
Simplified66.2%
Taylor expanded in k around inf 66.1%
associate-/l*72.5%
associate-/r/66.1%
unpow266.1%
unpow266.1%
Simplified66.1%
pow166.1%
associate-*l*66.0%
times-frac66.3%
Applied egg-rr66.3%
unpow166.3%
associate-*r*66.3%
*-commutative66.3%
associate-*r*69.4%
*-commutative69.4%
Simplified69.4%
add-cbrt-cube82.3%
Applied egg-rr82.3%
associate-*l*82.3%
associate-*r/82.2%
associate-/l*82.3%
associate-*l*82.3%
associate-*r/82.2%
associate-/l*82.3%
Simplified82.3%
if 3.2999999999999998e27 < k < 5.4999999999999999e170Initial program 51.7%
*-commutative51.7%
associate-*l*51.8%
associate-*r*51.7%
+-commutative51.7%
associate-+r+51.7%
metadata-eval51.7%
Simplified51.7%
Taylor expanded in k around inf 85.2%
associate-/l*78.5%
associate-/r/78.3%
unpow278.3%
unpow278.3%
Simplified78.3%
pow178.3%
associate-*l*78.3%
times-frac78.4%
Applied egg-rr78.4%
unpow178.4%
associate-*r*78.4%
*-commutative78.4%
associate-*r*88.9%
*-commutative88.9%
Simplified88.9%
associate-*r/96.1%
Applied egg-rr96.1%
if 5.4999999999999999e170 < k Initial program 50.2%
*-commutative50.2%
associate-*l*50.2%
associate-*r*50.2%
+-commutative50.2%
associate-+r+50.2%
metadata-eval50.2%
Simplified50.2%
Taylor expanded in k around inf 66.1%
associate-/l*65.7%
associate-/r/66.1%
unpow266.1%
unpow266.1%
Simplified66.1%
pow166.1%
associate-*l*66.1%
times-frac92.4%
Applied egg-rr92.4%
unpow192.4%
associate-*r*92.2%
*-commutative92.2%
associate-*r*99.5%
*-commutative99.5%
Simplified99.5%
Final simplification72.8%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (* (sin k) (tan k))))
(if (<= k 1.8e-130)
(/
2.0
(*
(+ 1.0 (+ 1.0 (pow (/ k t) 2.0)))
(* (tan k) (* (/ k l) (/ (pow t 3.0) l)))))
(if (<= k 1.75e-45)
(/ 2.0 (* (/ k (/ (cos k) k)) (/ (pow (sin k) 2.0) (/ l (/ t l)))))
(if (<= k 290.0)
(* (/ l (pow t 3.0)) (/ l (* k k)))
(if (<= k 1.8e+170)
(/ 2.0 (* t_1 (* (/ k l) (/ (* t k) l))))
(/ 2.0 (* t_1 (* (/ k l) (* t (/ k l)))))))))))
double code(double t, double l, double k) {
double t_1 = sin(k) * tan(k);
double tmp;
if (k <= 1.8e-130) {
tmp = 2.0 / ((1.0 + (1.0 + pow((k / t), 2.0))) * (tan(k) * ((k / l) * (pow(t, 3.0) / l))));
} else if (k <= 1.75e-45) {
tmp = 2.0 / ((k / (cos(k) / k)) * (pow(sin(k), 2.0) / (l / (t / l))));
} else if (k <= 290.0) {
tmp = (l / pow(t, 3.0)) * (l / (k * k));
} else if (k <= 1.8e+170) {
tmp = 2.0 / (t_1 * ((k / l) * ((t * k) / l)));
} else {
tmp = 2.0 / (t_1 * ((k / l) * (t * (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 = sin(k) * tan(k)
if (k <= 1.8d-130) then
tmp = 2.0d0 / ((1.0d0 + (1.0d0 + ((k / t) ** 2.0d0))) * (tan(k) * ((k / l) * ((t ** 3.0d0) / l))))
else if (k <= 1.75d-45) then
tmp = 2.0d0 / ((k / (cos(k) / k)) * ((sin(k) ** 2.0d0) / (l / (t / l))))
else if (k <= 290.0d0) then
tmp = (l / (t ** 3.0d0)) * (l / (k * k))
else if (k <= 1.8d+170) then
tmp = 2.0d0 / (t_1 * ((k / l) * ((t * k) / l)))
else
tmp = 2.0d0 / (t_1 * ((k / l) * (t * (k / l))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = Math.sin(k) * Math.tan(k);
double tmp;
if (k <= 1.8e-130) {
tmp = 2.0 / ((1.0 + (1.0 + Math.pow((k / t), 2.0))) * (Math.tan(k) * ((k / l) * (Math.pow(t, 3.0) / l))));
} else if (k <= 1.75e-45) {
tmp = 2.0 / ((k / (Math.cos(k) / k)) * (Math.pow(Math.sin(k), 2.0) / (l / (t / l))));
} else if (k <= 290.0) {
tmp = (l / Math.pow(t, 3.0)) * (l / (k * k));
} else if (k <= 1.8e+170) {
tmp = 2.0 / (t_1 * ((k / l) * ((t * k) / l)));
} else {
tmp = 2.0 / (t_1 * ((k / l) * (t * (k / l))));
}
return tmp;
}
def code(t, l, k): t_1 = math.sin(k) * math.tan(k) tmp = 0 if k <= 1.8e-130: tmp = 2.0 / ((1.0 + (1.0 + math.pow((k / t), 2.0))) * (math.tan(k) * ((k / l) * (math.pow(t, 3.0) / l)))) elif k <= 1.75e-45: tmp = 2.0 / ((k / (math.cos(k) / k)) * (math.pow(math.sin(k), 2.0) / (l / (t / l)))) elif k <= 290.0: tmp = (l / math.pow(t, 3.0)) * (l / (k * k)) elif k <= 1.8e+170: tmp = 2.0 / (t_1 * ((k / l) * ((t * k) / l))) else: tmp = 2.0 / (t_1 * ((k / l) * (t * (k / l)))) return tmp
function code(t, l, k) t_1 = Float64(sin(k) * tan(k)) tmp = 0.0 if (k <= 1.8e-130) tmp = Float64(2.0 / Float64(Float64(1.0 + Float64(1.0 + (Float64(k / t) ^ 2.0))) * Float64(tan(k) * Float64(Float64(k / l) * Float64((t ^ 3.0) / l))))); elseif (k <= 1.75e-45) tmp = Float64(2.0 / Float64(Float64(k / Float64(cos(k) / k)) * Float64((sin(k) ^ 2.0) / Float64(l / Float64(t / l))))); elseif (k <= 290.0) tmp = Float64(Float64(l / (t ^ 3.0)) * Float64(l / Float64(k * k))); elseif (k <= 1.8e+170) tmp = Float64(2.0 / Float64(t_1 * Float64(Float64(k / l) * Float64(Float64(t * k) / l)))); else tmp = Float64(2.0 / Float64(t_1 * Float64(Float64(k / l) * Float64(t * Float64(k / l))))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = sin(k) * tan(k); tmp = 0.0; if (k <= 1.8e-130) tmp = 2.0 / ((1.0 + (1.0 + ((k / t) ^ 2.0))) * (tan(k) * ((k / l) * ((t ^ 3.0) / l)))); elseif (k <= 1.75e-45) tmp = 2.0 / ((k / (cos(k) / k)) * ((sin(k) ^ 2.0) / (l / (t / l)))); elseif (k <= 290.0) tmp = (l / (t ^ 3.0)) * (l / (k * k)); elseif (k <= 1.8e+170) tmp = 2.0 / (t_1 * ((k / l) * ((t * k) / l))); else tmp = 2.0 / (t_1 * ((k / l) * (t * (k / l)))); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, 1.8e-130], 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[(N[(k / l), $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 1.75e-45], N[(2.0 / N[(N[(k / N[(N[Cos[k], $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] / N[(l / N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 290.0], N[(N[(l / N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision] * N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 1.8e+170], N[(2.0 / N[(t$95$1 * N[(N[(k / l), $MachinePrecision] * N[(N[(t * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(t$95$1 * N[(N[(k / l), $MachinePrecision] * N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sin k \cdot \tan k\\
\mathbf{if}\;k \leq 1.8 \cdot 10^{-130}:\\
\;\;\;\;\frac{2}{\left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right) \cdot \left(\tan k \cdot \left(\frac{k}{\ell} \cdot \frac{{t}^{3}}{\ell}\right)\right)}\\
\mathbf{elif}\;k \leq 1.75 \cdot 10^{-45}:\\
\;\;\;\;\frac{2}{\frac{k}{\frac{\cos k}{k}} \cdot \frac{{\sin k}^{2}}{\frac{\ell}{\frac{t}{\ell}}}}\\
\mathbf{elif}\;k \leq 290:\\
\;\;\;\;\frac{\ell}{{t}^{3}} \cdot \frac{\ell}{k \cdot k}\\
\mathbf{elif}\;k \leq 1.8 \cdot 10^{+170}:\\
\;\;\;\;\frac{2}{t_1 \cdot \left(\frac{k}{\ell} \cdot \frac{t \cdot k}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t_1 \cdot \left(\frac{k}{\ell} \cdot \left(t \cdot \frac{k}{\ell}\right)\right)}\\
\end{array}
\end{array}
if k < 1.8000000000000001e-130Initial program 53.1%
Taylor expanded in k around 0 52.8%
*-commutative52.8%
unpow252.8%
times-frac63.5%
Simplified63.5%
if 1.8000000000000001e-130 < k < 1.75e-45Initial program 51.9%
*-commutative51.9%
associate-*l*51.8%
associate-*r*51.8%
+-commutative51.8%
associate-+r+51.8%
metadata-eval51.8%
Simplified51.8%
Taylor expanded in k around inf 44.2%
times-frac59.6%
unpow259.6%
associate-/l*59.5%
associate-/l*75.6%
unpow275.6%
associate-/l*76.1%
Simplified76.1%
if 1.75e-45 < k < 290Initial program 72.4%
associate-*l*72.4%
associate-/l/72.4%
*-commutative72.4%
associate-*r/72.4%
associate-/l*72.4%
associate-/r/72.4%
Simplified72.3%
Taylor expanded in k around 0 81.8%
unpow281.8%
unpow281.8%
Simplified81.8%
times-frac82.0%
Applied egg-rr82.0%
if 290 < k < 1.8e170Initial program 57.5%
*-commutative57.5%
associate-*l*57.5%
associate-*r*57.4%
+-commutative57.4%
associate-+r+57.4%
metadata-eval57.4%
Simplified57.4%
Taylor expanded in k around inf 87.9%
associate-/l*82.4%
associate-/r/82.2%
unpow282.2%
unpow282.2%
Simplified82.2%
pow182.2%
associate-*l*82.2%
times-frac82.3%
Applied egg-rr82.3%
unpow182.3%
associate-*r*82.3%
*-commutative82.3%
associate-*r*90.9%
*-commutative90.9%
Simplified90.9%
associate-*r/96.8%
Applied egg-rr96.8%
if 1.8e170 < k Initial program 50.2%
*-commutative50.2%
associate-*l*50.2%
associate-*r*50.2%
+-commutative50.2%
associate-+r+50.2%
metadata-eval50.2%
Simplified50.2%
Taylor expanded in k around inf 66.1%
associate-/l*65.7%
associate-/r/66.1%
unpow266.1%
unpow266.1%
Simplified66.1%
pow166.1%
associate-*l*66.1%
times-frac92.4%
Applied egg-rr92.4%
unpow192.4%
associate-*r*92.2%
*-commutative92.2%
associate-*r*99.5%
*-commutative99.5%
Simplified99.5%
Final simplification72.9%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (* (sin k) (tan k))))
(if (<= k 1.7e-130)
(/
2.0
(*
(+ 1.0 (+ 1.0 (pow (/ k t) 2.0)))
(* (tan k) (/ (* (pow t 3.0) (/ k l)) l))))
(if (<= k 8.8e-45)
(/ 2.0 (* (/ k (/ (cos k) k)) (/ (pow (sin k) 2.0) (/ l (/ t l)))))
(if (<= k 210.0)
(* (/ l (pow t 3.0)) (/ l (* k k)))
(if (<= k 3.8e+170)
(/ 2.0 (* t_1 (* (/ k l) (/ (* t k) l))))
(/ 2.0 (* t_1 (* (/ k l) (* t (/ k l)))))))))))
double code(double t, double l, double k) {
double t_1 = sin(k) * tan(k);
double tmp;
if (k <= 1.7e-130) {
tmp = 2.0 / ((1.0 + (1.0 + pow((k / t), 2.0))) * (tan(k) * ((pow(t, 3.0) * (k / l)) / l)));
} else if (k <= 8.8e-45) {
tmp = 2.0 / ((k / (cos(k) / k)) * (pow(sin(k), 2.0) / (l / (t / l))));
} else if (k <= 210.0) {
tmp = (l / pow(t, 3.0)) * (l / (k * k));
} else if (k <= 3.8e+170) {
tmp = 2.0 / (t_1 * ((k / l) * ((t * k) / l)));
} else {
tmp = 2.0 / (t_1 * ((k / l) * (t * (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 = sin(k) * tan(k)
if (k <= 1.7d-130) then
tmp = 2.0d0 / ((1.0d0 + (1.0d0 + ((k / t) ** 2.0d0))) * (tan(k) * (((t ** 3.0d0) * (k / l)) / l)))
else if (k <= 8.8d-45) then
tmp = 2.0d0 / ((k / (cos(k) / k)) * ((sin(k) ** 2.0d0) / (l / (t / l))))
else if (k <= 210.0d0) then
tmp = (l / (t ** 3.0d0)) * (l / (k * k))
else if (k <= 3.8d+170) then
tmp = 2.0d0 / (t_1 * ((k / l) * ((t * k) / l)))
else
tmp = 2.0d0 / (t_1 * ((k / l) * (t * (k / l))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = Math.sin(k) * Math.tan(k);
double tmp;
if (k <= 1.7e-130) {
tmp = 2.0 / ((1.0 + (1.0 + Math.pow((k / t), 2.0))) * (Math.tan(k) * ((Math.pow(t, 3.0) * (k / l)) / l)));
} else if (k <= 8.8e-45) {
tmp = 2.0 / ((k / (Math.cos(k) / k)) * (Math.pow(Math.sin(k), 2.0) / (l / (t / l))));
} else if (k <= 210.0) {
tmp = (l / Math.pow(t, 3.0)) * (l / (k * k));
} else if (k <= 3.8e+170) {
tmp = 2.0 / (t_1 * ((k / l) * ((t * k) / l)));
} else {
tmp = 2.0 / (t_1 * ((k / l) * (t * (k / l))));
}
return tmp;
}
def code(t, l, k): t_1 = math.sin(k) * math.tan(k) tmp = 0 if k <= 1.7e-130: tmp = 2.0 / ((1.0 + (1.0 + math.pow((k / t), 2.0))) * (math.tan(k) * ((math.pow(t, 3.0) * (k / l)) / l))) elif k <= 8.8e-45: tmp = 2.0 / ((k / (math.cos(k) / k)) * (math.pow(math.sin(k), 2.0) / (l / (t / l)))) elif k <= 210.0: tmp = (l / math.pow(t, 3.0)) * (l / (k * k)) elif k <= 3.8e+170: tmp = 2.0 / (t_1 * ((k / l) * ((t * k) / l))) else: tmp = 2.0 / (t_1 * ((k / l) * (t * (k / l)))) return tmp
function code(t, l, k) t_1 = Float64(sin(k) * tan(k)) tmp = 0.0 if (k <= 1.7e-130) tmp = Float64(2.0 / Float64(Float64(1.0 + Float64(1.0 + (Float64(k / t) ^ 2.0))) * Float64(tan(k) * Float64(Float64((t ^ 3.0) * Float64(k / l)) / l)))); elseif (k <= 8.8e-45) tmp = Float64(2.0 / Float64(Float64(k / Float64(cos(k) / k)) * Float64((sin(k) ^ 2.0) / Float64(l / Float64(t / l))))); elseif (k <= 210.0) tmp = Float64(Float64(l / (t ^ 3.0)) * Float64(l / Float64(k * k))); elseif (k <= 3.8e+170) tmp = Float64(2.0 / Float64(t_1 * Float64(Float64(k / l) * Float64(Float64(t * k) / l)))); else tmp = Float64(2.0 / Float64(t_1 * Float64(Float64(k / l) * Float64(t * Float64(k / l))))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = sin(k) * tan(k); tmp = 0.0; if (k <= 1.7e-130) tmp = 2.0 / ((1.0 + (1.0 + ((k / t) ^ 2.0))) * (tan(k) * (((t ^ 3.0) * (k / l)) / l))); elseif (k <= 8.8e-45) tmp = 2.0 / ((k / (cos(k) / k)) * ((sin(k) ^ 2.0) / (l / (t / l)))); elseif (k <= 210.0) tmp = (l / (t ^ 3.0)) * (l / (k * k)); elseif (k <= 3.8e+170) tmp = 2.0 / (t_1 * ((k / l) * ((t * k) / l))); else tmp = 2.0 / (t_1 * ((k / l) * (t * (k / l)))); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, 1.7e-130], 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[(N[(N[Power[t, 3.0], $MachinePrecision] * N[(k / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 8.8e-45], N[(2.0 / N[(N[(k / N[(N[Cos[k], $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] / N[(l / N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 210.0], N[(N[(l / N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision] * N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 3.8e+170], N[(2.0 / N[(t$95$1 * N[(N[(k / l), $MachinePrecision] * N[(N[(t * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(t$95$1 * N[(N[(k / l), $MachinePrecision] * N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sin k \cdot \tan k\\
\mathbf{if}\;k \leq 1.7 \cdot 10^{-130}:\\
\;\;\;\;\frac{2}{\left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right) \cdot \left(\tan k \cdot \frac{{t}^{3} \cdot \frac{k}{\ell}}{\ell}\right)}\\
\mathbf{elif}\;k \leq 8.8 \cdot 10^{-45}:\\
\;\;\;\;\frac{2}{\frac{k}{\frac{\cos k}{k}} \cdot \frac{{\sin k}^{2}}{\frac{\ell}{\frac{t}{\ell}}}}\\
\mathbf{elif}\;k \leq 210:\\
\;\;\;\;\frac{\ell}{{t}^{3}} \cdot \frac{\ell}{k \cdot k}\\
\mathbf{elif}\;k \leq 3.8 \cdot 10^{+170}:\\
\;\;\;\;\frac{2}{t_1 \cdot \left(\frac{k}{\ell} \cdot \frac{t \cdot k}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t_1 \cdot \left(\frac{k}{\ell} \cdot \left(t \cdot \frac{k}{\ell}\right)\right)}\\
\end{array}
\end{array}
if k < 1.70000000000000003e-130Initial program 53.1%
Taylor expanded in k around 0 52.8%
*-commutative52.8%
unpow252.8%
times-frac63.5%
Simplified63.5%
associate-*l/63.5%
Applied egg-rr63.5%
if 1.70000000000000003e-130 < k < 8.79999999999999974e-45Initial program 51.9%
*-commutative51.9%
associate-*l*51.8%
associate-*r*51.8%
+-commutative51.8%
associate-+r+51.8%
metadata-eval51.8%
Simplified51.8%
Taylor expanded in k around inf 44.2%
times-frac59.6%
unpow259.6%
associate-/l*59.5%
associate-/l*75.6%
unpow275.6%
associate-/l*76.1%
Simplified76.1%
if 8.79999999999999974e-45 < k < 210Initial program 72.4%
associate-*l*72.4%
associate-/l/72.4%
*-commutative72.4%
associate-*r/72.4%
associate-/l*72.4%
associate-/r/72.4%
Simplified72.3%
Taylor expanded in k around 0 81.8%
unpow281.8%
unpow281.8%
Simplified81.8%
times-frac82.0%
Applied egg-rr82.0%
if 210 < k < 3.7999999999999998e170Initial program 57.5%
*-commutative57.5%
associate-*l*57.5%
associate-*r*57.4%
+-commutative57.4%
associate-+r+57.4%
metadata-eval57.4%
Simplified57.4%
Taylor expanded in k around inf 87.9%
associate-/l*82.4%
associate-/r/82.2%
unpow282.2%
unpow282.2%
Simplified82.2%
pow182.2%
associate-*l*82.2%
times-frac82.3%
Applied egg-rr82.3%
unpow182.3%
associate-*r*82.3%
*-commutative82.3%
associate-*r*90.9%
*-commutative90.9%
Simplified90.9%
associate-*r/96.8%
Applied egg-rr96.8%
if 3.7999999999999998e170 < k Initial program 50.2%
*-commutative50.2%
associate-*l*50.2%
associate-*r*50.2%
+-commutative50.2%
associate-+r+50.2%
metadata-eval50.2%
Simplified50.2%
Taylor expanded in k around inf 66.1%
associate-/l*65.7%
associate-/r/66.1%
unpow266.1%
unpow266.1%
Simplified66.1%
pow166.1%
associate-*l*66.1%
times-frac92.4%
Applied egg-rr92.4%
unpow192.4%
associate-*r*92.2%
*-commutative92.2%
associate-*r*99.5%
*-commutative99.5%
Simplified99.5%
Final simplification72.9%
(FPCore (t l k) :precision binary64 (if (<= k 2.26e-69) (/ 2.0 (* 2.0 (/ k (/ (/ l (/ (pow t 3.0) l)) k)))) (/ 2.0 (* (sin k) (* (tan k) (* (/ k l) (* k (/ t l))))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 2.26e-69) {
tmp = 2.0 / (2.0 * (k / ((l / (pow(t, 3.0) / l)) / k)));
} else {
tmp = 2.0 / (sin(k) * (tan(k) * ((k / l) * (k * (t / 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) :: tmp
if (k <= 2.26d-69) then
tmp = 2.0d0 / (2.0d0 * (k / ((l / ((t ** 3.0d0) / l)) / k)))
else
tmp = 2.0d0 / (sin(k) * (tan(k) * ((k / l) * (k * (t / l)))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 2.26e-69) {
tmp = 2.0 / (2.0 * (k / ((l / (Math.pow(t, 3.0) / l)) / k)));
} else {
tmp = 2.0 / (Math.sin(k) * (Math.tan(k) * ((k / l) * (k * (t / l)))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 2.26e-69: tmp = 2.0 / (2.0 * (k / ((l / (math.pow(t, 3.0) / l)) / k))) else: tmp = 2.0 / (math.sin(k) * (math.tan(k) * ((k / l) * (k * (t / l))))) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 2.26e-69) tmp = Float64(2.0 / Float64(2.0 * Float64(k / Float64(Float64(l / Float64((t ^ 3.0) / l)) / k)))); else tmp = Float64(2.0 / Float64(sin(k) * Float64(tan(k) * Float64(Float64(k / l) * Float64(k * Float64(t / l)))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 2.26e-69) tmp = 2.0 / (2.0 * (k / ((l / ((t ^ 3.0) / l)) / k))); else tmp = 2.0 / (sin(k) * (tan(k) * ((k / l) * (k * (t / l))))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 2.26e-69], N[(2.0 / N[(2.0 * N[(k / N[(N[(l / N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Sin[k], $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(N[(k / l), $MachinePrecision] * N[(k * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 2.26 \cdot 10^{-69}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{k}{\frac{\frac{\ell}{\frac{{t}^{3}}{\ell}}}{k}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\sin k \cdot \left(\tan k \cdot \left(\frac{k}{\ell} \cdot \left(k \cdot \frac{t}{\ell}\right)\right)\right)}\\
\end{array}
\end{array}
if k < 2.2599999999999999e-69Initial program 53.5%
*-commutative53.5%
associate-*l*46.8%
associate-*r*46.8%
+-commutative46.8%
associate-+r+46.8%
metadata-eval46.8%
Simplified46.8%
Taylor expanded in k around 0 47.7%
associate-/l*47.2%
unpow247.2%
unpow247.2%
associate-*r/56.2%
associate-/l*63.8%
associate-*r/54.8%
associate-/l*63.8%
Simplified63.8%
if 2.2599999999999999e-69 < k Initial program 55.9%
*-commutative55.9%
associate-*l*55.9%
associate-*r*55.9%
+-commutative55.9%
associate-+r+55.9%
metadata-eval55.9%
Simplified55.9%
Taylor expanded in k around inf 75.1%
associate-/l*72.5%
associate-/r/72.6%
unpow272.6%
unpow272.6%
Simplified72.6%
pow172.6%
associate-*l*72.6%
times-frac82.0%
Applied egg-rr82.0%
unpow182.0%
associate-*r*81.9%
*-commutative81.9%
associate-*r*89.7%
*-commutative89.7%
Simplified89.7%
pow189.7%
Applied egg-rr89.7%
unpow189.7%
*-commutative89.7%
associate-*l*89.8%
associate-*r/88.6%
associate-/l*89.8%
associate-/r/87.9%
Simplified87.9%
Final simplification70.7%
(FPCore (t l k) :precision binary64 (if (<= k 1.42e-68) (/ 2.0 (* 2.0 (/ k (/ (/ l (/ (pow t 3.0) l)) k)))) (/ 2.0 (* (* (sin k) (tan k)) (* (/ k l) (* t (/ k l)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 1.42e-68) {
tmp = 2.0 / (2.0 * (k / ((l / (pow(t, 3.0) / l)) / k)));
} else {
tmp = 2.0 / ((sin(k) * tan(k)) * ((k / l) * (t * (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) :: tmp
if (k <= 1.42d-68) then
tmp = 2.0d0 / (2.0d0 * (k / ((l / ((t ** 3.0d0) / l)) / k)))
else
tmp = 2.0d0 / ((sin(k) * tan(k)) * ((k / l) * (t * (k / l))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 1.42e-68) {
tmp = 2.0 / (2.0 * (k / ((l / (Math.pow(t, 3.0) / l)) / k)));
} else {
tmp = 2.0 / ((Math.sin(k) * Math.tan(k)) * ((k / l) * (t * (k / l))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 1.42e-68: tmp = 2.0 / (2.0 * (k / ((l / (math.pow(t, 3.0) / l)) / k))) else: tmp = 2.0 / ((math.sin(k) * math.tan(k)) * ((k / l) * (t * (k / l)))) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 1.42e-68) tmp = Float64(2.0 / Float64(2.0 * Float64(k / Float64(Float64(l / Float64((t ^ 3.0) / l)) / k)))); else tmp = Float64(2.0 / Float64(Float64(sin(k) * tan(k)) * Float64(Float64(k / l) * Float64(t * Float64(k / l))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 1.42e-68) tmp = 2.0 / (2.0 * (k / ((l / ((t ^ 3.0) / l)) / k))); else tmp = 2.0 / ((sin(k) * tan(k)) * ((k / l) * (t * (k / l)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 1.42e-68], N[(2.0 / N[(2.0 * N[(k / N[(N[(l / N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(k / l), $MachinePrecision] * N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.42 \cdot 10^{-68}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{k}{\frac{\frac{\ell}{\frac{{t}^{3}}{\ell}}}{k}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\sin k \cdot \tan k\right) \cdot \left(\frac{k}{\ell} \cdot \left(t \cdot \frac{k}{\ell}\right)\right)}\\
\end{array}
\end{array}
if k < 1.42e-68Initial program 53.5%
*-commutative53.5%
associate-*l*46.8%
associate-*r*46.8%
+-commutative46.8%
associate-+r+46.8%
metadata-eval46.8%
Simplified46.8%
Taylor expanded in k around 0 47.7%
associate-/l*47.2%
unpow247.2%
unpow247.2%
associate-*r/56.2%
associate-/l*63.8%
associate-*r/54.8%
associate-/l*63.8%
Simplified63.8%
if 1.42e-68 < k Initial program 55.9%
*-commutative55.9%
associate-*l*55.9%
associate-*r*55.9%
+-commutative55.9%
associate-+r+55.9%
metadata-eval55.9%
Simplified55.9%
Taylor expanded in k around inf 75.1%
associate-/l*72.5%
associate-/r/72.6%
unpow272.6%
unpow272.6%
Simplified72.6%
pow172.6%
associate-*l*72.6%
times-frac82.0%
Applied egg-rr82.0%
unpow182.0%
associate-*r*81.9%
*-commutative81.9%
associate-*r*89.7%
*-commutative89.7%
Simplified89.7%
Final simplification71.2%
(FPCore (t l k) :precision binary64 (if (<= k 9.6e-69) (/ 2.0 (* 2.0 (/ k (/ (/ l (/ (pow t 3.0) l)) k)))) (/ 2.0 (* (/ k l) (* (/ t (/ l k)) (* (sin k) (tan k)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 9.6e-69) {
tmp = 2.0 / (2.0 * (k / ((l / (pow(t, 3.0) / l)) / k)));
} else {
tmp = 2.0 / ((k / l) * ((t / (l / k)) * (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 (k <= 9.6d-69) then
tmp = 2.0d0 / (2.0d0 * (k / ((l / ((t ** 3.0d0) / l)) / k)))
else
tmp = 2.0d0 / ((k / l) * ((t / (l / k)) * (sin(k) * tan(k))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 9.6e-69) {
tmp = 2.0 / (2.0 * (k / ((l / (Math.pow(t, 3.0) / l)) / k)));
} else {
tmp = 2.0 / ((k / l) * ((t / (l / k)) * (Math.sin(k) * Math.tan(k))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 9.6e-69: tmp = 2.0 / (2.0 * (k / ((l / (math.pow(t, 3.0) / l)) / k))) else: tmp = 2.0 / ((k / l) * ((t / (l / k)) * (math.sin(k) * math.tan(k)))) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 9.6e-69) tmp = Float64(2.0 / Float64(2.0 * Float64(k / Float64(Float64(l / Float64((t ^ 3.0) / l)) / k)))); else tmp = Float64(2.0 / Float64(Float64(k / l) * Float64(Float64(t / Float64(l / k)) * Float64(sin(k) * tan(k))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 9.6e-69) tmp = 2.0 / (2.0 * (k / ((l / ((t ^ 3.0) / l)) / k))); else tmp = 2.0 / ((k / l) * ((t / (l / k)) * (sin(k) * tan(k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 9.6e-69], N[(2.0 / N[(2.0 * N[(k / N[(N[(l / N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k / l), $MachinePrecision] * N[(N[(t / N[(l / k), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 9.6 \cdot 10^{-69}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{k}{\frac{\frac{\ell}{\frac{{t}^{3}}{\ell}}}{k}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{k}{\ell} \cdot \left(\frac{t}{\frac{\ell}{k}} \cdot \left(\sin k \cdot \tan k\right)\right)}\\
\end{array}
\end{array}
if k < 9.6000000000000005e-69Initial program 53.5%
*-commutative53.5%
associate-*l*46.8%
associate-*r*46.8%
+-commutative46.8%
associate-+r+46.8%
metadata-eval46.8%
Simplified46.8%
Taylor expanded in k around 0 47.7%
associate-/l*47.2%
unpow247.2%
unpow247.2%
associate-*r/56.2%
associate-/l*63.8%
associate-*r/54.8%
associate-/l*63.8%
Simplified63.8%
if 9.6000000000000005e-69 < k Initial program 55.9%
*-commutative55.9%
associate-*l*55.9%
associate-*r*55.9%
+-commutative55.9%
associate-+r+55.9%
metadata-eval55.9%
Simplified55.9%
Taylor expanded in k around inf 75.1%
associate-/l*72.5%
associate-/r/72.6%
unpow272.6%
unpow272.6%
Simplified72.6%
pow172.6%
associate-*l*72.6%
times-frac82.0%
Applied egg-rr82.0%
unpow182.0%
associate-*r*81.9%
*-commutative81.9%
associate-*r*89.7%
*-commutative89.7%
Simplified89.7%
pow189.7%
Applied egg-rr89.7%
unpow189.7%
associate-*l*89.8%
associate-*r/88.5%
associate-/l*89.8%
Simplified89.8%
Final simplification71.2%
(FPCore (t l k)
:precision binary64
(if (<= t -1.5e-38)
(/ 2.0 (* 2.0 (/ k (/ (/ l (/ (pow t 3.0) l)) k))))
(if (<= t 4.6e-98)
(/ 2.0 (* (* (/ k l) (* t (/ k l))) (/ k (/ (cos k) k))))
(/ 2.0 (* 2.0 (* (pow t 3.0) (* (/ k l) (/ k l))))))))
double code(double t, double l, double k) {
double tmp;
if (t <= -1.5e-38) {
tmp = 2.0 / (2.0 * (k / ((l / (pow(t, 3.0) / l)) / k)));
} else if (t <= 4.6e-98) {
tmp = 2.0 / (((k / l) * (t * (k / l))) * (k / (cos(k) / k)));
} else {
tmp = 2.0 / (2.0 * (pow(t, 3.0) * ((k / l) * (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) :: tmp
if (t <= (-1.5d-38)) then
tmp = 2.0d0 / (2.0d0 * (k / ((l / ((t ** 3.0d0) / l)) / k)))
else if (t <= 4.6d-98) then
tmp = 2.0d0 / (((k / l) * (t * (k / l))) * (k / (cos(k) / k)))
else
tmp = 2.0d0 / (2.0d0 * ((t ** 3.0d0) * ((k / l) * (k / l))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (t <= -1.5e-38) {
tmp = 2.0 / (2.0 * (k / ((l / (Math.pow(t, 3.0) / l)) / k)));
} else if (t <= 4.6e-98) {
tmp = 2.0 / (((k / l) * (t * (k / l))) * (k / (Math.cos(k) / k)));
} else {
tmp = 2.0 / (2.0 * (Math.pow(t, 3.0) * ((k / l) * (k / l))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if t <= -1.5e-38: tmp = 2.0 / (2.0 * (k / ((l / (math.pow(t, 3.0) / l)) / k))) elif t <= 4.6e-98: tmp = 2.0 / (((k / l) * (t * (k / l))) * (k / (math.cos(k) / k))) else: tmp = 2.0 / (2.0 * (math.pow(t, 3.0) * ((k / l) * (k / l)))) return tmp
function code(t, l, k) tmp = 0.0 if (t <= -1.5e-38) tmp = Float64(2.0 / Float64(2.0 * Float64(k / Float64(Float64(l / Float64((t ^ 3.0) / l)) / k)))); elseif (t <= 4.6e-98) tmp = Float64(2.0 / Float64(Float64(Float64(k / l) * Float64(t * Float64(k / l))) * Float64(k / Float64(cos(k) / k)))); else tmp = Float64(2.0 / Float64(2.0 * Float64((t ^ 3.0) * Float64(Float64(k / l) * Float64(k / l))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (t <= -1.5e-38) tmp = 2.0 / (2.0 * (k / ((l / ((t ^ 3.0) / l)) / k))); elseif (t <= 4.6e-98) tmp = 2.0 / (((k / l) * (t * (k / l))) * (k / (cos(k) / k))); else tmp = 2.0 / (2.0 * ((t ^ 3.0) * ((k / l) * (k / l)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[t, -1.5e-38], N[(2.0 / N[(2.0 * N[(k / N[(N[(l / N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4.6e-98], N[(2.0 / N[(N[(N[(k / l), $MachinePrecision] * N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(k / N[(N[Cos[k], $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(2.0 * N[(N[Power[t, 3.0], $MachinePrecision] * N[(N[(k / l), $MachinePrecision] * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.5 \cdot 10^{-38}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{k}{\frac{\frac{\ell}{\frac{{t}^{3}}{\ell}}}{k}}}\\
\mathbf{elif}\;t \leq 4.6 \cdot 10^{-98}:\\
\;\;\;\;\frac{2}{\left(\frac{k}{\ell} \cdot \left(t \cdot \frac{k}{\ell}\right)\right) \cdot \frac{k}{\frac{\cos k}{k}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{2 \cdot \left({t}^{3} \cdot \left(\frac{k}{\ell} \cdot \frac{k}{\ell}\right)\right)}\\
\end{array}
\end{array}
if t < -1.49999999999999994e-38Initial program 69.0%
*-commutative69.0%
associate-*l*59.9%
associate-*r*59.9%
+-commutative59.9%
associate-+r+59.9%
metadata-eval59.9%
Simplified59.9%
Taylor expanded in k around 0 57.2%
associate-/l*57.1%
unpow257.1%
unpow257.1%
associate-*r/60.5%
associate-/l*69.7%
associate-*r/66.3%
associate-/l*69.6%
Simplified69.6%
if -1.49999999999999994e-38 < t < 4.60000000000000001e-98Initial program 26.5%
*-commutative26.5%
associate-*l*26.5%
associate-*r*26.5%
+-commutative26.5%
associate-+r+26.5%
metadata-eval26.5%
Simplified26.5%
Taylor expanded in k around inf 68.5%
times-frac69.1%
unpow269.1%
associate-/l*69.1%
associate-/l*69.1%
unpow269.1%
associate-/l*82.0%
Simplified82.0%
Taylor expanded in k around 0 60.6%
unpow260.6%
unpow260.6%
associate-*l*60.6%
times-frac74.0%
*-commutative74.0%
associate-*r/74.9%
Simplified74.9%
if 4.60000000000000001e-98 < t Initial program 73.1%
*-commutative73.1%
associate-*l*66.3%
associate-*r*66.3%
+-commutative66.3%
associate-+r+66.3%
metadata-eval66.3%
Simplified66.3%
Taylor expanded in l around 0 58.6%
associate-/l*57.4%
unpow257.4%
unpow257.4%
unpow257.4%
Simplified57.4%
Taylor expanded in k around 0 62.2%
associate-/l*61.2%
associate-/r/62.2%
unpow262.2%
unpow262.2%
times-frac73.2%
Simplified73.2%
Final simplification72.9%
(FPCore (t l k)
:precision binary64
(if (or (<= t -1.52e-37) (not (<= t 3e-100)))
(/ 2.0 (* 2.0 (* (pow t 3.0) (* (/ k l) (/ k l)))))
(*
2.0
(* (* (/ l (* k k)) (/ l t)) (+ -0.16666666666666666 (/ 1.0 (* k k)))))))
double code(double t, double l, double k) {
double tmp;
if ((t <= -1.52e-37) || !(t <= 3e-100)) {
tmp = 2.0 / (2.0 * (pow(t, 3.0) * ((k / l) * (k / l))));
} else {
tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (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 ((t <= (-1.52d-37)) .or. (.not. (t <= 3d-100))) then
tmp = 2.0d0 / (2.0d0 * ((t ** 3.0d0) * ((k / l) * (k / l))))
else
tmp = 2.0d0 * (((l / (k * k)) * (l / t)) * ((-0.16666666666666666d0) + (1.0d0 / (k * k))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((t <= -1.52e-37) || !(t <= 3e-100)) {
tmp = 2.0 / (2.0 * (Math.pow(t, 3.0) * ((k / l) * (k / l))));
} else {
tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (t <= -1.52e-37) or not (t <= 3e-100): tmp = 2.0 / (2.0 * (math.pow(t, 3.0) * ((k / l) * (k / l)))) else: tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))) return tmp
function code(t, l, k) tmp = 0.0 if ((t <= -1.52e-37) || !(t <= 3e-100)) tmp = Float64(2.0 / Float64(2.0 * Float64((t ^ 3.0) * Float64(Float64(k / l) * Float64(k / l))))); else tmp = Float64(2.0 * Float64(Float64(Float64(l / Float64(k * k)) * Float64(l / t)) * Float64(-0.16666666666666666 + Float64(1.0 / Float64(k * k))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((t <= -1.52e-37) || ~((t <= 3e-100))) tmp = 2.0 / (2.0 * ((t ^ 3.0) * ((k / l) * (k / l)))); else tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[Or[LessEqual[t, -1.52e-37], N[Not[LessEqual[t, 3e-100]], $MachinePrecision]], N[(2.0 / N[(2.0 * N[(N[Power[t, 3.0], $MachinePrecision] * N[(N[(k / l), $MachinePrecision] * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(-0.16666666666666666 + N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.52 \cdot 10^{-37} \lor \neg \left(t \leq 3 \cdot 10^{-100}\right):\\
\;\;\;\;\frac{2}{2 \cdot \left({t}^{3} \cdot \left(\frac{k}{\ell} \cdot \frac{k}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(\frac{\ell}{k \cdot k} \cdot \frac{\ell}{t}\right) \cdot \left(-0.16666666666666666 + \frac{1}{k \cdot k}\right)\right)\\
\end{array}
\end{array}
if t < -1.52e-37 or 3.0000000000000001e-100 < t Initial program 71.4%
*-commutative71.4%
associate-*l*63.6%
associate-*r*63.6%
+-commutative63.6%
associate-+r+63.6%
metadata-eval63.6%
Simplified63.6%
Taylor expanded in l around 0 56.7%
associate-/l*56.6%
unpow256.6%
unpow256.6%
unpow256.6%
Simplified56.6%
Taylor expanded in k around 0 60.1%
associate-/l*59.5%
associate-/r/58.8%
unpow258.8%
unpow258.8%
times-frac70.9%
Simplified70.9%
if -1.52e-37 < t < 3.0000000000000001e-100Initial program 26.5%
associate-/l/26.5%
associate-*l/25.8%
associate-*l/25.8%
associate-/r/25.8%
*-commutative25.8%
associate-/l/25.8%
associate-*r*25.8%
*-commutative25.8%
associate-*r*25.8%
*-commutative25.8%
Simplified25.8%
Taylor expanded in k around inf 68.2%
associate-/l*68.2%
unpow268.2%
associate-/r*68.2%
Simplified68.2%
Taylor expanded in k around 0 56.3%
Taylor expanded in l around 0 59.7%
unpow259.7%
associate-*l/59.8%
unpow259.8%
associate-/r*59.8%
associate-*r/59.8%
metadata-eval59.8%
div-sub59.8%
unpow259.8%
times-frac60.0%
*-commutative60.0%
associate-/l*56.5%
associate-/r/60.0%
unpow260.0%
times-frac74.1%
Simplified74.1%
Final simplification72.1%
(FPCore (t l k)
:precision binary64
(if (<= t -3.6e-39)
(/ 2.0 (* 2.0 (/ k (/ (/ l (/ (pow t 3.0) l)) k))))
(if (<= t 2.4e-95)
(*
2.0
(* (* (/ l (* k k)) (/ l t)) (+ -0.16666666666666666 (/ 1.0 (* k k)))))
(/ 2.0 (* 2.0 (* (pow t 3.0) (* (/ k l) (/ k l))))))))
double code(double t, double l, double k) {
double tmp;
if (t <= -3.6e-39) {
tmp = 2.0 / (2.0 * (k / ((l / (pow(t, 3.0) / l)) / k)));
} else if (t <= 2.4e-95) {
tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k))));
} else {
tmp = 2.0 / (2.0 * (pow(t, 3.0) * ((k / l) * (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) :: tmp
if (t <= (-3.6d-39)) then
tmp = 2.0d0 / (2.0d0 * (k / ((l / ((t ** 3.0d0) / l)) / k)))
else if (t <= 2.4d-95) then
tmp = 2.0d0 * (((l / (k * k)) * (l / t)) * ((-0.16666666666666666d0) + (1.0d0 / (k * k))))
else
tmp = 2.0d0 / (2.0d0 * ((t ** 3.0d0) * ((k / l) * (k / l))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (t <= -3.6e-39) {
tmp = 2.0 / (2.0 * (k / ((l / (Math.pow(t, 3.0) / l)) / k)));
} else if (t <= 2.4e-95) {
tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k))));
} else {
tmp = 2.0 / (2.0 * (Math.pow(t, 3.0) * ((k / l) * (k / l))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if t <= -3.6e-39: tmp = 2.0 / (2.0 * (k / ((l / (math.pow(t, 3.0) / l)) / k))) elif t <= 2.4e-95: tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))) else: tmp = 2.0 / (2.0 * (math.pow(t, 3.0) * ((k / l) * (k / l)))) return tmp
function code(t, l, k) tmp = 0.0 if (t <= -3.6e-39) tmp = Float64(2.0 / Float64(2.0 * Float64(k / Float64(Float64(l / Float64((t ^ 3.0) / l)) / k)))); elseif (t <= 2.4e-95) tmp = Float64(2.0 * Float64(Float64(Float64(l / Float64(k * k)) * Float64(l / t)) * Float64(-0.16666666666666666 + Float64(1.0 / Float64(k * k))))); else tmp = Float64(2.0 / Float64(2.0 * Float64((t ^ 3.0) * Float64(Float64(k / l) * Float64(k / l))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (t <= -3.6e-39) tmp = 2.0 / (2.0 * (k / ((l / ((t ^ 3.0) / l)) / k))); elseif (t <= 2.4e-95) tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))); else tmp = 2.0 / (2.0 * ((t ^ 3.0) * ((k / l) * (k / l)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[t, -3.6e-39], N[(2.0 / N[(2.0 * N[(k / N[(N[(l / N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.4e-95], N[(2.0 * N[(N[(N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(-0.16666666666666666 + N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(2.0 * N[(N[Power[t, 3.0], $MachinePrecision] * N[(N[(k / l), $MachinePrecision] * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -3.6 \cdot 10^{-39}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{k}{\frac{\frac{\ell}{\frac{{t}^{3}}{\ell}}}{k}}}\\
\mathbf{elif}\;t \leq 2.4 \cdot 10^{-95}:\\
\;\;\;\;2 \cdot \left(\left(\frac{\ell}{k \cdot k} \cdot \frac{\ell}{t}\right) \cdot \left(-0.16666666666666666 + \frac{1}{k \cdot k}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{2 \cdot \left({t}^{3} \cdot \left(\frac{k}{\ell} \cdot \frac{k}{\ell}\right)\right)}\\
\end{array}
\end{array}
if t < -3.6000000000000001e-39Initial program 69.0%
*-commutative69.0%
associate-*l*59.9%
associate-*r*59.9%
+-commutative59.9%
associate-+r+59.9%
metadata-eval59.9%
Simplified59.9%
Taylor expanded in k around 0 57.2%
associate-/l*57.1%
unpow257.1%
unpow257.1%
associate-*r/60.5%
associate-/l*69.7%
associate-*r/66.3%
associate-/l*69.6%
Simplified69.6%
if -3.6000000000000001e-39 < t < 2.4e-95Initial program 26.5%
associate-/l/26.5%
associate-*l/25.8%
associate-*l/25.8%
associate-/r/25.8%
*-commutative25.8%
associate-/l/25.8%
associate-*r*25.8%
*-commutative25.8%
associate-*r*25.8%
*-commutative25.8%
Simplified25.8%
Taylor expanded in k around inf 68.2%
associate-/l*68.2%
unpow268.2%
associate-/r*68.2%
Simplified68.2%
Taylor expanded in k around 0 56.3%
Taylor expanded in l around 0 59.7%
unpow259.7%
associate-*l/59.8%
unpow259.8%
associate-/r*59.8%
associate-*r/59.8%
metadata-eval59.8%
div-sub59.8%
unpow259.8%
times-frac60.0%
*-commutative60.0%
associate-/l*56.5%
associate-/r/60.0%
unpow260.0%
times-frac74.1%
Simplified74.1%
if 2.4e-95 < t Initial program 73.1%
*-commutative73.1%
associate-*l*66.3%
associate-*r*66.3%
+-commutative66.3%
associate-+r+66.3%
metadata-eval66.3%
Simplified66.3%
Taylor expanded in l around 0 58.6%
associate-/l*57.4%
unpow257.4%
unpow257.4%
unpow257.4%
Simplified57.4%
Taylor expanded in k around 0 62.2%
associate-/l*61.2%
associate-/r/62.2%
unpow262.2%
unpow262.2%
times-frac73.2%
Simplified73.2%
Final simplification72.6%
(FPCore (t l k)
:precision binary64
(if (or (<= t -4.6e-39) (not (<= t 2.6e-97)))
(* (/ l (pow t 3.0)) (/ (/ l k) k))
(*
2.0
(* (* (/ l (* k k)) (/ l t)) (+ -0.16666666666666666 (/ 1.0 (* k k)))))))
double code(double t, double l, double k) {
double tmp;
if ((t <= -4.6e-39) || !(t <= 2.6e-97)) {
tmp = (l / pow(t, 3.0)) * ((l / k) / k);
} else {
tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (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 ((t <= (-4.6d-39)) .or. (.not. (t <= 2.6d-97))) then
tmp = (l / (t ** 3.0d0)) * ((l / k) / k)
else
tmp = 2.0d0 * (((l / (k * k)) * (l / t)) * ((-0.16666666666666666d0) + (1.0d0 / (k * k))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((t <= -4.6e-39) || !(t <= 2.6e-97)) {
tmp = (l / Math.pow(t, 3.0)) * ((l / k) / k);
} else {
tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (t <= -4.6e-39) or not (t <= 2.6e-97): tmp = (l / math.pow(t, 3.0)) * ((l / k) / k) else: tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))) return tmp
function code(t, l, k) tmp = 0.0 if ((t <= -4.6e-39) || !(t <= 2.6e-97)) tmp = Float64(Float64(l / (t ^ 3.0)) * Float64(Float64(l / k) / k)); else tmp = Float64(2.0 * Float64(Float64(Float64(l / Float64(k * k)) * Float64(l / t)) * Float64(-0.16666666666666666 + Float64(1.0 / Float64(k * k))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((t <= -4.6e-39) || ~((t <= 2.6e-97))) tmp = (l / (t ^ 3.0)) * ((l / k) / k); else tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[Or[LessEqual[t, -4.6e-39], N[Not[LessEqual[t, 2.6e-97]], $MachinePrecision]], N[(N[(l / N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision] * N[(N[(l / k), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(-0.16666666666666666 + N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -4.6 \cdot 10^{-39} \lor \neg \left(t \leq 2.6 \cdot 10^{-97}\right):\\
\;\;\;\;\frac{\ell}{{t}^{3}} \cdot \frac{\frac{\ell}{k}}{k}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(\frac{\ell}{k \cdot k} \cdot \frac{\ell}{t}\right) \cdot \left(-0.16666666666666666 + \frac{1}{k \cdot k}\right)\right)\\
\end{array}
\end{array}
if t < -4.60000000000000016e-39 or 2.60000000000000007e-97 < t Initial program 71.4%
associate-*l*71.3%
associate-/l/71.3%
*-commutative71.3%
associate-*r/72.7%
associate-/l*72.3%
associate-/r/63.6%
Simplified67.1%
Taylor expanded in k around 0 60.1%
unpow260.1%
unpow260.1%
Simplified60.1%
Taylor expanded in l around 0 60.1%
unpow260.1%
times-frac62.3%
unpow262.3%
associate-/r*69.1%
Simplified69.1%
if -4.60000000000000016e-39 < t < 2.60000000000000007e-97Initial program 26.5%
associate-/l/26.5%
associate-*l/25.8%
associate-*l/25.8%
associate-/r/25.8%
*-commutative25.8%
associate-/l/25.8%
associate-*r*25.8%
*-commutative25.8%
associate-*r*25.8%
*-commutative25.8%
Simplified25.8%
Taylor expanded in k around inf 68.2%
associate-/l*68.2%
unpow268.2%
associate-/r*68.2%
Simplified68.2%
Taylor expanded in k around 0 56.3%
Taylor expanded in l around 0 59.7%
unpow259.7%
associate-*l/59.8%
unpow259.8%
associate-/r*59.8%
associate-*r/59.8%
metadata-eval59.8%
div-sub59.8%
unpow259.8%
times-frac60.0%
*-commutative60.0%
associate-/l*56.5%
associate-/r/60.0%
unpow260.0%
times-frac74.1%
Simplified74.1%
Final simplification71.0%
(FPCore (t l k)
:precision binary64
(if (<= t -1.25e-5)
(/ (* l l) (* k (* (pow t 3.0) k)))
(if (<= t 9e-95)
(*
2.0
(* (* (/ l (* k k)) (/ l t)) (+ -0.16666666666666666 (/ 1.0 (* k k)))))
(* (/ l (pow t 3.0)) (/ (/ l k) k)))))
double code(double t, double l, double k) {
double tmp;
if (t <= -1.25e-5) {
tmp = (l * l) / (k * (pow(t, 3.0) * k));
} else if (t <= 9e-95) {
tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k))));
} else {
tmp = (l / pow(t, 3.0)) * ((l / 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 (t <= (-1.25d-5)) then
tmp = (l * l) / (k * ((t ** 3.0d0) * k))
else if (t <= 9d-95) then
tmp = 2.0d0 * (((l / (k * k)) * (l / t)) * ((-0.16666666666666666d0) + (1.0d0 / (k * k))))
else
tmp = (l / (t ** 3.0d0)) * ((l / k) / k)
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (t <= -1.25e-5) {
tmp = (l * l) / (k * (Math.pow(t, 3.0) * k));
} else if (t <= 9e-95) {
tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k))));
} else {
tmp = (l / Math.pow(t, 3.0)) * ((l / k) / k);
}
return tmp;
}
def code(t, l, k): tmp = 0 if t <= -1.25e-5: tmp = (l * l) / (k * (math.pow(t, 3.0) * k)) elif t <= 9e-95: tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))) else: tmp = (l / math.pow(t, 3.0)) * ((l / k) / k) return tmp
function code(t, l, k) tmp = 0.0 if (t <= -1.25e-5) tmp = Float64(Float64(l * l) / Float64(k * Float64((t ^ 3.0) * k))); elseif (t <= 9e-95) tmp = Float64(2.0 * Float64(Float64(Float64(l / Float64(k * k)) * Float64(l / t)) * Float64(-0.16666666666666666 + Float64(1.0 / Float64(k * k))))); else tmp = Float64(Float64(l / (t ^ 3.0)) * Float64(Float64(l / k) / k)); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (t <= -1.25e-5) tmp = (l * l) / (k * ((t ^ 3.0) * k)); elseif (t <= 9e-95) tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))); else tmp = (l / (t ^ 3.0)) * ((l / k) / k); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[t, -1.25e-5], N[(N[(l * l), $MachinePrecision] / N[(k * N[(N[Power[t, 3.0], $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 9e-95], N[(2.0 * N[(N[(N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(-0.16666666666666666 + N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision] * N[(N[(l / k), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.25 \cdot 10^{-5}:\\
\;\;\;\;\frac{\ell \cdot \ell}{k \cdot \left({t}^{3} \cdot k\right)}\\
\mathbf{elif}\;t \leq 9 \cdot 10^{-95}:\\
\;\;\;\;2 \cdot \left(\left(\frac{\ell}{k \cdot k} \cdot \frac{\ell}{t}\right) \cdot \left(-0.16666666666666666 + \frac{1}{k \cdot k}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{{t}^{3}} \cdot \frac{\frac{\ell}{k}}{k}\\
\end{array}
\end{array}
if t < -1.25000000000000006e-5Initial program 67.6%
associate-*l*67.6%
associate-/l/67.6%
*-commutative67.6%
associate-*r/67.6%
associate-/l*67.6%
associate-/r/59.1%
Simplified62.8%
Taylor expanded in k around 0 57.5%
unpow257.5%
unpow257.5%
Simplified57.5%
Taylor expanded in l around 0 57.5%
unpow257.5%
unpow257.5%
associate-*r*65.9%
Simplified65.9%
if -1.25000000000000006e-5 < t < 9e-95Initial program 29.8%
associate-/l/29.7%
associate-*l/29.1%
associate-*l/28.2%
associate-/r/28.2%
*-commutative28.2%
associate-/l/28.2%
associate-*r*28.2%
*-commutative28.2%
associate-*r*28.2%
*-commutative28.2%
Simplified28.2%
Taylor expanded in k around inf 68.6%
associate-/l*68.6%
unpow268.6%
associate-/r*68.6%
Simplified68.6%
Taylor expanded in k around 0 56.4%
Taylor expanded in l around 0 59.7%
unpow259.7%
associate-*l/59.7%
unpow259.7%
associate-/r*59.7%
associate-*r/59.7%
metadata-eval59.7%
div-sub59.7%
unpow259.7%
times-frac59.9%
*-commutative59.9%
associate-/l*56.6%
associate-/r/59.9%
unpow259.9%
times-frac73.2%
Simplified73.2%
if 9e-95 < t Initial program 73.1%
associate-*l*73.1%
associate-/l/73.0%
*-commutative73.0%
associate-*r/75.5%
associate-/l*74.7%
associate-/r/66.4%
Simplified69.9%
Taylor expanded in k around 0 62.2%
unpow262.2%
unpow262.2%
Simplified62.2%
Taylor expanded in l around 0 62.2%
unpow262.2%
times-frac64.8%
unpow264.8%
associate-/r*72.3%
Simplified72.3%
Final simplification71.1%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ l (* k k))))
(if (<= t -1e-9)
(* (/ l (pow t 3.0)) t_1)
(* 2.0 (* (* t_1 (/ l t)) (+ -0.16666666666666666 (/ 1.0 (* k k))))))))
double code(double t, double l, double k) {
double t_1 = l / (k * k);
double tmp;
if (t <= -1e-9) {
tmp = (l / pow(t, 3.0)) * t_1;
} else {
tmp = 2.0 * ((t_1 * (l / t)) * (-0.16666666666666666 + (1.0 / (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) :: t_1
real(8) :: tmp
t_1 = l / (k * k)
if (t <= (-1d-9)) then
tmp = (l / (t ** 3.0d0)) * t_1
else
tmp = 2.0d0 * ((t_1 * (l / t)) * ((-0.16666666666666666d0) + (1.0d0 / (k * k))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = l / (k * k);
double tmp;
if (t <= -1e-9) {
tmp = (l / Math.pow(t, 3.0)) * t_1;
} else {
tmp = 2.0 * ((t_1 * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k))));
}
return tmp;
}
def code(t, l, k): t_1 = l / (k * k) tmp = 0 if t <= -1e-9: tmp = (l / math.pow(t, 3.0)) * t_1 else: tmp = 2.0 * ((t_1 * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))) return tmp
function code(t, l, k) t_1 = Float64(l / Float64(k * k)) tmp = 0.0 if (t <= -1e-9) tmp = Float64(Float64(l / (t ^ 3.0)) * t_1); else tmp = Float64(2.0 * Float64(Float64(t_1 * Float64(l / t)) * Float64(-0.16666666666666666 + Float64(1.0 / Float64(k * k))))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = l / (k * k); tmp = 0.0; if (t <= -1e-9) tmp = (l / (t ^ 3.0)) * t_1; else tmp = 2.0 * ((t_1 * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1e-9], N[(N[(l / N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision], N[(2.0 * N[(N[(t$95$1 * N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(-0.16666666666666666 + N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\ell}{k \cdot k}\\
\mathbf{if}\;t \leq -1 \cdot 10^{-9}:\\
\;\;\;\;\frac{\ell}{{t}^{3}} \cdot t_1\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(t_1 \cdot \frac{\ell}{t}\right) \cdot \left(-0.16666666666666666 + \frac{1}{k \cdot k}\right)\right)\\
\end{array}
\end{array}
if t < -1.00000000000000006e-9Initial program 67.6%
associate-*l*67.6%
associate-/l/67.6%
*-commutative67.6%
associate-*r/67.6%
associate-/l*67.6%
associate-/r/59.1%
Simplified62.8%
Taylor expanded in k around 0 57.5%
unpow257.5%
unpow257.5%
Simplified57.5%
times-frac59.5%
Applied egg-rr59.5%
if -1.00000000000000006e-9 < t Initial program 50.0%
associate-/l/50.0%
associate-*l/50.7%
associate-*l/49.1%
associate-/r/49.2%
*-commutative49.2%
associate-/l/49.2%
associate-*r*49.2%
*-commutative49.2%
associate-*r*49.2%
*-commutative49.2%
Simplified49.2%
Taylor expanded in k around inf 65.4%
associate-/l*65.4%
unpow265.4%
associate-/r*65.4%
Simplified65.4%
Taylor expanded in k around 0 57.4%
Taylor expanded in l around 0 60.2%
unpow260.2%
associate-*l/59.8%
unpow259.8%
associate-/r*59.8%
associate-*r/59.8%
metadata-eval59.8%
div-sub59.8%
unpow259.8%
times-frac60.3%
*-commutative60.3%
associate-/l*58.0%
associate-/r/60.8%
unpow260.8%
times-frac68.1%
Simplified68.1%
Final simplification66.0%
(FPCore (t l k) :precision binary64 (if (<= (* l l) 1e-240) (* (* l l) (/ 2.0 (* t (+ (* (* k k) -6.0) -36.0)))) (* (* l l) (* (/ 2.0 (* k k)) (/ (/ (/ 1.0 t) k) k)))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-240) {
tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0)));
} else {
tmp = (l * l) * ((2.0 / (k * k)) * (((1.0 / 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-240) then
tmp = (l * l) * (2.0d0 / (t * (((k * k) * (-6.0d0)) + (-36.0d0))))
else
tmp = (l * l) * ((2.0d0 / (k * k)) * (((1.0d0 / 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-240) {
tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0)));
} else {
tmp = (l * l) * ((2.0 / (k * k)) * (((1.0 / t) / k) / k));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (l * l) <= 1e-240: tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0))) else: tmp = (l * l) * ((2.0 / (k * k)) * (((1.0 / t) / k) / k)) return tmp
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-240) tmp = Float64(Float64(l * l) * Float64(2.0 / Float64(t * Float64(Float64(Float64(k * k) * -6.0) + -36.0)))); else tmp = Float64(Float64(l * l) * Float64(Float64(2.0 / Float64(k * k)) * Float64(Float64(Float64(1.0 / t) / k) / k))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((l * l) <= 1e-240) tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0))); else tmp = (l * l) * ((2.0 / (k * k)) * (((1.0 / t) / k) / k)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e-240], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(t * N[(N[(N[(k * k), $MachinePrecision] * -6.0), $MachinePrecision] + -36.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(1.0 / t), $MachinePrecision] / k), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-240}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{t \cdot \left(\left(k \cdot k\right) \cdot -6 + -36\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \left(\frac{2}{k \cdot k} \cdot \frac{\frac{\frac{1}{t}}{k}}{k}\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 9.9999999999999997e-241Initial program 54.1%
associate-/l/54.0%
associate-*l/54.6%
associate-*l/53.2%
associate-/r/53.2%
*-commutative53.2%
associate-/l/53.2%
associate-*r*53.2%
*-commutative53.2%
associate-*r*53.2%
*-commutative53.2%
Simplified53.2%
Taylor expanded in k around inf 51.4%
associate-/l*51.4%
unpow251.4%
associate-/r*51.4%
Simplified51.4%
Taylor expanded in k around 0 48.4%
Taylor expanded in k around inf 64.3%
+-commutative64.3%
unpow264.3%
associate-*r*64.3%
distribute-rgt-out64.3%
Simplified64.3%
if 9.9999999999999997e-241 < (*.f64 l l) Initial program 54.2%
associate-/l/54.2%
associate-*l/54.8%
associate-*l/53.7%
associate-/r/53.2%
*-commutative53.2%
associate-/l/53.2%
associate-*r*53.2%
*-commutative53.2%
associate-*r*53.2%
*-commutative53.2%
Simplified53.2%
Taylor expanded in k around inf 63.2%
associate-/l*63.2%
unpow263.2%
associate-/r*63.2%
Simplified63.2%
Taylor expanded in k around 0 52.2%
*-commutative52.2%
associate-/r*52.2%
unpow252.2%
Simplified52.2%
associate-/r/52.2%
associate-/r*52.2%
Applied egg-rr52.2%
Final simplification56.4%
(FPCore (t l k) :precision binary64 (if (<= (* l l) 1e-240) (* (* l l) (/ 2.0 (* t (+ (* (* k k) -6.0) -36.0)))) (* (* l l) (/ 2.0 (/ (* k k) (/ (/ 1.0 t) (* k k)))))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-240) {
tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0)));
} else {
tmp = (l * l) * (2.0 / ((k * k) / ((1.0 / 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-240) then
tmp = (l * l) * (2.0d0 / (t * (((k * k) * (-6.0d0)) + (-36.0d0))))
else
tmp = (l * l) * (2.0d0 / ((k * k) / ((1.0d0 / 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-240) {
tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0)));
} else {
tmp = (l * l) * (2.0 / ((k * k) / ((1.0 / t) / (k * k))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (l * l) <= 1e-240: tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0))) else: tmp = (l * l) * (2.0 / ((k * k) / ((1.0 / t) / (k * k)))) return tmp
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-240) tmp = Float64(Float64(l * l) * Float64(2.0 / Float64(t * Float64(Float64(Float64(k * k) * -6.0) + -36.0)))); else tmp = Float64(Float64(l * l) * Float64(2.0 / Float64(Float64(k * k) / Float64(Float64(1.0 / t) / Float64(k * k))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((l * l) <= 1e-240) tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0))); else tmp = (l * l) * (2.0 / ((k * k) / ((1.0 / t) / (k * k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e-240], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(t * N[(N[(N[(k * k), $MachinePrecision] * -6.0), $MachinePrecision] + -36.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[(k * k), $MachinePrecision] / N[(N[(1.0 / t), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-240}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{t \cdot \left(\left(k \cdot k\right) \cdot -6 + -36\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{\frac{k \cdot k}{\frac{\frac{1}{t}}{k \cdot k}}}\\
\end{array}
\end{array}
if (*.f64 l l) < 9.9999999999999997e-241Initial program 54.1%
associate-/l/54.0%
associate-*l/54.6%
associate-*l/53.2%
associate-/r/53.2%
*-commutative53.2%
associate-/l/53.2%
associate-*r*53.2%
*-commutative53.2%
associate-*r*53.2%
*-commutative53.2%
Simplified53.2%
Taylor expanded in k around inf 51.4%
associate-/l*51.4%
unpow251.4%
associate-/r*51.4%
Simplified51.4%
Taylor expanded in k around 0 48.4%
Taylor expanded in k around inf 64.3%
+-commutative64.3%
unpow264.3%
associate-*r*64.3%
distribute-rgt-out64.3%
Simplified64.3%
if 9.9999999999999997e-241 < (*.f64 l l) Initial program 54.2%
associate-/l/54.2%
associate-*l/54.8%
associate-*l/53.7%
associate-/r/53.2%
*-commutative53.2%
associate-/l/53.2%
associate-*r*53.2%
*-commutative53.2%
associate-*r*53.2%
*-commutative53.2%
Simplified53.2%
Taylor expanded in k around inf 63.2%
associate-/l*63.2%
unpow263.2%
associate-/r*63.2%
Simplified63.2%
Taylor expanded in k around 0 52.2%
*-commutative52.2%
associate-/r*52.2%
unpow252.2%
Simplified52.2%
Final simplification56.4%
(FPCore (t l k) :precision binary64 (if (<= (* l l) 1e-240) (* (* l l) (/ 2.0 (* t (+ (* (* k k) -6.0) -36.0)))) (/ (* 2.0 (* l l)) (* (* k k) (/ (* k k) (/ 1.0 t))))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-240) {
tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0)));
} else {
tmp = (2.0 * (l * l)) / ((k * k) * ((k * k) / (1.0 / 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 ((l * l) <= 1d-240) then
tmp = (l * l) * (2.0d0 / (t * (((k * k) * (-6.0d0)) + (-36.0d0))))
else
tmp = (2.0d0 * (l * l)) / ((k * k) * ((k * k) / (1.0d0 / t)))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-240) {
tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0)));
} else {
tmp = (2.0 * (l * l)) / ((k * k) * ((k * k) / (1.0 / t)));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (l * l) <= 1e-240: tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0))) else: tmp = (2.0 * (l * l)) / ((k * k) * ((k * k) / (1.0 / t))) return tmp
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-240) tmp = Float64(Float64(l * l) * Float64(2.0 / Float64(t * Float64(Float64(Float64(k * k) * -6.0) + -36.0)))); else tmp = Float64(Float64(2.0 * Float64(l * l)) / Float64(Float64(k * k) * Float64(Float64(k * k) / Float64(1.0 / t)))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((l * l) <= 1e-240) tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0))); else tmp = (2.0 * (l * l)) / ((k * k) * ((k * k) / (1.0 / t))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e-240], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(t * N[(N[(N[(k * k), $MachinePrecision] * -6.0), $MachinePrecision] + -36.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(l * l), $MachinePrecision]), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] / N[(1.0 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-240}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{t \cdot \left(\left(k \cdot k\right) \cdot -6 + -36\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \left(\ell \cdot \ell\right)}{\left(k \cdot k\right) \cdot \frac{k \cdot k}{\frac{1}{t}}}\\
\end{array}
\end{array}
if (*.f64 l l) < 9.9999999999999997e-241Initial program 54.1%
associate-/l/54.0%
associate-*l/54.6%
associate-*l/53.2%
associate-/r/53.2%
*-commutative53.2%
associate-/l/53.2%
associate-*r*53.2%
*-commutative53.2%
associate-*r*53.2%
*-commutative53.2%
Simplified53.2%
Taylor expanded in k around inf 51.4%
associate-/l*51.4%
unpow251.4%
associate-/r*51.4%
Simplified51.4%
Taylor expanded in k around 0 48.4%
Taylor expanded in k around inf 64.3%
+-commutative64.3%
unpow264.3%
associate-*r*64.3%
distribute-rgt-out64.3%
Simplified64.3%
if 9.9999999999999997e-241 < (*.f64 l l) Initial program 54.2%
associate-/l/54.2%
associate-*l/54.8%
associate-*l/53.7%
associate-/r/53.2%
*-commutative53.2%
associate-/l/53.2%
associate-*r*53.2%
*-commutative53.2%
associate-*r*53.2%
*-commutative53.2%
Simplified53.2%
Taylor expanded in k around inf 63.2%
associate-/l*63.2%
unpow263.2%
associate-/r*63.2%
Simplified63.2%
Taylor expanded in k around 0 52.2%
*-commutative52.2%
associate-/r*52.2%
unpow252.2%
Simplified52.2%
associate-*r/52.4%
associate-/r/52.4%
Applied egg-rr52.4%
Final simplification56.5%
(FPCore (t l k)
:precision binary64
(if (<= l 1.7e-146)
(* (* l l) (/ 2.0 (* t (+ (* (* k k) -6.0) -36.0))))
(*
-2.0
(* (/ (* l l) (* k k)) (/ (+ 0.16666666666666666 (/ -1.0 (* k k))) t)))))
double code(double t, double l, double k) {
double tmp;
if (l <= 1.7e-146) {
tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0)));
} else {
tmp = -2.0 * (((l * l) / (k * k)) * ((0.16666666666666666 + (-1.0 / (k * 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 (l <= 1.7d-146) then
tmp = (l * l) * (2.0d0 / (t * (((k * k) * (-6.0d0)) + (-36.0d0))))
else
tmp = (-2.0d0) * (((l * l) / (k * k)) * ((0.16666666666666666d0 + ((-1.0d0) / (k * k))) / t))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (l <= 1.7e-146) {
tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0)));
} else {
tmp = -2.0 * (((l * l) / (k * k)) * ((0.16666666666666666 + (-1.0 / (k * k))) / t));
}
return tmp;
}
def code(t, l, k): tmp = 0 if l <= 1.7e-146: tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0))) else: tmp = -2.0 * (((l * l) / (k * k)) * ((0.16666666666666666 + (-1.0 / (k * k))) / t)) return tmp
function code(t, l, k) tmp = 0.0 if (l <= 1.7e-146) tmp = Float64(Float64(l * l) * Float64(2.0 / Float64(t * Float64(Float64(Float64(k * k) * -6.0) + -36.0)))); else tmp = Float64(-2.0 * Float64(Float64(Float64(l * l) / Float64(k * k)) * Float64(Float64(0.16666666666666666 + Float64(-1.0 / Float64(k * k))) / t))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (l <= 1.7e-146) tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0))); else tmp = -2.0 * (((l * l) / (k * k)) * ((0.16666666666666666 + (-1.0 / (k * k))) / t)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[l, 1.7e-146], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(t * N[(N[(N[(k * k), $MachinePrecision] * -6.0), $MachinePrecision] + -36.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-2.0 * N[(N[(N[(l * l), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(0.16666666666666666 + N[(-1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \leq 1.7 \cdot 10^{-146}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{t \cdot \left(\left(k \cdot k\right) \cdot -6 + -36\right)}\\
\mathbf{else}:\\
\;\;\;\;-2 \cdot \left(\frac{\ell \cdot \ell}{k \cdot k} \cdot \frac{0.16666666666666666 + \frac{-1}{k \cdot k}}{t}\right)\\
\end{array}
\end{array}
if l < 1.7e-146Initial program 56.7%
associate-/l/56.7%
associate-*l/57.2%
associate-*l/55.6%
associate-/r/55.6%
*-commutative55.6%
associate-/l/55.6%
associate-*r*55.6%
*-commutative55.6%
associate-*r*55.6%
*-commutative55.6%
Simplified55.6%
Taylor expanded in k around inf 58.2%
associate-/l*58.2%
unpow258.2%
associate-/r*58.2%
Simplified58.2%
Taylor expanded in k around 0 54.4%
Taylor expanded in k around inf 46.1%
+-commutative46.1%
unpow246.1%
associate-*r*46.1%
distribute-rgt-out46.1%
Simplified46.1%
if 1.7e-146 < l Initial program 48.8%
associate-/l/48.8%
associate-*l/49.5%
associate-*l/49.2%
associate-/r/48.0%
*-commutative48.0%
associate-/l/48.0%
associate-*r*48.0%
*-commutative48.0%
associate-*r*48.0%
*-commutative48.0%
Simplified48.0%
Taylor expanded in k around inf 61.1%
associate-/l*61.1%
unpow261.1%
associate-/r*61.1%
Simplified61.1%
Taylor expanded in k around 0 50.0%
Taylor expanded in t around -inf 53.4%
unpow253.4%
times-frac52.1%
unpow252.1%
unpow252.1%
Simplified52.1%
Final simplification48.0%
(FPCore (t l k) :precision binary64 (* 2.0 (* (* (/ l (* k k)) (/ l t)) (+ -0.16666666666666666 (/ 1.0 (* k k))))))
double code(double t, double l, double k) {
return 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (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 = 2.0d0 * (((l / (k * k)) * (l / t)) * ((-0.16666666666666666d0) + (1.0d0 / (k * k))))
end function
public static double code(double t, double l, double k) {
return 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k))));
}
def code(t, l, k): return 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k))))
function code(t, l, k) return Float64(2.0 * Float64(Float64(Float64(l / Float64(k * k)) * Float64(l / t)) * Float64(-0.16666666666666666 + Float64(1.0 / Float64(k * k))))) end
function tmp = code(t, l, k) tmp = 2.0 * (((l / (k * k)) * (l / t)) * (-0.16666666666666666 + (1.0 / (k * k)))); end
code[t_, l_, k_] := N[(2.0 * N[(N[(N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(-0.16666666666666666 + N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \left(\left(\frac{\ell}{k \cdot k} \cdot \frac{\ell}{t}\right) \cdot \left(-0.16666666666666666 + \frac{1}{k \cdot k}\right)\right)
\end{array}
Initial program 54.2%
associate-/l/54.2%
associate-*l/54.7%
associate-*l/53.5%
associate-/r/53.2%
*-commutative53.2%
associate-/l/53.2%
associate-*r*53.2%
*-commutative53.2%
associate-*r*53.2%
*-commutative53.2%
Simplified53.2%
Taylor expanded in k around inf 59.1%
associate-/l*59.1%
unpow259.1%
associate-/r*59.1%
Simplified59.1%
Taylor expanded in k around 0 53.0%
Taylor expanded in l around 0 55.5%
unpow255.5%
associate-*l/55.2%
unpow255.2%
associate-/r*55.2%
associate-*r/55.2%
metadata-eval55.2%
div-sub55.2%
unpow255.2%
times-frac55.6%
*-commutative55.6%
associate-/l*53.5%
associate-/r/56.0%
unpow256.0%
times-frac63.2%
Simplified63.2%
Final simplification63.2%
(FPCore (t l k) :precision binary64 (* (* l l) (/ 2.0 (* t (+ (* (* k k) -6.0) -36.0)))))
double code(double t, double l, double k) {
return (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.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 * l) * (2.0d0 / (t * (((k * k) * (-6.0d0)) + (-36.0d0))))
end function
public static double code(double t, double l, double k) {
return (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0)));
}
def code(t, l, k): return (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0)))
function code(t, l, k) return Float64(Float64(l * l) * Float64(2.0 / Float64(t * Float64(Float64(Float64(k * k) * -6.0) + -36.0)))) end
function tmp = code(t, l, k) tmp = (l * l) * (2.0 / (t * (((k * k) * -6.0) + -36.0))); end
code[t_, l_, k_] := N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(t * N[(N[(N[(k * k), $MachinePrecision] * -6.0), $MachinePrecision] + -36.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\ell \cdot \ell\right) \cdot \frac{2}{t \cdot \left(\left(k \cdot k\right) \cdot -6 + -36\right)}
\end{array}
Initial program 54.2%
associate-/l/54.2%
associate-*l/54.7%
associate-*l/53.5%
associate-/r/53.2%
*-commutative53.2%
associate-/l/53.2%
associate-*r*53.2%
*-commutative53.2%
associate-*r*53.2%
*-commutative53.2%
Simplified53.2%
Taylor expanded in k around inf 59.1%
associate-/l*59.1%
unpow259.1%
associate-/r*59.1%
Simplified59.1%
Taylor expanded in k around 0 53.0%
Taylor expanded in k around inf 40.3%
+-commutative40.3%
unpow240.3%
associate-*r*40.3%
distribute-rgt-out40.3%
Simplified40.3%
Final simplification40.3%
(FPCore (t l k) :precision binary64 (* (* (/ l (* k k)) (/ l t)) -0.3333333333333333))
double code(double t, double l, double k) {
return ((l / (k * k)) * (l / t)) * -0.3333333333333333;
}
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 * k)) * (l / t)) * (-0.3333333333333333d0)
end function
public static double code(double t, double l, double k) {
return ((l / (k * k)) * (l / t)) * -0.3333333333333333;
}
def code(t, l, k): return ((l / (k * k)) * (l / t)) * -0.3333333333333333
function code(t, l, k) return Float64(Float64(Float64(l / Float64(k * k)) * Float64(l / t)) * -0.3333333333333333) end
function tmp = code(t, l, k) tmp = ((l / (k * k)) * (l / t)) * -0.3333333333333333; end
code[t_, l_, k_] := N[(N[(N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision] * -0.3333333333333333), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{\ell}{k \cdot k} \cdot \frac{\ell}{t}\right) \cdot -0.3333333333333333
\end{array}
Initial program 54.2%
associate-/l/54.2%
associate-*l/54.7%
associate-*l/53.5%
associate-/r/53.2%
*-commutative53.2%
associate-/l/53.2%
associate-*r*53.2%
*-commutative53.2%
associate-*r*53.2%
*-commutative53.2%
Simplified53.2%
Taylor expanded in k around inf 59.1%
associate-/l*59.1%
unpow259.1%
associate-/r*59.1%
Simplified59.1%
Taylor expanded in k around 0 53.0%
Taylor expanded in k around inf 35.3%
unpow235.3%
unpow235.3%
*-commutative35.3%
times-frac36.9%
Simplified36.9%
Final simplification36.9%
(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(Float64(l * Float64(l / 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[(N[(l * N[(l / t), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.3333333333333333 \cdot \frac{\ell \cdot \frac{\ell}{t}}{k \cdot k}
\end{array}
Initial program 54.2%
associate-/l/54.2%
associate-*l/54.7%
associate-*l/53.5%
associate-/r/53.2%
*-commutative53.2%
associate-/l/53.2%
associate-*r*53.2%
*-commutative53.2%
associate-*r*53.2%
*-commutative53.2%
Simplified53.2%
Taylor expanded in k around inf 59.1%
associate-/l*59.1%
unpow259.1%
associate-/r*59.1%
Simplified59.1%
Taylor expanded in k around 0 53.0%
Taylor expanded in k around inf 35.3%
unpow235.3%
unpow235.3%
*-commutative35.3%
times-frac36.9%
Simplified36.9%
associate-*r/37.4%
*-commutative37.4%
Applied egg-rr37.4%
Final simplification37.4%
herbie shell --seed 2023216
(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))))