
(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 19 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}
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.25e-60)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ (* (cos k_m) l) k_m)))
(/
2.0
(* (/ (/ k_m (cos k_m)) l) (* (pow (sin k_m) 2.0) (* t (/ k_m l)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.25e-60) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / ((cos(k_m) * l) / k_m));
} else {
tmp = 2.0 / (((k_m / cos(k_m)) / l) * (pow(sin(k_m), 2.0) * (t * (k_m / l))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.25e-60) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(Float64(cos(k_m) * l) / k_m))); else tmp = Float64(2.0 / Float64(Float64(Float64(k_m / cos(k_m)) / l) * Float64((sin(k_m) ^ 2.0) * Float64(t * Float64(k_m / l))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.25e-60], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(k$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * N[(t * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.25 \cdot 10^{-60}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{\cos k\_m \cdot \ell}{k\_m}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\frac{k\_m}{\cos k\_m}}{\ell} \cdot \left({\sin k\_m}^{2} \cdot \left(t \cdot \frac{k\_m}{\ell}\right)\right)}\\
\end{array}
\end{array}
if k < 1.25e-60Initial program 38.3%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6493.9
Applied rewrites93.9%
Taylor expanded in k around 0
Applied rewrites81.7%
Applied rewrites84.8%
if 1.25e-60 < k Initial program 28.8%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6489.0
Applied rewrites89.0%
Applied rewrites97.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= (* l l) 2e-221)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ (* (cos k_m) l) k_m)))
(if (<= (* l l) 5e+306)
(/ 2.0 (* (pow (* (sin k_m) k_m) 2.0) (/ t (* (* l (cos k_m)) l))))
(/
2.0
(*
(/ t l)
(*
(* (* (* (tan k_m) (sin k_m)) (* (/ t l) t)) (/ k_m t))
(/ k_m t)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if ((l * l) <= 2e-221) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / ((cos(k_m) * l) / k_m));
} else if ((l * l) <= 5e+306) {
tmp = 2.0 / (pow((sin(k_m) * k_m), 2.0) * (t / ((l * cos(k_m)) * l)));
} else {
tmp = 2.0 / ((t / l) * ((((tan(k_m) * sin(k_m)) * ((t / l) * t)) * (k_m / t)) * (k_m / t)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (Float64(l * l) <= 2e-221) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(Float64(cos(k_m) * l) / k_m))); elseif (Float64(l * l) <= 5e+306) tmp = Float64(2.0 / Float64((Float64(sin(k_m) * k_m) ^ 2.0) * Float64(t / Float64(Float64(l * cos(k_m)) * l)))); else tmp = Float64(2.0 / Float64(Float64(t / l) * Float64(Float64(Float64(Float64(tan(k_m) * sin(k_m)) * Float64(Float64(t / l) * t)) * Float64(k_m / t)) * Float64(k_m / t)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[N[(l * l), $MachinePrecision], 2e-221], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 5e+306], N[(2.0 / N[(N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * k$95$m), $MachinePrecision], 2.0], $MachinePrecision] * N[(t / N[(N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t / l), $MachinePrecision] * N[(N[(N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 2 \cdot 10^{-221}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{\cos k\_m \cdot \ell}{k\_m}}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 5 \cdot 10^{+306}:\\
\;\;\;\;\frac{2}{{\left(\sin k\_m \cdot k\_m\right)}^{2} \cdot \frac{t}{\left(\ell \cdot \cos k\_m\right) \cdot \ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{t}{\ell} \cdot \left(\left(\left(\left(\tan k\_m \cdot \sin k\_m\right) \cdot \left(\frac{t}{\ell} \cdot t\right)\right) \cdot \frac{k\_m}{t}\right) \cdot \frac{k\_m}{t}\right)}\\
\end{array}
\end{array}
if (*.f64 l l) < 2.00000000000000003e-221Initial program 24.8%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6489.2
Applied rewrites89.2%
Taylor expanded in k around 0
Applied rewrites84.6%
Applied rewrites91.5%
if 2.00000000000000003e-221 < (*.f64 l l) < 4.99999999999999993e306Initial program 41.8%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6496.0
Applied rewrites96.0%
Applied rewrites98.5%
Applied rewrites87.2%
Applied rewrites87.6%
if 4.99999999999999993e306 < (*.f64 l l) Initial program 37.8%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
associate-*l*N/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites57.8%
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-pow.f64N/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6463.5
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
lift-/.f64N/A
*-commutativeN/A
lower-*.f6480.1
Applied rewrites80.1%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (* (cos k_m) l)))
(if (<= k_m 3e-66)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ t_1 k_m)))
(if (<= k_m 1.5e+178)
(/ 2.0 (* (/ (pow (* (sin k_m) k_m) 2.0) l) (/ t (* l (cos k_m)))))
(/ 2.0 (* (pow (sin k_m) 2.0) (* (* k_m t) (/ k_m (* t_1 l)))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = cos(k_m) * l;
double tmp;
if (k_m <= 3e-66) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / (t_1 / k_m));
} else if (k_m <= 1.5e+178) {
tmp = 2.0 / ((pow((sin(k_m) * k_m), 2.0) / l) * (t / (l * cos(k_m))));
} else {
tmp = 2.0 / (pow(sin(k_m), 2.0) * ((k_m * t) * (k_m / (t_1 * l))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(cos(k_m) * l) tmp = 0.0 if (k_m <= 3e-66) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(t_1 / k_m))); elseif (k_m <= 1.5e+178) tmp = Float64(2.0 / Float64(Float64((Float64(sin(k_m) * k_m) ^ 2.0) / l) * Float64(t / Float64(l * cos(k_m))))); else tmp = Float64(2.0 / Float64((sin(k_m) ^ 2.0) * Float64(Float64(k_m * t) * Float64(k_m / Float64(t_1 * l))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision]}, If[LessEqual[k$95$m, 3e-66], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(t$95$1 / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.5e+178], N[(2.0 / N[(N[(N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t / N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * N[(N[(k$95$m * t), $MachinePrecision] * N[(k$95$m / N[(t$95$1 * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := \cos k\_m \cdot \ell\\
\mathbf{if}\;k\_m \leq 3 \cdot 10^{-66}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{t\_1}{k\_m}}}\\
\mathbf{elif}\;k\_m \leq 1.5 \cdot 10^{+178}:\\
\;\;\;\;\frac{2}{\frac{{\left(\sin k\_m \cdot k\_m\right)}^{2}}{\ell} \cdot \frac{t}{\ell \cdot \cos k\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\sin k\_m}^{2} \cdot \left(\left(k\_m \cdot t\right) \cdot \frac{k\_m}{t\_1 \cdot \ell}\right)}\\
\end{array}
\end{array}
if k < 3.0000000000000002e-66Initial program 38.5%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6494.1
Applied rewrites94.1%
Taylor expanded in k around 0
Applied rewrites81.7%
Applied rewrites84.9%
if 3.0000000000000002e-66 < k < 1.50000000000000008e178Initial program 21.4%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6487.7
Applied rewrites87.7%
Applied rewrites97.5%
Applied rewrites73.8%
Applied rewrites93.3%
if 1.50000000000000008e178 < k Initial program 40.4%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6490.5
Applied rewrites90.5%
Applied rewrites97.1%
Applied rewrites83.1%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (pow (sin k_m) 2.0)) (t_2 (* (cos k_m) l)))
(if (<= k_m 1.05e-59)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ t_2 k_m)))
(if (<= k_m 2e+118)
(* (/ (* 2.0 (cos k_m)) (* (* k_m k_m) t_1)) (* l (/ l t)))
(/ 2.0 (* t_1 (* (* k_m t) (/ k_m (* t_2 l)))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = pow(sin(k_m), 2.0);
double t_2 = cos(k_m) * l;
double tmp;
if (k_m <= 1.05e-59) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / (t_2 / k_m));
} else if (k_m <= 2e+118) {
tmp = ((2.0 * cos(k_m)) / ((k_m * k_m) * t_1)) * (l * (l / t));
} else {
tmp = 2.0 / (t_1 * ((k_m * t) * (k_m / (t_2 * l))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = sin(k_m) ^ 2.0 t_2 = Float64(cos(k_m) * l) tmp = 0.0 if (k_m <= 1.05e-59) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(t_2 / k_m))); elseif (k_m <= 2e+118) tmp = Float64(Float64(Float64(2.0 * cos(k_m)) / Float64(Float64(k_m * k_m) * t_1)) * Float64(l * Float64(l / t))); else tmp = Float64(2.0 / Float64(t_1 * Float64(Float64(k_m * t) * Float64(k_m / Float64(t_2 * l))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision]}, If[LessEqual[k$95$m, 1.05e-59], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(t$95$2 / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2e+118], N[(N[(N[(2.0 * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] * N[(l * N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(t$95$1 * N[(N[(k$95$m * t), $MachinePrecision] * N[(k$95$m / N[(t$95$2 * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := {\sin k\_m}^{2}\\
t_2 := \cos k\_m \cdot \ell\\
\mathbf{if}\;k\_m \leq 1.05 \cdot 10^{-59}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{t\_2}{k\_m}}}\\
\mathbf{elif}\;k\_m \leq 2 \cdot 10^{+118}:\\
\;\;\;\;\frac{2 \cdot \cos k\_m}{\left(k\_m \cdot k\_m\right) \cdot t\_1} \cdot \left(\ell \cdot \frac{\ell}{t}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t\_1 \cdot \left(\left(k\_m \cdot t\right) \cdot \frac{k\_m}{t\_2 \cdot \ell}\right)}\\
\end{array}
\end{array}
if k < 1.04999999999999998e-59Initial program 38.3%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6493.9
Applied rewrites93.9%
Taylor expanded in k around 0
Applied rewrites81.7%
Applied rewrites84.8%
if 1.04999999999999998e-59 < k < 1.99999999999999993e118Initial program 27.3%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6488.5
Applied rewrites88.5%
Applied rewrites99.6%
Taylor expanded in t around 0
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
unpow2N/A
associate-/l*N/A
Applied rewrites94.3%
if 1.99999999999999993e118 < k Initial program 30.0%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6489.4
Applied rewrites89.4%
Applied rewrites96.2%
Applied rewrites79.6%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (pow (sin k_m) 2.0)) (t_2 (* (cos k_m) l)))
(if (<= k_m 1.05e-59)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ t_2 k_m)))
(if (<= k_m 2e+118)
(* (/ (* 2.0 (cos k_m)) (* (* k_m k_m) t_1)) (* l (/ l t)))
(/ 2.0 (* k_m (/ (* (* t_1 t) k_m) (* t_2 l))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = pow(sin(k_m), 2.0);
double t_2 = cos(k_m) * l;
double tmp;
if (k_m <= 1.05e-59) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / (t_2 / k_m));
} else if (k_m <= 2e+118) {
tmp = ((2.0 * cos(k_m)) / ((k_m * k_m) * t_1)) * (l * (l / t));
} else {
tmp = 2.0 / (k_m * (((t_1 * t) * k_m) / (t_2 * l)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = sin(k_m) ^ 2.0 t_2 = Float64(cos(k_m) * l) tmp = 0.0 if (k_m <= 1.05e-59) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(t_2 / k_m))); elseif (k_m <= 2e+118) tmp = Float64(Float64(Float64(2.0 * cos(k_m)) / Float64(Float64(k_m * k_m) * t_1)) * Float64(l * Float64(l / t))); else tmp = Float64(2.0 / Float64(k_m * Float64(Float64(Float64(t_1 * t) * k_m) / Float64(t_2 * l)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision]}, If[LessEqual[k$95$m, 1.05e-59], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(t$95$2 / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2e+118], N[(N[(N[(2.0 * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] * N[(l * N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(k$95$m * N[(N[(N[(t$95$1 * t), $MachinePrecision] * k$95$m), $MachinePrecision] / N[(t$95$2 * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := {\sin k\_m}^{2}\\
t_2 := \cos k\_m \cdot \ell\\
\mathbf{if}\;k\_m \leq 1.05 \cdot 10^{-59}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{t\_2}{k\_m}}}\\
\mathbf{elif}\;k\_m \leq 2 \cdot 10^{+118}:\\
\;\;\;\;\frac{2 \cdot \cos k\_m}{\left(k\_m \cdot k\_m\right) \cdot t\_1} \cdot \left(\ell \cdot \frac{\ell}{t}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{k\_m \cdot \frac{\left(t\_1 \cdot t\right) \cdot k\_m}{t\_2 \cdot \ell}}\\
\end{array}
\end{array}
if k < 1.04999999999999998e-59Initial program 38.3%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6493.9
Applied rewrites93.9%
Taylor expanded in k around 0
Applied rewrites81.7%
Applied rewrites84.8%
if 1.04999999999999998e-59 < k < 1.99999999999999993e118Initial program 27.3%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6488.5
Applied rewrites88.5%
Applied rewrites99.6%
Taylor expanded in t around 0
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
unpow2N/A
associate-/l*N/A
Applied rewrites94.3%
if 1.99999999999999993e118 < k Initial program 30.0%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6489.4
Applied rewrites89.4%
Applied rewrites96.2%
Applied rewrites77.3%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (pow (sin k_m) 2.0)) (t_2 (* (cos k_m) l)))
(if (<= k_m 1.05e-59)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ t_2 k_m)))
(if (<= k_m 1.85e+81)
(* (/ (* 2.0 (cos k_m)) (* (* k_m k_m) t_1)) (* l (/ l t)))
(/ 2.0 (* k_m (* (* t_1 t) (/ k_m (* t_2 l)))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = pow(sin(k_m), 2.0);
double t_2 = cos(k_m) * l;
double tmp;
if (k_m <= 1.05e-59) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / (t_2 / k_m));
} else if (k_m <= 1.85e+81) {
tmp = ((2.0 * cos(k_m)) / ((k_m * k_m) * t_1)) * (l * (l / t));
} else {
tmp = 2.0 / (k_m * ((t_1 * t) * (k_m / (t_2 * l))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = sin(k_m) ^ 2.0 t_2 = Float64(cos(k_m) * l) tmp = 0.0 if (k_m <= 1.05e-59) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(t_2 / k_m))); elseif (k_m <= 1.85e+81) tmp = Float64(Float64(Float64(2.0 * cos(k_m)) / Float64(Float64(k_m * k_m) * t_1)) * Float64(l * Float64(l / t))); else tmp = Float64(2.0 / Float64(k_m * Float64(Float64(t_1 * t) * Float64(k_m / Float64(t_2 * l))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision]}, If[LessEqual[k$95$m, 1.05e-59], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(t$95$2 / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.85e+81], N[(N[(N[(2.0 * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] * N[(l * N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(k$95$m * N[(N[(t$95$1 * t), $MachinePrecision] * N[(k$95$m / N[(t$95$2 * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := {\sin k\_m}^{2}\\
t_2 := \cos k\_m \cdot \ell\\
\mathbf{if}\;k\_m \leq 1.05 \cdot 10^{-59}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{t\_2}{k\_m}}}\\
\mathbf{elif}\;k\_m \leq 1.85 \cdot 10^{+81}:\\
\;\;\;\;\frac{2 \cdot \cos k\_m}{\left(k\_m \cdot k\_m\right) \cdot t\_1} \cdot \left(\ell \cdot \frac{\ell}{t}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{k\_m \cdot \left(\left(t\_1 \cdot t\right) \cdot \frac{k\_m}{t\_2 \cdot \ell}\right)}\\
\end{array}
\end{array}
if k < 1.04999999999999998e-59Initial program 38.3%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6493.9
Applied rewrites93.9%
Taylor expanded in k around 0
Applied rewrites81.7%
Applied rewrites84.8%
if 1.04999999999999998e-59 < k < 1.85e81Initial program 27.5%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6487.0
Applied rewrites87.0%
Applied rewrites99.6%
Taylor expanded in t around 0
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
unpow2N/A
associate-/l*N/A
Applied rewrites93.6%
if 1.85e81 < k Initial program 29.7%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6490.3
Applied rewrites90.3%
Applied rewrites96.5%
Applied rewrites79.4%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= l 9e-161)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ (* (cos k_m) l) k_m)))
(if (<= l 1.02e+154)
(*
(/ (* 2.0 (cos k_m)) (* (* k_m k_m) t))
(/ (* l l) (pow (sin k_m) 2.0)))
(/
2.0
(*
(/ t l)
(*
(* (* (* (tan k_m) (sin k_m)) (* (/ t l) t)) (/ k_m t))
(/ k_m t)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (l <= 9e-161) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / ((cos(k_m) * l) / k_m));
} else if (l <= 1.02e+154) {
tmp = ((2.0 * cos(k_m)) / ((k_m * k_m) * t)) * ((l * l) / pow(sin(k_m), 2.0));
} else {
tmp = 2.0 / ((t / l) * ((((tan(k_m) * sin(k_m)) * ((t / l) * t)) * (k_m / t)) * (k_m / t)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (l <= 9e-161) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(Float64(cos(k_m) * l) / k_m))); elseif (l <= 1.02e+154) tmp = Float64(Float64(Float64(2.0 * cos(k_m)) / Float64(Float64(k_m * k_m) * t)) * Float64(Float64(l * l) / (sin(k_m) ^ 2.0))); else tmp = Float64(2.0 / Float64(Float64(t / l) * Float64(Float64(Float64(Float64(tan(k_m) * sin(k_m)) * Float64(Float64(t / l) * t)) * Float64(k_m / t)) * Float64(k_m / t)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[l, 9e-161], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[l, 1.02e+154], N[(N[(N[(2.0 * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t / l), $MachinePrecision] * N[(N[(N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;\ell \leq 9 \cdot 10^{-161}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{\cos k\_m \cdot \ell}{k\_m}}}\\
\mathbf{elif}\;\ell \leq 1.02 \cdot 10^{+154}:\\
\;\;\;\;\frac{2 \cdot \cos k\_m}{\left(k\_m \cdot k\_m\right) \cdot t} \cdot \frac{\ell \cdot \ell}{{\sin k\_m}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{t}{\ell} \cdot \left(\left(\left(\left(\tan k\_m \cdot \sin k\_m\right) \cdot \left(\frac{t}{\ell} \cdot t\right)\right) \cdot \frac{k\_m}{t}\right) \cdot \frac{k\_m}{t}\right)}\\
\end{array}
\end{array}
if l < 8.9999999999999993e-161Initial program 33.9%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6491.4
Applied rewrites91.4%
Taylor expanded in k around 0
Applied rewrites75.7%
Applied rewrites78.7%
if 8.9999999999999993e-161 < l < 1.02000000000000007e154Initial program 34.0%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6495.9
Applied rewrites95.9%
Taylor expanded in t around 0
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6489.7
Applied rewrites89.7%
if 1.02000000000000007e154 < l Initial program 45.8%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
associate-*l*N/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites60.7%
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-pow.f64N/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6461.7
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
lift-/.f64N/A
*-commutativeN/A
lower-*.f6476.7
Applied rewrites76.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (* (cos k_m) l)))
(if (<= k_m 1.2e-60)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ t_1 k_m)))
(/ 2.0 (* (/ k_m t_1) (* (* (pow (sin k_m) 2.0) (/ k_m l)) t))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = cos(k_m) * l;
double tmp;
if (k_m <= 1.2e-60) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / (t_1 / k_m));
} else {
tmp = 2.0 / ((k_m / t_1) * ((pow(sin(k_m), 2.0) * (k_m / l)) * t));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(cos(k_m) * l) tmp = 0.0 if (k_m <= 1.2e-60) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(t_1 / k_m))); else tmp = Float64(2.0 / Float64(Float64(k_m / t_1) * Float64(Float64((sin(k_m) ^ 2.0) * Float64(k_m / l)) * t))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision]}, If[LessEqual[k$95$m, 1.2e-60], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(t$95$1 / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k$95$m / t$95$1), $MachinePrecision] * N[(N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := \cos k\_m \cdot \ell\\
\mathbf{if}\;k\_m \leq 1.2 \cdot 10^{-60}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{t\_1}{k\_m}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{k\_m}{t\_1} \cdot \left(\left({\sin k\_m}^{2} \cdot \frac{k\_m}{\ell}\right) \cdot t\right)}\\
\end{array}
\end{array}
if k < 1.20000000000000005e-60Initial program 38.3%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6493.9
Applied rewrites93.9%
Taylor expanded in k around 0
Applied rewrites81.7%
Applied rewrites84.8%
if 1.20000000000000005e-60 < k Initial program 28.8%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6489.0
Applied rewrites89.0%
Applied rewrites97.7%
Applied rewrites97.6%
Applied rewrites97.6%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.05e-59)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ (* (cos k_m) l) k_m)))
(*
(/ (* 2.0 (cos k_m)) (* (* k_m k_m) (pow (sin k_m) 2.0)))
(* l (/ l t)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.05e-59) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / ((cos(k_m) * l) / k_m));
} else {
tmp = ((2.0 * cos(k_m)) / ((k_m * k_m) * pow(sin(k_m), 2.0))) * (l * (l / t));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.05e-59) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(Float64(cos(k_m) * l) / k_m))); else tmp = Float64(Float64(Float64(2.0 * cos(k_m)) / Float64(Float64(k_m * k_m) * (sin(k_m) ^ 2.0))) * Float64(l * Float64(l / t))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.05e-59], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l * N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.05 \cdot 10^{-59}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{\cos k\_m \cdot \ell}{k\_m}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \cos k\_m}{\left(k\_m \cdot k\_m\right) \cdot {\sin k\_m}^{2}} \cdot \left(\ell \cdot \frac{\ell}{t}\right)\\
\end{array}
\end{array}
if k < 1.04999999999999998e-59Initial program 38.3%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6493.9
Applied rewrites93.9%
Taylor expanded in k around 0
Applied rewrites81.7%
Applied rewrites84.8%
if 1.04999999999999998e-59 < k Initial program 28.8%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6489.0
Applied rewrites89.0%
Applied rewrites97.7%
Taylor expanded in t around 0
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
unpow2N/A
associate-/l*N/A
Applied rewrites77.1%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 14600000000.0)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ (* (cos k_m) l) k_m)))
(if (<= k_m 4.8e+210)
(/
2.0
(*
(/ t l)
(* (* (* (* (tan k_m) (sin k_m)) (* (/ t l) t)) (/ k_m t)) (/ k_m t))))
(/ 2.0 (* (/ k_m l) (/ (* (* (pow (sin k_m) 2.0) t) k_m) l))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 14600000000.0) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / ((cos(k_m) * l) / k_m));
} else if (k_m <= 4.8e+210) {
tmp = 2.0 / ((t / l) * ((((tan(k_m) * sin(k_m)) * ((t / l) * t)) * (k_m / t)) * (k_m / t)));
} else {
tmp = 2.0 / ((k_m / l) * (((pow(sin(k_m), 2.0) * t) * k_m) / l));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 14600000000.0) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(Float64(cos(k_m) * l) / k_m))); elseif (k_m <= 4.8e+210) tmp = Float64(2.0 / Float64(Float64(t / l) * Float64(Float64(Float64(Float64(tan(k_m) * sin(k_m)) * Float64(Float64(t / l) * t)) * Float64(k_m / t)) * Float64(k_m / t)))); else tmp = Float64(2.0 / Float64(Float64(k_m / l) * Float64(Float64(Float64((sin(k_m) ^ 2.0) * t) * k_m) / l))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 14600000000.0], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 4.8e+210], N[(2.0 / N[(N[(t / l), $MachinePrecision] * N[(N[(N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k$95$m / l), $MachinePrecision] * N[(N[(N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * t), $MachinePrecision] * k$95$m), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 14600000000:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{\cos k\_m \cdot \ell}{k\_m}}}\\
\mathbf{elif}\;k\_m \leq 4.8 \cdot 10^{+210}:\\
\;\;\;\;\frac{2}{\frac{t}{\ell} \cdot \left(\left(\left(\left(\tan k\_m \cdot \sin k\_m\right) \cdot \left(\frac{t}{\ell} \cdot t\right)\right) \cdot \frac{k\_m}{t}\right) \cdot \frac{k\_m}{t}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{k\_m}{\ell} \cdot \frac{\left({\sin k\_m}^{2} \cdot t\right) \cdot k\_m}{\ell}}\\
\end{array}
\end{array}
if k < 1.46e10Initial program 38.0%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6493.0
Applied rewrites93.0%
Taylor expanded in k around 0
Applied rewrites81.8%
Applied rewrites84.7%
if 1.46e10 < k < 4.79999999999999977e210Initial program 20.1%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
associate-*l*N/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites51.3%
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-pow.f64N/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6473.6
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
lift-/.f64N/A
*-commutativeN/A
lower-*.f6477.1
Applied rewrites77.1%
if 4.79999999999999977e210 < k Initial program 37.6%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6492.3
Applied rewrites92.3%
Taylor expanded in k around 0
Applied rewrites69.1%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.25e-60)
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ (* (cos k_m) l) k_m)))
(/ 2.0 (* (/ k_m l) (* (pow (sin k_m) 2.0) (* t (/ k_m l)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.25e-60) {
tmp = 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / ((cos(k_m) * l) / k_m));
} else {
tmp = 2.0 / ((k_m / l) * (pow(sin(k_m), 2.0) * (t * (k_m / l))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.25e-60) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(Float64(cos(k_m) * l) / k_m))); else tmp = Float64(2.0 / Float64(Float64(k_m / l) * Float64((sin(k_m) ^ 2.0) * Float64(t * Float64(k_m / l))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.25e-60], N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k$95$m / l), $MachinePrecision] * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * N[(t * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.25 \cdot 10^{-60}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{\cos k\_m \cdot \ell}{k\_m}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{k\_m}{\ell} \cdot \left({\sin k\_m}^{2} \cdot \left(t \cdot \frac{k\_m}{\ell}\right)\right)}\\
\end{array}
\end{array}
if k < 1.25e-60Initial program 38.3%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6493.9
Applied rewrites93.9%
Taylor expanded in k around 0
Applied rewrites81.7%
Applied rewrites84.8%
if 1.25e-60 < k Initial program 28.8%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6489.0
Applied rewrites89.0%
Applied rewrites97.7%
Taylor expanded in k around 0
Applied rewrites63.3%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(/
2.0
(/
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(/ k_m l))
(/ (* (cos k_m) l) k_m))))k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / ((((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * (k_m / l)) / ((cos(k_m) * l) / k_m));
}
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64(Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(k_m / l)) / Float64(Float64(cos(k_m) * l) / k_m))) end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{\frac{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{k\_m}{\ell}}{\frac{\cos k\_m \cdot \ell}{k\_m}}}
\end{array}
Initial program 35.5%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6492.5
Applied rewrites92.5%
Taylor expanded in k around 0
Applied rewrites74.2%
Applied rewrites76.4%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(/
2.0
(*
(/ (/ k_m (cos k_m)) l)
(*
(/ k_m l)
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / (((k_m / cos(k_m)) / l) * ((k_m / l) * ((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m)));
}
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64(Float64(Float64(k_m / cos(k_m)) / l) * Float64(Float64(k_m / l) * Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m)))) end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(N[(N[(k$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(N[(k$95$m / l), $MachinePrecision] * N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{\frac{\frac{k\_m}{\cos k\_m}}{\ell} \cdot \left(\frac{k\_m}{\ell} \cdot \left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right)\right)}
\end{array}
Initial program 35.5%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6492.5
Applied rewrites92.5%
Taylor expanded in k around 0
Applied rewrites74.2%
Applied rewrites76.3%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 3.2e-164)
(/
2.0
(*
(*
(*
(fma
(* (fma 0.044444444444444446 (* k_m k_m) -0.3333333333333333) t)
(* k_m k_m)
t)
k_m)
k_m)
(* (/ k_m l) (/ k_m l))))
(pow (/ (* (* k_m (/ k_m l)) t) (* 2.0 (/ l (* k_m k_m)))) -1.0)))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 3.2e-164) {
tmp = 2.0 / (((fma((fma(0.044444444444444446, (k_m * k_m), -0.3333333333333333) * t), (k_m * k_m), t) * k_m) * k_m) * ((k_m / l) * (k_m / l)));
} else {
tmp = pow((((k_m * (k_m / l)) * t) / (2.0 * (l / (k_m * k_m)))), -1.0);
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 3.2e-164) tmp = Float64(2.0 / Float64(Float64(Float64(fma(Float64(fma(0.044444444444444446, Float64(k_m * k_m), -0.3333333333333333) * t), Float64(k_m * k_m), t) * k_m) * k_m) * Float64(Float64(k_m / l) * Float64(k_m / l)))); else tmp = Float64(Float64(Float64(k_m * Float64(k_m / l)) * t) / Float64(2.0 * Float64(l / Float64(k_m * k_m)))) ^ -1.0; end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 3.2e-164], N[(2.0 / N[(N[(N[(N[(N[(N[(0.044444444444444446 * N[(k$95$m * k$95$m), $MachinePrecision] + -0.3333333333333333), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(N[(k$95$m * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision] / N[(2.0 * N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 3.2 \cdot 10^{-164}:\\
\;\;\;\;\frac{2}{\left(\left(\mathsf{fma}\left(\mathsf{fma}\left(0.044444444444444446, k\_m \cdot k\_m, -0.3333333333333333\right) \cdot t, k\_m \cdot k\_m, t\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \left(\frac{k\_m}{\ell} \cdot \frac{k\_m}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{\left(k\_m \cdot \frac{k\_m}{\ell}\right) \cdot t}{2 \cdot \frac{\ell}{k\_m \cdot k\_m}}\right)}^{-1}\\
\end{array}
\end{array}
if k < 3.2e-164Initial program 35.5%
Taylor expanded in t around 0
unpow2N/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6494.6
Applied rewrites94.6%
Taylor expanded in k around 0
Applied rewrites81.1%
Taylor expanded in k around 0
Applied rewrites80.9%
Applied rewrites84.0%
if 3.2e-164 < k Initial program 35.5%
Taylor expanded in k around 0
associate-*r/N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f6464.0
Applied rewrites64.0%
Applied rewrites65.1%
Applied rewrites65.8%
Final simplification77.5%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (pow (/ (* (* k_m (/ k_m l)) t) (* 2.0 (/ l (* k_m k_m)))) -1.0))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return pow((((k_m * (k_m / l)) * t) / (2.0 * (l / (k_m * k_m)))), -1.0);
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = (((k_m * (k_m / l)) * t) / (2.0d0 * (l / (k_m * k_m)))) ** (-1.0d0)
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return Math.pow((((k_m * (k_m / l)) * t) / (2.0 * (l / (k_m * k_m)))), -1.0);
}
k_m = math.fabs(k) def code(t, l, k_m): return math.pow((((k_m * (k_m / l)) * t) / (2.0 * (l / (k_m * k_m)))), -1.0)
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(Float64(k_m * Float64(k_m / l)) * t) / Float64(2.0 * Float64(l / Float64(k_m * k_m)))) ^ -1.0 end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = (((k_m * (k_m / l)) * t) / (2.0 * (l / (k_m * k_m)))) ^ -1.0; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[Power[N[(N[(N[(k$95$m * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision] / N[(2.0 * N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
{\left(\frac{\left(k\_m \cdot \frac{k\_m}{\ell}\right) \cdot t}{2 \cdot \frac{\ell}{k\_m \cdot k\_m}}\right)}^{-1}
\end{array}
Initial program 35.5%
Taylor expanded in k around 0
associate-*r/N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f6470.9
Applied rewrites70.9%
Applied rewrites71.3%
Applied rewrites75.9%
Final simplification75.9%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* (* (/ l (* (* k_m k_m) t)) 2.0) (/ (/ l k_m) k_m)))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return ((l / ((k_m * k_m) * t)) * 2.0) * ((l / k_m) / k_m);
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = ((l / ((k_m * k_m) * t)) * 2.0d0) * ((l / k_m) / k_m)
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return ((l / ((k_m * k_m) * t)) * 2.0) * ((l / k_m) / k_m);
}
k_m = math.fabs(k) def code(t, l, k_m): return ((l / ((k_m * k_m) * t)) * 2.0) * ((l / k_m) / k_m)
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(Float64(l / Float64(Float64(k_m * k_m) * t)) * 2.0) * Float64(Float64(l / k_m) / k_m)) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = ((l / ((k_m * k_m) * t)) * 2.0) * ((l / k_m) / k_m); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(N[(l / N[(N[(k$95$m * k$95$m), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision] * N[(N[(l / k$95$m), $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\left(\frac{\ell}{\left(k\_m \cdot k\_m\right) \cdot t} \cdot 2\right) \cdot \frac{\frac{\ell}{k\_m}}{k\_m}
\end{array}
Initial program 35.5%
Taylor expanded in k around 0
associate-*r/N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f6470.9
Applied rewrites70.9%
Applied rewrites72.8%
Taylor expanded in t around 0
Applied rewrites74.8%
Applied rewrites74.9%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* (* (/ l (* (* t k_m) k_m)) 2.0) (/ l (* k_m k_m))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return ((l / ((t * k_m) * k_m)) * 2.0) * (l / (k_m * k_m));
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = ((l / ((t * k_m) * k_m)) * 2.0d0) * (l / (k_m * k_m))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return ((l / ((t * k_m) * k_m)) * 2.0) * (l / (k_m * k_m));
}
k_m = math.fabs(k) def code(t, l, k_m): return ((l / ((t * k_m) * k_m)) * 2.0) * (l / (k_m * k_m))
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(Float64(l / Float64(Float64(t * k_m) * k_m)) * 2.0) * Float64(l / Float64(k_m * k_m))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = ((l / ((t * k_m) * k_m)) * 2.0) * (l / (k_m * k_m)); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(N[(l / N[(N[(t * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision] * N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\left(\frac{\ell}{\left(t \cdot k\_m\right) \cdot k\_m} \cdot 2\right) \cdot \frac{\ell}{k\_m \cdot k\_m}
\end{array}
Initial program 35.5%
Taylor expanded in k around 0
associate-*r/N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f6470.9
Applied rewrites70.9%
Applied rewrites72.8%
Taylor expanded in t around 0
Applied rewrites74.8%
Applied rewrites74.9%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* (* (/ l (* (* k_m k_m) t)) 2.0) (/ l (* k_m k_m))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return ((l / ((k_m * k_m) * t)) * 2.0) * (l / (k_m * k_m));
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = ((l / ((k_m * k_m) * t)) * 2.0d0) * (l / (k_m * k_m))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return ((l / ((k_m * k_m) * t)) * 2.0) * (l / (k_m * k_m));
}
k_m = math.fabs(k) def code(t, l, k_m): return ((l / ((k_m * k_m) * t)) * 2.0) * (l / (k_m * k_m))
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(Float64(l / Float64(Float64(k_m * k_m) * t)) * 2.0) * Float64(l / Float64(k_m * k_m))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = ((l / ((k_m * k_m) * t)) * 2.0) * (l / (k_m * k_m)); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(N[(l / N[(N[(k$95$m * k$95$m), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision] * N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\left(\frac{\ell}{\left(k\_m \cdot k\_m\right) \cdot t} \cdot 2\right) \cdot \frac{\ell}{k\_m \cdot k\_m}
\end{array}
Initial program 35.5%
Taylor expanded in k around 0
associate-*r/N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f6470.9
Applied rewrites70.9%
Applied rewrites72.8%
Taylor expanded in t around 0
Applied rewrites74.8%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* (* l (/ 2.0 (* (* k_m k_m) t))) (/ l (* k_m k_m))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return (l * (2.0 / ((k_m * k_m) * t))) * (l / (k_m * k_m));
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = (l * (2.0d0 / ((k_m * k_m) * t))) * (l / (k_m * k_m))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return (l * (2.0 / ((k_m * k_m) * t))) * (l / (k_m * k_m));
}
k_m = math.fabs(k) def code(t, l, k_m): return (l * (2.0 / ((k_m * k_m) * t))) * (l / (k_m * k_m))
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(l * Float64(2.0 / Float64(Float64(k_m * k_m) * t))) * Float64(l / Float64(k_m * k_m))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = (l * (2.0 / ((k_m * k_m) * t))) * (l / (k_m * k_m)); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(l * N[(2.0 / N[(N[(k$95$m * k$95$m), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\left(\ell \cdot \frac{2}{\left(k\_m \cdot k\_m\right) \cdot t}\right) \cdot \frac{\ell}{k\_m \cdot k\_m}
\end{array}
Initial program 35.5%
Taylor expanded in k around 0
associate-*r/N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f6470.9
Applied rewrites70.9%
Applied rewrites72.8%
Taylor expanded in t around 0
Applied rewrites74.8%
Applied rewrites74.8%
herbie shell --seed 2024321
(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))))