
(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 13 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
(if (<= (* l l) 4e+275)
(/ 2.0 (* (/ k l) (* (sin k) (* (sin k) (/ (* k t) (* l (cos k)))))))
(/
2.0
(* (/ k l) (* (/ (* k (fma (cos (+ k k)) -0.5 0.5)) l) (/ t (cos k)))))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 4e+275) {
tmp = 2.0 / ((k / l) * (sin(k) * (sin(k) * ((k * t) / (l * cos(k))))));
} else {
tmp = 2.0 / ((k / l) * (((k * fma(cos((k + k)), -0.5, 0.5)) / l) * (t / cos(k))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 4e+275) tmp = Float64(2.0 / Float64(Float64(k / l) * Float64(sin(k) * Float64(sin(k) * Float64(Float64(k * t) / Float64(l * cos(k))))))); else tmp = Float64(2.0 / Float64(Float64(k / l) * Float64(Float64(Float64(k * fma(cos(Float64(k + k)), -0.5, 0.5)) / l) * Float64(t / cos(k))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 4e+275], N[(2.0 / N[(N[(k / l), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(k * t), $MachinePrecision] / N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k / l), $MachinePrecision] * N[(N[(N[(k * N[(N[Cos[N[(k + k), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(t / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 4 \cdot 10^{+275}:\\
\;\;\;\;\frac{2}{\frac{k}{\ell} \cdot \left(\sin k \cdot \left(\sin k \cdot \frac{k \cdot t}{\ell \cdot \cos k}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{k}{\ell} \cdot \left(\frac{k \cdot \mathsf{fma}\left(\cos \left(k + k\right), -0.5, 0.5\right)}{\ell} \cdot \frac{t}{\cos k}\right)}\\
\end{array}
\end{array}
if (*.f64 l l) < 3.99999999999999984e275Initial program 39.9%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6485.0
Simplified85.0%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-*l*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
Applied egg-rr82.0%
lift-+.f64N/A
lift-cos.f64N/A
lift-cos.f64N/A
lift-+.f64N/A
count-2N/A
sqr-sin-aN/A
pow2N/A
lower-pow.f64N/A
lower-sin.f6494.8
Applied egg-rr94.8%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-/l*N/A
lift-pow.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6498.8
Applied egg-rr98.8%
if 3.99999999999999984e275 < (*.f64 l l) Initial program 36.8%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6467.9
Simplified67.9%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-*l*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
Applied egg-rr89.9%
lift-+.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift--.f64N/A
associate-*r*N/A
lift-cos.f64N/A
times-fracN/A
lower-*.f64N/A
Applied egg-rr99.3%
(FPCore (t l k)
:precision binary64
(if (<= k 1.15e-19)
(/ 2.0 (* (* (/ k l) (* k k)) (/ (* k t) l)))
(/
2.0
(* (/ k l) (* (/ (* k (fma (cos (+ k k)) -0.5 0.5)) l) (/ t (cos k)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 1.15e-19) {
tmp = 2.0 / (((k / l) * (k * k)) * ((k * t) / l));
} else {
tmp = 2.0 / ((k / l) * (((k * fma(cos((k + k)), -0.5, 0.5)) / l) * (t / cos(k))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (k <= 1.15e-19) tmp = Float64(2.0 / Float64(Float64(Float64(k / l) * Float64(k * k)) * Float64(Float64(k * t) / l))); else tmp = Float64(2.0 / Float64(Float64(k / l) * Float64(Float64(Float64(k * fma(cos(Float64(k + k)), -0.5, 0.5)) / l) * Float64(t / cos(k))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[k, 1.15e-19], N[(2.0 / N[(N[(N[(k / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k / l), $MachinePrecision] * N[(N[(N[(k * N[(N[Cos[N[(k + k), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(t / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.15 \cdot 10^{-19}:\\
\;\;\;\;\frac{2}{\left(\frac{k}{\ell} \cdot \left(k \cdot k\right)\right) \cdot \frac{k \cdot t}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{k}{\ell} \cdot \left(\frac{k \cdot \mathsf{fma}\left(\cos \left(k + k\right), -0.5, 0.5\right)}{\ell} \cdot \frac{t}{\cos k}\right)}\\
\end{array}
\end{array}
if k < 1.1499999999999999e-19Initial program 40.1%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6472.1
Simplified72.1%
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6481.9
Applied egg-rr81.9%
if 1.1499999999999999e-19 < k Initial program 36.0%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6479.8
Simplified79.8%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-*l*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
Applied egg-rr93.1%
lift-+.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift--.f64N/A
associate-*r*N/A
lift-cos.f64N/A
times-fracN/A
lower-*.f64N/A
Applied egg-rr98.4%
Final simplification86.4%
(FPCore (t l k)
:precision binary64
(if (<= (* l l) 1e-179)
(/ (/ (/ 2.0 (* k k)) (/ k l)) (/ (* k t) l))
(if (<= (* l l) 2e+305)
(/ 2.0 (* k (* (* t (/ k (* l l))) (* (sin k) (tan k)))))
(*
l
(/
(* l (* 2.0 (cos k)))
(* (fma (cos (+ k k)) -0.5 0.5) (* k (* k t))))))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-179) {
tmp = ((2.0 / (k * k)) / (k / l)) / ((k * t) / l);
} else if ((l * l) <= 2e+305) {
tmp = 2.0 / (k * ((t * (k / (l * l))) * (sin(k) * tan(k))));
} else {
tmp = l * ((l * (2.0 * cos(k))) / (fma(cos((k + k)), -0.5, 0.5) * (k * (k * t))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-179) tmp = Float64(Float64(Float64(2.0 / Float64(k * k)) / Float64(k / l)) / Float64(Float64(k * t) / l)); elseif (Float64(l * l) <= 2e+305) tmp = Float64(2.0 / Float64(k * Float64(Float64(t * Float64(k / Float64(l * l))) * Float64(sin(k) * tan(k))))); else tmp = Float64(l * Float64(Float64(l * Float64(2.0 * cos(k))) / Float64(fma(cos(Float64(k + k)), -0.5, 0.5) * Float64(k * Float64(k * t))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e-179], N[(N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] / N[(k / l), $MachinePrecision]), $MachinePrecision] / N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+305], N[(2.0 / N[(k * N[(N[(t * N[(k / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l * N[(N[(l * N[(2.0 * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[N[(k + k), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision] * N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-179}:\\
\;\;\;\;\frac{\frac{\frac{2}{k \cdot k}}{\frac{k}{\ell}}}{\frac{k \cdot t}{\ell}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+305}:\\
\;\;\;\;\frac{2}{k \cdot \left(\left(t \cdot \frac{k}{\ell \cdot \ell}\right) \cdot \left(\sin k \cdot \tan k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \frac{\ell \cdot \left(2 \cdot \cos k\right)}{\mathsf{fma}\left(\cos \left(k + k\right), -0.5, 0.5\right) \cdot \left(k \cdot \left(k \cdot t\right)\right)}\\
\end{array}
\end{array}
if (*.f64 l l) < 1e-179Initial program 31.9%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.0
Simplified71.0%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-/r*N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f6494.6
Applied egg-rr94.6%
if 1e-179 < (*.f64 l l) < 1.9999999999999999e305Initial program 46.8%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6494.3
Simplified94.3%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
Applied egg-rr99.7%
if 1.9999999999999999e305 < (*.f64 l l) Initial program 35.5%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6466.9
Simplified66.9%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-*l*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
Applied egg-rr88.3%
Applied egg-rr79.9%
Final simplification92.8%
(FPCore (t l k)
:precision binary64
(if (<= (* l l) 1e-179)
(/ (/ (/ 2.0 (* k k)) (/ k l)) (/ (* k t) l))
(if (<= (* l l) 2e+305)
(/ 2.0 (* k (* (* t (/ k (* l l))) (* (sin k) (tan k)))))
(*
(* l (cos k))
(* l (/ 2.0 (* (fma (cos (+ k k)) -0.5 0.5) (* k (* k t)))))))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-179) {
tmp = ((2.0 / (k * k)) / (k / l)) / ((k * t) / l);
} else if ((l * l) <= 2e+305) {
tmp = 2.0 / (k * ((t * (k / (l * l))) * (sin(k) * tan(k))));
} else {
tmp = (l * cos(k)) * (l * (2.0 / (fma(cos((k + k)), -0.5, 0.5) * (k * (k * t)))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-179) tmp = Float64(Float64(Float64(2.0 / Float64(k * k)) / Float64(k / l)) / Float64(Float64(k * t) / l)); elseif (Float64(l * l) <= 2e+305) tmp = Float64(2.0 / Float64(k * Float64(Float64(t * Float64(k / Float64(l * l))) * Float64(sin(k) * tan(k))))); else tmp = Float64(Float64(l * cos(k)) * Float64(l * Float64(2.0 / Float64(fma(cos(Float64(k + k)), -0.5, 0.5) * Float64(k * Float64(k * t)))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e-179], N[(N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] / N[(k / l), $MachinePrecision]), $MachinePrecision] / N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+305], N[(2.0 / N[(k * N[(N[(t * N[(k / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision] * N[(l * N[(2.0 / N[(N[(N[Cos[N[(k + k), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision] * N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-179}:\\
\;\;\;\;\frac{\frac{\frac{2}{k \cdot k}}{\frac{k}{\ell}}}{\frac{k \cdot t}{\ell}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+305}:\\
\;\;\;\;\frac{2}{k \cdot \left(\left(t \cdot \frac{k}{\ell \cdot \ell}\right) \cdot \left(\sin k \cdot \tan k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \cos k\right) \cdot \left(\ell \cdot \frac{2}{\mathsf{fma}\left(\cos \left(k + k\right), -0.5, 0.5\right) \cdot \left(k \cdot \left(k \cdot t\right)\right)}\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 1e-179Initial program 31.9%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.0
Simplified71.0%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-/r*N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f6494.6
Applied egg-rr94.6%
if 1e-179 < (*.f64 l l) < 1.9999999999999999e305Initial program 46.8%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6494.3
Simplified94.3%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
Applied egg-rr99.7%
if 1.9999999999999999e305 < (*.f64 l l) Initial program 35.5%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6466.9
Simplified66.9%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-*l*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
Applied egg-rr88.3%
Applied egg-rr79.9%
Final simplification92.8%
(FPCore (t l k)
:precision binary64
(if (<= k 1.15e-19)
(/ 2.0 (* (* (/ k l) (* k k)) (/ (* k t) l)))
(*
(* 2.0 (/ l k))
(* (/ (cos k) k) (/ l (* t (fma (cos (+ k k)) -0.5 0.5)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 1.15e-19) {
tmp = 2.0 / (((k / l) * (k * k)) * ((k * t) / l));
} else {
tmp = (2.0 * (l / k)) * ((cos(k) / k) * (l / (t * fma(cos((k + k)), -0.5, 0.5))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (k <= 1.15e-19) tmp = Float64(2.0 / Float64(Float64(Float64(k / l) * Float64(k * k)) * Float64(Float64(k * t) / l))); else tmp = Float64(Float64(2.0 * Float64(l / k)) * Float64(Float64(cos(k) / k) * Float64(l / Float64(t * fma(cos(Float64(k + k)), -0.5, 0.5))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[k, 1.15e-19], N[(2.0 / N[(N[(N[(k / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(l / k), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Cos[k], $MachinePrecision] / k), $MachinePrecision] * N[(l / N[(t * N[(N[Cos[N[(k + k), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.15 \cdot 10^{-19}:\\
\;\;\;\;\frac{2}{\left(\frac{k}{\ell} \cdot \left(k \cdot k\right)\right) \cdot \frac{k \cdot t}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{k}\right) \cdot \left(\frac{\cos k}{k} \cdot \frac{\ell}{t \cdot \mathsf{fma}\left(\cos \left(k + k\right), -0.5, 0.5\right)}\right)\\
\end{array}
\end{array}
if k < 1.1499999999999999e-19Initial program 40.1%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6472.1
Simplified72.1%
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6481.9
Applied egg-rr81.9%
if 1.1499999999999999e-19 < k Initial program 36.0%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6479.8
Simplified79.8%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-*l*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
Applied egg-rr93.1%
Applied egg-rr93.9%
lift-cos.f64N/A
*-commutativeN/A
lift-+.f64N/A
lift-cos.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6497.0
Applied egg-rr97.0%
Final simplification86.1%
(FPCore (t l k)
:precision binary64
(if (<= k 1.15e-19)
(/ 2.0 (* (* (/ k l) (* k k)) (/ (* k t) l)))
(*
(* 2.0 (/ l k))
(* l (/ (cos k) (* (* k t) (fma (cos (+ k k)) -0.5 0.5)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 1.15e-19) {
tmp = 2.0 / (((k / l) * (k * k)) * ((k * t) / l));
} else {
tmp = (2.0 * (l / k)) * (l * (cos(k) / ((k * t) * fma(cos((k + k)), -0.5, 0.5))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (k <= 1.15e-19) tmp = Float64(2.0 / Float64(Float64(Float64(k / l) * Float64(k * k)) * Float64(Float64(k * t) / l))); else tmp = Float64(Float64(2.0 * Float64(l / k)) * Float64(l * Float64(cos(k) / Float64(Float64(k * t) * fma(cos(Float64(k + k)), -0.5, 0.5))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[k, 1.15e-19], N[(2.0 / N[(N[(N[(k / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(l / k), $MachinePrecision]), $MachinePrecision] * N[(l * N[(N[Cos[k], $MachinePrecision] / N[(N[(k * t), $MachinePrecision] * N[(N[Cos[N[(k + k), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.15 \cdot 10^{-19}:\\
\;\;\;\;\frac{2}{\left(\frac{k}{\ell} \cdot \left(k \cdot k\right)\right) \cdot \frac{k \cdot t}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{k}\right) \cdot \left(\ell \cdot \frac{\cos k}{\left(k \cdot t\right) \cdot \mathsf{fma}\left(\cos \left(k + k\right), -0.5, 0.5\right)}\right)\\
\end{array}
\end{array}
if k < 1.1499999999999999e-19Initial program 40.1%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6472.1
Simplified72.1%
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6481.9
Applied egg-rr81.9%
if 1.1499999999999999e-19 < k Initial program 36.0%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6479.8
Simplified79.8%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-*l*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f64N/A
Applied egg-rr93.1%
Applied egg-rr93.9%
lift-cos.f64N/A
lift-+.f64N/A
lift-cos.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
Applied egg-rr94.0%
Final simplification85.2%
(FPCore (t l k) :precision binary64 (if (<= (* l l) 1e-179) (/ (/ (/ 2.0 (* k k)) (/ k l)) (/ (* k t) l)) (/ 2.0 (* k (* (* t (/ k (* l l))) (* (sin k) (tan k)))))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-179) {
tmp = ((2.0 / (k * k)) / (k / l)) / ((k * t) / l);
} else {
tmp = 2.0 / (k * ((t * (k / (l * 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 ((l * l) <= 1d-179) then
tmp = ((2.0d0 / (k * k)) / (k / l)) / ((k * t) / l)
else
tmp = 2.0d0 / (k * ((t * (k / (l * l))) * (sin(k) * tan(k))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-179) {
tmp = ((2.0 / (k * k)) / (k / l)) / ((k * t) / l);
} else {
tmp = 2.0 / (k * ((t * (k / (l * l))) * (Math.sin(k) * Math.tan(k))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (l * l) <= 1e-179: tmp = ((2.0 / (k * k)) / (k / l)) / ((k * t) / l) else: tmp = 2.0 / (k * ((t * (k / (l * l))) * (math.sin(k) * math.tan(k)))) return tmp
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-179) tmp = Float64(Float64(Float64(2.0 / Float64(k * k)) / Float64(k / l)) / Float64(Float64(k * t) / l)); else tmp = Float64(2.0 / Float64(k * Float64(Float64(t * Float64(k / Float64(l * l))) * Float64(sin(k) * tan(k))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((l * l) <= 1e-179) tmp = ((2.0 / (k * k)) / (k / l)) / ((k * t) / l); else tmp = 2.0 / (k * ((t * (k / (l * l))) * (sin(k) * tan(k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e-179], N[(N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] / N[(k / l), $MachinePrecision]), $MachinePrecision] / N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(k * N[(N[(t * N[(k / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-179}:\\
\;\;\;\;\frac{\frac{\frac{2}{k \cdot k}}{\frac{k}{\ell}}}{\frac{k \cdot t}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{k \cdot \left(\left(t \cdot \frac{k}{\ell \cdot \ell}\right) \cdot \left(\sin k \cdot \tan k\right)\right)}\\
\end{array}
\end{array}
if (*.f64 l l) < 1e-179Initial program 31.9%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.0
Simplified71.0%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-/r*N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f6494.6
Applied egg-rr94.6%
if 1e-179 < (*.f64 l l) Initial program 42.3%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6483.5
Simplified83.5%
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
Applied egg-rr87.0%
Final simplification89.5%
(FPCore (t l k) :precision binary64 (if (<= (* l l) 1e-179) (/ (/ (/ 2.0 (* k k)) (/ k l)) (/ (* k t) l)) (/ 2.0 (* k (/ (* k (* k (* k t))) (* (* l l) (cos k)))))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-179) {
tmp = ((2.0 / (k * k)) / (k / l)) / ((k * t) / l);
} else {
tmp = 2.0 / (k * ((k * (k * (k * t))) / ((l * l) * cos(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-179) then
tmp = ((2.0d0 / (k * k)) / (k / l)) / ((k * t) / l)
else
tmp = 2.0d0 / (k * ((k * (k * (k * t))) / ((l * l) * cos(k))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-179) {
tmp = ((2.0 / (k * k)) / (k / l)) / ((k * t) / l);
} else {
tmp = 2.0 / (k * ((k * (k * (k * t))) / ((l * l) * Math.cos(k))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (l * l) <= 1e-179: tmp = ((2.0 / (k * k)) / (k / l)) / ((k * t) / l) else: tmp = 2.0 / (k * ((k * (k * (k * t))) / ((l * l) * math.cos(k)))) return tmp
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-179) tmp = Float64(Float64(Float64(2.0 / Float64(k * k)) / Float64(k / l)) / Float64(Float64(k * t) / l)); else tmp = Float64(2.0 / Float64(k * Float64(Float64(k * Float64(k * Float64(k * t))) / Float64(Float64(l * l) * cos(k))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((l * l) <= 1e-179) tmp = ((2.0 / (k * k)) / (k / l)) / ((k * t) / l); else tmp = 2.0 / (k * ((k * (k * (k * t))) / ((l * l) * cos(k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e-179], N[(N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] / N[(k / l), $MachinePrecision]), $MachinePrecision] / N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(k * N[(N[(k * N[(k * N[(k * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(l * l), $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-179}:\\
\;\;\;\;\frac{\frac{\frac{2}{k \cdot k}}{\frac{k}{\ell}}}{\frac{k \cdot t}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{k \cdot \frac{k \cdot \left(k \cdot \left(k \cdot t\right)\right)}{\left(\ell \cdot \ell\right) \cdot \cos k}}\\
\end{array}
\end{array}
if (*.f64 l l) < 1e-179Initial program 31.9%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.0
Simplified71.0%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-/r*N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f6494.6
Applied egg-rr94.6%
if 1e-179 < (*.f64 l l) Initial program 42.3%
Taylor expanded in t around 0
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-cos.f6483.5
Simplified83.5%
Taylor expanded in k around 0
unpow2N/A
lower-*.f6468.2
Simplified68.2%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6469.4
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-*.f6470.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f6470.0
Applied egg-rr70.0%
Final simplification78.0%
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (/ k l) (* k k)) (/ (* k t) l))))
double code(double t, double l, double k) {
return 2.0 / (((k / l) * (k * k)) * ((k * t) / l));
}
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 / (((k / l) * (k * k)) * ((k * t) / l))
end function
public static double code(double t, double l, double k) {
return 2.0 / (((k / l) * (k * k)) * ((k * t) / l));
}
def code(t, l, k): return 2.0 / (((k / l) * (k * k)) * ((k * t) / l))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(k / l) * Float64(k * k)) * Float64(Float64(k * t) / l))) end
function tmp = code(t, l, k) tmp = 2.0 / (((k / l) * (k * k)) * ((k * t) / l)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(k / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(\frac{k}{\ell} \cdot \left(k \cdot k\right)\right) \cdot \frac{k \cdot t}{\ell}}
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.9
Simplified68.9%
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6475.7
Applied egg-rr75.7%
Final simplification75.7%
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* k k) (* (/ k l) (/ (* k t) l)))))
double code(double t, double l, double k) {
return 2.0 / ((k * k) * ((k / l) * ((k * t) / l)));
}
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 / ((k * k) * ((k / l) * ((k * t) / l)))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((k * k) * ((k / l) * ((k * t) / l)));
}
def code(t, l, k): return 2.0 / ((k * k) * ((k / l) * ((k * t) / l)))
function code(t, l, k) return Float64(2.0 / Float64(Float64(k * k) * Float64(Float64(k / l) * Float64(Float64(k * t) / l)))) end
function tmp = code(t, l, k) tmp = 2.0 / ((k * k) * ((k / l) * ((k * t) / l))); end
code[t_, l_, k_] := N[(2.0 / N[(N[(k * k), $MachinePrecision] * N[(N[(k / l), $MachinePrecision] * N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(k \cdot k\right) \cdot \left(\frac{k}{\ell} \cdot \frac{k \cdot t}{\ell}\right)}
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.9
Simplified68.9%
lift-*.f64N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6475.3
Applied egg-rr75.3%
Final simplification75.3%
(FPCore (t l k) :precision binary64 (* (/ l (* k t)) (* 2.0 (/ l (* k (* k k))))))
double code(double t, double l, double k) {
return (l / (k * t)) * (2.0 * (l / (k * (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 = (l / (k * t)) * (2.0d0 * (l / (k * (k * k))))
end function
public static double code(double t, double l, double k) {
return (l / (k * t)) * (2.0 * (l / (k * (k * k))));
}
def code(t, l, k): return (l / (k * t)) * (2.0 * (l / (k * (k * k))))
function code(t, l, k) return Float64(Float64(l / Float64(k * t)) * Float64(2.0 * Float64(l / Float64(k * Float64(k * k))))) end
function tmp = code(t, l, k) tmp = (l / (k * t)) * (2.0 * (l / (k * (k * k)))); end
code[t_, l_, k_] := N[(N[(l / N[(k * t), $MachinePrecision]), $MachinePrecision] * N[(2.0 * N[(l / N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell}{k \cdot t} \cdot \left(2 \cdot \frac{\ell}{k \cdot \left(k \cdot k\right)}\right)
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.9
Simplified68.9%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-/.f64N/A
associate-*r/N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6468.0
Applied egg-rr68.0%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l/N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6468.0
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lift-*.f64N/A
cube-multN/A
lower-*.f64N/A
cube-multN/A
lift-*.f64N/A
lower-*.f6467.6
Applied egg-rr67.6%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f6474.9
Applied egg-rr74.9%
Final simplification74.9%
(FPCore (t l k) :precision binary64 (* 2.0 (* (/ l (* t (* k k))) (/ l (* k k)))))
double code(double t, double l, double k) {
return 2.0 * ((l / (t * (k * k))) * (l / (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 / (t * (k * k))) * (l / (k * k)))
end function
public static double code(double t, double l, double k) {
return 2.0 * ((l / (t * (k * k))) * (l / (k * k)));
}
def code(t, l, k): return 2.0 * ((l / (t * (k * k))) * (l / (k * k)))
function code(t, l, k) return Float64(2.0 * Float64(Float64(l / Float64(t * Float64(k * k))) * Float64(l / Float64(k * k)))) end
function tmp = code(t, l, k) tmp = 2.0 * ((l / (t * (k * k))) * (l / (k * k))); end
code[t_, l_, k_] := N[(2.0 * N[(N[(l / N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \left(\frac{\ell}{t \cdot \left(k \cdot k\right)} \cdot \frac{\ell}{k \cdot k}\right)
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.9
Simplified68.9%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-/.f64N/A
associate-*r/N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6468.0
Applied egg-rr68.0%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l/N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6468.0
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lift-*.f64N/A
cube-multN/A
lower-*.f64N/A
cube-multN/A
lift-*.f64N/A
lower-*.f6467.6
Applied egg-rr67.6%
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6474.9
Applied egg-rr74.9%
(FPCore (t l k) :precision binary64 (* 2.0 (* l (/ l (* k (* k (* t (* k k))))))))
double code(double t, double l, double k) {
return 2.0 * (l * (l / (k * (k * (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 = 2.0d0 * (l * (l / (k * (k * (t * (k * k))))))
end function
public static double code(double t, double l, double k) {
return 2.0 * (l * (l / (k * (k * (t * (k * k))))));
}
def code(t, l, k): return 2.0 * (l * (l / (k * (k * (t * (k * k))))))
function code(t, l, k) return Float64(2.0 * Float64(l * Float64(l / Float64(k * Float64(k * Float64(t * Float64(k * k))))))) end
function tmp = code(t, l, k) tmp = 2.0 * (l * (l / (k * (k * (t * (k * k)))))); end
code[t_, l_, k_] := N[(2.0 * N[(l * N[(l / N[(k * N[(k * N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \left(\ell \cdot \frac{\ell}{k \cdot \left(k \cdot \left(t \cdot \left(k \cdot k\right)\right)\right)}\right)
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
associate-/l*N/A
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.9
Simplified68.9%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-/.f64N/A
associate-*r/N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6468.0
Applied egg-rr68.0%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l/N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6468.0
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lift-*.f64N/A
cube-multN/A
lower-*.f64N/A
cube-multN/A
lift-*.f64N/A
lower-*.f6467.6
Applied egg-rr67.6%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lower-*.f64N/A
Applied egg-rr72.7%
Final simplification72.7%
herbie shell --seed 2024219
(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))))