
(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 20 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.15e-19)
(* (pow (/ (* t (* k_m k_m)) (* 2.0 l)) -1.0) (pow (/ (* k_m k_m) l) -1.0))
(*
(* (/ l (* k_m (fma (cos (+ k_m k_m)) -0.5 0.5))) (/ 2.0 t))
(/ (* l (cos k_m)) k_m))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.15e-19) {
tmp = pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * pow(((k_m * k_m) / l), -1.0);
} else {
tmp = ((l / (k_m * fma(cos((k_m + k_m)), -0.5, 0.5))) * (2.0 / t)) * ((l * cos(k_m)) / k_m);
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.15e-19) tmp = Float64((Float64(Float64(t * Float64(k_m * k_m)) / Float64(2.0 * l)) ^ -1.0) * (Float64(Float64(k_m * k_m) / l) ^ -1.0)); else tmp = Float64(Float64(Float64(l / Float64(k_m * fma(cos(Float64(k_m + k_m)), -0.5, 0.5))) * Float64(2.0 / t)) * Float64(Float64(l * cos(k_m)) / k_m)); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.15e-19], N[(N[Power[N[(N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(2.0 * l), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * N[Power[N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(l / N[(k$95$m * N[(N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(2.0 / t), $MachinePrecision]), $MachinePrecision] * N[(N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.15 \cdot 10^{-19}:\\
\;\;\;\;{\left(\frac{t \cdot \left(k\_m \cdot k\_m\right)}{2 \cdot \ell}\right)}^{-1} \cdot {\left(\frac{k\_m \cdot k\_m}{\ell}\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\ell}{k\_m \cdot \mathsf{fma}\left(\cos \left(k\_m + k\_m\right), -0.5, 0.5\right)} \cdot \frac{2}{t}\right) \cdot \frac{\ell \cdot \cos k\_m}{k\_m}\\
\end{array}
\end{array}
if k < 1.1499999999999999e-19Initial program 40.1%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.3
Applied rewrites68.3%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
clear-numN/A
inv-powN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
metadata-evalN/A
unpow-prod-downN/A
lower-*.f64N/A
Applied rewrites81.2%
if 1.1499999999999999e-19 < k Initial program 36.0%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6480.6
Applied rewrites80.6%
lift-cos.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
Applied rewrites93.9%
*-commutativeN/A
lift-+.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
Applied rewrites99.4%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.15e-19)
(* (pow (/ (* t (* k_m k_m)) (* 2.0 l)) -1.0) (pow (/ (* k_m k_m) l) -1.0))
(*
(/ (* l (cos k_m)) k_m)
(* (/ 2.0 (* k_m (fma (cos (+ k_m k_m)) -0.5 0.5))) (/ l t)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.15e-19) {
tmp = pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * pow(((k_m * k_m) / l), -1.0);
} else {
tmp = ((l * cos(k_m)) / k_m) * ((2.0 / (k_m * fma(cos((k_m + k_m)), -0.5, 0.5))) * (l / t));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.15e-19) tmp = Float64((Float64(Float64(t * Float64(k_m * k_m)) / Float64(2.0 * l)) ^ -1.0) * (Float64(Float64(k_m * k_m) / l) ^ -1.0)); else tmp = Float64(Float64(Float64(l * cos(k_m)) / k_m) * Float64(Float64(2.0 / Float64(k_m * fma(cos(Float64(k_m + k_m)), -0.5, 0.5))) * 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.15e-19], N[(N[Power[N[(N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(2.0 * l), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * N[Power[N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision] * N[(N[(2.0 / N[(k$95$m * N[(N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.15 \cdot 10^{-19}:\\
\;\;\;\;{\left(\frac{t \cdot \left(k\_m \cdot k\_m\right)}{2 \cdot \ell}\right)}^{-1} \cdot {\left(\frac{k\_m \cdot k\_m}{\ell}\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \cos k\_m}{k\_m} \cdot \left(\frac{2}{k\_m \cdot \mathsf{fma}\left(\cos \left(k\_m + k\_m\right), -0.5, 0.5\right)} \cdot \frac{\ell}{t}\right)\\
\end{array}
\end{array}
if k < 1.1499999999999999e-19Initial program 40.1%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.3
Applied rewrites68.3%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
clear-numN/A
inv-powN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
metadata-evalN/A
unpow-prod-downN/A
lower-*.f64N/A
Applied rewrites81.2%
if 1.1499999999999999e-19 < k Initial program 36.0%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6480.6
Applied rewrites80.6%
lift-cos.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
Applied rewrites93.9%
lift-+.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
Applied rewrites96.9%
Final simplification85.6%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.15e-19)
(* (pow (/ (* t (* k_m k_m)) (* 2.0 l)) -1.0) (pow (/ (* k_m k_m) l) -1.0))
(*
(/ (* l (cos k_m)) k_m)
(/ (* 2.0 l) (* (- 0.5 (* (cos (+ k_m k_m)) 0.5)) (* k_m t))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.15e-19) {
tmp = pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * pow(((k_m * k_m) / l), -1.0);
} else {
tmp = ((l * cos(k_m)) / k_m) * ((2.0 * l) / ((0.5 - (cos((k_m + k_m)) * 0.5)) * (k_m * t)));
}
return tmp;
}
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
real(8) :: tmp
if (k_m <= 1.15d-19) then
tmp = (((t * (k_m * k_m)) / (2.0d0 * l)) ** (-1.0d0)) * (((k_m * k_m) / l) ** (-1.0d0))
else
tmp = ((l * cos(k_m)) / k_m) * ((2.0d0 * l) / ((0.5d0 - (cos((k_m + k_m)) * 0.5d0)) * (k_m * t)))
end if
code = tmp
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.15e-19) {
tmp = Math.pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * Math.pow(((k_m * k_m) / l), -1.0);
} else {
tmp = ((l * Math.cos(k_m)) / k_m) * ((2.0 * l) / ((0.5 - (Math.cos((k_m + k_m)) * 0.5)) * (k_m * t)));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 1.15e-19: tmp = math.pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * math.pow(((k_m * k_m) / l), -1.0) else: tmp = ((l * math.cos(k_m)) / k_m) * ((2.0 * l) / ((0.5 - (math.cos((k_m + k_m)) * 0.5)) * (k_m * t))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.15e-19) tmp = Float64((Float64(Float64(t * Float64(k_m * k_m)) / Float64(2.0 * l)) ^ -1.0) * (Float64(Float64(k_m * k_m) / l) ^ -1.0)); else tmp = Float64(Float64(Float64(l * cos(k_m)) / k_m) * Float64(Float64(2.0 * l) / Float64(Float64(0.5 - Float64(cos(Float64(k_m + k_m)) * 0.5)) * Float64(k_m * t)))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 1.15e-19) tmp = (((t * (k_m * k_m)) / (2.0 * l)) ^ -1.0) * (((k_m * k_m) / l) ^ -1.0); else tmp = ((l * cos(k_m)) / k_m) * ((2.0 * l) / ((0.5 - (cos((k_m + k_m)) * 0.5)) * (k_m * t))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.15e-19], N[(N[Power[N[(N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(2.0 * l), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * N[Power[N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision] * N[(N[(2.0 * l), $MachinePrecision] / N[(N[(0.5 - N[(N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.15 \cdot 10^{-19}:\\
\;\;\;\;{\left(\frac{t \cdot \left(k\_m \cdot k\_m\right)}{2 \cdot \ell}\right)}^{-1} \cdot {\left(\frac{k\_m \cdot k\_m}{\ell}\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \cos k\_m}{k\_m} \cdot \frac{2 \cdot \ell}{\left(0.5 - \cos \left(k\_m + k\_m\right) \cdot 0.5\right) \cdot \left(k\_m \cdot t\right)}\\
\end{array}
\end{array}
if k < 1.1499999999999999e-19Initial program 40.1%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.3
Applied rewrites68.3%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
clear-numN/A
inv-powN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
metadata-evalN/A
unpow-prod-downN/A
lower-*.f64N/A
Applied rewrites81.2%
if 1.1499999999999999e-19 < k Initial program 36.0%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6480.6
Applied rewrites80.6%
lift-cos.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
Applied rewrites93.9%
Final simplification84.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.15e-19)
(* (pow (/ (* t (* k_m k_m)) (* 2.0 l)) -1.0) (pow (/ (* k_m k_m) l) -1.0))
(*
(/ (* 2.0 l) k_m)
(/ (* l (cos k_m)) (* (- 0.5 (* (cos (+ k_m k_m)) 0.5)) (* k_m t))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.15e-19) {
tmp = pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * pow(((k_m * k_m) / l), -1.0);
} else {
tmp = ((2.0 * l) / k_m) * ((l * cos(k_m)) / ((0.5 - (cos((k_m + k_m)) * 0.5)) * (k_m * t)));
}
return tmp;
}
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
real(8) :: tmp
if (k_m <= 1.15d-19) then
tmp = (((t * (k_m * k_m)) / (2.0d0 * l)) ** (-1.0d0)) * (((k_m * k_m) / l) ** (-1.0d0))
else
tmp = ((2.0d0 * l) / k_m) * ((l * cos(k_m)) / ((0.5d0 - (cos((k_m + k_m)) * 0.5d0)) * (k_m * t)))
end if
code = tmp
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.15e-19) {
tmp = Math.pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * Math.pow(((k_m * k_m) / l), -1.0);
} else {
tmp = ((2.0 * l) / k_m) * ((l * Math.cos(k_m)) / ((0.5 - (Math.cos((k_m + k_m)) * 0.5)) * (k_m * t)));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 1.15e-19: tmp = math.pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * math.pow(((k_m * k_m) / l), -1.0) else: tmp = ((2.0 * l) / k_m) * ((l * math.cos(k_m)) / ((0.5 - (math.cos((k_m + k_m)) * 0.5)) * (k_m * t))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.15e-19) tmp = Float64((Float64(Float64(t * Float64(k_m * k_m)) / Float64(2.0 * l)) ^ -1.0) * (Float64(Float64(k_m * k_m) / l) ^ -1.0)); else tmp = Float64(Float64(Float64(2.0 * l) / k_m) * Float64(Float64(l * cos(k_m)) / Float64(Float64(0.5 - Float64(cos(Float64(k_m + k_m)) * 0.5)) * Float64(k_m * t)))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 1.15e-19) tmp = (((t * (k_m * k_m)) / (2.0 * l)) ^ -1.0) * (((k_m * k_m) / l) ^ -1.0); else tmp = ((2.0 * l) / k_m) * ((l * cos(k_m)) / ((0.5 - (cos((k_m + k_m)) * 0.5)) * (k_m * t))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.15e-19], N[(N[Power[N[(N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(2.0 * l), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * N[Power[N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 * l), $MachinePrecision] / k$95$m), $MachinePrecision] * N[(N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(N[(0.5 - N[(N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.15 \cdot 10^{-19}:\\
\;\;\;\;{\left(\frac{t \cdot \left(k\_m \cdot k\_m\right)}{2 \cdot \ell}\right)}^{-1} \cdot {\left(\frac{k\_m \cdot k\_m}{\ell}\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \ell}{k\_m} \cdot \frac{\ell \cdot \cos k\_m}{\left(0.5 - \cos \left(k\_m + k\_m\right) \cdot 0.5\right) \cdot \left(k\_m \cdot t\right)}\\
\end{array}
\end{array}
if k < 1.1499999999999999e-19Initial program 40.1%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.3
Applied rewrites68.3%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
clear-numN/A
inv-powN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
metadata-evalN/A
unpow-prod-downN/A
lower-*.f64N/A
Applied rewrites81.2%
if 1.1499999999999999e-19 < k Initial program 36.0%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6480.6
Applied rewrites80.6%
lift-cos.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f6494.0
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites93.9%
Final simplification84.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.15e-19)
(* (pow (/ (* t (* k_m k_m)) (* 2.0 l)) -1.0) (pow (/ (* k_m k_m) l) -1.0))
(*
(/ (* l (cos k_m)) k_m)
(* l (/ 2.0 (* (fma (cos (+ k_m k_m)) -0.5 0.5) (* k_m t)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.15e-19) {
tmp = pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * pow(((k_m * k_m) / l), -1.0);
} else {
tmp = ((l * cos(k_m)) / k_m) * (l * (2.0 / (fma(cos((k_m + k_m)), -0.5, 0.5) * (k_m * t))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.15e-19) tmp = Float64((Float64(Float64(t * Float64(k_m * k_m)) / Float64(2.0 * l)) ^ -1.0) * (Float64(Float64(k_m * k_m) / l) ^ -1.0)); else tmp = Float64(Float64(Float64(l * cos(k_m)) / k_m) * Float64(l * Float64(2.0 / Float64(fma(cos(Float64(k_m + k_m)), -0.5, 0.5) * Float64(k_m * t))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.15e-19], N[(N[Power[N[(N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(2.0 * l), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * N[Power[N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision] * N[(l * N[(2.0 / N[(N[(N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.15 \cdot 10^{-19}:\\
\;\;\;\;{\left(\frac{t \cdot \left(k\_m \cdot k\_m\right)}{2 \cdot \ell}\right)}^{-1} \cdot {\left(\frac{k\_m \cdot k\_m}{\ell}\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \cos k\_m}{k\_m} \cdot \left(\ell \cdot \frac{2}{\mathsf{fma}\left(\cos \left(k\_m + k\_m\right), -0.5, 0.5\right) \cdot \left(k\_m \cdot t\right)}\right)\\
\end{array}
\end{array}
if k < 1.1499999999999999e-19Initial program 40.1%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.3
Applied rewrites68.3%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
clear-numN/A
inv-powN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
metadata-evalN/A
unpow-prod-downN/A
lower-*.f64N/A
Applied rewrites81.2%
if 1.1499999999999999e-19 < k Initial program 36.0%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6480.6
Applied rewrites80.6%
lift-cos.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
Applied rewrites93.9%
Applied rewrites93.9%
Final simplification84.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.15e-19)
(* (pow (/ (* t (* k_m k_m)) (* 2.0 l)) -1.0) (pow (/ (* k_m k_m) l) -1.0))
(*
l
(*
(/ (cos k_m) k_m)
(/ (* 2.0 l) (* (fma (cos (+ k_m k_m)) -0.5 0.5) (* k_m t)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.15e-19) {
tmp = pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * pow(((k_m * k_m) / l), -1.0);
} else {
tmp = l * ((cos(k_m) / k_m) * ((2.0 * l) / (fma(cos((k_m + k_m)), -0.5, 0.5) * (k_m * t))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.15e-19) tmp = Float64((Float64(Float64(t * Float64(k_m * k_m)) / Float64(2.0 * l)) ^ -1.0) * (Float64(Float64(k_m * k_m) / l) ^ -1.0)); else tmp = Float64(l * Float64(Float64(cos(k_m) / k_m) * Float64(Float64(2.0 * l) / Float64(fma(cos(Float64(k_m + k_m)), -0.5, 0.5) * Float64(k_m * t))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.15e-19], N[(N[Power[N[(N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(2.0 * l), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * N[Power[N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision], N[(l * N[(N[(N[Cos[k$95$m], $MachinePrecision] / k$95$m), $MachinePrecision] * N[(N[(2.0 * l), $MachinePrecision] / N[(N[(N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.15 \cdot 10^{-19}:\\
\;\;\;\;{\left(\frac{t \cdot \left(k\_m \cdot k\_m\right)}{2 \cdot \ell}\right)}^{-1} \cdot {\left(\frac{k\_m \cdot k\_m}{\ell}\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \left(\frac{\cos k\_m}{k\_m} \cdot \frac{2 \cdot \ell}{\mathsf{fma}\left(\cos \left(k\_m + k\_m\right), -0.5, 0.5\right) \cdot \left(k\_m \cdot t\right)}\right)\\
\end{array}
\end{array}
if k < 1.1499999999999999e-19Initial program 40.1%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.3
Applied rewrites68.3%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
clear-numN/A
inv-powN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
metadata-evalN/A
unpow-prod-downN/A
lower-*.f64N/A
Applied rewrites81.2%
if 1.1499999999999999e-19 < k Initial program 36.0%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6480.6
Applied rewrites80.6%
lift-cos.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
Applied rewrites93.9%
lift-+.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift-/.f64N/A
*-commutativeN/A
Applied rewrites90.1%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.15e-19)
(* (pow (/ (* t (* k_m k_m)) (* 2.0 l)) -1.0) (pow (/ (* k_m k_m) l) -1.0))
(*
(* 2.0 l)
(/
(* l (cos k_m))
(* k_m (* (fma (cos (+ k_m k_m)) -0.5 0.5) (* k_m t)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.15e-19) {
tmp = pow(((t * (k_m * k_m)) / (2.0 * l)), -1.0) * pow(((k_m * k_m) / l), -1.0);
} else {
tmp = (2.0 * l) * ((l * cos(k_m)) / (k_m * (fma(cos((k_m + k_m)), -0.5, 0.5) * (k_m * t))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.15e-19) tmp = Float64((Float64(Float64(t * Float64(k_m * k_m)) / Float64(2.0 * l)) ^ -1.0) * (Float64(Float64(k_m * k_m) / l) ^ -1.0)); else tmp = Float64(Float64(2.0 * l) * Float64(Float64(l * cos(k_m)) / Float64(k_m * Float64(fma(cos(Float64(k_m + k_m)), -0.5, 0.5) * Float64(k_m * t))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.15e-19], N[(N[Power[N[(N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(2.0 * l), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * N[Power[N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * l), $MachinePrecision] * N[(N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(k$95$m * N[(N[(N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.15 \cdot 10^{-19}:\\
\;\;\;\;{\left(\frac{t \cdot \left(k\_m \cdot k\_m\right)}{2 \cdot \ell}\right)}^{-1} \cdot {\left(\frac{k\_m \cdot k\_m}{\ell}\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \ell\right) \cdot \frac{\ell \cdot \cos k\_m}{k\_m \cdot \left(\mathsf{fma}\left(\cos \left(k\_m + k\_m\right), -0.5, 0.5\right) \cdot \left(k\_m \cdot t\right)\right)}\\
\end{array}
\end{array}
if k < 1.1499999999999999e-19Initial program 40.1%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.3
Applied rewrites68.3%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
clear-numN/A
inv-powN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
metadata-evalN/A
unpow-prod-downN/A
lower-*.f64N/A
Applied rewrites81.2%
if 1.1499999999999999e-19 < k Initial program 36.0%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6480.6
Applied rewrites80.6%
lift-cos.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
Applied rewrites93.9%
Applied rewrites84.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (/ l (* k_m k_m))))
(if (<= k_m 2.1e-73)
(* (/ 2.0 t) (* t_1 t_1))
(if (<= k_m 2.8e+78)
(*
(/ (* l (cos k_m)) k_m)
(/
(* (/ l t) (fma 0.6666666666666666 (* k_m k_m) 2.0))
(* k_m (* k_m k_m))))
(*
(/ (* 2.0 l) (* (fma (cos (+ k_m k_m)) -0.5 0.5) (* k_m t)))
(/ l k_m))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = l / (k_m * k_m);
double tmp;
if (k_m <= 2.1e-73) {
tmp = (2.0 / t) * (t_1 * t_1);
} else if (k_m <= 2.8e+78) {
tmp = ((l * cos(k_m)) / k_m) * (((l / t) * fma(0.6666666666666666, (k_m * k_m), 2.0)) / (k_m * (k_m * k_m)));
} else {
tmp = ((2.0 * l) / (fma(cos((k_m + k_m)), -0.5, 0.5) * (k_m * t))) * (l / k_m);
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(l / Float64(k_m * k_m)) tmp = 0.0 if (k_m <= 2.1e-73) tmp = Float64(Float64(2.0 / t) * Float64(t_1 * t_1)); elseif (k_m <= 2.8e+78) tmp = Float64(Float64(Float64(l * cos(k_m)) / k_m) * Float64(Float64(Float64(l / t) * fma(0.6666666666666666, Float64(k_m * k_m), 2.0)) / Float64(k_m * Float64(k_m * k_m)))); else tmp = Float64(Float64(Float64(2.0 * l) / Float64(fma(cos(Float64(k_m + k_m)), -0.5, 0.5) * Float64(k_m * t))) * Float64(l / k_m)); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k$95$m, 2.1e-73], N[(N[(2.0 / t), $MachinePrecision] * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.8e+78], N[(N[(N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision] * N[(N[(N[(l / t), $MachinePrecision] * N[(0.6666666666666666 * N[(k$95$m * k$95$m), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 * l), $MachinePrecision] / N[(N[(N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / k$95$m), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := \frac{\ell}{k\_m \cdot k\_m}\\
\mathbf{if}\;k\_m \leq 2.1 \cdot 10^{-73}:\\
\;\;\;\;\frac{2}{t} \cdot \left(t\_1 \cdot t\_1\right)\\
\mathbf{elif}\;k\_m \leq 2.8 \cdot 10^{+78}:\\
\;\;\;\;\frac{\ell \cdot \cos k\_m}{k\_m} \cdot \frac{\frac{\ell}{t} \cdot \mathsf{fma}\left(0.6666666666666666, k\_m \cdot k\_m, 2\right)}{k\_m \cdot \left(k\_m \cdot k\_m\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \ell}{\mathsf{fma}\left(\cos \left(k\_m + k\_m\right), -0.5, 0.5\right) \cdot \left(k\_m \cdot t\right)} \cdot \frac{\ell}{k\_m}\\
\end{array}
\end{array}
if k < 2.0999999999999999e-73Initial program 39.5%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6467.2
Applied rewrites67.2%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
cube-unmultN/A
lower-*.f64N/A
cube-unmultN/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6466.7
Applied rewrites66.7%
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6479.1
Applied rewrites79.1%
if 2.0999999999999999e-73 < k < 2.8000000000000001e78Initial program 40.9%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6488.2
Applied rewrites88.2%
lift-cos.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
Applied rewrites73.2%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
lower-/.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lower-/.f64N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
Applied rewrites80.5%
if 2.8000000000000001e78 < k Initial program 36.1%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6475.9
Applied rewrites75.9%
Taylor expanded in k around 0
lower-*.f64N/A
unpow2N/A
lower-*.f6463.4
Applied rewrites63.4%
Applied rewrites66.6%
Final simplification76.7%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 30.0) (* (/ l (* k_m k_m)) (/ (* 2.0 l) (* t (* k_m k_m)))) (* (/ (* 2.0 l) (* (fma (cos (+ k_m k_m)) -0.5 0.5) (* k_m t))) (/ l k_m))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 30.0) {
tmp = (l / (k_m * k_m)) * ((2.0 * l) / (t * (k_m * k_m)));
} else {
tmp = ((2.0 * l) / (fma(cos((k_m + k_m)), -0.5, 0.5) * (k_m * t))) * (l / k_m);
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 30.0) tmp = Float64(Float64(l / Float64(k_m * k_m)) * Float64(Float64(2.0 * l) / Float64(t * Float64(k_m * k_m)))); else tmp = Float64(Float64(Float64(2.0 * l) / Float64(fma(cos(Float64(k_m + k_m)), -0.5, 0.5) * Float64(k_m * t))) * Float64(l / k_m)); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 30.0], N[(N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * l), $MachinePrecision] / N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 * l), $MachinePrecision] / N[(N[(N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / k$95$m), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 30:\\
\;\;\;\;\frac{\ell}{k\_m \cdot k\_m} \cdot \frac{2 \cdot \ell}{t \cdot \left(k\_m \cdot k\_m\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \ell}{\mathsf{fma}\left(\cos \left(k\_m + k\_m\right), -0.5, 0.5\right) \cdot \left(k\_m \cdot t\right)} \cdot \frac{\ell}{k\_m}\\
\end{array}
\end{array}
if k < 30Initial program 40.0%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6467.7
Applied rewrites67.7%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6480.5
Applied rewrites80.5%
if 30 < k Initial program 36.1%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6479.8
Applied rewrites79.8%
Taylor expanded in k around 0
lower-*.f64N/A
unpow2N/A
lower-*.f6462.8
Applied rewrites62.8%
Applied rewrites65.6%
Final simplification76.5%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 30.0) (* (/ l (* k_m k_m)) (/ (* 2.0 l) (* t (* k_m k_m)))) (* (* 2.0 l) (/ l (* k_m (* (fma (cos (+ k_m k_m)) -0.5 0.5) (* k_m t)))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 30.0) {
tmp = (l / (k_m * k_m)) * ((2.0 * l) / (t * (k_m * k_m)));
} else {
tmp = (2.0 * l) * (l / (k_m * (fma(cos((k_m + k_m)), -0.5, 0.5) * (k_m * t))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 30.0) tmp = Float64(Float64(l / Float64(k_m * k_m)) * Float64(Float64(2.0 * l) / Float64(t * Float64(k_m * k_m)))); else tmp = Float64(Float64(2.0 * l) * Float64(l / Float64(k_m * Float64(fma(cos(Float64(k_m + k_m)), -0.5, 0.5) * Float64(k_m * t))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 30.0], N[(N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * l), $MachinePrecision] / N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * l), $MachinePrecision] * N[(l / N[(k$95$m * N[(N[(N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision] * -0.5 + 0.5), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 30:\\
\;\;\;\;\frac{\ell}{k\_m \cdot k\_m} \cdot \frac{2 \cdot \ell}{t \cdot \left(k\_m \cdot k\_m\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \ell\right) \cdot \frac{\ell}{k\_m \cdot \left(\mathsf{fma}\left(\cos \left(k\_m + k\_m\right), -0.5, 0.5\right) \cdot \left(k\_m \cdot t\right)\right)}\\
\end{array}
\end{array}
if k < 30Initial program 40.0%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6467.7
Applied rewrites67.7%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6480.5
Applied rewrites80.5%
if 30 < k Initial program 36.1%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f6479.8
Applied rewrites79.8%
Taylor expanded in k around 0
lower-*.f64N/A
unpow2N/A
lower-*.f6462.8
Applied rewrites62.8%
associate-*r*N/A
lift-*.f64N/A
lift-sin.f64N/A
lift-pow.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6464.8
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-pow.f64N/A
unpow2N/A
lift-sin.f64N/A
lift-sin.f64N/A
sqr-sin-aN/A
count-2N/A
lift-+.f64N/A
Applied rewrites64.8%
Final simplification76.3%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= l 1e-148) (* (* 2.0 l) (/ l (* k_m (* t (* k_m (* k_m k_m)))))) (* (/ 2.0 (* t (* k_m k_m))) (/ (* l l) (* k_m k_m)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (l <= 1e-148) {
tmp = (2.0 * l) * (l / (k_m * (t * (k_m * (k_m * k_m)))));
} else {
tmp = (2.0 / (t * (k_m * k_m))) * ((l * l) / (k_m * k_m));
}
return tmp;
}
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
real(8) :: tmp
if (l <= 1d-148) then
tmp = (2.0d0 * l) * (l / (k_m * (t * (k_m * (k_m * k_m)))))
else
tmp = (2.0d0 / (t * (k_m * k_m))) * ((l * l) / (k_m * k_m))
end if
code = tmp
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (l <= 1e-148) {
tmp = (2.0 * l) * (l / (k_m * (t * (k_m * (k_m * k_m)))));
} else {
tmp = (2.0 / (t * (k_m * k_m))) * ((l * l) / (k_m * k_m));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if l <= 1e-148: tmp = (2.0 * l) * (l / (k_m * (t * (k_m * (k_m * k_m))))) else: tmp = (2.0 / (t * (k_m * k_m))) * ((l * l) / (k_m * k_m)) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (l <= 1e-148) tmp = Float64(Float64(2.0 * l) * Float64(l / Float64(k_m * Float64(t * Float64(k_m * Float64(k_m * k_m)))))); else tmp = Float64(Float64(2.0 / Float64(t * Float64(k_m * k_m))) * Float64(Float64(l * l) / Float64(k_m * k_m))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (l <= 1e-148) tmp = (2.0 * l) * (l / (k_m * (t * (k_m * (k_m * k_m))))); else tmp = (2.0 / (t * (k_m * k_m))) * ((l * l) / (k_m * k_m)); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[l, 1e-148], N[(N[(2.0 * l), $MachinePrecision] * N[(l / N[(k$95$m * N[(t * N[(k$95$m * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;\ell \leq 10^{-148}:\\
\;\;\;\;\left(2 \cdot \ell\right) \cdot \frac{\ell}{k\_m \cdot \left(t \cdot \left(k\_m \cdot \left(k\_m \cdot k\_m\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t \cdot \left(k\_m \cdot k\_m\right)} \cdot \frac{\ell \cdot \ell}{k\_m \cdot k\_m}\\
\end{array}
\end{array}
if l < 9.99999999999999936e-149Initial program 38.9%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6465.1
Applied rewrites65.1%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6473.2
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
cube-unmultN/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
cube-unmultN/A
lift-*.f64N/A
lower-*.f6475.2
Applied rewrites75.2%
if 9.99999999999999936e-149 < l Initial program 39.0%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6465.7
Applied rewrites65.7%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f6470.3
Applied rewrites70.3%
Final simplification73.2%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* (* (/ 2.0 t) (/ l k_m)) (/ l (* k_m (* k_m k_m)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return ((2.0 / t) * (l / k_m)) * (l / (k_m * (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 = ((2.0d0 / t) * (l / k_m)) * (l / (k_m * (k_m * k_m)))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return ((2.0 / t) * (l / k_m)) * (l / (k_m * (k_m * k_m)));
}
k_m = math.fabs(k) def code(t, l, k_m): return ((2.0 / t) * (l / k_m)) * (l / (k_m * (k_m * k_m)))
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(Float64(2.0 / t) * Float64(l / k_m)) * Float64(l / Float64(k_m * Float64(k_m * k_m)))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = ((2.0 / t) * (l / k_m)) * (l / (k_m * (k_m * k_m))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(N[(2.0 / t), $MachinePrecision] * N[(l / k$95$m), $MachinePrecision]), $MachinePrecision] * N[(l / N[(k$95$m * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\left(\frac{2}{t} \cdot \frac{\ell}{k\_m}\right) \cdot \frac{\ell}{k\_m \cdot \left(k\_m \cdot k\_m\right)}
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6465.4
Applied rewrites65.4%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
times-fracN/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
cube-unmultN/A
lower-*.f64N/A
cube-unmultN/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6465.1
Applied rewrites65.1%
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-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6474.9
Applied rewrites74.9%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* (/ l (* k_m k_m)) (/ (* 2.0 l) (* t (* k_m k_m)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return (l / (k_m * k_m)) * ((2.0 * l) / (t * (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)) * ((2.0d0 * l) / (t * (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)) * ((2.0 * l) / (t * (k_m * k_m)));
}
k_m = math.fabs(k) def code(t, l, k_m): return (l / (k_m * k_m)) * ((2.0 * l) / (t * (k_m * k_m)))
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(l / Float64(k_m * k_m)) * Float64(Float64(2.0 * l) / Float64(t * Float64(k_m * k_m)))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = (l / (k_m * k_m)) * ((2.0 * l) / (t * (k_m * k_m))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * l), $MachinePrecision] / N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{\ell}{k\_m \cdot k\_m} \cdot \frac{2 \cdot \ell}{t \cdot \left(k\_m \cdot k\_m\right)}
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6465.4
Applied rewrites65.4%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6474.9
Applied rewrites74.9%
Final simplification74.9%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* (* 2.0 l) (/ l (* k_m (* t (* k_m (* k_m k_m)))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return (2.0 * l) * (l / (k_m * (t * (k_m * (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 = (2.0d0 * l) * (l / (k_m * (t * (k_m * (k_m * k_m)))))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return (2.0 * l) * (l / (k_m * (t * (k_m * (k_m * k_m)))));
}
k_m = math.fabs(k) def code(t, l, k_m): return (2.0 * l) * (l / (k_m * (t * (k_m * (k_m * k_m)))))
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(2.0 * l) * Float64(l / Float64(k_m * Float64(t * Float64(k_m * Float64(k_m * k_m)))))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = (2.0 * l) * (l / (k_m * (t * (k_m * (k_m * k_m))))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(2.0 * l), $MachinePrecision] * N[(l / N[(k$95$m * N[(t * N[(k$95$m * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\left(2 \cdot \ell\right) \cdot \frac{\ell}{k\_m \cdot \left(t \cdot \left(k\_m \cdot \left(k\_m \cdot k\_m\right)\right)\right)}
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-evalN/A
pow-sqrN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6465.4
Applied rewrites65.4%
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6470.3
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
cube-unmultN/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
cube-unmultN/A
lift-*.f64N/A
lower-*.f6472.3
Applied rewrites72.3%
Final simplification72.3%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* k_m (/ k_m (/ t (* (* l l) -0.0205026455026455)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return k_m * (k_m / (t / ((l * l) * -0.0205026455026455)));
}
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 / (t / ((l * l) * (-0.0205026455026455d0))))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return k_m * (k_m / (t / ((l * l) * -0.0205026455026455)));
}
k_m = math.fabs(k) def code(t, l, k_m): return k_m * (k_m / (t / ((l * l) * -0.0205026455026455)))
k_m = abs(k) function code(t, l, k_m) return Float64(k_m * Float64(k_m / Float64(t / Float64(Float64(l * l) * -0.0205026455026455)))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = k_m * (k_m / (t / ((l * l) * -0.0205026455026455))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(k$95$m * N[(k$95$m / N[(t / N[(N[(l * l), $MachinePrecision] * -0.0205026455026455), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
k\_m \cdot \frac{k\_m}{\frac{t}{\left(\ell \cdot \ell\right) \cdot -0.0205026455026455}}
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
Applied rewrites23.0%
Taylor expanded in k around inf
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6420.2
Applied rewrites20.2%
lift-*.f64N/A
lift-*.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6420.7
Applied rewrites20.7%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* k_m (* k_m (/ (* (* l l) -0.0205026455026455) t))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return k_m * (k_m * (((l * l) * -0.0205026455026455) / t));
}
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 * l) * (-0.0205026455026455d0)) / t))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return k_m * (k_m * (((l * l) * -0.0205026455026455) / t));
}
k_m = math.fabs(k) def code(t, l, k_m): return k_m * (k_m * (((l * l) * -0.0205026455026455) / t))
k_m = abs(k) function code(t, l, k_m) return Float64(k_m * Float64(k_m * Float64(Float64(Float64(l * l) * -0.0205026455026455) / t))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = k_m * (k_m * (((l * l) * -0.0205026455026455) / t)); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(k$95$m * N[(k$95$m * N[(N[(N[(l * l), $MachinePrecision] * -0.0205026455026455), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
k\_m \cdot \left(k\_m \cdot \frac{\left(\ell \cdot \ell\right) \cdot -0.0205026455026455}{t}\right)
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
Applied rewrites23.0%
Taylor expanded in k around inf
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6420.2
Applied rewrites20.2%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* k_m (* -0.0205026455026455 (/ (* k_m (* l l)) t))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return k_m * (-0.0205026455026455 * ((k_m * (l * l)) / t));
}
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 * ((-0.0205026455026455d0) * ((k_m * (l * l)) / t))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return k_m * (-0.0205026455026455 * ((k_m * (l * l)) / t));
}
k_m = math.fabs(k) def code(t, l, k_m): return k_m * (-0.0205026455026455 * ((k_m * (l * l)) / t))
k_m = abs(k) function code(t, l, k_m) return Float64(k_m * Float64(-0.0205026455026455 * Float64(Float64(k_m * Float64(l * l)) / t))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = k_m * (-0.0205026455026455 * ((k_m * (l * l)) / t)); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(k$95$m * N[(-0.0205026455026455 * N[(N[(k$95$m * N[(l * l), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
k\_m \cdot \left(-0.0205026455026455 \cdot \frac{k\_m \cdot \left(\ell \cdot \ell\right)}{t}\right)
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
Applied rewrites23.0%
Taylor expanded in k around inf
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6420.2
Applied rewrites20.2%
Taylor expanded in k around 0
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6419.1
Applied rewrites19.1%
Final simplification19.1%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* k_m (* -0.0205026455026455 (* l (/ (* k_m l) t)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return k_m * (-0.0205026455026455 * (l * ((k_m * l) / t)));
}
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 * ((-0.0205026455026455d0) * (l * ((k_m * l) / t)))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return k_m * (-0.0205026455026455 * (l * ((k_m * l) / t)));
}
k_m = math.fabs(k) def code(t, l, k_m): return k_m * (-0.0205026455026455 * (l * ((k_m * l) / t)))
k_m = abs(k) function code(t, l, k_m) return Float64(k_m * Float64(-0.0205026455026455 * Float64(l * Float64(Float64(k_m * l) / t)))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = k_m * (-0.0205026455026455 * (l * ((k_m * l) / t))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(k$95$m * N[(-0.0205026455026455 * N[(l * N[(N[(k$95$m * l), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
k\_m \cdot \left(-0.0205026455026455 \cdot \left(\ell \cdot \frac{k\_m \cdot \ell}{t}\right)\right)
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
Applied rewrites23.0%
Taylor expanded in k around inf
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6420.2
Applied rewrites20.2%
Taylor expanded in k around 0
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6419.1
Applied rewrites19.1%
associate-*l*N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f6410.8
Applied rewrites10.8%
Final simplification10.8%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* k_m (* -0.0205026455026455 (* l (* l (/ k_m t))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return k_m * (-0.0205026455026455 * (l * (l * (k_m / t))));
}
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 * ((-0.0205026455026455d0) * (l * (l * (k_m / t))))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return k_m * (-0.0205026455026455 * (l * (l * (k_m / t))));
}
k_m = math.fabs(k) def code(t, l, k_m): return k_m * (-0.0205026455026455 * (l * (l * (k_m / t))))
k_m = abs(k) function code(t, l, k_m) return Float64(k_m * Float64(-0.0205026455026455 * Float64(l * Float64(l * Float64(k_m / t))))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = k_m * (-0.0205026455026455 * (l * (l * (k_m / t)))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(k$95$m * N[(-0.0205026455026455 * N[(l * N[(l * N[(k$95$m / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
k\_m \cdot \left(-0.0205026455026455 \cdot \left(\ell \cdot \left(\ell \cdot \frac{k\_m}{t}\right)\right)\right)
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
Applied rewrites23.0%
Taylor expanded in k around inf
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6420.2
Applied rewrites20.2%
Taylor expanded in k around 0
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6419.1
Applied rewrites19.1%
lift-*.f64N/A
associate-/l*N/A
lift-*.f64N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6410.8
Applied rewrites10.8%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* -0.0205026455026455 (* l (/ (* (* k_m k_m) l) t))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return -0.0205026455026455 * (l * (((k_m * k_m) * l) / t));
}
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 = (-0.0205026455026455d0) * (l * (((k_m * k_m) * l) / t))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return -0.0205026455026455 * (l * (((k_m * k_m) * l) / t));
}
k_m = math.fabs(k) def code(t, l, k_m): return -0.0205026455026455 * (l * (((k_m * k_m) * l) / t))
k_m = abs(k) function code(t, l, k_m) return Float64(-0.0205026455026455 * Float64(l * Float64(Float64(Float64(k_m * k_m) * l) / t))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = -0.0205026455026455 * (l * (((k_m * k_m) * l) / t)); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(-0.0205026455026455 * N[(l * N[(N[(N[(k$95$m * k$95$m), $MachinePrecision] * l), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
-0.0205026455026455 \cdot \left(\ell \cdot \frac{\left(k\_m \cdot k\_m\right) \cdot \ell}{t}\right)
\end{array}
Initial program 39.0%
Taylor expanded in k around 0
Applied rewrites23.0%
Taylor expanded in k around inf
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6420.2
Applied rewrites20.2%
Taylor expanded in k around 0
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6419.1
Applied rewrites19.1%
Taylor expanded in k around 0
lower-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f647.0
Applied rewrites7.0%
Final simplification7.0%
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))))