
(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 18 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
(let* ((t_1 (* (sin k_m) (* t (tan k_m)))))
(if (<= k_m 1.36e+29)
(/ 2.0 (* 2.0 (* (* t (/ (* t (sin k_m)) l)) (* (/ t l) (tan k_m)))))
(if (<= k_m 3e+150)
(/
(/ 2.0 t_1)
(/ (* k_m (* k_m (fma (* t t) (/ 2.0 (* k_m k_m)) 1.0))) (* l l)))
(/
2.0
(*
t
(* (/ 1.0 l) (* t_1 (* (/ t l) (fma (/ k_m t) (/ k_m t) 2.0))))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = sin(k_m) * (t * tan(k_m));
double tmp;
if (k_m <= 1.36e+29) {
tmp = 2.0 / (2.0 * ((t * ((t * sin(k_m)) / l)) * ((t / l) * tan(k_m))));
} else if (k_m <= 3e+150) {
tmp = (2.0 / t_1) / ((k_m * (k_m * fma((t * t), (2.0 / (k_m * k_m)), 1.0))) / (l * l));
} else {
tmp = 2.0 / (t * ((1.0 / l) * (t_1 * ((t / l) * fma((k_m / t), (k_m / t), 2.0)))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(sin(k_m) * Float64(t * tan(k_m))) tmp = 0.0 if (k_m <= 1.36e+29) tmp = Float64(2.0 / Float64(2.0 * Float64(Float64(t * Float64(Float64(t * sin(k_m)) / l)) * Float64(Float64(t / l) * tan(k_m))))); elseif (k_m <= 3e+150) tmp = Float64(Float64(2.0 / t_1) / Float64(Float64(k_m * Float64(k_m * fma(Float64(t * t), Float64(2.0 / Float64(k_m * k_m)), 1.0))) / Float64(l * l))); else tmp = Float64(2.0 / Float64(t * Float64(Float64(1.0 / l) * Float64(t_1 * Float64(Float64(t / l) * fma(Float64(k_m / t), Float64(k_m / t), 2.0)))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(N[Sin[k$95$m], $MachinePrecision] * N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k$95$m, 1.36e+29], N[(2.0 / N[(2.0 * N[(N[(t * N[(N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 3e+150], N[(N[(2.0 / t$95$1), $MachinePrecision] / N[(N[(k$95$m * N[(k$95$m * N[(N[(t * t), $MachinePrecision] * N[(2.0 / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(t * N[(N[(1.0 / l), $MachinePrecision] * N[(t$95$1 * N[(N[(t / l), $MachinePrecision] * N[(N[(k$95$m / t), $MachinePrecision] * N[(k$95$m / t), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := \sin k\_m \cdot \left(t \cdot \tan k\_m\right)\\
\mathbf{if}\;k\_m \leq 1.36 \cdot 10^{+29}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\left(t \cdot \frac{t \cdot \sin k\_m}{\ell}\right) \cdot \left(\frac{t}{\ell} \cdot \tan k\_m\right)\right)}\\
\mathbf{elif}\;k\_m \leq 3 \cdot 10^{+150}:\\
\;\;\;\;\frac{\frac{2}{t\_1}}{\frac{k\_m \cdot \left(k\_m \cdot \mathsf{fma}\left(t \cdot t, \frac{2}{k\_m \cdot k\_m}, 1\right)\right)}{\ell \cdot \ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t \cdot \left(\frac{1}{\ell} \cdot \left(t\_1 \cdot \left(\frac{t}{\ell} \cdot \mathsf{fma}\left(\frac{k\_m}{t}, \frac{k\_m}{t}, 2\right)\right)\right)\right)}\\
\end{array}
\end{array}
if k < 1.36e29Initial program 57.6%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6418.9
Applied egg-rr18.9%
Applied egg-rr75.5%
Taylor expanded in k around 0
Simplified75.8%
lift-/.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*l*N/A
associate-*r*N/A
lower-*.f64N/A
Applied egg-rr82.2%
if 1.36e29 < k < 3.00000000000000012e150Initial program 52.1%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6410.5
Applied egg-rr10.5%
Applied egg-rr52.3%
Taylor expanded in k around inf
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6485.8
Simplified85.8%
if 3.00000000000000012e150 < k Initial program 32.9%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6415.6
Applied egg-rr15.6%
Applied egg-rr57.1%
Applied egg-rr58.7%
associate-/r*N/A
lift-/.f64N/A
associate-*r/N/A
associate-*l/N/A
lift-/.f64N/A
lower-fma.f6468.5
Applied egg-rr68.5%
Final simplification81.1%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<=
(/
2.0
(*
(* (tan k_m) (* (sin k_m) (/ (pow t 3.0) (* l l))))
(+ 1.0 (+ 1.0 (pow (/ k_m t) 2.0)))))
4e+305)
(/ (* l l) (* k_m (* k_m (* t (* t t)))))
(/ (/ (/ (* l l) k_m) (* k_m t)) t)))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if ((2.0 / ((tan(k_m) * (sin(k_m) * (pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + pow((k_m / t), 2.0))))) <= 4e+305) {
tmp = (l * l) / (k_m * (k_m * (t * (t * t))));
} else {
tmp = (((l * l) / k_m) / (k_m * t)) / 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 ((2.0d0 / ((tan(k_m) * (sin(k_m) * ((t ** 3.0d0) / (l * l)))) * (1.0d0 + (1.0d0 + ((k_m / t) ** 2.0d0))))) <= 4d+305) then
tmp = (l * l) / (k_m * (k_m * (t * (t * t))))
else
tmp = (((l * l) / k_m) / (k_m * t)) / 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 ((2.0 / ((Math.tan(k_m) * (Math.sin(k_m) * (Math.pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + Math.pow((k_m / t), 2.0))))) <= 4e+305) {
tmp = (l * l) / (k_m * (k_m * (t * (t * t))));
} else {
tmp = (((l * l) / k_m) / (k_m * t)) / t;
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if (2.0 / ((math.tan(k_m) * (math.sin(k_m) * (math.pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + math.pow((k_m / t), 2.0))))) <= 4e+305: tmp = (l * l) / (k_m * (k_m * (t * (t * t)))) else: tmp = (((l * l) / k_m) / (k_m * t)) / t return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (Float64(2.0 / Float64(Float64(tan(k_m) * Float64(sin(k_m) * Float64((t ^ 3.0) / Float64(l * l)))) * Float64(1.0 + Float64(1.0 + (Float64(k_m / t) ^ 2.0))))) <= 4e+305) tmp = Float64(Float64(l * l) / Float64(k_m * Float64(k_m * Float64(t * Float64(t * t))))); else tmp = Float64(Float64(Float64(Float64(l * l) / k_m) / Float64(k_m * t)) / t); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if ((2.0 / ((tan(k_m) * (sin(k_m) * ((t ^ 3.0) / (l * l)))) * (1.0 + (1.0 + ((k_m / t) ^ 2.0))))) <= 4e+305) tmp = (l * l) / (k_m * (k_m * (t * (t * t)))); else tmp = (((l * l) / k_m) / (k_m * t)) / t; end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(1.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 4e+305], N[(N[(l * l), $MachinePrecision] / N[(k$95$m * N[(k$95$m * N[(t * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(l * l), $MachinePrecision] / k$95$m), $MachinePrecision] / N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{2}{\left(\tan k\_m \cdot \left(\sin k\_m \cdot \frac{{t}^{3}}{\ell \cdot \ell}\right)\right) \cdot \left(1 + \left(1 + {\left(\frac{k\_m}{t}\right)}^{2}\right)\right)} \leq 4 \cdot 10^{+305}:\\
\;\;\;\;\frac{\ell \cdot \ell}{k\_m \cdot \left(k\_m \cdot \left(t \cdot \left(t \cdot t\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{\ell \cdot \ell}{k\_m}}{k\_m \cdot t}}{t}\\
\end{array}
\end{array}
if (/.f64 #s(literal 2 binary64) (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t #s(literal 3 binary64)) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 k t) #s(literal 2 binary64))) #s(literal 1 binary64)))) < 3.9999999999999998e305Initial program 80.7%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6456.0
Simplified56.0%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6456.6
Simplified56.6%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6447.7
Simplified47.7%
Taylor expanded in k around -inf
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6478.0
Simplified78.0%
if 3.9999999999999998e305 < (/.f64 #s(literal 2 binary64) (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t #s(literal 3 binary64)) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 k t) #s(literal 2 binary64))) #s(literal 1 binary64)))) Initial program 16.3%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6417.6
Simplified17.6%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6422.4
Simplified22.4%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6432.3
Simplified32.3%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f6440.6
Applied egg-rr40.6%
Final simplification62.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<=
(/
2.0
(*
(* (tan k_m) (* (sin k_m) (/ (pow t 3.0) (* l l))))
(+ 1.0 (+ 1.0 (pow (/ k_m t) 2.0)))))
4e+305)
(/ (* l l) (* k_m (* k_m (* t (* t t)))))
(* (/ l (* k_m k_m)) (/ l (* t t)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if ((2.0 / ((tan(k_m) * (sin(k_m) * (pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + pow((k_m / t), 2.0))))) <= 4e+305) {
tmp = (l * l) / (k_m * (k_m * (t * (t * t))));
} else {
tmp = (l / (k_m * k_m)) * (l / (t * 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 ((2.0d0 / ((tan(k_m) * (sin(k_m) * ((t ** 3.0d0) / (l * l)))) * (1.0d0 + (1.0d0 + ((k_m / t) ** 2.0d0))))) <= 4d+305) then
tmp = (l * l) / (k_m * (k_m * (t * (t * t))))
else
tmp = (l / (k_m * k_m)) * (l / (t * 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 ((2.0 / ((Math.tan(k_m) * (Math.sin(k_m) * (Math.pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + Math.pow((k_m / t), 2.0))))) <= 4e+305) {
tmp = (l * l) / (k_m * (k_m * (t * (t * t))));
} else {
tmp = (l / (k_m * k_m)) * (l / (t * t));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if (2.0 / ((math.tan(k_m) * (math.sin(k_m) * (math.pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + math.pow((k_m / t), 2.0))))) <= 4e+305: tmp = (l * l) / (k_m * (k_m * (t * (t * t)))) else: tmp = (l / (k_m * k_m)) * (l / (t * t)) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (Float64(2.0 / Float64(Float64(tan(k_m) * Float64(sin(k_m) * Float64((t ^ 3.0) / Float64(l * l)))) * Float64(1.0 + Float64(1.0 + (Float64(k_m / t) ^ 2.0))))) <= 4e+305) tmp = Float64(Float64(l * l) / Float64(k_m * Float64(k_m * Float64(t * Float64(t * t))))); else tmp = Float64(Float64(l / Float64(k_m * k_m)) * Float64(l / Float64(t * t))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if ((2.0 / ((tan(k_m) * (sin(k_m) * ((t ^ 3.0) / (l * l)))) * (1.0 + (1.0 + ((k_m / t) ^ 2.0))))) <= 4e+305) tmp = (l * l) / (k_m * (k_m * (t * (t * t)))); else tmp = (l / (k_m * k_m)) * (l / (t * t)); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(1.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 4e+305], N[(N[(l * l), $MachinePrecision] / N[(k$95$m * N[(k$95$m * N[(t * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(l / N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{2}{\left(\tan k\_m \cdot \left(\sin k\_m \cdot \frac{{t}^{3}}{\ell \cdot \ell}\right)\right) \cdot \left(1 + \left(1 + {\left(\frac{k\_m}{t}\right)}^{2}\right)\right)} \leq 4 \cdot 10^{+305}:\\
\;\;\;\;\frac{\ell \cdot \ell}{k\_m \cdot \left(k\_m \cdot \left(t \cdot \left(t \cdot t\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k\_m \cdot k\_m} \cdot \frac{\ell}{t \cdot t}\\
\end{array}
\end{array}
if (/.f64 #s(literal 2 binary64) (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t #s(literal 3 binary64)) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 k t) #s(literal 2 binary64))) #s(literal 1 binary64)))) < 3.9999999999999998e305Initial program 80.7%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6456.0
Simplified56.0%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6456.6
Simplified56.6%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6447.7
Simplified47.7%
Taylor expanded in k around -inf
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6478.0
Simplified78.0%
if 3.9999999999999998e305 < (/.f64 #s(literal 2 binary64) (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t #s(literal 3 binary64)) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 k t) #s(literal 2 binary64))) #s(literal 1 binary64)))) Initial program 16.3%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6417.6
Simplified17.6%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6422.4
Simplified22.4%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6432.3
Simplified32.3%
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6440.1
Applied egg-rr40.1%
Final simplification62.4%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<=
(/
2.0
(*
(* (tan k_m) (* (sin k_m) (/ (pow t 3.0) (* l l))))
(+ 1.0 (+ 1.0 (pow (/ k_m t) 2.0)))))
4e+305)
(/ (* l l) (* k_m (* k_m (* t (* t t)))))
(* l (/ l (* k_m (* k_m (* t t)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if ((2.0 / ((tan(k_m) * (sin(k_m) * (pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + pow((k_m / t), 2.0))))) <= 4e+305) {
tmp = (l * l) / (k_m * (k_m * (t * (t * t))));
} else {
tmp = l * (l / (k_m * (k_m * (t * 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 ((2.0d0 / ((tan(k_m) * (sin(k_m) * ((t ** 3.0d0) / (l * l)))) * (1.0d0 + (1.0d0 + ((k_m / t) ** 2.0d0))))) <= 4d+305) then
tmp = (l * l) / (k_m * (k_m * (t * (t * t))))
else
tmp = l * (l / (k_m * (k_m * (t * 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 ((2.0 / ((Math.tan(k_m) * (Math.sin(k_m) * (Math.pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + Math.pow((k_m / t), 2.0))))) <= 4e+305) {
tmp = (l * l) / (k_m * (k_m * (t * (t * t))));
} else {
tmp = l * (l / (k_m * (k_m * (t * t))));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if (2.0 / ((math.tan(k_m) * (math.sin(k_m) * (math.pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + math.pow((k_m / t), 2.0))))) <= 4e+305: tmp = (l * l) / (k_m * (k_m * (t * (t * t)))) else: tmp = l * (l / (k_m * (k_m * (t * t)))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (Float64(2.0 / Float64(Float64(tan(k_m) * Float64(sin(k_m) * Float64((t ^ 3.0) / Float64(l * l)))) * Float64(1.0 + Float64(1.0 + (Float64(k_m / t) ^ 2.0))))) <= 4e+305) tmp = Float64(Float64(l * l) / Float64(k_m * Float64(k_m * Float64(t * Float64(t * t))))); else tmp = Float64(l * Float64(l / Float64(k_m * Float64(k_m * Float64(t * t))))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if ((2.0 / ((tan(k_m) * (sin(k_m) * ((t ^ 3.0) / (l * l)))) * (1.0 + (1.0 + ((k_m / t) ^ 2.0))))) <= 4e+305) tmp = (l * l) / (k_m * (k_m * (t * (t * t)))); else tmp = l * (l / (k_m * (k_m * (t * t)))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(1.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 4e+305], N[(N[(l * l), $MachinePrecision] / N[(k$95$m * N[(k$95$m * N[(t * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l * N[(l / N[(k$95$m * N[(k$95$m * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{2}{\left(\tan k\_m \cdot \left(\sin k\_m \cdot \frac{{t}^{3}}{\ell \cdot \ell}\right)\right) \cdot \left(1 + \left(1 + {\left(\frac{k\_m}{t}\right)}^{2}\right)\right)} \leq 4 \cdot 10^{+305}:\\
\;\;\;\;\frac{\ell \cdot \ell}{k\_m \cdot \left(k\_m \cdot \left(t \cdot \left(t \cdot t\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \frac{\ell}{k\_m \cdot \left(k\_m \cdot \left(t \cdot t\right)\right)}\\
\end{array}
\end{array}
if (/.f64 #s(literal 2 binary64) (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t #s(literal 3 binary64)) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 k t) #s(literal 2 binary64))) #s(literal 1 binary64)))) < 3.9999999999999998e305Initial program 80.7%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6456.0
Simplified56.0%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6456.6
Simplified56.6%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6447.7
Simplified47.7%
Taylor expanded in k around -inf
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6478.0
Simplified78.0%
if 3.9999999999999998e305 < (/.f64 #s(literal 2 binary64) (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t #s(literal 3 binary64)) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 k t) #s(literal 2 binary64))) #s(literal 1 binary64)))) Initial program 16.3%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6417.6
Simplified17.6%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6422.4
Simplified22.4%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6432.3
Simplified32.3%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6439.5
Applied egg-rr39.5%
Final simplification62.2%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (* (sin k_m) (* t (tan k_m)))))
(if (<= k_m 1.36e+29)
(/ 2.0 (* 2.0 (* (* t (/ (* t (sin k_m)) l)) (* (/ t l) (tan k_m)))))
(if (<= k_m 3e+150)
(/
(/ 2.0 t_1)
(/ (* k_m (* k_m (fma (* t t) (/ 2.0 (* k_m k_m)) 1.0))) (* l l)))
(/
(/ 2.0 (/ t l))
(* t_1 (* (/ t l) (fma k_m (/ k_m (* t t)) 2.0))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = sin(k_m) * (t * tan(k_m));
double tmp;
if (k_m <= 1.36e+29) {
tmp = 2.0 / (2.0 * ((t * ((t * sin(k_m)) / l)) * ((t / l) * tan(k_m))));
} else if (k_m <= 3e+150) {
tmp = (2.0 / t_1) / ((k_m * (k_m * fma((t * t), (2.0 / (k_m * k_m)), 1.0))) / (l * l));
} else {
tmp = (2.0 / (t / l)) / (t_1 * ((t / l) * fma(k_m, (k_m / (t * t)), 2.0)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(sin(k_m) * Float64(t * tan(k_m))) tmp = 0.0 if (k_m <= 1.36e+29) tmp = Float64(2.0 / Float64(2.0 * Float64(Float64(t * Float64(Float64(t * sin(k_m)) / l)) * Float64(Float64(t / l) * tan(k_m))))); elseif (k_m <= 3e+150) tmp = Float64(Float64(2.0 / t_1) / Float64(Float64(k_m * Float64(k_m * fma(Float64(t * t), Float64(2.0 / Float64(k_m * k_m)), 1.0))) / Float64(l * l))); else tmp = Float64(Float64(2.0 / Float64(t / l)) / Float64(t_1 * Float64(Float64(t / l) * fma(k_m, Float64(k_m / Float64(t * t)), 2.0)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(N[Sin[k$95$m], $MachinePrecision] * N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k$95$m, 1.36e+29], N[(2.0 / N[(2.0 * N[(N[(t * N[(N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 3e+150], N[(N[(2.0 / t$95$1), $MachinePrecision] / N[(N[(k$95$m * N[(k$95$m * N[(N[(t * t), $MachinePrecision] * N[(2.0 / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(t / l), $MachinePrecision]), $MachinePrecision] / N[(t$95$1 * N[(N[(t / l), $MachinePrecision] * N[(k$95$m * N[(k$95$m / N[(t * t), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := \sin k\_m \cdot \left(t \cdot \tan k\_m\right)\\
\mathbf{if}\;k\_m \leq 1.36 \cdot 10^{+29}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\left(t \cdot \frac{t \cdot \sin k\_m}{\ell}\right) \cdot \left(\frac{t}{\ell} \cdot \tan k\_m\right)\right)}\\
\mathbf{elif}\;k\_m \leq 3 \cdot 10^{+150}:\\
\;\;\;\;\frac{\frac{2}{t\_1}}{\frac{k\_m \cdot \left(k\_m \cdot \mathsf{fma}\left(t \cdot t, \frac{2}{k\_m \cdot k\_m}, 1\right)\right)}{\ell \cdot \ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\frac{t}{\ell}}}{t\_1 \cdot \left(\frac{t}{\ell} \cdot \mathsf{fma}\left(k\_m, \frac{k\_m}{t \cdot t}, 2\right)\right)}\\
\end{array}
\end{array}
if k < 1.36e29Initial program 57.6%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6418.9
Applied egg-rr18.9%
Applied egg-rr75.5%
Taylor expanded in k around 0
Simplified75.8%
lift-/.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*l*N/A
associate-*r*N/A
lower-*.f64N/A
Applied egg-rr82.2%
if 1.36e29 < k < 3.00000000000000012e150Initial program 52.1%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6410.5
Applied egg-rr10.5%
Applied egg-rr52.3%
Taylor expanded in k around inf
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6485.8
Simplified85.8%
if 3.00000000000000012e150 < k Initial program 32.9%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6415.6
Applied egg-rr15.6%
Applied egg-rr57.1%
Applied egg-rr58.7%
Final simplification80.0%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (/ 2.0 (* (sin k_m) (* t (tan k_m))))))
(if (<= k_m 1.36e+29)
(/ 2.0 (* 2.0 (* (* t (/ (* t (sin k_m)) l)) (* (/ t l) (tan k_m)))))
(if (<= k_m 3e+150)
(/
t_1
(/ (* k_m (* k_m (fma (* t t) (/ 2.0 (* k_m k_m)) 1.0))) (* l l)))
(/ t_1 (* (/ t l) (/ (* t (fma k_m (/ k_m (* t t)) 2.0)) l)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = 2.0 / (sin(k_m) * (t * tan(k_m)));
double tmp;
if (k_m <= 1.36e+29) {
tmp = 2.0 / (2.0 * ((t * ((t * sin(k_m)) / l)) * ((t / l) * tan(k_m))));
} else if (k_m <= 3e+150) {
tmp = t_1 / ((k_m * (k_m * fma((t * t), (2.0 / (k_m * k_m)), 1.0))) / (l * l));
} else {
tmp = t_1 / ((t / l) * ((t * fma(k_m, (k_m / (t * t)), 2.0)) / l));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(2.0 / Float64(sin(k_m) * Float64(t * tan(k_m)))) tmp = 0.0 if (k_m <= 1.36e+29) tmp = Float64(2.0 / Float64(2.0 * Float64(Float64(t * Float64(Float64(t * sin(k_m)) / l)) * Float64(Float64(t / l) * tan(k_m))))); elseif (k_m <= 3e+150) tmp = Float64(t_1 / Float64(Float64(k_m * Float64(k_m * fma(Float64(t * t), Float64(2.0 / Float64(k_m * k_m)), 1.0))) / Float64(l * l))); else tmp = Float64(t_1 / Float64(Float64(t / l) * Float64(Float64(t * fma(k_m, Float64(k_m / Float64(t * t)), 2.0)) / l))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k$95$m, 1.36e+29], N[(2.0 / N[(2.0 * N[(N[(t * N[(N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 3e+150], N[(t$95$1 / N[(N[(k$95$m * N[(k$95$m * N[(N[(t * t), $MachinePrecision] * N[(2.0 / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 / N[(N[(t / l), $MachinePrecision] * N[(N[(t * N[(k$95$m * N[(k$95$m / N[(t * t), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := \frac{2}{\sin k\_m \cdot \left(t \cdot \tan k\_m\right)}\\
\mathbf{if}\;k\_m \leq 1.36 \cdot 10^{+29}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\left(t \cdot \frac{t \cdot \sin k\_m}{\ell}\right) \cdot \left(\frac{t}{\ell} \cdot \tan k\_m\right)\right)}\\
\mathbf{elif}\;k\_m \leq 3 \cdot 10^{+150}:\\
\;\;\;\;\frac{t\_1}{\frac{k\_m \cdot \left(k\_m \cdot \mathsf{fma}\left(t \cdot t, \frac{2}{k\_m \cdot k\_m}, 1\right)\right)}{\ell \cdot \ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_1}{\frac{t}{\ell} \cdot \frac{t \cdot \mathsf{fma}\left(k\_m, \frac{k\_m}{t \cdot t}, 2\right)}{\ell}}\\
\end{array}
\end{array}
if k < 1.36e29Initial program 57.6%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6418.9
Applied egg-rr18.9%
Applied egg-rr75.5%
Taylor expanded in k around 0
Simplified75.8%
lift-/.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*l*N/A
associate-*r*N/A
lower-*.f64N/A
Applied egg-rr82.2%
if 1.36e29 < k < 3.00000000000000012e150Initial program 52.1%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6410.5
Applied egg-rr10.5%
Applied egg-rr52.3%
Taylor expanded in k around inf
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6485.8
Simplified85.8%
if 3.00000000000000012e150 < k Initial program 32.9%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6415.6
Applied egg-rr15.6%
Applied egg-rr33.3%
lift-*.f64N/A
lift-/.f64N/A
lift-fma.f64N/A
associate-*l*N/A
times-fracN/A
lift-/.f64N/A
lower-*.f64N/A
lower-/.f64N/A
Applied egg-rr58.6%
Final simplification80.0%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (* t (tan k_m))))
(if (<= k_m 1.36e+29)
(/ 2.0 (* 2.0 (* (* t (/ (* t (sin k_m)) l)) (* (/ t l) (tan k_m)))))
(if (<= k_m 3e+150)
(/
(/ 2.0 (* (sin k_m) t_1))
(/ (* k_m (* k_m (fma (* t t) (/ 2.0 (* k_m k_m)) 1.0))) (* l l)))
(/
2.0
(*
(sin k_m)
(* (/ t l) (* (/ t l) (* t_1 (fma k_m (/ k_m (* t t)) 2.0))))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = t * tan(k_m);
double tmp;
if (k_m <= 1.36e+29) {
tmp = 2.0 / (2.0 * ((t * ((t * sin(k_m)) / l)) * ((t / l) * tan(k_m))));
} else if (k_m <= 3e+150) {
tmp = (2.0 / (sin(k_m) * t_1)) / ((k_m * (k_m * fma((t * t), (2.0 / (k_m * k_m)), 1.0))) / (l * l));
} else {
tmp = 2.0 / (sin(k_m) * ((t / l) * ((t / l) * (t_1 * fma(k_m, (k_m / (t * t)), 2.0)))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(t * tan(k_m)) tmp = 0.0 if (k_m <= 1.36e+29) tmp = Float64(2.0 / Float64(2.0 * Float64(Float64(t * Float64(Float64(t * sin(k_m)) / l)) * Float64(Float64(t / l) * tan(k_m))))); elseif (k_m <= 3e+150) tmp = Float64(Float64(2.0 / Float64(sin(k_m) * t_1)) / Float64(Float64(k_m * Float64(k_m * fma(Float64(t * t), Float64(2.0 / Float64(k_m * k_m)), 1.0))) / Float64(l * l))); else tmp = Float64(2.0 / Float64(sin(k_m) * Float64(Float64(t / l) * Float64(Float64(t / l) * Float64(t_1 * fma(k_m, Float64(k_m / Float64(t * t)), 2.0)))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k$95$m, 1.36e+29], N[(2.0 / N[(2.0 * N[(N[(t * N[(N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 3e+150], N[(N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * N[(k$95$m * N[(N[(t * t), $MachinePrecision] * N[(2.0 / N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(t$95$1 * N[(k$95$m * N[(k$95$m / N[(t * t), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := t \cdot \tan k\_m\\
\mathbf{if}\;k\_m \leq 1.36 \cdot 10^{+29}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\left(t \cdot \frac{t \cdot \sin k\_m}{\ell}\right) \cdot \left(\frac{t}{\ell} \cdot \tan k\_m\right)\right)}\\
\mathbf{elif}\;k\_m \leq 3 \cdot 10^{+150}:\\
\;\;\;\;\frac{\frac{2}{\sin k\_m \cdot t\_1}}{\frac{k\_m \cdot \left(k\_m \cdot \mathsf{fma}\left(t \cdot t, \frac{2}{k\_m \cdot k\_m}, 1\right)\right)}{\ell \cdot \ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left(\frac{t}{\ell} \cdot \left(\frac{t}{\ell} \cdot \left(t\_1 \cdot \mathsf{fma}\left(k\_m, \frac{k\_m}{t \cdot t}, 2\right)\right)\right)\right)}\\
\end{array}
\end{array}
if k < 1.36e29Initial program 57.6%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6418.9
Applied egg-rr18.9%
Applied egg-rr75.5%
Taylor expanded in k around 0
Simplified75.8%
lift-/.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*l*N/A
associate-*r*N/A
lower-*.f64N/A
Applied egg-rr82.2%
if 1.36e29 < k < 3.00000000000000012e150Initial program 52.1%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6410.5
Applied egg-rr10.5%
Applied egg-rr52.3%
Taylor expanded in k around inf
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6485.8
Simplified85.8%
if 3.00000000000000012e150 < k Initial program 32.9%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6415.6
Applied egg-rr15.6%
Applied egg-rr57.1%
lift-/.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f6457.2
Applied egg-rr57.2%
Applied egg-rr58.7%
Final simplification80.0%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (* t (tan k_m))))
(if (<= k_m 1.36e+29)
(/ 2.0 (* 2.0 (* (* t (/ (* t (sin k_m)) l)) (* (/ t l) (tan k_m)))))
(if (<= k_m 3e+150)
(/ (/ 2.0 (* (sin k_m) t_1)) (/ (* k_m k_m) (* l l)))
(/
2.0
(*
(sin k_m)
(* (/ t l) (* (/ t l) (* t_1 (fma k_m (/ k_m (* t t)) 2.0))))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = t * tan(k_m);
double tmp;
if (k_m <= 1.36e+29) {
tmp = 2.0 / (2.0 * ((t * ((t * sin(k_m)) / l)) * ((t / l) * tan(k_m))));
} else if (k_m <= 3e+150) {
tmp = (2.0 / (sin(k_m) * t_1)) / ((k_m * k_m) / (l * l));
} else {
tmp = 2.0 / (sin(k_m) * ((t / l) * ((t / l) * (t_1 * fma(k_m, (k_m / (t * t)), 2.0)))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(t * tan(k_m)) tmp = 0.0 if (k_m <= 1.36e+29) tmp = Float64(2.0 / Float64(2.0 * Float64(Float64(t * Float64(Float64(t * sin(k_m)) / l)) * Float64(Float64(t / l) * tan(k_m))))); elseif (k_m <= 3e+150) tmp = Float64(Float64(2.0 / Float64(sin(k_m) * t_1)) / Float64(Float64(k_m * k_m) / Float64(l * l))); else tmp = Float64(2.0 / Float64(sin(k_m) * Float64(Float64(t / l) * Float64(Float64(t / l) * Float64(t_1 * fma(k_m, Float64(k_m / Float64(t * t)), 2.0)))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k$95$m, 1.36e+29], N[(2.0 / N[(2.0 * N[(N[(t * N[(N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 3e+150], N[(N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(t$95$1 * N[(k$95$m * N[(k$95$m / N[(t * t), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := t \cdot \tan k\_m\\
\mathbf{if}\;k\_m \leq 1.36 \cdot 10^{+29}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\left(t \cdot \frac{t \cdot \sin k\_m}{\ell}\right) \cdot \left(\frac{t}{\ell} \cdot \tan k\_m\right)\right)}\\
\mathbf{elif}\;k\_m \leq 3 \cdot 10^{+150}:\\
\;\;\;\;\frac{\frac{2}{\sin k\_m \cdot t\_1}}{\frac{k\_m \cdot k\_m}{\ell \cdot \ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left(\frac{t}{\ell} \cdot \left(\frac{t}{\ell} \cdot \left(t\_1 \cdot \mathsf{fma}\left(k\_m, \frac{k\_m}{t \cdot t}, 2\right)\right)\right)\right)}\\
\end{array}
\end{array}
if k < 1.36e29Initial program 57.6%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6418.9
Applied egg-rr18.9%
Applied egg-rr75.5%
Taylor expanded in k around 0
Simplified75.8%
lift-/.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*l*N/A
associate-*r*N/A
lower-*.f64N/A
Applied egg-rr82.2%
if 1.36e29 < k < 3.00000000000000012e150Initial program 52.1%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6410.5
Applied egg-rr10.5%
Applied egg-rr52.3%
Taylor expanded in t around 0
unpow2N/A
lower-*.f6475.6
Simplified75.6%
if 3.00000000000000012e150 < k Initial program 32.9%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6415.6
Applied egg-rr15.6%
Applied egg-rr57.1%
lift-/.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f6457.2
Applied egg-rr57.2%
Applied egg-rr58.7%
Final simplification79.0%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (* t (tan k_m))))
(if (<= t 9.8e-47)
(/ (/ 2.0 (* (sin k_m) t_1)) (/ (* k_m k_m) (* l l)))
(/
(/ 2.0 (/ (* t (sin k_m)) l))
(* (/ t l) (* t_1 (fma k_m (/ k_m (* t t)) 2.0)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = t * tan(k_m);
double tmp;
if (t <= 9.8e-47) {
tmp = (2.0 / (sin(k_m) * t_1)) / ((k_m * k_m) / (l * l));
} else {
tmp = (2.0 / ((t * sin(k_m)) / l)) / ((t / l) * (t_1 * fma(k_m, (k_m / (t * t)), 2.0)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(t * tan(k_m)) tmp = 0.0 if (t <= 9.8e-47) tmp = Float64(Float64(2.0 / Float64(sin(k_m) * t_1)) / Float64(Float64(k_m * k_m) / Float64(l * l))); else tmp = Float64(Float64(2.0 / Float64(Float64(t * sin(k_m)) / l)) / Float64(Float64(t / l) * Float64(t_1 * fma(k_m, Float64(k_m / Float64(t * t)), 2.0)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, 9.8e-47], N[(N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / N[(N[(t / l), $MachinePrecision] * N[(t$95$1 * N[(k$95$m * N[(k$95$m / N[(t * t), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := t \cdot \tan k\_m\\
\mathbf{if}\;t \leq 9.8 \cdot 10^{-47}:\\
\;\;\;\;\frac{\frac{2}{\sin k\_m \cdot t\_1}}{\frac{k\_m \cdot k\_m}{\ell \cdot \ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\frac{t \cdot \sin k\_m}{\ell}}}{\frac{t}{\ell} \cdot \left(t\_1 \cdot \mathsf{fma}\left(k\_m, \frac{k\_m}{t \cdot t}, 2\right)\right)}\\
\end{array}
\end{array}
if t < 9.800000000000001e-47Initial program 54.4%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6410.4
Applied egg-rr10.4%
Applied egg-rr45.5%
Taylor expanded in t around 0
unpow2N/A
lower-*.f6465.8
Simplified65.8%
if 9.800000000000001e-47 < t Initial program 54.1%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6434.5
Applied egg-rr34.5%
Applied egg-rr83.6%
lift-/.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f6493.3
Applied egg-rr93.3%
Applied egg-rr89.8%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.95e-25)
(/
2.0
(*
2.0
(*
(/ t l)
(*
(/ t l)
(* (sin k_m) (* t (fma k_m (* (* k_m k_m) 0.3333333333333333) k_m)))))))
(if (<= k_m 1.5e+89)
(/ 2.0 (* 2.0 (* (tan k_m) (* t (* t (/ (* t (sin k_m)) (* l l)))))))
(/ (/ 2.0 (* (sin k_m) (* t (tan k_m)))) (/ (* k_m k_m) (* l l))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.95e-25) {
tmp = 2.0 / (2.0 * ((t / l) * ((t / l) * (sin(k_m) * (t * fma(k_m, ((k_m * k_m) * 0.3333333333333333), k_m))))));
} else if (k_m <= 1.5e+89) {
tmp = 2.0 / (2.0 * (tan(k_m) * (t * (t * ((t * sin(k_m)) / (l * l))))));
} else {
tmp = (2.0 / (sin(k_m) * (t * tan(k_m)))) / ((k_m * k_m) / (l * l));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.95e-25) tmp = Float64(2.0 / Float64(2.0 * Float64(Float64(t / l) * Float64(Float64(t / l) * Float64(sin(k_m) * Float64(t * fma(k_m, Float64(Float64(k_m * k_m) * 0.3333333333333333), k_m))))))); elseif (k_m <= 1.5e+89) tmp = Float64(2.0 / Float64(2.0 * Float64(tan(k_m) * Float64(t * Float64(t * Float64(Float64(t * sin(k_m)) / Float64(l * l))))))); else tmp = Float64(Float64(2.0 / Float64(sin(k_m) * Float64(t * tan(k_m)))) / Float64(Float64(k_m * k_m) / Float64(l * l))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.95e-25], N[(2.0 / N[(2.0 * N[(N[(t / l), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(t * N[(k$95$m * N[(N[(k$95$m * k$95$m), $MachinePrecision] * 0.3333333333333333), $MachinePrecision] + k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.5e+89], N[(2.0 / N[(2.0 * N[(N[Tan[k$95$m], $MachinePrecision] * N[(t * N[(t * N[(N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.95 \cdot 10^{-25}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\frac{t}{\ell} \cdot \left(\frac{t}{\ell} \cdot \left(\sin k\_m \cdot \left(t \cdot \mathsf{fma}\left(k\_m, \left(k\_m \cdot k\_m\right) \cdot 0.3333333333333333, k\_m\right)\right)\right)\right)\right)}\\
\mathbf{elif}\;k\_m \leq 1.5 \cdot 10^{+89}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\tan k\_m \cdot \left(t \cdot \left(t \cdot \frac{t \cdot \sin k\_m}{\ell \cdot \ell}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\sin k\_m \cdot \left(t \cdot \tan k\_m\right)}}{\frac{k\_m \cdot k\_m}{\ell \cdot \ell}}\\
\end{array}
\end{array}
if k < 1.95e-25Initial program 58.6%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6419.4
Applied egg-rr19.4%
Applied egg-rr76.3%
Taylor expanded in k around 0
Simplified75.9%
Taylor expanded in k around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6475.3
Simplified75.3%
if 1.95e-25 < k < 1.50000000000000006e89Initial program 47.3%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6412.4
Applied egg-rr12.4%
lift-log.f64N/A
lift-*.f64N/A
lift-log.f64N/A
lift-*.f64N/A
div-expN/A
lift-*.f64N/A
lift-log.f64N/A
pow-to-expN/A
pow2N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-log.f64N/A
exp-to-powN/A
pow3N/A
lift-*.f64N/A
associate-*r/N/A
lift-/.f64N/A
lift-sin.f64N/A
Applied egg-rr54.6%
Taylor expanded in k around 0
Simplified71.2%
if 1.50000000000000006e89 < k Initial program 39.1%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6413.8
Applied egg-rr13.8%
Applied egg-rr39.5%
Taylor expanded in t around 0
unpow2N/A
lower-*.f6455.1
Simplified55.1%
Final simplification71.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 1.95e-25)
(/
2.0
(*
2.0
(*
(/ t l)
(*
(/ t l)
(* (sin k_m) (* t (fma k_m (* (* k_m k_m) 0.3333333333333333) k_m)))))))
(if (<= k_m 3.2e+145)
(/ 2.0 (* 2.0 (* (tan k_m) (* t (* t (/ (* t (sin k_m)) (* l l)))))))
(/
2.0
(*
t
(*
(/ 1.0 l)
(*
(* (/ t l) (fma k_m (/ k_m (* t t)) 2.0))
(*
(* t (tan k_m))
(fma k_m (* k_m (* k_m -0.16666666666666666)) k_m)))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.95e-25) {
tmp = 2.0 / (2.0 * ((t / l) * ((t / l) * (sin(k_m) * (t * fma(k_m, ((k_m * k_m) * 0.3333333333333333), k_m))))));
} else if (k_m <= 3.2e+145) {
tmp = 2.0 / (2.0 * (tan(k_m) * (t * (t * ((t * sin(k_m)) / (l * l))))));
} else {
tmp = 2.0 / (t * ((1.0 / l) * (((t / l) * fma(k_m, (k_m / (t * t)), 2.0)) * ((t * tan(k_m)) * fma(k_m, (k_m * (k_m * -0.16666666666666666)), k_m)))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.95e-25) tmp = Float64(2.0 / Float64(2.0 * Float64(Float64(t / l) * Float64(Float64(t / l) * Float64(sin(k_m) * Float64(t * fma(k_m, Float64(Float64(k_m * k_m) * 0.3333333333333333), k_m))))))); elseif (k_m <= 3.2e+145) tmp = Float64(2.0 / Float64(2.0 * Float64(tan(k_m) * Float64(t * Float64(t * Float64(Float64(t * sin(k_m)) / Float64(l * l))))))); else tmp = Float64(2.0 / Float64(t * Float64(Float64(1.0 / l) * Float64(Float64(Float64(t / l) * fma(k_m, Float64(k_m / Float64(t * t)), 2.0)) * Float64(Float64(t * tan(k_m)) * fma(k_m, Float64(k_m * Float64(k_m * -0.16666666666666666)), k_m)))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.95e-25], N[(2.0 / N[(2.0 * N[(N[(t / l), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(t * N[(k$95$m * N[(N[(k$95$m * k$95$m), $MachinePrecision] * 0.3333333333333333), $MachinePrecision] + k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 3.2e+145], N[(2.0 / N[(2.0 * N[(N[Tan[k$95$m], $MachinePrecision] * N[(t * N[(t * N[(N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(t * N[(N[(1.0 / l), $MachinePrecision] * N[(N[(N[(t / l), $MachinePrecision] * N[(k$95$m * N[(k$95$m / N[(t * t), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] * N[(N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(k$95$m * N[(k$95$m * N[(k$95$m * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.95 \cdot 10^{-25}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\frac{t}{\ell} \cdot \left(\frac{t}{\ell} \cdot \left(\sin k\_m \cdot \left(t \cdot \mathsf{fma}\left(k\_m, \left(k\_m \cdot k\_m\right) \cdot 0.3333333333333333, k\_m\right)\right)\right)\right)\right)}\\
\mathbf{elif}\;k\_m \leq 3.2 \cdot 10^{+145}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\tan k\_m \cdot \left(t \cdot \left(t \cdot \frac{t \cdot \sin k\_m}{\ell \cdot \ell}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t \cdot \left(\frac{1}{\ell} \cdot \left(\left(\frac{t}{\ell} \cdot \mathsf{fma}\left(k\_m, \frac{k\_m}{t \cdot t}, 2\right)\right) \cdot \left(\left(t \cdot \tan k\_m\right) \cdot \mathsf{fma}\left(k\_m, k\_m \cdot \left(k\_m \cdot -0.16666666666666666\right), k\_m\right)\right)\right)\right)}\\
\end{array}
\end{array}
if k < 1.95e-25Initial program 58.6%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6419.4
Applied egg-rr19.4%
Applied egg-rr76.3%
Taylor expanded in k around 0
Simplified75.9%
Taylor expanded in k around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6475.3
Simplified75.3%
if 1.95e-25 < k < 3.20000000000000008e145Initial program 49.3%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6412.1
Applied egg-rr12.1%
lift-log.f64N/A
lift-*.f64N/A
lift-log.f64N/A
lift-*.f64N/A
div-expN/A
lift-*.f64N/A
lift-log.f64N/A
pow-to-expN/A
pow2N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-log.f64N/A
exp-to-powN/A
pow3N/A
lift-*.f64N/A
associate-*r/N/A
lift-/.f64N/A
lift-sin.f64N/A
Applied egg-rr59.9%
Taylor expanded in k around 0
Simplified75.3%
if 3.20000000000000008e145 < k Initial program 34.0%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6414.6
Applied egg-rr14.6%
Applied egg-rr56.7%
Applied egg-rr58.2%
Taylor expanded in k around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f6442.5
Simplified42.5%
Final simplification71.4%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 1.36e+29) (/ 2.0 (* 2.0 (* (* t (/ (* t (sin k_m)) l)) (* (/ t l) (tan k_m))))) (/ (/ 2.0 (* (sin k_m) (* t (tan k_m)))) (/ (* k_m k_m) (* l l)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.36e+29) {
tmp = 2.0 / (2.0 * ((t * ((t * sin(k_m)) / l)) * ((t / l) * tan(k_m))));
} else {
tmp = (2.0 / (sin(k_m) * (t * tan(k_m)))) / ((k_m * k_m) / (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 <= 1.36d+29) then
tmp = 2.0d0 / (2.0d0 * ((t * ((t * sin(k_m)) / l)) * ((t / l) * tan(k_m))))
else
tmp = (2.0d0 / (sin(k_m) * (t * tan(k_m)))) / ((k_m * k_m) / (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 <= 1.36e+29) {
tmp = 2.0 / (2.0 * ((t * ((t * Math.sin(k_m)) / l)) * ((t / l) * Math.tan(k_m))));
} else {
tmp = (2.0 / (Math.sin(k_m) * (t * Math.tan(k_m)))) / ((k_m * k_m) / (l * l));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 1.36e+29: tmp = 2.0 / (2.0 * ((t * ((t * math.sin(k_m)) / l)) * ((t / l) * math.tan(k_m)))) else: tmp = (2.0 / (math.sin(k_m) * (t * math.tan(k_m)))) / ((k_m * k_m) / (l * l)) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.36e+29) tmp = Float64(2.0 / Float64(2.0 * Float64(Float64(t * Float64(Float64(t * sin(k_m)) / l)) * Float64(Float64(t / l) * tan(k_m))))); else tmp = Float64(Float64(2.0 / Float64(sin(k_m) * Float64(t * tan(k_m)))) / Float64(Float64(k_m * k_m) / 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 <= 1.36e+29) tmp = 2.0 / (2.0 * ((t * ((t * sin(k_m)) / l)) * ((t / l) * tan(k_m)))); else tmp = (2.0 / (sin(k_m) * (t * tan(k_m)))) / ((k_m * k_m) / (l * l)); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.36e+29], N[(2.0 / N[(2.0 * N[(N[(t * N[(N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.36 \cdot 10^{+29}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\left(t \cdot \frac{t \cdot \sin k\_m}{\ell}\right) \cdot \left(\frac{t}{\ell} \cdot \tan k\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\sin k\_m \cdot \left(t \cdot \tan k\_m\right)}}{\frac{k\_m \cdot k\_m}{\ell \cdot \ell}}\\
\end{array}
\end{array}
if k < 1.36e29Initial program 57.6%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6418.9
Applied egg-rr18.9%
Applied egg-rr75.5%
Taylor expanded in k around 0
Simplified75.8%
lift-/.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
associate-*l*N/A
associate-*r*N/A
lower-*.f64N/A
Applied egg-rr82.2%
if 1.36e29 < k Initial program 42.3%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6413.1
Applied egg-rr13.1%
Applied egg-rr42.7%
Taylor expanded in t around 0
unpow2N/A
lower-*.f6457.9
Simplified57.9%
Final simplification77.0%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 1.36e+29) (/ 2.0 (* 2.0 (* t (* (/ (* t (sin k_m)) l) (* (/ t l) (tan k_m)))))) (/ (/ 2.0 (* (sin k_m) (* t (tan k_m)))) (/ (* k_m k_m) (* l l)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.36e+29) {
tmp = 2.0 / (2.0 * (t * (((t * sin(k_m)) / l) * ((t / l) * tan(k_m)))));
} else {
tmp = (2.0 / (sin(k_m) * (t * tan(k_m)))) / ((k_m * k_m) / (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 <= 1.36d+29) then
tmp = 2.0d0 / (2.0d0 * (t * (((t * sin(k_m)) / l) * ((t / l) * tan(k_m)))))
else
tmp = (2.0d0 / (sin(k_m) * (t * tan(k_m)))) / ((k_m * k_m) / (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 <= 1.36e+29) {
tmp = 2.0 / (2.0 * (t * (((t * Math.sin(k_m)) / l) * ((t / l) * Math.tan(k_m)))));
} else {
tmp = (2.0 / (Math.sin(k_m) * (t * Math.tan(k_m)))) / ((k_m * k_m) / (l * l));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 1.36e+29: tmp = 2.0 / (2.0 * (t * (((t * math.sin(k_m)) / l) * ((t / l) * math.tan(k_m))))) else: tmp = (2.0 / (math.sin(k_m) * (t * math.tan(k_m)))) / ((k_m * k_m) / (l * l)) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.36e+29) tmp = Float64(2.0 / Float64(2.0 * Float64(t * Float64(Float64(Float64(t * sin(k_m)) / l) * Float64(Float64(t / l) * tan(k_m)))))); else tmp = Float64(Float64(2.0 / Float64(sin(k_m) * Float64(t * tan(k_m)))) / Float64(Float64(k_m * k_m) / 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 <= 1.36e+29) tmp = 2.0 / (2.0 * (t * (((t * sin(k_m)) / l) * ((t / l) * tan(k_m))))); else tmp = (2.0 / (sin(k_m) * (t * tan(k_m)))) / ((k_m * k_m) / (l * l)); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.36e+29], N[(2.0 / N[(2.0 * N[(t * N[(N[(N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.36 \cdot 10^{+29}:\\
\;\;\;\;\frac{2}{2 \cdot \left(t \cdot \left(\frac{t \cdot \sin k\_m}{\ell} \cdot \left(\frac{t}{\ell} \cdot \tan k\_m\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\sin k\_m \cdot \left(t \cdot \tan k\_m\right)}}{\frac{k\_m \cdot k\_m}{\ell \cdot \ell}}\\
\end{array}
\end{array}
if k < 1.36e29Initial program 57.6%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6418.9
Applied egg-rr18.9%
Applied egg-rr75.5%
Taylor expanded in k around 0
Simplified75.8%
lift-/.f64N/A
lift-/.f64N/A
lift-sin.f64N/A
lift-tan.f64N/A
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
associate-*l*N/A
Applied egg-rr81.4%
if 1.36e29 < k Initial program 42.3%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6413.1
Applied egg-rr13.1%
Applied egg-rr42.7%
Taylor expanded in t around 0
unpow2N/A
lower-*.f6457.9
Simplified57.9%
Final simplification76.3%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(/
2.0
(*
2.0
(*
(/ t l)
(*
(/ t l)
(*
(* t (tan k_m))
(fma k_m (* k_m (* k_m -0.16666666666666666)) k_m)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / (2.0 * ((t / l) * ((t / l) * ((t * tan(k_m)) * fma(k_m, (k_m * (k_m * -0.16666666666666666)), k_m)))));
}
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64(2.0 * Float64(Float64(t / l) * Float64(Float64(t / l) * Float64(Float64(t * tan(k_m)) * fma(k_m, Float64(k_m * Float64(k_m * -0.16666666666666666)), k_m)))))) end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(2.0 * N[(N[(t / l), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(N[(t * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(k$95$m * N[(k$95$m * N[(k$95$m * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{2 \cdot \left(\frac{t}{\ell} \cdot \left(\frac{t}{\ell} \cdot \left(\left(t \cdot \tan k\_m\right) \cdot \mathsf{fma}\left(k\_m, k\_m \cdot \left(k\_m \cdot -0.16666666666666666\right), k\_m\right)\right)\right)\right)}
\end{array}
Initial program 54.3%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6417.7
Applied egg-rr17.7%
Applied egg-rr71.8%
Taylor expanded in k around 0
Simplified70.0%
Taylor expanded in k around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f6469.2
Simplified69.2%
Final simplification69.2%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(/
2.0
(*
2.0
(*
(/ t l)
(*
(/ t l)
(* (sin k_m) (* t (fma k_m (* (* k_m k_m) 0.3333333333333333) k_m))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / (2.0 * ((t / l) * ((t / l) * (sin(k_m) * (t * fma(k_m, ((k_m * k_m) * 0.3333333333333333), k_m))))));
}
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64(2.0 * Float64(Float64(t / l) * Float64(Float64(t / l) * Float64(sin(k_m) * Float64(t * fma(k_m, Float64(Float64(k_m * k_m) * 0.3333333333333333), k_m))))))) end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(2.0 * N[(N[(t / l), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(t * N[(k$95$m * N[(N[(k$95$m * k$95$m), $MachinePrecision] * 0.3333333333333333), $MachinePrecision] + k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{2 \cdot \left(\frac{t}{\ell} \cdot \left(\frac{t}{\ell} \cdot \left(\sin k\_m \cdot \left(t \cdot \mathsf{fma}\left(k\_m, \left(k\_m \cdot k\_m\right) \cdot 0.3333333333333333, k\_m\right)\right)\right)\right)\right)}
\end{array}
Initial program 54.3%
pow-to-expN/A
pow2N/A
pow-to-expN/A
div-expN/A
lower-exp.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-log.f64N/A
lower-*.f64N/A
lower-log.f6417.7
Applied egg-rr17.7%
Applied egg-rr71.8%
Taylor expanded in k around 0
Simplified70.0%
Taylor expanded in k around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.8
Simplified68.8%
Final simplification68.8%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 1.75e-24) (* l (/ l (* k_m (* k_m (* t t))))) (/ (* l l) (* t (* t (* k_m k_m))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.75e-24) {
tmp = l * (l / (k_m * (k_m * (t * t))));
} else {
tmp = (l * l) / (t * (t * (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 (k_m <= 1.75d-24) then
tmp = l * (l / (k_m * (k_m * (t * t))))
else
tmp = (l * l) / (t * (t * (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 (k_m <= 1.75e-24) {
tmp = l * (l / (k_m * (k_m * (t * t))));
} else {
tmp = (l * l) / (t * (t * (k_m * k_m)));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 1.75e-24: tmp = l * (l / (k_m * (k_m * (t * t)))) else: tmp = (l * l) / (t * (t * (k_m * k_m))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.75e-24) tmp = Float64(l * Float64(l / Float64(k_m * Float64(k_m * Float64(t * t))))); else tmp = Float64(Float64(l * l) / Float64(t * Float64(t * 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 (k_m <= 1.75e-24) tmp = l * (l / (k_m * (k_m * (t * t)))); else tmp = (l * l) / (t * (t * (k_m * k_m))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.75e-24], N[(l * N[(l / N[(k$95$m * N[(k$95$m * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] / N[(t * N[(t * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.75 \cdot 10^{-24}:\\
\;\;\;\;\ell \cdot \frac{\ell}{k\_m \cdot \left(k\_m \cdot \left(t \cdot t\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \ell}{t \cdot \left(t \cdot \left(k\_m \cdot k\_m\right)\right)}\\
\end{array}
\end{array}
if k < 1.7499999999999998e-24Initial program 58.6%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6440.0
Simplified40.0%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6441.8
Simplified41.8%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6440.0
Simplified40.0%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6443.3
Applied egg-rr43.3%
if 1.7499999999999998e-24 < k Initial program 42.8%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6440.8
Simplified40.8%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6444.6
Simplified44.6%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6445.1
Simplified45.1%
lift-*.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6449.5
Applied egg-rr49.5%
Final simplification45.0%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 1.75e-24) (* l (/ l (* k_m (* k_m (* t t))))) (/ (* l l) (* k_m (* t (* k_m t))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.75e-24) {
tmp = l * (l / (k_m * (k_m * (t * t))));
} else {
tmp = (l * l) / (k_m * (t * (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.75d-24) then
tmp = l * (l / (k_m * (k_m * (t * t))))
else
tmp = (l * l) / (k_m * (t * (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.75e-24) {
tmp = l * (l / (k_m * (k_m * (t * t))));
} else {
tmp = (l * l) / (k_m * (t * (k_m * t)));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 1.75e-24: tmp = l * (l / (k_m * (k_m * (t * t)))) else: tmp = (l * l) / (k_m * (t * (k_m * t))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.75e-24) tmp = Float64(l * Float64(l / Float64(k_m * Float64(k_m * Float64(t * t))))); else tmp = Float64(Float64(l * l) / Float64(k_m * Float64(t * 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.75e-24) tmp = l * (l / (k_m * (k_m * (t * t)))); else tmp = (l * l) / (k_m * (t * (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.75e-24], N[(l * N[(l / N[(k$95$m * N[(k$95$m * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] / N[(k$95$m * N[(t * 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.75 \cdot 10^{-24}:\\
\;\;\;\;\ell \cdot \frac{\ell}{k\_m \cdot \left(k\_m \cdot \left(t \cdot t\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \ell}{k\_m \cdot \left(t \cdot \left(k\_m \cdot t\right)\right)}\\
\end{array}
\end{array}
if k < 1.7499999999999998e-24Initial program 58.6%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6440.0
Simplified40.0%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6441.8
Simplified41.8%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6440.0
Simplified40.0%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6443.3
Applied egg-rr43.3%
if 1.7499999999999998e-24 < k Initial program 42.8%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6440.8
Simplified40.8%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6444.6
Simplified44.6%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6445.1
Simplified45.1%
associate-*r*N/A
lower-*.f64N/A
lower-*.f6448.0
Applied egg-rr48.0%
Final simplification44.6%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (* l (/ l (* k_m (* k_m (* t t))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return l * (l / (k_m * (k_m * (t * 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 = l * (l / (k_m * (k_m * (t * t))))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return l * (l / (k_m * (k_m * (t * t))));
}
k_m = math.fabs(k) def code(t, l, k_m): return l * (l / (k_m * (k_m * (t * t))))
k_m = abs(k) function code(t, l, k_m) return Float64(l * Float64(l / Float64(k_m * Float64(k_m * Float64(t * t))))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = l * (l / (k_m * (k_m * (t * t)))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(l * N[(l / N[(k$95$m * N[(k$95$m * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\ell \cdot \frac{\ell}{k\_m \cdot \left(k\_m \cdot \left(t \cdot t\right)\right)}
\end{array}
Initial program 54.3%
Taylor expanded in k around 0
lower-/.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-pow.f6440.3
Simplified40.3%
Taylor expanded in t around inf
unpow2N/A
lower-*.f6442.6
Simplified42.6%
Taylor expanded in l around inf
unpow2N/A
lower-*.f6441.4
Simplified41.4%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6443.4
Applied egg-rr43.4%
Final simplification43.4%
herbie shell --seed 2024214
(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))))