
(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 17 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 6.2e-5)
(* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m)))
(*
(/ (/ (cos k_m) (/ k_m l)) t)
(/ l (* k_m (+ 0.25 (* (cos (* k_m 2.0)) -0.25)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 6.2e-5) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = ((cos(k_m) / (k_m / l)) / t) * (l / (k_m * (0.25 + (cos((k_m * 2.0)) * -0.25))));
}
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 <= 6.2d-5) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else
tmp = ((cos(k_m) / (k_m / l)) / t) * (l / (k_m * (0.25d0 + (cos((k_m * 2.0d0)) * (-0.25d0)))))
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 <= 6.2e-5) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = ((Math.cos(k_m) / (k_m / l)) / t) * (l / (k_m * (0.25 + (Math.cos((k_m * 2.0)) * -0.25))));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 6.2e-5: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) else: tmp = ((math.cos(k_m) / (k_m / l)) / t) * (l / (k_m * (0.25 + (math.cos((k_m * 2.0)) * -0.25)))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 6.2e-5) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); else tmp = Float64(Float64(Float64(cos(k_m) / Float64(k_m / l)) / t) * Float64(l / Float64(k_m * Float64(0.25 + Float64(cos(Float64(k_m * 2.0)) * -0.25))))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 6.2e-5) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); else tmp = ((cos(k_m) / (k_m / l)) / t) * (l / (k_m * (0.25 + (cos((k_m * 2.0)) * -0.25)))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 6.2e-5], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] * N[(l / N[(k$95$m * N[(0.25 + N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 6.2 \cdot 10^{-5}:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\cos k\_m}{\frac{k\_m}{\ell}}}{t} \cdot \frac{\ell}{k\_m \cdot \left(0.25 + \cos \left(k\_m \cdot 2\right) \cdot -0.25\right)}\\
\end{array}
\end{array}
if k < 6.20000000000000027e-5Initial program 37.1%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6470.6%
Simplified70.6%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.7%
Applied egg-rr70.7%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6478.5%
Applied egg-rr78.5%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6481.6%
Applied egg-rr81.6%
if 6.20000000000000027e-5 < k Initial program 24.0%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6471.9%
Simplified71.9%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr75.3%
associate-/l/N/A
associate-/r*N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
Applied egg-rr82.8%
associate-/r/N/A
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f6498.9%
Applied egg-rr98.9%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 6.2e-5)
(* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m)))
(*
(/ (cos k_m) (/ k_m l))
(/ (/ l k_m) (* t (+ 0.25 (* (cos (* k_m 2.0)) -0.25)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 6.2e-5) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = (cos(k_m) / (k_m / l)) * ((l / k_m) / (t * (0.25 + (cos((k_m * 2.0)) * -0.25))));
}
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 <= 6.2d-5) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else
tmp = (cos(k_m) / (k_m / l)) * ((l / k_m) / (t * (0.25d0 + (cos((k_m * 2.0d0)) * (-0.25d0)))))
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 <= 6.2e-5) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = (Math.cos(k_m) / (k_m / l)) * ((l / k_m) / (t * (0.25 + (Math.cos((k_m * 2.0)) * -0.25))));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 6.2e-5: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) else: tmp = (math.cos(k_m) / (k_m / l)) * ((l / k_m) / (t * (0.25 + (math.cos((k_m * 2.0)) * -0.25)))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 6.2e-5) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); else tmp = Float64(Float64(cos(k_m) / Float64(k_m / l)) * Float64(Float64(l / k_m) / Float64(t * Float64(0.25 + Float64(cos(Float64(k_m * 2.0)) * -0.25))))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 6.2e-5) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); else tmp = (cos(k_m) / (k_m / l)) * ((l / k_m) / (t * (0.25 + (cos((k_m * 2.0)) * -0.25)))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 6.2e-5], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] * N[(N[(l / k$95$m), $MachinePrecision] / N[(t * N[(0.25 + N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 6.2 \cdot 10^{-5}:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{else}:\\
\;\;\;\;\frac{\cos k\_m}{\frac{k\_m}{\ell}} \cdot \frac{\frac{\ell}{k\_m}}{t \cdot \left(0.25 + \cos \left(k\_m \cdot 2\right) \cdot -0.25\right)}\\
\end{array}
\end{array}
if k < 6.20000000000000027e-5Initial program 37.1%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6470.6%
Simplified70.6%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.7%
Applied egg-rr70.7%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6478.5%
Applied egg-rr78.5%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6481.6%
Applied egg-rr81.6%
if 6.20000000000000027e-5 < k Initial program 24.0%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6471.9%
Simplified71.9%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr75.3%
associate-/l/N/A
associate-/r*N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
Applied egg-rr82.8%
div-invN/A
associate-*l*N/A
times-fracN/A
*-lowering-*.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
/-lowering-/.f64N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f6498.9%
Applied egg-rr98.9%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 6.2e-5)
(* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m)))
(*
(/ l k_m)
(/ (/ (* l (cos k_m)) t) (* k_m (+ 0.25 (* (cos (* k_m 2.0)) -0.25)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 6.2e-5) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = (l / k_m) * (((l * cos(k_m)) / t) / (k_m * (0.25 + (cos((k_m * 2.0)) * -0.25))));
}
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 <= 6.2d-5) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else
tmp = (l / k_m) * (((l * cos(k_m)) / t) / (k_m * (0.25d0 + (cos((k_m * 2.0d0)) * (-0.25d0)))))
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 <= 6.2e-5) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = (l / k_m) * (((l * Math.cos(k_m)) / t) / (k_m * (0.25 + (Math.cos((k_m * 2.0)) * -0.25))));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 6.2e-5: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) else: tmp = (l / k_m) * (((l * math.cos(k_m)) / t) / (k_m * (0.25 + (math.cos((k_m * 2.0)) * -0.25)))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 6.2e-5) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); else tmp = Float64(Float64(l / k_m) * Float64(Float64(Float64(l * cos(k_m)) / t) / Float64(k_m * Float64(0.25 + Float64(cos(Float64(k_m * 2.0)) * -0.25))))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 6.2e-5) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); else tmp = (l / k_m) * (((l * cos(k_m)) / t) / (k_m * (0.25 + (cos((k_m * 2.0)) * -0.25)))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 6.2e-5], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k$95$m), $MachinePrecision] * N[(N[(N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] / N[(k$95$m * N[(0.25 + N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 6.2 \cdot 10^{-5}:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k\_m} \cdot \frac{\frac{\ell \cdot \cos k\_m}{t}}{k\_m \cdot \left(0.25 + \cos \left(k\_m \cdot 2\right) \cdot -0.25\right)}\\
\end{array}
\end{array}
if k < 6.20000000000000027e-5Initial program 37.1%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6470.6%
Simplified70.6%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.7%
Applied egg-rr70.7%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6478.5%
Applied egg-rr78.5%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6481.6%
Applied egg-rr81.6%
if 6.20000000000000027e-5 < k Initial program 24.0%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6471.9%
Simplified71.9%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr75.3%
associate-/l/N/A
times-fracN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
div-invN/A
metadata-evalN/A
associate-*l*N/A
/-rgt-identityN/A
clear-numN/A
div-invN/A
Applied egg-rr93.2%
Final simplification84.4%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 6.2e-5)
(* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m)))
(*
l
(/
(/ (cos k_m) (/ (* k_m t) l))
(* k_m (+ 0.25 (* (cos (* k_m 2.0)) -0.25)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 6.2e-5) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = l * ((cos(k_m) / ((k_m * t) / l)) / (k_m * (0.25 + (cos((k_m * 2.0)) * -0.25))));
}
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 <= 6.2d-5) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else
tmp = l * ((cos(k_m) / ((k_m * t) / l)) / (k_m * (0.25d0 + (cos((k_m * 2.0d0)) * (-0.25d0)))))
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 <= 6.2e-5) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = l * ((Math.cos(k_m) / ((k_m * t) / l)) / (k_m * (0.25 + (Math.cos((k_m * 2.0)) * -0.25))));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 6.2e-5: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) else: tmp = l * ((math.cos(k_m) / ((k_m * t) / l)) / (k_m * (0.25 + (math.cos((k_m * 2.0)) * -0.25)))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 6.2e-5) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); else tmp = Float64(l * Float64(Float64(cos(k_m) / Float64(Float64(k_m * t) / l)) / Float64(k_m * Float64(0.25 + Float64(cos(Float64(k_m * 2.0)) * -0.25))))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 6.2e-5) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); else tmp = l * ((cos(k_m) / ((k_m * t) / l)) / (k_m * (0.25 + (cos((k_m * 2.0)) * -0.25)))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 6.2e-5], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l * N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[(N[(k$95$m * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * N[(0.25 + N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 6.2 \cdot 10^{-5}:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \frac{\frac{\cos k\_m}{\frac{k\_m \cdot t}{\ell}}}{k\_m \cdot \left(0.25 + \cos \left(k\_m \cdot 2\right) \cdot -0.25\right)}\\
\end{array}
\end{array}
if k < 6.20000000000000027e-5Initial program 37.1%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6470.6%
Simplified70.6%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.7%
Applied egg-rr70.7%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6478.5%
Applied egg-rr78.5%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6481.6%
Applied egg-rr81.6%
if 6.20000000000000027e-5 < k Initial program 24.0%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6471.9%
Simplified71.9%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr75.3%
associate-/l/N/A
associate-/l*N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
div-invN/A
metadata-evalN/A
associate-*l*N/A
Applied egg-rr90.2%
Final simplification83.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 6.2e-5)
(* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m)))
(*
l
(/
(/ (cos k_m) (/ k_m l))
(* k_m (* t (+ 0.25 (* (cos (* k_m 2.0)) -0.25))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 6.2e-5) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = l * ((cos(k_m) / (k_m / l)) / (k_m * (t * (0.25 + (cos((k_m * 2.0)) * -0.25)))));
}
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 <= 6.2d-5) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else
tmp = l * ((cos(k_m) / (k_m / l)) / (k_m * (t * (0.25d0 + (cos((k_m * 2.0d0)) * (-0.25d0))))))
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 <= 6.2e-5) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = l * ((Math.cos(k_m) / (k_m / l)) / (k_m * (t * (0.25 + (Math.cos((k_m * 2.0)) * -0.25)))));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 6.2e-5: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) else: tmp = l * ((math.cos(k_m) / (k_m / l)) / (k_m * (t * (0.25 + (math.cos((k_m * 2.0)) * -0.25))))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 6.2e-5) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); else tmp = Float64(l * Float64(Float64(cos(k_m) / Float64(k_m / l)) / Float64(k_m * Float64(t * Float64(0.25 + Float64(cos(Float64(k_m * 2.0)) * -0.25)))))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 6.2e-5) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); else tmp = l * ((cos(k_m) / (k_m / l)) / (k_m * (t * (0.25 + (cos((k_m * 2.0)) * -0.25))))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 6.2e-5], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l * N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * N[(t * N[(0.25 + N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 6.2 \cdot 10^{-5}:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \frac{\frac{\cos k\_m}{\frac{k\_m}{\ell}}}{k\_m \cdot \left(t \cdot \left(0.25 + \cos \left(k\_m \cdot 2\right) \cdot -0.25\right)\right)}\\
\end{array}
\end{array}
if k < 6.20000000000000027e-5Initial program 37.1%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6470.6%
Simplified70.6%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.7%
Applied egg-rr70.7%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6478.5%
Applied egg-rr78.5%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6481.6%
Applied egg-rr81.6%
if 6.20000000000000027e-5 < k Initial program 24.0%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6471.9%
Simplified71.9%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr75.3%
associate-/l/N/A
associate-/r*N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
Applied egg-rr82.8%
associate-/r/N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr90.1%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 3.4e-105)
(/ (/ (/ (* l (/ l (* k_m k_m))) k_m) t) (/ k_m 2.0))
(if (<= k_m 4.8e+106)
(*
(/ (* l (cos k_m)) k_m)
(/
(* (/ l t) (+ 2.0 (* (* k_m k_m) 0.6666666666666666)))
(* k_m (* k_m k_m))))
(/
(/ 1.0 (/ t (/ (* l 2.0) (* k_m (+ 0.5 (* (cos (* k_m 2.0)) -0.5))))))
(/ k_m l)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 3.4e-105) {
tmp = (((l * (l / (k_m * k_m))) / k_m) / t) / (k_m / 2.0);
} else if (k_m <= 4.8e+106) {
tmp = ((l * cos(k_m)) / k_m) * (((l / t) * (2.0 + ((k_m * k_m) * 0.6666666666666666))) / (k_m * (k_m * k_m)));
} else {
tmp = (1.0 / (t / ((l * 2.0) / (k_m * (0.5 + (cos((k_m * 2.0)) * -0.5)))))) / (k_m / l);
}
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 <= 3.4d-105) then
tmp = (((l * (l / (k_m * k_m))) / k_m) / t) / (k_m / 2.0d0)
else if (k_m <= 4.8d+106) then
tmp = ((l * cos(k_m)) / k_m) * (((l / t) * (2.0d0 + ((k_m * k_m) * 0.6666666666666666d0))) / (k_m * (k_m * k_m)))
else
tmp = (1.0d0 / (t / ((l * 2.0d0) / (k_m * (0.5d0 + (cos((k_m * 2.0d0)) * (-0.5d0))))))) / (k_m / l)
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 <= 3.4e-105) {
tmp = (((l * (l / (k_m * k_m))) / k_m) / t) / (k_m / 2.0);
} else if (k_m <= 4.8e+106) {
tmp = ((l * Math.cos(k_m)) / k_m) * (((l / t) * (2.0 + ((k_m * k_m) * 0.6666666666666666))) / (k_m * (k_m * k_m)));
} else {
tmp = (1.0 / (t / ((l * 2.0) / (k_m * (0.5 + (Math.cos((k_m * 2.0)) * -0.5)))))) / (k_m / l);
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 3.4e-105: tmp = (((l * (l / (k_m * k_m))) / k_m) / t) / (k_m / 2.0) elif k_m <= 4.8e+106: tmp = ((l * math.cos(k_m)) / k_m) * (((l / t) * (2.0 + ((k_m * k_m) * 0.6666666666666666))) / (k_m * (k_m * k_m))) else: tmp = (1.0 / (t / ((l * 2.0) / (k_m * (0.5 + (math.cos((k_m * 2.0)) * -0.5)))))) / (k_m / l) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 3.4e-105) tmp = Float64(Float64(Float64(Float64(l * Float64(l / Float64(k_m * k_m))) / k_m) / t) / Float64(k_m / 2.0)); elseif (k_m <= 4.8e+106) tmp = Float64(Float64(Float64(l * cos(k_m)) / k_m) * Float64(Float64(Float64(l / t) * Float64(2.0 + Float64(Float64(k_m * k_m) * 0.6666666666666666))) / Float64(k_m * Float64(k_m * k_m)))); else tmp = Float64(Float64(1.0 / Float64(t / Float64(Float64(l * 2.0) / Float64(k_m * Float64(0.5 + Float64(cos(Float64(k_m * 2.0)) * -0.5)))))) / Float64(k_m / l)); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 3.4e-105) tmp = (((l * (l / (k_m * k_m))) / k_m) / t) / (k_m / 2.0); elseif (k_m <= 4.8e+106) tmp = ((l * cos(k_m)) / k_m) * (((l / t) * (2.0 + ((k_m * k_m) * 0.6666666666666666))) / (k_m * (k_m * k_m))); else tmp = (1.0 / (t / ((l * 2.0) / (k_m * (0.5 + (cos((k_m * 2.0)) * -0.5)))))) / (k_m / l); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 3.4e-105], N[(N[(N[(N[(l * N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision] / t), $MachinePrecision] / N[(k$95$m / 2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 4.8e+106], N[(N[(N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision] * N[(N[(N[(l / t), $MachinePrecision] * N[(2.0 + N[(N[(k$95$m * k$95$m), $MachinePrecision] * 0.6666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 / N[(t / N[(N[(l * 2.0), $MachinePrecision] / N[(k$95$m * N[(0.5 + N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m / l), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 3.4 \cdot 10^{-105}:\\
\;\;\;\;\frac{\frac{\frac{\ell \cdot \frac{\ell}{k\_m \cdot k\_m}}{k\_m}}{t}}{\frac{k\_m}{2}}\\
\mathbf{elif}\;k\_m \leq 4.8 \cdot 10^{+106}:\\
\;\;\;\;\frac{\ell \cdot \cos k\_m}{k\_m} \cdot \frac{\frac{\ell}{t} \cdot \left(2 + \left(k\_m \cdot k\_m\right) \cdot 0.6666666666666666\right)}{k\_m \cdot \left(k\_m \cdot k\_m\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{\frac{t}{\frac{\ell \cdot 2}{k\_m \cdot \left(0.5 + \cos \left(k\_m \cdot 2\right) \cdot -0.5\right)}}}}{\frac{k\_m}{\ell}}\\
\end{array}
\end{array}
if k < 3.39999999999999992e-105Initial program 37.0%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6476.3%
Simplified76.3%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr73.6%
Taylor expanded in k around 0
associate-/r*N/A
/-lowering-/.f64N/A
unpow3N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
unpow2N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6480.7%
Simplified80.7%
if 3.39999999999999992e-105 < k < 4.8000000000000001e106Initial program 33.7%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6486.9%
Simplified86.9%
associate-*r*N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr79.9%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
/-lowering-/.f64N/A
Simplified77.7%
if 4.8000000000000001e106 < k Initial program 21.8%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6465.1%
Simplified65.1%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.5%
Simplified52.5%
associate-*r*N/A
*-commutativeN/A
times-fracN/A
clear-numN/A
un-div-invN/A
unpow2N/A
sqr-sin-aN/A
cancel-sign-sub-invN/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
/-lowering-/.f64N/A
Applied egg-rr57.5%
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
cancel-sign-sub-invN/A
*-commutativeN/A
sqr-sin-aN/A
unpow2N/A
associate-/r*N/A
*-commutativeN/A
associate-*r*N/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr58.8%
Final simplification76.6%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 136000000000.0)
(* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m)))
(if (<= k_m 2.9e+104)
(/ (/ (/ (* l (* l (cos k_m))) (* k_m t)) (* k_m k_m)) (/ k_m 2.0))
(/
(/ 1.0 (/ t (/ (* l 2.0) (* k_m (+ 0.5 (* (cos (* k_m 2.0)) -0.5))))))
(/ k_m l)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 136000000000.0) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else if (k_m <= 2.9e+104) {
tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0);
} else {
tmp = (1.0 / (t / ((l * 2.0) / (k_m * (0.5 + (cos((k_m * 2.0)) * -0.5)))))) / (k_m / l);
}
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 <= 136000000000.0d0) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else if (k_m <= 2.9d+104) then
tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0d0)
else
tmp = (1.0d0 / (t / ((l * 2.0d0) / (k_m * (0.5d0 + (cos((k_m * 2.0d0)) * (-0.5d0))))))) / (k_m / l)
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 <= 136000000000.0) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else if (k_m <= 2.9e+104) {
tmp = (((l * (l * Math.cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0);
} else {
tmp = (1.0 / (t / ((l * 2.0) / (k_m * (0.5 + (Math.cos((k_m * 2.0)) * -0.5)))))) / (k_m / l);
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 136000000000.0: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) elif k_m <= 2.9e+104: tmp = (((l * (l * math.cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0) else: tmp = (1.0 / (t / ((l * 2.0) / (k_m * (0.5 + (math.cos((k_m * 2.0)) * -0.5)))))) / (k_m / l) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 136000000000.0) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); elseif (k_m <= 2.9e+104) tmp = Float64(Float64(Float64(Float64(l * Float64(l * cos(k_m))) / Float64(k_m * t)) / Float64(k_m * k_m)) / Float64(k_m / 2.0)); else tmp = Float64(Float64(1.0 / Float64(t / Float64(Float64(l * 2.0) / Float64(k_m * Float64(0.5 + Float64(cos(Float64(k_m * 2.0)) * -0.5)))))) / Float64(k_m / l)); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 136000000000.0) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); elseif (k_m <= 2.9e+104) tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0); else tmp = (1.0 / (t / ((l * 2.0) / (k_m * (0.5 + (cos((k_m * 2.0)) * -0.5)))))) / (k_m / l); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 136000000000.0], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.9e+104], N[(N[(N[(N[(l * N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(k$95$m / 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 / N[(t / N[(N[(l * 2.0), $MachinePrecision] / N[(k$95$m * N[(0.5 + N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m / l), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 136000000000:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{elif}\;k\_m \leq 2.9 \cdot 10^{+104}:\\
\;\;\;\;\frac{\frac{\frac{\ell \cdot \left(\ell \cdot \cos k\_m\right)}{k\_m \cdot t}}{k\_m \cdot k\_m}}{\frac{k\_m}{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{\frac{t}{\frac{\ell \cdot 2}{k\_m \cdot \left(0.5 + \cos \left(k\_m \cdot 2\right) \cdot -0.5\right)}}}}{\frac{k\_m}{\ell}}\\
\end{array}
\end{array}
if k < 1.36e11Initial program 37.3%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6469.9%
Simplified69.9%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.0%
Applied egg-rr70.0%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6477.8%
Applied egg-rr77.8%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6480.8%
Applied egg-rr80.8%
if 1.36e11 < k < 2.8999999999999998e104Initial program 26.2%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6484.1%
Simplified84.1%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr84.6%
Taylor expanded in k around 0
unpow2N/A
*-lowering-*.f6450.9%
Simplified50.9%
if 2.8999999999999998e104 < k Initial program 21.8%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6465.1%
Simplified65.1%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.5%
Simplified52.5%
associate-*r*N/A
*-commutativeN/A
times-fracN/A
clear-numN/A
un-div-invN/A
unpow2N/A
sqr-sin-aN/A
cancel-sign-sub-invN/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
/-lowering-/.f64N/A
Applied egg-rr57.5%
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
cancel-sign-sub-invN/A
*-commutativeN/A
sqr-sin-aN/A
unpow2N/A
associate-/r*N/A
*-commutativeN/A
associate-*r*N/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr58.8%
Final simplification75.0%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 136000000000.0)
(* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m)))
(if (<= k_m 2.9e+104)
(/ (/ (/ (* l (* l (cos k_m))) (* k_m t)) (* k_m k_m)) (/ k_m 2.0))
(/
(/ (/ 2.0 (+ 0.5 (* (cos (* k_m 2.0)) -0.5))) (/ t (/ l k_m)))
(/ k_m l)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 136000000000.0) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else if (k_m <= 2.9e+104) {
tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0);
} else {
tmp = ((2.0 / (0.5 + (cos((k_m * 2.0)) * -0.5))) / (t / (l / k_m))) / (k_m / l);
}
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 <= 136000000000.0d0) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else if (k_m <= 2.9d+104) then
tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0d0)
else
tmp = ((2.0d0 / (0.5d0 + (cos((k_m * 2.0d0)) * (-0.5d0)))) / (t / (l / k_m))) / (k_m / l)
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 <= 136000000000.0) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else if (k_m <= 2.9e+104) {
tmp = (((l * (l * Math.cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0);
} else {
tmp = ((2.0 / (0.5 + (Math.cos((k_m * 2.0)) * -0.5))) / (t / (l / k_m))) / (k_m / l);
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 136000000000.0: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) elif k_m <= 2.9e+104: tmp = (((l * (l * math.cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0) else: tmp = ((2.0 / (0.5 + (math.cos((k_m * 2.0)) * -0.5))) / (t / (l / k_m))) / (k_m / l) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 136000000000.0) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); elseif (k_m <= 2.9e+104) tmp = Float64(Float64(Float64(Float64(l * Float64(l * cos(k_m))) / Float64(k_m * t)) / Float64(k_m * k_m)) / Float64(k_m / 2.0)); else tmp = Float64(Float64(Float64(2.0 / Float64(0.5 + Float64(cos(Float64(k_m * 2.0)) * -0.5))) / Float64(t / Float64(l / k_m))) / Float64(k_m / l)); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 136000000000.0) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); elseif (k_m <= 2.9e+104) tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0); else tmp = ((2.0 / (0.5 + (cos((k_m * 2.0)) * -0.5))) / (t / (l / k_m))) / (k_m / l); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 136000000000.0], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.9e+104], N[(N[(N[(N[(l * N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(k$95$m / 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / N[(0.5 + N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t / N[(l / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m / l), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 136000000000:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{elif}\;k\_m \leq 2.9 \cdot 10^{+104}:\\
\;\;\;\;\frac{\frac{\frac{\ell \cdot \left(\ell \cdot \cos k\_m\right)}{k\_m \cdot t}}{k\_m \cdot k\_m}}{\frac{k\_m}{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{2}{0.5 + \cos \left(k\_m \cdot 2\right) \cdot -0.5}}{\frac{t}{\frac{\ell}{k\_m}}}}{\frac{k\_m}{\ell}}\\
\end{array}
\end{array}
if k < 1.36e11Initial program 37.3%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6469.9%
Simplified69.9%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.0%
Applied egg-rr70.0%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6477.8%
Applied egg-rr77.8%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6480.8%
Applied egg-rr80.8%
if 1.36e11 < k < 2.8999999999999998e104Initial program 26.2%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6484.1%
Simplified84.1%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr84.6%
Taylor expanded in k around 0
unpow2N/A
*-lowering-*.f6450.9%
Simplified50.9%
if 2.8999999999999998e104 < k Initial program 21.8%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6465.1%
Simplified65.1%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.5%
Simplified52.5%
associate-*r*N/A
*-commutativeN/A
times-fracN/A
clear-numN/A
un-div-invN/A
unpow2N/A
sqr-sin-aN/A
cancel-sign-sub-invN/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
/-lowering-/.f64N/A
Applied egg-rr57.5%
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
cancel-sign-sub-invN/A
*-commutativeN/A
sqr-sin-aN/A
unpow2N/A
associate-/r*N/A
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr58.8%
Final simplification75.0%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 70000000000.0)
(* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m)))
(if (<= k_m 2.9e+104)
(/ (/ (/ (* l (* l (cos k_m))) (* k_m t)) (* k_m k_m)) (/ k_m 2.0))
(/
(* (/ l t) (/ 2.0 (* k_m (+ 0.5 (* (cos (* k_m 2.0)) -0.5)))))
(/ k_m l)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 70000000000.0) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else if (k_m <= 2.9e+104) {
tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0);
} else {
tmp = ((l / t) * (2.0 / (k_m * (0.5 + (cos((k_m * 2.0)) * -0.5))))) / (k_m / l);
}
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 <= 70000000000.0d0) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else if (k_m <= 2.9d+104) then
tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0d0)
else
tmp = ((l / t) * (2.0d0 / (k_m * (0.5d0 + (cos((k_m * 2.0d0)) * (-0.5d0)))))) / (k_m / l)
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 <= 70000000000.0) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else if (k_m <= 2.9e+104) {
tmp = (((l * (l * Math.cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0);
} else {
tmp = ((l / t) * (2.0 / (k_m * (0.5 + (Math.cos((k_m * 2.0)) * -0.5))))) / (k_m / l);
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 70000000000.0: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) elif k_m <= 2.9e+104: tmp = (((l * (l * math.cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0) else: tmp = ((l / t) * (2.0 / (k_m * (0.5 + (math.cos((k_m * 2.0)) * -0.5))))) / (k_m / l) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 70000000000.0) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); elseif (k_m <= 2.9e+104) tmp = Float64(Float64(Float64(Float64(l * Float64(l * cos(k_m))) / Float64(k_m * t)) / Float64(k_m * k_m)) / Float64(k_m / 2.0)); else tmp = Float64(Float64(Float64(l / t) * Float64(2.0 / Float64(k_m * Float64(0.5 + Float64(cos(Float64(k_m * 2.0)) * -0.5))))) / Float64(k_m / l)); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 70000000000.0) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); elseif (k_m <= 2.9e+104) tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0); else tmp = ((l / t) * (2.0 / (k_m * (0.5 + (cos((k_m * 2.0)) * -0.5))))) / (k_m / l); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 70000000000.0], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.9e+104], N[(N[(N[(N[(l * N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(k$95$m / 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(l / t), $MachinePrecision] * N[(2.0 / N[(k$95$m * N[(0.5 + N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m / l), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 70000000000:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{elif}\;k\_m \leq 2.9 \cdot 10^{+104}:\\
\;\;\;\;\frac{\frac{\frac{\ell \cdot \left(\ell \cdot \cos k\_m\right)}{k\_m \cdot t}}{k\_m \cdot k\_m}}{\frac{k\_m}{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell}{t} \cdot \frac{2}{k\_m \cdot \left(0.5 + \cos \left(k\_m \cdot 2\right) \cdot -0.5\right)}}{\frac{k\_m}{\ell}}\\
\end{array}
\end{array}
if k < 7e10Initial program 37.3%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6469.9%
Simplified69.9%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.0%
Applied egg-rr70.0%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6477.8%
Applied egg-rr77.8%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6480.8%
Applied egg-rr80.8%
if 7e10 < k < 2.8999999999999998e104Initial program 26.2%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6484.1%
Simplified84.1%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr84.6%
Taylor expanded in k around 0
unpow2N/A
*-lowering-*.f6450.9%
Simplified50.9%
if 2.8999999999999998e104 < k Initial program 21.8%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6465.1%
Simplified65.1%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.5%
Simplified52.5%
associate-*r*N/A
*-commutativeN/A
times-fracN/A
clear-numN/A
un-div-invN/A
unpow2N/A
sqr-sin-aN/A
cancel-sign-sub-invN/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
/-lowering-/.f64N/A
Applied egg-rr57.5%
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
Applied egg-rr57.9%
Final simplification74.9%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 136000000000.0)
(* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m)))
(if (<= k_m 2.9e+104)
(/ (/ (/ (* l (* l (cos k_m))) (* k_m t)) (* k_m k_m)) (/ k_m 2.0))
(/
1.0
(*
(/ k_m l)
(/ (* k_m (+ 0.25 (* (cos (* k_m 2.0)) -0.25))) (/ l t)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 136000000000.0) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else if (k_m <= 2.9e+104) {
tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0);
} else {
tmp = 1.0 / ((k_m / l) * ((k_m * (0.25 + (cos((k_m * 2.0)) * -0.25))) / (l / 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 <= 136000000000.0d0) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else if (k_m <= 2.9d+104) then
tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0d0)
else
tmp = 1.0d0 / ((k_m / l) * ((k_m * (0.25d0 + (cos((k_m * 2.0d0)) * (-0.25d0)))) / (l / 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 <= 136000000000.0) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else if (k_m <= 2.9e+104) {
tmp = (((l * (l * Math.cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0);
} else {
tmp = 1.0 / ((k_m / l) * ((k_m * (0.25 + (Math.cos((k_m * 2.0)) * -0.25))) / (l / t)));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 136000000000.0: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) elif k_m <= 2.9e+104: tmp = (((l * (l * math.cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0) else: tmp = 1.0 / ((k_m / l) * ((k_m * (0.25 + (math.cos((k_m * 2.0)) * -0.25))) / (l / t))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 136000000000.0) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); elseif (k_m <= 2.9e+104) tmp = Float64(Float64(Float64(Float64(l * Float64(l * cos(k_m))) / Float64(k_m * t)) / Float64(k_m * k_m)) / Float64(k_m / 2.0)); else tmp = Float64(1.0 / Float64(Float64(k_m / l) * Float64(Float64(k_m * Float64(0.25 + Float64(cos(Float64(k_m * 2.0)) * -0.25))) / Float64(l / t)))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 136000000000.0) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); elseif (k_m <= 2.9e+104) tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0); else tmp = 1.0 / ((k_m / l) * ((k_m * (0.25 + (cos((k_m * 2.0)) * -0.25))) / (l / t))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 136000000000.0], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.9e+104], N[(N[(N[(N[(l * N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(k$95$m / 2.0), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(k$95$m / l), $MachinePrecision] * N[(N[(k$95$m * N[(0.25 + N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 136000000000:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{elif}\;k\_m \leq 2.9 \cdot 10^{+104}:\\
\;\;\;\;\frac{\frac{\frac{\ell \cdot \left(\ell \cdot \cos k\_m\right)}{k\_m \cdot t}}{k\_m \cdot k\_m}}{\frac{k\_m}{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{k\_m}{\ell} \cdot \frac{k\_m \cdot \left(0.25 + \cos \left(k\_m \cdot 2\right) \cdot -0.25\right)}{\frac{\ell}{t}}}\\
\end{array}
\end{array}
if k < 1.36e11Initial program 37.3%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6469.9%
Simplified69.9%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.0%
Applied egg-rr70.0%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6477.8%
Applied egg-rr77.8%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6480.8%
Applied egg-rr80.8%
if 1.36e11 < k < 2.8999999999999998e104Initial program 26.2%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6484.1%
Simplified84.1%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr84.6%
Taylor expanded in k around 0
unpow2N/A
*-lowering-*.f6450.9%
Simplified50.9%
if 2.8999999999999998e104 < k Initial program 21.8%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6465.1%
Simplified65.1%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.5%
Simplified52.5%
Applied egg-rr57.9%
Final simplification74.9%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= (* l l) 0.0) (* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m))) (/ (/ (/ (* l (* l (cos k_m))) (* k_m t)) (* k_m k_m)) (/ k_m 2.0))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0);
}
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 * l) <= 0.0d0) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else
tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0d0)
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 * l) <= 0.0) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = (((l * (l * Math.cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0);
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if (l * l) <= 0.0: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) else: tmp = (((l * (l * math.cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (Float64(l * l) <= 0.0) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); else tmp = Float64(Float64(Float64(Float64(l * Float64(l * cos(k_m))) / Float64(k_m * t)) / Float64(k_m * k_m)) / Float64(k_m / 2.0)); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if ((l * l) <= 0.0) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); else tmp = (((l * (l * cos(k_m))) / (k_m * t)) / (k_m * k_m)) / (k_m / 2.0); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[N[(l * l), $MachinePrecision], 0.0], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(l * N[(l * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(k$95$m / 2.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 0:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{\ell \cdot \left(\ell \cdot \cos k\_m\right)}{k\_m \cdot t}}{k\_m \cdot k\_m}}{\frac{k\_m}{2}}\\
\end{array}
\end{array}
if (*.f64 l l) < 0.0Initial program 12.5%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6458.5%
Simplified58.5%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6458.5%
Applied egg-rr58.5%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6483.0%
Applied egg-rr83.0%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6488.1%
Applied egg-rr88.1%
if 0.0 < (*.f64 l l) Initial program 40.1%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6481.4%
Simplified81.4%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr78.6%
Taylor expanded in k around 0
unpow2N/A
*-lowering-*.f6473.5%
Simplified73.5%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 5.1e-6) (* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m))) (/ 1.0 (/ (* k_m (* k_m k_m)) (/ 2.0 (/ k_m (/ l (/ t l))))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 5.1e-6) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = 1.0 / ((k_m * (k_m * k_m)) / (2.0 / (k_m / (l / (t / l)))));
}
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 <= 5.1d-6) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else
tmp = 1.0d0 / ((k_m * (k_m * k_m)) / (2.0d0 / (k_m / (l / (t / l)))))
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 <= 5.1e-6) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = 1.0 / ((k_m * (k_m * k_m)) / (2.0 / (k_m / (l / (t / l)))));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 5.1e-6: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) else: tmp = 1.0 / ((k_m * (k_m * k_m)) / (2.0 / (k_m / (l / (t / l))))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 5.1e-6) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); else tmp = Float64(1.0 / Float64(Float64(k_m * Float64(k_m * k_m)) / Float64(2.0 / Float64(k_m / Float64(l / Float64(t / l)))))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 5.1e-6) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); else tmp = 1.0 / ((k_m * (k_m * k_m)) / (2.0 / (k_m / (l / (t / l))))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 5.1e-6], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(k$95$m * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] / N[(2.0 / N[(k$95$m / N[(l / N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 5.1 \cdot 10^{-6}:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{k\_m \cdot \left(k\_m \cdot k\_m\right)}{\frac{2}{\frac{k\_m}{\frac{\ell}{\frac{t}{\ell}}}}}}\\
\end{array}
\end{array}
if k < 5.1000000000000003e-6Initial program 37.3%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6470.9%
Simplified70.9%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6471.0%
Applied egg-rr71.0%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6478.5%
Applied egg-rr78.5%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6481.6%
Applied egg-rr81.6%
if 5.1000000000000003e-6 < k Initial program 23.7%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6442.6%
Simplified42.6%
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6442.6%
Applied egg-rr42.6%
frac-2negN/A
/-lowering-/.f64N/A
metadata-evalN/A
neg-sub0N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6442.6%
Applied egg-rr42.6%
clear-numN/A
frac-2negN/A
metadata-evalN/A
/-lowering-/.f64N/A
neg-sub0N/A
--lowering--.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
sub0-negN/A
frac-2negN/A
associate-/l/N/A
/-lowering-/.f64N/A
Applied egg-rr48.5%
Final simplification72.6%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 14.2) (* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m))) (/ (/ (/ 2.0 k_m) (* k_m (* k_m k_m))) (/ t (* l l)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 14.2) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = ((2.0 / k_m) / (k_m * (k_m * k_m))) / (t / (l * l));
}
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 <= 14.2d0) then
tmp = (l * 2.0d0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
else
tmp = ((2.0d0 / k_m) / (k_m * (k_m * k_m))) / (t / (l * l))
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 <= 14.2) {
tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
} else {
tmp = ((2.0 / k_m) / (k_m * (k_m * k_m))) / (t / (l * l));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 14.2: tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)) else: tmp = ((2.0 / k_m) / (k_m * (k_m * k_m))) / (t / (l * l)) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 14.2) tmp = Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))); else tmp = Float64(Float64(Float64(2.0 / k_m) / Float64(k_m * Float64(k_m * k_m))) / Float64(t / Float64(l * l))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 14.2) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); else tmp = ((2.0 / k_m) / (k_m * (k_m * k_m))) / (t / (l * l)); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 14.2], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / k$95$m), $MachinePrecision] / N[(k$95$m * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 14.2:\\
\;\;\;\;\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{2}{k\_m}}{k\_m \cdot \left(k\_m \cdot k\_m\right)}}{\frac{t}{\ell \cdot \ell}}\\
\end{array}
\end{array}
if k < 14.199999999999999Initial program 37.1%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6470.6%
Simplified70.6%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.7%
Applied egg-rr70.7%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6478.5%
Applied egg-rr78.5%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6481.6%
Applied egg-rr81.6%
if 14.199999999999999 < k Initial program 24.0%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6443.2%
Simplified43.2%
associate-/r*N/A
frac-2negN/A
distribute-frac-negN/A
neg-lowering-neg.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
neg-sub0N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6444.7%
Applied egg-rr44.7%
Final simplification72.1%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= t 4.9e+167) (/ (/ (/ (* l (/ l (* k_m k_m))) k_m) t) (/ k_m 2.0)) (/ (/ (* l 2.0) (* k_m (* k_m (* k_m t)))) (/ k_m l))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 4.9e+167) {
tmp = (((l * (l / (k_m * k_m))) / k_m) / t) / (k_m / 2.0);
} else {
tmp = ((l * 2.0) / (k_m * (k_m * (k_m * t)))) / (k_m / l);
}
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 (t <= 4.9d+167) then
tmp = (((l * (l / (k_m * k_m))) / k_m) / t) / (k_m / 2.0d0)
else
tmp = ((l * 2.0d0) / (k_m * (k_m * (k_m * t)))) / (k_m / l)
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 (t <= 4.9e+167) {
tmp = (((l * (l / (k_m * k_m))) / k_m) / t) / (k_m / 2.0);
} else {
tmp = ((l * 2.0) / (k_m * (k_m * (k_m * t)))) / (k_m / l);
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if t <= 4.9e+167: tmp = (((l * (l / (k_m * k_m))) / k_m) / t) / (k_m / 2.0) else: tmp = ((l * 2.0) / (k_m * (k_m * (k_m * t)))) / (k_m / l) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 4.9e+167) tmp = Float64(Float64(Float64(Float64(l * Float64(l / Float64(k_m * k_m))) / k_m) / t) / Float64(k_m / 2.0)); else tmp = Float64(Float64(Float64(l * 2.0) / Float64(k_m * Float64(k_m * Float64(k_m * t)))) / Float64(k_m / l)); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (t <= 4.9e+167) tmp = (((l * (l / (k_m * k_m))) / k_m) / t) / (k_m / 2.0); else tmp = ((l * 2.0) / (k_m * (k_m * (k_m * t)))) / (k_m / l); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 4.9e+167], N[(N[(N[(N[(l * N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision] / t), $MachinePrecision] / N[(k$95$m / 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(l * 2.0), $MachinePrecision] / N[(k$95$m * N[(k$95$m * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k$95$m / l), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 4.9 \cdot 10^{+167}:\\
\;\;\;\;\frac{\frac{\frac{\ell \cdot \frac{\ell}{k\_m \cdot k\_m}}{k\_m}}{t}}{\frac{k\_m}{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell \cdot 2}{k\_m \cdot \left(k\_m \cdot \left(k\_m \cdot t\right)\right)}}{\frac{k\_m}{\ell}}\\
\end{array}
\end{array}
if t < 4.9000000000000003e167Initial program 36.0%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6476.5%
Simplified76.5%
times-fracN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr73.9%
Taylor expanded in k around 0
associate-/r*N/A
/-lowering-/.f64N/A
unpow3N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
unpow2N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6473.0%
Simplified73.0%
if 4.9000000000000003e167 < t Initial program 13.0%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f6474.2%
Simplified74.2%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6470.7%
Simplified70.7%
associate-*r*N/A
*-commutativeN/A
times-fracN/A
clear-numN/A
un-div-invN/A
unpow2N/A
sqr-sin-aN/A
cancel-sign-sub-invN/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
/-lowering-/.f64N/A
Applied egg-rr59.4%
Taylor expanded in k around 0
cube-multN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6478.8%
Simplified78.8%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* (* l 2.0) (/ (/ (/ l k_m) (* k_m t)) (* k_m k_m))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return (l * 2.0) * (((l / k_m) / (k_m * 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 * 2.0d0) * (((l / k_m) / (k_m * 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 * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m));
}
k_m = math.fabs(k) def code(t, l, k_m): return (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m))
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(l * 2.0) * Float64(Float64(Float64(l / k_m) / Float64(k_m * t)) / Float64(k_m * k_m))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = (l * 2.0) * (((l / k_m) / (k_m * t)) / (k_m * k_m)); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(l / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\left(\ell \cdot 2\right) \cdot \frac{\frac{\frac{\ell}{k\_m}}{k\_m \cdot t}}{k\_m \cdot k\_m}
\end{array}
Initial program 33.9%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6463.8%
Simplified63.8%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.9%
Applied egg-rr63.9%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.3%
Applied egg-rr70.3%
associate-/r*N/A
associate-*r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6472.7%
Applied egg-rr72.7%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* (* l 2.0) (/ l (* k_m (* k_m (* t (* k_m k_m)))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return (l * 2.0) * (l / (k_m * (k_m * (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 * 2.0d0) * (l / (k_m * (k_m * (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 * 2.0) * (l / (k_m * (k_m * (t * (k_m * k_m)))));
}
k_m = math.fabs(k) def code(t, l, k_m): return (l * 2.0) * (l / (k_m * (k_m * (t * (k_m * k_m)))))
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(l * 2.0) * Float64(l / Float64(k_m * Float64(k_m * Float64(t * Float64(k_m * k_m)))))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = (l * 2.0) * (l / (k_m * (k_m * (t * (k_m * k_m))))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(l * 2.0), $MachinePrecision] * N[(l / N[(k$95$m * N[(k$95$m * N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\left(\ell \cdot 2\right) \cdot \frac{\ell}{k\_m \cdot \left(k\_m \cdot \left(t \cdot \left(k\_m \cdot k\_m\right)\right)\right)}
\end{array}
Initial program 33.9%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6463.8%
Simplified63.8%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.9%
Applied egg-rr63.9%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.3%
Applied egg-rr70.3%
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.8%
Applied egg-rr70.8%
Final simplification70.8%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* (* l 2.0) (/ l (* k_m (* t (* k_m (* k_m k_m)))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return (l * 2.0) * (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 = (l * 2.0d0) * (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 (l * 2.0) * (l / (k_m * (t * (k_m * (k_m * k_m)))));
}
k_m = math.fabs(k) def code(t, l, k_m): return (l * 2.0) * (l / (k_m * (t * (k_m * (k_m * k_m)))))
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(l * 2.0) * 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 = (l * 2.0) * (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[(l * 2.0), $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(\ell \cdot 2\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 33.9%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6463.8%
Simplified63.8%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.9%
Applied egg-rr63.9%
associate-*l/N/A
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.3%
Applied egg-rr70.3%
herbie shell --seed 2024192
(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))))