
(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 4.4e-130)
(/ 2.0 (* (/ (* (* k_m 2.0) t) (/ l t)) (/ k_m (/ l t))))
(if (<= k_m 3.4e+78)
(/
2.0
(*
(/ t l)
(* (* (tan k_m) (/ (sin k_m) l)) (fma (* t t) 2.0 (* k_m k_m)))))
(*
(/ (* (cos k_m) l) k_m)
(/ (* l 2.0) (* (* (pow (sin k_m) 2.0) t) k_m))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 4.4e-130) {
tmp = 2.0 / ((((k_m * 2.0) * t) / (l / t)) * (k_m / (l / t)));
} else if (k_m <= 3.4e+78) {
tmp = 2.0 / ((t / l) * ((tan(k_m) * (sin(k_m) / l)) * fma((t * t), 2.0, (k_m * k_m))));
} else {
tmp = ((cos(k_m) * l) / k_m) * ((l * 2.0) / ((pow(sin(k_m), 2.0) * t) * k_m));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 4.4e-130) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(k_m * 2.0) * t) / Float64(l / t)) * Float64(k_m / Float64(l / t)))); elseif (k_m <= 3.4e+78) tmp = Float64(2.0 / Float64(Float64(t / l) * Float64(Float64(tan(k_m) * Float64(sin(k_m) / l)) * fma(Float64(t * t), 2.0, Float64(k_m * k_m))))); else tmp = Float64(Float64(Float64(cos(k_m) * l) / k_m) * Float64(Float64(l * 2.0) / Float64(Float64((sin(k_m) ^ 2.0) * t) * k_m))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 4.4e-130], N[(2.0 / N[(N[(N[(N[(k$95$m * 2.0), $MachinePrecision] * t), $MachinePrecision] / N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 3.4e+78], N[(2.0 / N[(N[(t / l), $MachinePrecision] * N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] / k$95$m), $MachinePrecision] * N[(N[(l * 2.0), $MachinePrecision] / N[(N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * t), $MachinePrecision] * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 4.4 \cdot 10^{-130}:\\
\;\;\;\;\frac{2}{\frac{\left(k\_m \cdot 2\right) \cdot t}{\frac{\ell}{t}} \cdot \frac{k\_m}{\frac{\ell}{t}}}\\
\mathbf{elif}\;k\_m \leq 3.4 \cdot 10^{+78}:\\
\;\;\;\;\frac{2}{\frac{t}{\ell} \cdot \left(\left(\tan k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\cos k\_m \cdot \ell}{k\_m} \cdot \frac{\ell \cdot 2}{\left({\sin k\_m}^{2} \cdot t\right) \cdot k\_m}\\
\end{array}
\end{array}
if k < 4.3999999999999997e-130Initial program 54.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6454.8
Applied rewrites54.8%
Applied rewrites72.6%
Applied rewrites75.4%
if 4.3999999999999997e-130 < k < 3.40000000000000007e78Initial program 58.6%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites83.9%
Applied rewrites94.2%
Applied rewrites96.0%
if 3.40000000000000007e78 < k Initial program 60.3%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites69.9%
Applied rewrites71.7%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites77.0%
Applied rewrites94.5%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<=
(/
2.0
(*
(* (* (/ (pow t 3.0) (* l l)) (sin k_m)) (tan k_m))
(+ (+ 1.0 (pow (/ k_m t) 2.0)) 1.0)))
INFINITY)
(/ 2.0 (* (* (* (/ (* k_m 2.0) l) (* t t)) (/ k_m l)) t))
(/ 2.0 (* (* (* k_m k_m) 2.0) (* (/ t l) (* (/ t l) t))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if ((2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k_m)) * tan(k_m)) * ((1.0 + pow((k_m / t), 2.0)) + 1.0))) <= ((double) INFINITY)) {
tmp = 2.0 / (((((k_m * 2.0) / l) * (t * t)) * (k_m / l)) * t);
} else {
tmp = 2.0 / (((k_m * k_m) * 2.0) * ((t / l) * ((t / l) * t)));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if ((2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k_m)) * Math.tan(k_m)) * ((1.0 + Math.pow((k_m / t), 2.0)) + 1.0))) <= Double.POSITIVE_INFINITY) {
tmp = 2.0 / (((((k_m * 2.0) / l) * (t * t)) * (k_m / l)) * t);
} else {
tmp = 2.0 / (((k_m * k_m) * 2.0) * ((t / l) * ((t / l) * t)));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if (2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k_m)) * math.tan(k_m)) * ((1.0 + math.pow((k_m / t), 2.0)) + 1.0))) <= math.inf: tmp = 2.0 / (((((k_m * 2.0) / l) * (t * t)) * (k_m / l)) * t) else: tmp = 2.0 / (((k_m * k_m) * 2.0) * ((t / l) * ((t / l) * t))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k_m)) * tan(k_m)) * Float64(Float64(1.0 + (Float64(k_m / t) ^ 2.0)) + 1.0))) <= Inf) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(Float64(k_m * 2.0) / l) * Float64(t * t)) * Float64(k_m / l)) * t)); else tmp = Float64(2.0 / Float64(Float64(Float64(k_m * k_m) * 2.0) * Float64(Float64(t / l) * Float64(Float64(t / l) * t)))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if ((2.0 / (((((t ^ 3.0) / (l * l)) * sin(k_m)) * tan(k_m)) * ((1.0 + ((k_m / t) ^ 2.0)) + 1.0))) <= Inf) tmp = 2.0 / (((((k_m * 2.0) / l) * (t * t)) * (k_m / l)) * t); else tmp = 2.0 / (((k_m * k_m) * 2.0) * ((t / l) * ((t / l) * 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[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(2.0 / N[(N[(N[(N[(N[(k$95$m * 2.0), $MachinePrecision] / l), $MachinePrecision] * N[(t * t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(k$95$m * k$95$m), $MachinePrecision] * 2.0), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\_m\right) \cdot \tan k\_m\right) \cdot \left(\left(1 + {\left(\frac{k\_m}{t}\right)}^{2}\right) + 1\right)} \leq \infty:\\
\;\;\;\;\frac{2}{\left(\left(\frac{k\_m \cdot 2}{\ell} \cdot \left(t \cdot t\right)\right) \cdot \frac{k\_m}{\ell}\right) \cdot t}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(k\_m \cdot k\_m\right) \cdot 2\right) \cdot \left(\frac{t}{\ell} \cdot \left(\frac{t}{\ell} \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)))) < +inf.0Initial program 84.1%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6474.7
Applied rewrites74.7%
Applied rewrites74.7%
Applied rewrites84.2%
if +inf.0 < (/.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 0.0%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6423.1
Applied rewrites23.1%
Applied rewrites19.7%
Applied rewrites40.8%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<=
(*
(* (* (/ (pow t 3.0) (* l l)) (sin k_m)) (tan k_m))
(+ (+ 1.0 (pow (/ k_m t) 2.0)) 1.0))
INFINITY)
(/ 2.0 (* k_m (* (/ t l) (* (/ (* k_m 2.0) l) (* t t)))))
(/ 2.0 (* (* (* k_m k_m) 2.0) (* (/ t l) (* (/ t l) t))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (((((pow(t, 3.0) / (l * l)) * sin(k_m)) * tan(k_m)) * ((1.0 + pow((k_m / t), 2.0)) + 1.0)) <= ((double) INFINITY)) {
tmp = 2.0 / (k_m * ((t / l) * (((k_m * 2.0) / l) * (t * t))));
} else {
tmp = 2.0 / (((k_m * k_m) * 2.0) * ((t / l) * ((t / l) * t)));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k_m)) * Math.tan(k_m)) * ((1.0 + Math.pow((k_m / t), 2.0)) + 1.0)) <= Double.POSITIVE_INFINITY) {
tmp = 2.0 / (k_m * ((t / l) * (((k_m * 2.0) / l) * (t * t))));
} else {
tmp = 2.0 / (((k_m * k_m) * 2.0) * ((t / l) * ((t / l) * t)));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if ((((math.pow(t, 3.0) / (l * l)) * math.sin(k_m)) * math.tan(k_m)) * ((1.0 + math.pow((k_m / t), 2.0)) + 1.0)) <= math.inf: tmp = 2.0 / (k_m * ((t / l) * (((k_m * 2.0) / l) * (t * t)))) else: tmp = 2.0 / (((k_m * k_m) * 2.0) * ((t / l) * ((t / l) * t))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k_m)) * tan(k_m)) * Float64(Float64(1.0 + (Float64(k_m / t) ^ 2.0)) + 1.0)) <= Inf) tmp = Float64(2.0 / Float64(k_m * Float64(Float64(t / l) * Float64(Float64(Float64(k_m * 2.0) / l) * Float64(t * t))))); else tmp = Float64(2.0 / Float64(Float64(Float64(k_m * k_m) * 2.0) * Float64(Float64(t / l) * Float64(Float64(t / l) * t)))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if ((((((t ^ 3.0) / (l * l)) * sin(k_m)) * tan(k_m)) * ((1.0 + ((k_m / t) ^ 2.0)) + 1.0)) <= Inf) tmp = 2.0 / (k_m * ((t / l) * (((k_m * 2.0) / l) * (t * t)))); else tmp = 2.0 / (((k_m * k_m) * 2.0) * ((t / l) * ((t / l) * t))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], Infinity], N[(2.0 / N[(k$95$m * N[(N[(t / l), $MachinePrecision] * N[(N[(N[(k$95$m * 2.0), $MachinePrecision] / l), $MachinePrecision] * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(k$95$m * k$95$m), $MachinePrecision] * 2.0), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\_m\right) \cdot \tan k\_m\right) \cdot \left(\left(1 + {\left(\frac{k\_m}{t}\right)}^{2}\right) + 1\right) \leq \infty:\\
\;\;\;\;\frac{2}{k\_m \cdot \left(\frac{t}{\ell} \cdot \left(\frac{k\_m \cdot 2}{\ell} \cdot \left(t \cdot t\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(k\_m \cdot k\_m\right) \cdot 2\right) \cdot \left(\frac{t}{\ell} \cdot \left(\frac{t}{\ell} \cdot t\right)\right)}\\
\end{array}
\end{array}
if (*.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))) < +inf.0Initial program 84.1%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6474.7
Applied rewrites74.7%
Applied rewrites74.7%
Applied rewrites82.1%
if +inf.0 < (*.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 0.0%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6423.1
Applied rewrites23.1%
Applied rewrites19.7%
Applied rewrites40.8%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 4.4e-130)
(/ 2.0 (* (/ (* (* k_m 2.0) t) (/ l t)) (/ k_m (/ l t))))
(if (<= k_m 9.6e+126)
(/
2.0
(*
(/ t l)
(* (* (tan k_m) (/ (sin k_m) l)) (fma (* t t) 2.0 (* k_m k_m)))))
(/
(* (* (* l l) (cos k_m)) 2.0)
(* (* (* (pow (sin k_m) 2.0) t) k_m) k_m)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 4.4e-130) {
tmp = 2.0 / ((((k_m * 2.0) * t) / (l / t)) * (k_m / (l / t)));
} else if (k_m <= 9.6e+126) {
tmp = 2.0 / ((t / l) * ((tan(k_m) * (sin(k_m) / l)) * fma((t * t), 2.0, (k_m * k_m))));
} else {
tmp = (((l * l) * cos(k_m)) * 2.0) / (((pow(sin(k_m), 2.0) * t) * k_m) * k_m);
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 4.4e-130) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(k_m * 2.0) * t) / Float64(l / t)) * Float64(k_m / Float64(l / t)))); elseif (k_m <= 9.6e+126) tmp = Float64(2.0 / Float64(Float64(t / l) * Float64(Float64(tan(k_m) * Float64(sin(k_m) / l)) * fma(Float64(t * t), 2.0, Float64(k_m * k_m))))); else tmp = Float64(Float64(Float64(Float64(l * l) * cos(k_m)) * 2.0) / Float64(Float64(Float64((sin(k_m) ^ 2.0) * t) * k_m) * k_m)); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 4.4e-130], N[(2.0 / N[(N[(N[(N[(k$95$m * 2.0), $MachinePrecision] * t), $MachinePrecision] / N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 9.6e+126], N[(2.0 / N[(N[(t / l), $MachinePrecision] * N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(l * l), $MachinePrecision] * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision] / N[(N[(N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 4.4 \cdot 10^{-130}:\\
\;\;\;\;\frac{2}{\frac{\left(k\_m \cdot 2\right) \cdot t}{\frac{\ell}{t}} \cdot \frac{k\_m}{\frac{\ell}{t}}}\\
\mathbf{elif}\;k\_m \leq 9.6 \cdot 10^{+126}:\\
\;\;\;\;\frac{2}{\frac{t}{\ell} \cdot \left(\left(\tan k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\left(\ell \cdot \ell\right) \cdot \cos k\_m\right) \cdot 2}{\left(\left({\sin k\_m}^{2} \cdot t\right) \cdot k\_m\right) \cdot k\_m}\\
\end{array}
\end{array}
if k < 4.3999999999999997e-130Initial program 54.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6454.8
Applied rewrites54.8%
Applied rewrites72.6%
Applied rewrites75.4%
if 4.3999999999999997e-130 < k < 9.60000000000000048e126Initial program 60.4%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites82.5%
Applied rewrites92.5%
Applied rewrites95.3%
if 9.60000000000000048e126 < k Initial program 58.2%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites68.1%
Applied rewrites68.1%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites77.2%
Applied rewrites77.2%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (fma (* t t) 2.0 (* k_m k_m))))
(if (<= k_m 4.4e-130)
(/ 2.0 (* (/ (* (* k_m 2.0) t) (/ l t)) (/ k_m (/ l t))))
(if (<= k_m 0.125)
(/
2.0
(*
t_1
(*
(*
(fma
(fma
(*
(fma
0.03432539682539683
(/ (* k_m k_m) l)
(/ 0.08611111111111111 l))
k_m)
k_m
(/ 0.16666666666666666 l))
(* k_m k_m)
(pow l -1.0))
(* k_m k_m))
(/ t l))))
(if (<= k_m 3.8e+147)
(/ 2.0 (* t_1 (/ (* (* (tan k_m) (sin k_m)) t) (* l l))))
(/
(* (* (* (cos k_m) l) l) 2.0)
(* (* (* (- 0.5 (* 0.5 (cos (+ k_m k_m)))) t) k_m) k_m)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = fma((t * t), 2.0, (k_m * k_m));
double tmp;
if (k_m <= 4.4e-130) {
tmp = 2.0 / ((((k_m * 2.0) * t) / (l / t)) * (k_m / (l / t)));
} else if (k_m <= 0.125) {
tmp = 2.0 / (t_1 * ((fma(fma((fma(0.03432539682539683, ((k_m * k_m) / l), (0.08611111111111111 / l)) * k_m), k_m, (0.16666666666666666 / l)), (k_m * k_m), pow(l, -1.0)) * (k_m * k_m)) * (t / l)));
} else if (k_m <= 3.8e+147) {
tmp = 2.0 / (t_1 * (((tan(k_m) * sin(k_m)) * t) / (l * l)));
} else {
tmp = (((cos(k_m) * l) * l) * 2.0) / ((((0.5 - (0.5 * cos((k_m + k_m)))) * t) * k_m) * k_m);
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = fma(Float64(t * t), 2.0, Float64(k_m * k_m)) tmp = 0.0 if (k_m <= 4.4e-130) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(k_m * 2.0) * t) / Float64(l / t)) * Float64(k_m / Float64(l / t)))); elseif (k_m <= 0.125) tmp = Float64(2.0 / Float64(t_1 * Float64(Float64(fma(fma(Float64(fma(0.03432539682539683, Float64(Float64(k_m * k_m) / l), Float64(0.08611111111111111 / l)) * k_m), k_m, Float64(0.16666666666666666 / l)), Float64(k_m * k_m), (l ^ -1.0)) * Float64(k_m * k_m)) * Float64(t / l)))); elseif (k_m <= 3.8e+147) tmp = Float64(2.0 / Float64(t_1 * Float64(Float64(Float64(tan(k_m) * sin(k_m)) * t) / Float64(l * l)))); else tmp = Float64(Float64(Float64(Float64(cos(k_m) * l) * l) * 2.0) / Float64(Float64(Float64(Float64(0.5 - Float64(0.5 * cos(Float64(k_m + k_m)))) * t) * k_m) * k_m)); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k$95$m, 4.4e-130], N[(2.0 / N[(N[(N[(N[(k$95$m * 2.0), $MachinePrecision] * t), $MachinePrecision] / N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 0.125], N[(2.0 / N[(t$95$1 * N[(N[(N[(N[(N[(N[(0.03432539682539683 * N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision] + N[(0.08611111111111111 / l), $MachinePrecision]), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m + N[(0.16666666666666666 / l), $MachinePrecision]), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + N[Power[l, -1.0], $MachinePrecision]), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 3.8e+147], N[(2.0 / N[(t$95$1 * N[(N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] * l), $MachinePrecision] * 2.0), $MachinePrecision] / N[(N[(N[(N[(0.5 - N[(0.5 * N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right)\\
\mathbf{if}\;k\_m \leq 4.4 \cdot 10^{-130}:\\
\;\;\;\;\frac{2}{\frac{\left(k\_m \cdot 2\right) \cdot t}{\frac{\ell}{t}} \cdot \frac{k\_m}{\frac{\ell}{t}}}\\
\mathbf{elif}\;k\_m \leq 0.125:\\
\;\;\;\;\frac{2}{t\_1 \cdot \left(\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.03432539682539683, \frac{k\_m \cdot k\_m}{\ell}, \frac{0.08611111111111111}{\ell}\right) \cdot k\_m, k\_m, \frac{0.16666666666666666}{\ell}\right), k\_m \cdot k\_m, {\ell}^{-1}\right) \cdot \left(k\_m \cdot k\_m\right)\right) \cdot \frac{t}{\ell}\right)}\\
\mathbf{elif}\;k\_m \leq 3.8 \cdot 10^{+147}:\\
\;\;\;\;\frac{2}{t\_1 \cdot \frac{\left(\tan k\_m \cdot \sin k\_m\right) \cdot t}{\ell \cdot \ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\left(\cos k\_m \cdot \ell\right) \cdot \ell\right) \cdot 2}{\left(\left(\left(0.5 - 0.5 \cdot \cos \left(k\_m + k\_m\right)\right) \cdot t\right) \cdot k\_m\right) \cdot k\_m}\\
\end{array}
\end{array}
if k < 4.3999999999999997e-130Initial program 54.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6454.8
Applied rewrites54.8%
Applied rewrites72.6%
Applied rewrites75.4%
if 4.3999999999999997e-130 < k < 0.125Initial program 58.3%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites84.4%
Applied rewrites98.7%
Taylor expanded in k around 0
Applied rewrites98.7%
if 0.125 < k < 3.7999999999999997e147Initial program 62.9%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites83.7%
Applied rewrites88.8%
Applied rewrites83.9%
if 3.7999999999999997e147 < k Initial program 56.8%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites62.9%
Applied rewrites63.0%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites73.5%
Applied rewrites73.4%
Final simplification79.3%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 4.4e-130)
(/ 2.0 (* (/ (* (* k_m 2.0) t) (/ l t)) (/ k_m (/ l t))))
(if (<= k_m 9.6e+126)
(/
2.0
(*
(/ t l)
(* (* (tan k_m) (/ (sin k_m) l)) (fma (* t t) 2.0 (* k_m k_m)))))
(/
(* (* (* (cos k_m) l) l) 2.0)
(* (* (* (- 0.5 (* 0.5 (cos (+ k_m k_m)))) t) k_m) k_m)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 4.4e-130) {
tmp = 2.0 / ((((k_m * 2.0) * t) / (l / t)) * (k_m / (l / t)));
} else if (k_m <= 9.6e+126) {
tmp = 2.0 / ((t / l) * ((tan(k_m) * (sin(k_m) / l)) * fma((t * t), 2.0, (k_m * k_m))));
} else {
tmp = (((cos(k_m) * l) * l) * 2.0) / ((((0.5 - (0.5 * cos((k_m + k_m)))) * t) * k_m) * k_m);
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 4.4e-130) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(k_m * 2.0) * t) / Float64(l / t)) * Float64(k_m / Float64(l / t)))); elseif (k_m <= 9.6e+126) tmp = Float64(2.0 / Float64(Float64(t / l) * Float64(Float64(tan(k_m) * Float64(sin(k_m) / l)) * fma(Float64(t * t), 2.0, Float64(k_m * k_m))))); else tmp = Float64(Float64(Float64(Float64(cos(k_m) * l) * l) * 2.0) / Float64(Float64(Float64(Float64(0.5 - Float64(0.5 * cos(Float64(k_m + k_m)))) * t) * k_m) * k_m)); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 4.4e-130], N[(2.0 / N[(N[(N[(N[(k$95$m * 2.0), $MachinePrecision] * t), $MachinePrecision] / N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 9.6e+126], N[(2.0 / N[(N[(t / l), $MachinePrecision] * N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] * l), $MachinePrecision] * 2.0), $MachinePrecision] / N[(N[(N[(N[(0.5 - N[(0.5 * N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 4.4 \cdot 10^{-130}:\\
\;\;\;\;\frac{2}{\frac{\left(k\_m \cdot 2\right) \cdot t}{\frac{\ell}{t}} \cdot \frac{k\_m}{\frac{\ell}{t}}}\\
\mathbf{elif}\;k\_m \leq 9.6 \cdot 10^{+126}:\\
\;\;\;\;\frac{2}{\frac{t}{\ell} \cdot \left(\left(\tan k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\left(\cos k\_m \cdot \ell\right) \cdot \ell\right) \cdot 2}{\left(\left(\left(0.5 - 0.5 \cdot \cos \left(k\_m + k\_m\right)\right) \cdot t\right) \cdot k\_m\right) \cdot k\_m}\\
\end{array}
\end{array}
if k < 4.3999999999999997e-130Initial program 54.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6454.8
Applied rewrites54.8%
Applied rewrites72.6%
Applied rewrites75.4%
if 4.3999999999999997e-130 < k < 9.60000000000000048e126Initial program 60.4%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites82.5%
Applied rewrites92.5%
Applied rewrites95.3%
if 9.60000000000000048e126 < k Initial program 58.2%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites68.1%
Applied rewrites68.1%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites77.2%
Applied rewrites77.1%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 4.4e-130)
(/ 2.0 (* (/ (* (* k_m 2.0) t) (/ l t)) (/ k_m (/ l t))))
(if (<= k_m 15000000000.0)
(/
2.0
(*
(fma (* t t) 2.0 (* k_m k_m))
(*
(*
(fma
(fma
(*
(fma
0.03432539682539683
(/ (* k_m k_m) l)
(/ 0.08611111111111111 l))
k_m)
k_m
(/ 0.16666666666666666 l))
(* k_m k_m)
(pow l -1.0))
(* k_m k_m))
(/ t l))))
(/
(* (* (* (cos k_m) l) l) 2.0)
(* (* (* (- 0.5 (* 0.5 (cos (+ k_m k_m)))) t) k_m) k_m)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 4.4e-130) {
tmp = 2.0 / ((((k_m * 2.0) * t) / (l / t)) * (k_m / (l / t)));
} else if (k_m <= 15000000000.0) {
tmp = 2.0 / (fma((t * t), 2.0, (k_m * k_m)) * ((fma(fma((fma(0.03432539682539683, ((k_m * k_m) / l), (0.08611111111111111 / l)) * k_m), k_m, (0.16666666666666666 / l)), (k_m * k_m), pow(l, -1.0)) * (k_m * k_m)) * (t / l)));
} else {
tmp = (((cos(k_m) * l) * l) * 2.0) / ((((0.5 - (0.5 * cos((k_m + k_m)))) * t) * k_m) * k_m);
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 4.4e-130) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(k_m * 2.0) * t) / Float64(l / t)) * Float64(k_m / Float64(l / t)))); elseif (k_m <= 15000000000.0) tmp = Float64(2.0 / Float64(fma(Float64(t * t), 2.0, Float64(k_m * k_m)) * Float64(Float64(fma(fma(Float64(fma(0.03432539682539683, Float64(Float64(k_m * k_m) / l), Float64(0.08611111111111111 / l)) * k_m), k_m, Float64(0.16666666666666666 / l)), Float64(k_m * k_m), (l ^ -1.0)) * Float64(k_m * k_m)) * Float64(t / l)))); else tmp = Float64(Float64(Float64(Float64(cos(k_m) * l) * l) * 2.0) / Float64(Float64(Float64(Float64(0.5 - Float64(0.5 * cos(Float64(k_m + k_m)))) * t) * k_m) * k_m)); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 4.4e-130], N[(2.0 / N[(N[(N[(N[(k$95$m * 2.0), $MachinePrecision] * t), $MachinePrecision] / N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 15000000000.0], N[(2.0 / N[(N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(N[(N[(0.03432539682539683 * N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision] + N[(0.08611111111111111 / l), $MachinePrecision]), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m + N[(0.16666666666666666 / l), $MachinePrecision]), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision] + N[Power[l, -1.0], $MachinePrecision]), $MachinePrecision] * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[Cos[k$95$m], $MachinePrecision] * l), $MachinePrecision] * l), $MachinePrecision] * 2.0), $MachinePrecision] / N[(N[(N[(N[(0.5 - N[(0.5 * N[Cos[N[(k$95$m + k$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 4.4 \cdot 10^{-130}:\\
\;\;\;\;\frac{2}{\frac{\left(k\_m \cdot 2\right) \cdot t}{\frac{\ell}{t}} \cdot \frac{k\_m}{\frac{\ell}{t}}}\\
\mathbf{elif}\;k\_m \leq 15000000000:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right) \cdot \left(\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.03432539682539683, \frac{k\_m \cdot k\_m}{\ell}, \frac{0.08611111111111111}{\ell}\right) \cdot k\_m, k\_m, \frac{0.16666666666666666}{\ell}\right), k\_m \cdot k\_m, {\ell}^{-1}\right) \cdot \left(k\_m \cdot k\_m\right)\right) \cdot \frac{t}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\left(\cos k\_m \cdot \ell\right) \cdot \ell\right) \cdot 2}{\left(\left(\left(0.5 - 0.5 \cdot \cos \left(k\_m + k\_m\right)\right) \cdot t\right) \cdot k\_m\right) \cdot k\_m}\\
\end{array}
\end{array}
if k < 4.3999999999999997e-130Initial program 54.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6454.8
Applied rewrites54.8%
Applied rewrites72.6%
Applied rewrites75.4%
if 4.3999999999999997e-130 < k < 1.5e10Initial program 59.1%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites85.7%
Applied rewrites98.8%
Taylor expanded in k around 0
Applied rewrites96.1%
if 1.5e10 < k Initial program 59.7%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites72.7%
Applied rewrites75.4%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites75.4%
Applied rewrites75.1%
Final simplification78.1%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= t 2.2e-54)
(/
2.0
(*
(fma (* t t) 2.0 (* k_m k_m))
(*
(* (* (fma (/ (* k_m k_m) l) 0.16666666666666666 (pow l -1.0)) k_m) k_m)
(/ t l))))
(/ 2.0 (* (/ (* (* k_m 2.0) t) (/ l t)) (/ k_m (/ l t))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 2.2e-54) {
tmp = 2.0 / (fma((t * t), 2.0, (k_m * k_m)) * (((fma(((k_m * k_m) / l), 0.16666666666666666, pow(l, -1.0)) * k_m) * k_m) * (t / l)));
} else {
tmp = 2.0 / ((((k_m * 2.0) * t) / (l / t)) * (k_m / (l / t)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 2.2e-54) tmp = Float64(2.0 / Float64(fma(Float64(t * t), 2.0, Float64(k_m * k_m)) * Float64(Float64(Float64(fma(Float64(Float64(k_m * k_m) / l), 0.16666666666666666, (l ^ -1.0)) * k_m) * k_m) * Float64(t / l)))); else tmp = Float64(2.0 / Float64(Float64(Float64(Float64(k_m * 2.0) * t) / Float64(l / t)) * Float64(k_m / Float64(l / t)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 2.2e-54], N[(2.0 / N[(N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision] * 0.16666666666666666 + N[Power[l, -1.0], $MachinePrecision]), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[(k$95$m * 2.0), $MachinePrecision] * t), $MachinePrecision] / N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 2.2 \cdot 10^{-54}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right) \cdot \left(\left(\left(\mathsf{fma}\left(\frac{k\_m \cdot k\_m}{\ell}, 0.16666666666666666, {\ell}^{-1}\right) \cdot k\_m\right) \cdot k\_m\right) \cdot \frac{t}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\left(k\_m \cdot 2\right) \cdot t}{\frac{\ell}{t}} \cdot \frac{k\_m}{\frac{\ell}{t}}}\\
\end{array}
\end{array}
if t < 2.2e-54Initial program 50.7%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites74.6%
Applied rewrites81.6%
Taylor expanded in k around 0
Applied rewrites71.4%
if 2.2e-54 < t Initial program 73.4%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6462.6
Applied rewrites62.6%
Applied rewrites81.5%
Applied rewrites85.1%
Final simplification74.9%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= t 5.8e-63) (/ 2.0 (* (fma (* t t) 2.0 (* k_m k_m)) (/ (/ (* (* k_m t) k_m) l) l))) (/ 2.0 (* (/ (* (* k_m 2.0) t) (/ l t)) (/ k_m (/ l t))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 5.8e-63) {
tmp = 2.0 / (fma((t * t), 2.0, (k_m * k_m)) * ((((k_m * t) * k_m) / l) / l));
} else {
tmp = 2.0 / ((((k_m * 2.0) * t) / (l / t)) * (k_m / (l / t)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 5.8e-63) tmp = Float64(2.0 / Float64(fma(Float64(t * t), 2.0, Float64(k_m * k_m)) * Float64(Float64(Float64(Float64(k_m * t) * k_m) / l) / l))); else tmp = Float64(2.0 / Float64(Float64(Float64(Float64(k_m * 2.0) * t) / Float64(l / t)) * Float64(k_m / Float64(l / t)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 5.8e-63], N[(2.0 / N[(N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(k$95$m * t), $MachinePrecision] * k$95$m), $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[(k$95$m * 2.0), $MachinePrecision] * t), $MachinePrecision] / N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 5.8 \cdot 10^{-63}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right) \cdot \frac{\frac{\left(k\_m \cdot t\right) \cdot k\_m}{\ell}}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\left(k\_m \cdot 2\right) \cdot t}{\frac{\ell}{t}} \cdot \frac{k\_m}{\frac{\ell}{t}}}\\
\end{array}
\end{array}
if t < 5.7999999999999995e-63Initial program 50.7%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites74.6%
Applied rewrites81.6%
Applied rewrites82.7%
Taylor expanded in k around 0
Applied rewrites70.1%
if 5.7999999999999995e-63 < t Initial program 73.4%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6462.6
Applied rewrites62.6%
Applied rewrites81.5%
Applied rewrites85.1%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= t 5.8e-63) (/ 2.0 (* (fma (* t t) 2.0 (* k_m k_m)) (/ (/ (* (* k_m t) k_m) l) l))) (/ 2.0 (* (* (* 2.0 (/ k_m (/ l t))) (* k_m t)) (/ t l)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 5.8e-63) {
tmp = 2.0 / (fma((t * t), 2.0, (k_m * k_m)) * ((((k_m * t) * k_m) / l) / l));
} else {
tmp = 2.0 / (((2.0 * (k_m / (l / t))) * (k_m * t)) * (t / l));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 5.8e-63) tmp = Float64(2.0 / Float64(fma(Float64(t * t), 2.0, Float64(k_m * k_m)) * Float64(Float64(Float64(Float64(k_m * t) * k_m) / l) / l))); else tmp = Float64(2.0 / Float64(Float64(Float64(2.0 * Float64(k_m / Float64(l / t))) * Float64(k_m * t)) * Float64(t / l))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 5.8e-63], N[(2.0 / N[(N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(k$95$m * t), $MachinePrecision] * k$95$m), $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(2.0 * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 5.8 \cdot 10^{-63}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right) \cdot \frac{\frac{\left(k\_m \cdot t\right) \cdot k\_m}{\ell}}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(2 \cdot \frac{k\_m}{\frac{\ell}{t}}\right) \cdot \left(k\_m \cdot t\right)\right) \cdot \frac{t}{\ell}}\\
\end{array}
\end{array}
if t < 5.7999999999999995e-63Initial program 50.7%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites74.6%
Applied rewrites81.6%
Applied rewrites82.7%
Taylor expanded in k around 0
Applied rewrites70.1%
if 5.7999999999999995e-63 < t Initial program 73.4%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6462.6
Applied rewrites62.6%
Applied rewrites61.4%
Applied rewrites82.5%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= t 4.2e-7) (/ 2.0 (* (fma (* t t) 2.0 (* k_m k_m)) (* t (/ (/ (* k_m k_m) l) l)))) (/ 2.0 (* (* (* 2.0 (/ k_m (/ l t))) (* k_m t)) (/ t l)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 4.2e-7) {
tmp = 2.0 / (fma((t * t), 2.0, (k_m * k_m)) * (t * (((k_m * k_m) / l) / l)));
} else {
tmp = 2.0 / (((2.0 * (k_m / (l / t))) * (k_m * t)) * (t / l));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 4.2e-7) tmp = Float64(2.0 / Float64(fma(Float64(t * t), 2.0, Float64(k_m * k_m)) * Float64(t * Float64(Float64(Float64(k_m * k_m) / l) / l)))); else tmp = Float64(2.0 / Float64(Float64(Float64(2.0 * Float64(k_m / Float64(l / t))) * Float64(k_m * t)) * Float64(t / l))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 4.2e-7], N[(2.0 / N[(N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(t * N[(N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(2.0 * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 4.2 \cdot 10^{-7}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right) \cdot \left(t \cdot \frac{\frac{k\_m \cdot k\_m}{\ell}}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(2 \cdot \frac{k\_m}{\frac{\ell}{t}}\right) \cdot \left(k\_m \cdot t\right)\right) \cdot \frac{t}{\ell}}\\
\end{array}
\end{array}
if t < 4.2e-7Initial program 52.1%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites75.0%
Applied rewrites82.1%
Taylor expanded in k around 0
Applied rewrites67.2%
if 4.2e-7 < t Initial program 70.8%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6463.4
Applied rewrites63.4%
Applied rewrites62.1%
Applied rewrites85.3%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= t 4.2e-7) (/ 2.0 (* (/ (/ (* k_m k_m) l) l) (* t (fma (* t t) 2.0 (* k_m k_m))))) (/ 2.0 (* (* (* 2.0 (/ k_m (/ l t))) (* k_m t)) (/ t l)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 4.2e-7) {
tmp = 2.0 / ((((k_m * k_m) / l) / l) * (t * fma((t * t), 2.0, (k_m * k_m))));
} else {
tmp = 2.0 / (((2.0 * (k_m / (l / t))) * (k_m * t)) * (t / l));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 4.2e-7) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(k_m * k_m) / l) / l) * Float64(t * fma(Float64(t * t), 2.0, Float64(k_m * k_m))))); else tmp = Float64(2.0 / Float64(Float64(Float64(2.0 * Float64(k_m / Float64(l / t))) * Float64(k_m * t)) * Float64(t / l))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 4.2e-7], N[(2.0 / N[(N[(N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision] * N[(t * N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(2.0 * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(k$95$m * t), $MachinePrecision]), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 4.2 \cdot 10^{-7}:\\
\;\;\;\;\frac{2}{\frac{\frac{k\_m \cdot k\_m}{\ell}}{\ell} \cdot \left(t \cdot \mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(2 \cdot \frac{k\_m}{\frac{\ell}{t}}\right) \cdot \left(k\_m \cdot t\right)\right) \cdot \frac{t}{\ell}}\\
\end{array}
\end{array}
if t < 4.2e-7Initial program 52.1%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites75.0%
Taylor expanded in k around 0
Applied rewrites66.3%
if 4.2e-7 < t Initial program 70.8%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6463.4
Applied rewrites63.4%
Applied rewrites62.1%
Applied rewrites85.3%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= t 8.2e-61) (/ 2.0 (* (/ (/ (* k_m k_m) l) l) (* t (fma (* t t) 2.0 (* k_m k_m))))) (/ 2.0 (* (* k_m 2.0) (* (* (/ t l) t) (/ k_m (/ l t)))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 8.2e-61) {
tmp = 2.0 / ((((k_m * k_m) / l) / l) * (t * fma((t * t), 2.0, (k_m * k_m))));
} else {
tmp = 2.0 / ((k_m * 2.0) * (((t / l) * t) * (k_m / (l / t))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 8.2e-61) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(k_m * k_m) / l) / l) * Float64(t * fma(Float64(t * t), 2.0, Float64(k_m * k_m))))); else tmp = Float64(2.0 / Float64(Float64(k_m * 2.0) * Float64(Float64(Float64(t / l) * t) * Float64(k_m / Float64(l / t))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 8.2e-61], N[(2.0 / N[(N[(N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision] * N[(t * N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k$95$m * 2.0), $MachinePrecision] * N[(N[(N[(t / l), $MachinePrecision] * t), $MachinePrecision] * N[(k$95$m / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 8.2 \cdot 10^{-61}:\\
\;\;\;\;\frac{2}{\frac{\frac{k\_m \cdot k\_m}{\ell}}{\ell} \cdot \left(t \cdot \mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(k\_m \cdot 2\right) \cdot \left(\left(\frac{t}{\ell} \cdot t\right) \cdot \frac{k\_m}{\frac{\ell}{t}}\right)}\\
\end{array}
\end{array}
if t < 8.19999999999999998e-61Initial program 50.7%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites74.6%
Taylor expanded in k around 0
Applied rewrites66.7%
if 8.19999999999999998e-61 < t Initial program 73.4%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6462.6
Applied rewrites62.6%
Applied rewrites61.4%
Applied rewrites79.5%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= t 2.25e-66) (/ 2.0 (* (/ (/ (* k_m k_m) l) l) (* t (fma (* t t) 2.0 (* k_m k_m))))) (/ 2.0 (* (* (* (/ (* k_m 2.0) l) (* t t)) (/ k_m l)) t))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 2.25e-66) {
tmp = 2.0 / ((((k_m * k_m) / l) / l) * (t * fma((t * t), 2.0, (k_m * k_m))));
} else {
tmp = 2.0 / (((((k_m * 2.0) / l) * (t * t)) * (k_m / l)) * t);
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 2.25e-66) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(k_m * k_m) / l) / l) * Float64(t * fma(Float64(t * t), 2.0, Float64(k_m * k_m))))); else tmp = Float64(2.0 / Float64(Float64(Float64(Float64(Float64(k_m * 2.0) / l) * Float64(t * t)) * Float64(k_m / l)) * t)); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 2.25e-66], N[(2.0 / N[(N[(N[(N[(k$95$m * k$95$m), $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision] * N[(t * N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[(N[(k$95$m * 2.0), $MachinePrecision] / l), $MachinePrecision] * N[(t * t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 2.25 \cdot 10^{-66}:\\
\;\;\;\;\frac{2}{\frac{\frac{k\_m \cdot k\_m}{\ell}}{\ell} \cdot \left(t \cdot \mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{k\_m \cdot 2}{\ell} \cdot \left(t \cdot t\right)\right) \cdot \frac{k\_m}{\ell}\right) \cdot t}\\
\end{array}
\end{array}
if t < 2.2499999999999999e-66Initial program 50.7%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites74.6%
Taylor expanded in k around 0
Applied rewrites66.7%
if 2.2499999999999999e-66 < t Initial program 73.4%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6462.6
Applied rewrites62.6%
Applied rewrites61.4%
Applied rewrites80.7%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 4.8e-130) (/ 2.0 (* (* (* (/ (* k_m 2.0) l) (* t t)) (/ k_m l)) t)) (/ 2.0 (* (fma (* t t) 2.0 (* k_m k_m)) (* (* (/ t (* l l)) k_m) k_m)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 4.8e-130) {
tmp = 2.0 / (((((k_m * 2.0) / l) * (t * t)) * (k_m / l)) * t);
} else {
tmp = 2.0 / (fma((t * t), 2.0, (k_m * k_m)) * (((t / (l * l)) * k_m) * k_m));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 4.8e-130) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(Float64(k_m * 2.0) / l) * Float64(t * t)) * Float64(k_m / l)) * t)); else tmp = Float64(2.0 / Float64(fma(Float64(t * t), 2.0, Float64(k_m * k_m)) * Float64(Float64(Float64(t / Float64(l * l)) * k_m) * k_m))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 4.8e-130], N[(2.0 / N[(N[(N[(N[(N[(k$95$m * 2.0), $MachinePrecision] / l), $MachinePrecision] * N[(t * t), $MachinePrecision]), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(t * t), $MachinePrecision] * 2.0 + N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(t / N[(l * l), $MachinePrecision]), $MachinePrecision] * k$95$m), $MachinePrecision] * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 4.8 \cdot 10^{-130}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{k\_m \cdot 2}{\ell} \cdot \left(t \cdot t\right)\right) \cdot \frac{k\_m}{\ell}\right) \cdot t}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(t \cdot t, 2, k\_m \cdot k\_m\right) \cdot \left(\left(\frac{t}{\ell \cdot \ell} \cdot k\_m\right) \cdot k\_m\right)}\\
\end{array}
\end{array}
if k < 4.79999999999999993e-130Initial program 54.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6454.8
Applied rewrites54.8%
Applied rewrites51.3%
Applied rewrites69.4%
if 4.79999999999999993e-130 < k Initial program 59.5%
Taylor expanded in t around 0
distribute-rgt-inN/A
associate-*l*N/A
associate-*l/N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
Applied rewrites76.8%
Applied rewrites82.8%
Applied rewrites82.8%
Taylor expanded in k around 0
Applied rewrites67.9%
Final simplification68.8%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (/ 2.0 (* k_m (* (/ t l) (* (/ (* k_m 2.0) l) (* t t))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / (k_m * ((t / l) * (((k_m * 2.0) / l) * (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 = 2.0d0 / (k_m * ((t / l) * (((k_m * 2.0d0) / l) * (t * t))))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return 2.0 / (k_m * ((t / l) * (((k_m * 2.0) / l) * (t * t))));
}
k_m = math.fabs(k) def code(t, l, k_m): return 2.0 / (k_m * ((t / l) * (((k_m * 2.0) / l) * (t * t))))
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64(k_m * Float64(Float64(t / l) * Float64(Float64(Float64(k_m * 2.0) / l) * Float64(t * t))))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = 2.0 / (k_m * ((t / l) * (((k_m * 2.0) / l) * (t * t)))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(k$95$m * N[(N[(t / l), $MachinePrecision] * N[(N[(N[(k$95$m * 2.0), $MachinePrecision] / l), $MachinePrecision] * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{k\_m \cdot \left(\frac{t}{\ell} \cdot \left(\frac{k\_m \cdot 2}{\ell} \cdot \left(t \cdot t\right)\right)\right)}
\end{array}
Initial program 56.5%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6457.8
Applied rewrites57.8%
Applied rewrites54.6%
Applied rewrites65.3%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (/ 2.0 (* (* (* k_m k_m) 2.0) (* t (/ (* t t) (* l l))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / (((k_m * k_m) * 2.0) * (t * ((t * t) / (l * l))));
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = 2.0d0 / (((k_m * k_m) * 2.0d0) * (t * ((t * t) / (l * l))))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return 2.0 / (((k_m * k_m) * 2.0) * (t * ((t * t) / (l * l))));
}
k_m = math.fabs(k) def code(t, l, k_m): return 2.0 / (((k_m * k_m) * 2.0) * (t * ((t * t) / (l * l))))
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64(Float64(Float64(k_m * k_m) * 2.0) * Float64(t * Float64(Float64(t * t) / Float64(l * l))))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = 2.0 / (((k_m * k_m) * 2.0) * (t * ((t * t) / (l * l)))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(N[(N[(k$95$m * k$95$m), $MachinePrecision] * 2.0), $MachinePrecision] * N[(t * N[(N[(t * t), $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{\left(\left(k\_m \cdot k\_m\right) \cdot 2\right) \cdot \left(t \cdot \frac{t \cdot t}{\ell \cdot \ell}\right)}
\end{array}
Initial program 56.5%
Taylor expanded in k around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6457.8
Applied rewrites57.8%
Applied rewrites54.6%
Applied rewrites57.8%
herbie shell --seed 2024302
(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))))