
(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}
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(let* ((t_2 (pow (sin k) 2.0)) (t_3 (/ (sqrt 2.0) k)))
(*
t_s
(if (<= t_m 1.05e-107)
(/
2.0
(*
(/ (* t_m (pow k 2.0)) (cos k))
(exp (+ (log t_2) (* -2.0 (log l_m))))))
(*
(/
(* t_m t_3)
(pow (* t_m (* (pow (cbrt l_m) -2.0) (cbrt (/ t_2 (cos k))))) 2.0))
(/ (* t_3 (pow (cbrt l_m) 2.0)) (cbrt (* (sin k) (tan k)))))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double t_2 = pow(sin(k), 2.0);
double t_3 = sqrt(2.0) / k;
double tmp;
if (t_m <= 1.05e-107) {
tmp = 2.0 / (((t_m * pow(k, 2.0)) / cos(k)) * exp((log(t_2) + (-2.0 * log(l_m)))));
} else {
tmp = ((t_m * t_3) / pow((t_m * (pow(cbrt(l_m), -2.0) * cbrt((t_2 / cos(k))))), 2.0)) * ((t_3 * pow(cbrt(l_m), 2.0)) / cbrt((sin(k) * tan(k))));
}
return t_s * tmp;
}
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double t_2 = Math.pow(Math.sin(k), 2.0);
double t_3 = Math.sqrt(2.0) / k;
double tmp;
if (t_m <= 1.05e-107) {
tmp = 2.0 / (((t_m * Math.pow(k, 2.0)) / Math.cos(k)) * Math.exp((Math.log(t_2) + (-2.0 * Math.log(l_m)))));
} else {
tmp = ((t_m * t_3) / Math.pow((t_m * (Math.pow(Math.cbrt(l_m), -2.0) * Math.cbrt((t_2 / Math.cos(k))))), 2.0)) * ((t_3 * Math.pow(Math.cbrt(l_m), 2.0)) / Math.cbrt((Math.sin(k) * Math.tan(k))));
}
return t_s * tmp;
}
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) t_2 = sin(k) ^ 2.0 t_3 = Float64(sqrt(2.0) / k) tmp = 0.0 if (t_m <= 1.05e-107) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) / cos(k)) * exp(Float64(log(t_2) + Float64(-2.0 * log(l_m)))))); else tmp = Float64(Float64(Float64(t_m * t_3) / (Float64(t_m * Float64((cbrt(l_m) ^ -2.0) * cbrt(Float64(t_2 / cos(k))))) ^ 2.0)) * Float64(Float64(t_3 * (cbrt(l_m) ^ 2.0)) / cbrt(Float64(sin(k) * tan(k))))); end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := Block[{t$95$2 = N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$3 = N[(N[Sqrt[2.0], $MachinePrecision] / k), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 1.05e-107], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * N[Exp[N[(N[Log[t$95$2], $MachinePrecision] + N[(-2.0 * N[Log[l$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t$95$m * t$95$3), $MachinePrecision] / N[Power[N[(t$95$m * N[(N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], -2.0], $MachinePrecision] * N[Power[N[(t$95$2 / N[Cos[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(t$95$3 * N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[Power[N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := {\sin k}^{2}\\
t_3 := \frac{\sqrt{2}}{k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.05 \cdot 10^{-107}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{2}}{\cos k} \cdot e^{\log t\_2 + -2 \cdot \log l\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_m \cdot t\_3}{{\left(t\_m \cdot \left({\left(\sqrt[3]{l\_m}\right)}^{-2} \cdot \sqrt[3]{\frac{t\_2}{\cos k}}\right)\right)}^{2}} \cdot \frac{t\_3 \cdot {\left(\sqrt[3]{l\_m}\right)}^{2}}{\sqrt[3]{\sin k \cdot \tan k}}\\
\end{array}
\end{array}
\end{array}
if t < 1.05e-107Initial program 30.2%
Simplified30.2%
Taylor expanded in t around 0 75.2%
associate-*r*75.2%
*-commutative75.2%
times-frac78.0%
Simplified78.0%
pow278.0%
add-log-exp59.9%
div-inv59.8%
exp-prod64.1%
pow264.1%
pow-flip64.1%
metadata-eval64.1%
Applied egg-rr64.1%
pow-exp59.8%
rem-log-exp78.0%
pow-to-exp35.6%
pow-to-exp17.3%
prod-exp18.9%
rem-log-exp18.9%
pow-to-exp44.3%
rem-log-exp40.9%
pow-to-exp76.6%
log-pow44.3%
Applied egg-rr44.3%
if 1.05e-107 < t Initial program 26.2%
Simplified26.2%
Applied egg-rr88.4%
associate-/r/88.4%
associate-/r*88.4%
associate-/r/88.4%
Simplified88.4%
associate-*r/87.5%
Applied egg-rr87.5%
associate-/l*88.4%
associate-*r/88.4%
associate-*l*88.5%
associate-/l/85.2%
*-commutative85.2%
times-frac88.5%
*-inverses88.5%
*-lft-identity88.5%
Simplified88.5%
Taylor expanded in k around inf 88.6%
Final simplification59.0%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(let* ((t_2 (cbrt (* (sin k) (tan k)))) (t_3 (/ (sqrt 2.0) k)))
(*
t_s
(if (<= t_m 8e-108)
(/
2.0
(*
(/ (* t_m (pow k 2.0)) (cos k))
(exp (+ (log (pow (sin k) 2.0)) (* -2.0 (log l_m))))))
(*
(/ (* t_3 (pow (cbrt l_m) 2.0)) t_2)
(/ (* t_m t_3) (pow (* t_m (* (pow (cbrt l_m) -2.0) t_2)) 2.0)))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double t_2 = cbrt((sin(k) * tan(k)));
double t_3 = sqrt(2.0) / k;
double tmp;
if (t_m <= 8e-108) {
tmp = 2.0 / (((t_m * pow(k, 2.0)) / cos(k)) * exp((log(pow(sin(k), 2.0)) + (-2.0 * log(l_m)))));
} else {
tmp = ((t_3 * pow(cbrt(l_m), 2.0)) / t_2) * ((t_m * t_3) / pow((t_m * (pow(cbrt(l_m), -2.0) * t_2)), 2.0));
}
return t_s * tmp;
}
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double t_2 = Math.cbrt((Math.sin(k) * Math.tan(k)));
double t_3 = Math.sqrt(2.0) / k;
double tmp;
if (t_m <= 8e-108) {
tmp = 2.0 / (((t_m * Math.pow(k, 2.0)) / Math.cos(k)) * Math.exp((Math.log(Math.pow(Math.sin(k), 2.0)) + (-2.0 * Math.log(l_m)))));
} else {
tmp = ((t_3 * Math.pow(Math.cbrt(l_m), 2.0)) / t_2) * ((t_m * t_3) / Math.pow((t_m * (Math.pow(Math.cbrt(l_m), -2.0) * t_2)), 2.0));
}
return t_s * tmp;
}
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) t_2 = cbrt(Float64(sin(k) * tan(k))) t_3 = Float64(sqrt(2.0) / k) tmp = 0.0 if (t_m <= 8e-108) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) / cos(k)) * exp(Float64(log((sin(k) ^ 2.0)) + Float64(-2.0 * log(l_m)))))); else tmp = Float64(Float64(Float64(t_3 * (cbrt(l_m) ^ 2.0)) / t_2) * Float64(Float64(t_m * t_3) / (Float64(t_m * Float64((cbrt(l_m) ^ -2.0) * t_2)) ^ 2.0))); end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := Block[{t$95$2 = N[Power[N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]}, Block[{t$95$3 = N[(N[Sqrt[2.0], $MachinePrecision] / k), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 8e-108], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * N[Exp[N[(N[Log[N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision] + N[(-2.0 * N[Log[l$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t$95$3 * N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision] * N[(N[(t$95$m * t$95$3), $MachinePrecision] / N[Power[N[(t$95$m * N[(N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], -2.0], $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sqrt[3]{\sin k \cdot \tan k}\\
t_3 := \frac{\sqrt{2}}{k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 8 \cdot 10^{-108}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{2}}{\cos k} \cdot e^{\log \left({\sin k}^{2}\right) + -2 \cdot \log l\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_3 \cdot {\left(\sqrt[3]{l\_m}\right)}^{2}}{t\_2} \cdot \frac{t\_m \cdot t\_3}{{\left(t\_m \cdot \left({\left(\sqrt[3]{l\_m}\right)}^{-2} \cdot t\_2\right)\right)}^{2}}\\
\end{array}
\end{array}
\end{array}
if t < 8.00000000000000032e-108Initial program 30.2%
Simplified30.2%
Taylor expanded in t around 0 75.2%
associate-*r*75.2%
*-commutative75.2%
times-frac78.0%
Simplified78.0%
pow278.0%
add-log-exp59.9%
div-inv59.8%
exp-prod64.1%
pow264.1%
pow-flip64.1%
metadata-eval64.1%
Applied egg-rr64.1%
pow-exp59.8%
rem-log-exp78.0%
pow-to-exp35.6%
pow-to-exp17.3%
prod-exp18.9%
rem-log-exp18.9%
pow-to-exp44.3%
rem-log-exp40.9%
pow-to-exp76.6%
log-pow44.3%
Applied egg-rr44.3%
if 8.00000000000000032e-108 < t Initial program 26.2%
Simplified26.2%
Applied egg-rr88.4%
associate-/r/88.4%
associate-/r*88.4%
associate-/r/88.4%
Simplified88.4%
associate-*r/87.5%
Applied egg-rr87.5%
associate-/l*88.4%
associate-*r/88.4%
associate-*l*88.5%
associate-/l/85.2%
*-commutative85.2%
times-frac88.5%
*-inverses88.5%
*-lft-identity88.5%
Simplified88.5%
Final simplification59.0%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(let* ((t_2 (cbrt (* (sin k) (tan k)))) (t_3 (/ (sqrt 2.0) k)))
(*
t_s
(if (<= t_m 7.5e-108)
(/
2.0
(*
(/ (* t_m (pow k 2.0)) (cos k))
(exp (+ (log (pow (sin k) 2.0)) (* -2.0 (log l_m))))))
(*
(/ (* t_m t_3) (pow (* t_m (* (pow (cbrt l_m) -2.0) t_2)) 2.0))
(* t_3 (/ (pow (cbrt l_m) 2.0) t_2)))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double t_2 = cbrt((sin(k) * tan(k)));
double t_3 = sqrt(2.0) / k;
double tmp;
if (t_m <= 7.5e-108) {
tmp = 2.0 / (((t_m * pow(k, 2.0)) / cos(k)) * exp((log(pow(sin(k), 2.0)) + (-2.0 * log(l_m)))));
} else {
tmp = ((t_m * t_3) / pow((t_m * (pow(cbrt(l_m), -2.0) * t_2)), 2.0)) * (t_3 * (pow(cbrt(l_m), 2.0) / t_2));
}
return t_s * tmp;
}
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double t_2 = Math.cbrt((Math.sin(k) * Math.tan(k)));
double t_3 = Math.sqrt(2.0) / k;
double tmp;
if (t_m <= 7.5e-108) {
tmp = 2.0 / (((t_m * Math.pow(k, 2.0)) / Math.cos(k)) * Math.exp((Math.log(Math.pow(Math.sin(k), 2.0)) + (-2.0 * Math.log(l_m)))));
} else {
tmp = ((t_m * t_3) / Math.pow((t_m * (Math.pow(Math.cbrt(l_m), -2.0) * t_2)), 2.0)) * (t_3 * (Math.pow(Math.cbrt(l_m), 2.0) / t_2));
}
return t_s * tmp;
}
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) t_2 = cbrt(Float64(sin(k) * tan(k))) t_3 = Float64(sqrt(2.0) / k) tmp = 0.0 if (t_m <= 7.5e-108) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) / cos(k)) * exp(Float64(log((sin(k) ^ 2.0)) + Float64(-2.0 * log(l_m)))))); else tmp = Float64(Float64(Float64(t_m * t_3) / (Float64(t_m * Float64((cbrt(l_m) ^ -2.0) * t_2)) ^ 2.0)) * Float64(t_3 * Float64((cbrt(l_m) ^ 2.0) / t_2))); end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := Block[{t$95$2 = N[Power[N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]}, Block[{t$95$3 = N[(N[Sqrt[2.0], $MachinePrecision] / k), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 7.5e-108], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * N[Exp[N[(N[Log[N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision] + N[(-2.0 * N[Log[l$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t$95$m * t$95$3), $MachinePrecision] / N[Power[N[(t$95$m * N[(N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], -2.0], $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$3 * N[(N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], 2.0], $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sqrt[3]{\sin k \cdot \tan k}\\
t_3 := \frac{\sqrt{2}}{k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 7.5 \cdot 10^{-108}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{2}}{\cos k} \cdot e^{\log \left({\sin k}^{2}\right) + -2 \cdot \log l\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_m \cdot t\_3}{{\left(t\_m \cdot \left({\left(\sqrt[3]{l\_m}\right)}^{-2} \cdot t\_2\right)\right)}^{2}} \cdot \left(t\_3 \cdot \frac{{\left(\sqrt[3]{l\_m}\right)}^{2}}{t\_2}\right)\\
\end{array}
\end{array}
\end{array}
if t < 7.4999999999999993e-108Initial program 30.2%
Simplified30.2%
Taylor expanded in t around 0 75.2%
associate-*r*75.2%
*-commutative75.2%
times-frac78.0%
Simplified78.0%
pow278.0%
add-log-exp59.9%
div-inv59.8%
exp-prod64.1%
pow264.1%
pow-flip64.1%
metadata-eval64.1%
Applied egg-rr64.1%
pow-exp59.8%
rem-log-exp78.0%
pow-to-exp35.6%
pow-to-exp17.3%
prod-exp18.9%
rem-log-exp18.9%
pow-to-exp44.3%
rem-log-exp40.9%
pow-to-exp76.6%
log-pow44.3%
Applied egg-rr44.3%
if 7.4999999999999993e-108 < t Initial program 26.2%
Simplified26.2%
Applied egg-rr88.4%
associate-/r/88.4%
associate-/r*88.4%
associate-/r/88.4%
Simplified88.4%
associate-*r/87.5%
Applied egg-rr87.5%
associate-/l*88.4%
associate-*r/88.4%
associate-*l*88.5%
associate-/l/85.2%
*-commutative85.2%
times-frac88.5%
*-inverses88.5%
*-lft-identity88.5%
Simplified88.5%
associate-*r/87.6%
associate-/l*89.5%
associate-*r*89.4%
Applied egg-rr89.4%
associate-/l*90.3%
associate-*r/88.4%
associate-*l*88.5%
associate-*r/88.4%
Simplified88.4%
Final simplification59.0%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(let* ((t_2 (/ (sqrt 2.0) k)))
(*
t_s
(if (<= l_m 2.3e-158)
(*
(/
(* t_m t_2)
(pow (* t_m (* (pow (cbrt l_m) -2.0) (cbrt (* (sin k) (tan k))))) 2.0))
(/ (* t_2 (pow (cbrt l_m) 2.0)) (cbrt (pow k 2.0))))
(/
2.0
(*
(/ (exp (+ (log t_m) (* 2.0 (log k)))) (cos k))
(/ (pow (sin k) 2.0) (pow l_m 2.0))))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double t_2 = sqrt(2.0) / k;
double tmp;
if (l_m <= 2.3e-158) {
tmp = ((t_m * t_2) / pow((t_m * (pow(cbrt(l_m), -2.0) * cbrt((sin(k) * tan(k))))), 2.0)) * ((t_2 * pow(cbrt(l_m), 2.0)) / cbrt(pow(k, 2.0)));
} else {
tmp = 2.0 / ((exp((log(t_m) + (2.0 * log(k)))) / cos(k)) * (pow(sin(k), 2.0) / pow(l_m, 2.0)));
}
return t_s * tmp;
}
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double t_2 = Math.sqrt(2.0) / k;
double tmp;
if (l_m <= 2.3e-158) {
tmp = ((t_m * t_2) / Math.pow((t_m * (Math.pow(Math.cbrt(l_m), -2.0) * Math.cbrt((Math.sin(k) * Math.tan(k))))), 2.0)) * ((t_2 * Math.pow(Math.cbrt(l_m), 2.0)) / Math.cbrt(Math.pow(k, 2.0)));
} else {
tmp = 2.0 / ((Math.exp((Math.log(t_m) + (2.0 * Math.log(k)))) / Math.cos(k)) * (Math.pow(Math.sin(k), 2.0) / Math.pow(l_m, 2.0)));
}
return t_s * tmp;
}
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) t_2 = Float64(sqrt(2.0) / k) tmp = 0.0 if (l_m <= 2.3e-158) tmp = Float64(Float64(Float64(t_m * t_2) / (Float64(t_m * Float64((cbrt(l_m) ^ -2.0) * cbrt(Float64(sin(k) * tan(k))))) ^ 2.0)) * Float64(Float64(t_2 * (cbrt(l_m) ^ 2.0)) / cbrt((k ^ 2.0)))); else tmp = Float64(2.0 / Float64(Float64(exp(Float64(log(t_m) + Float64(2.0 * log(k)))) / cos(k)) * Float64((sin(k) ^ 2.0) / (l_m ^ 2.0)))); end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := Block[{t$95$2 = N[(N[Sqrt[2.0], $MachinePrecision] / k), $MachinePrecision]}, N[(t$95$s * If[LessEqual[l$95$m, 2.3e-158], N[(N[(N[(t$95$m * t$95$2), $MachinePrecision] / N[Power[N[(t$95$m * N[(N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], -2.0], $MachinePrecision] * N[Power[N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(t$95$2 * N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[Power[N[Power[k, 2.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Exp[N[(N[Log[t$95$m], $MachinePrecision] + N[(2.0 * N[Log[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] / N[Power[l$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{\sqrt{2}}{k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;l\_m \leq 2.3 \cdot 10^{-158}:\\
\;\;\;\;\frac{t\_m \cdot t\_2}{{\left(t\_m \cdot \left({\left(\sqrt[3]{l\_m}\right)}^{-2} \cdot \sqrt[3]{\sin k \cdot \tan k}\right)\right)}^{2}} \cdot \frac{t\_2 \cdot {\left(\sqrt[3]{l\_m}\right)}^{2}}{\sqrt[3]{{k}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{e^{\log t\_m + 2 \cdot \log k}}{\cos k} \cdot \frac{{\sin k}^{2}}{{l\_m}^{2}}}\\
\end{array}
\end{array}
\end{array}
if l < 2.2999999999999999e-158Initial program 29.8%
Simplified29.8%
Applied egg-rr86.3%
associate-/r/86.3%
associate-/r*86.4%
associate-/r/86.4%
Simplified86.4%
associate-*r/85.9%
Applied egg-rr85.9%
associate-/l*86.4%
associate-*r/86.4%
associate-*l*86.5%
associate-/l/84.7%
*-commutative84.7%
times-frac86.5%
*-inverses86.5%
*-lft-identity86.5%
Simplified86.5%
Taylor expanded in k around 0 78.0%
if 2.2999999999999999e-158 < l Initial program 27.4%
Simplified27.4%
Taylor expanded in t around 0 79.8%
associate-*r*79.9%
*-commutative79.9%
times-frac83.6%
Simplified83.6%
add-log-exp37.8%
*-commutative37.8%
exp-prod28.1%
Applied egg-rr28.1%
pow-exp37.8%
rem-log-exp83.6%
add-exp-log34.0%
pow-to-exp15.6%
prod-exp16.4%
rem-log-exp15.6%
pow-to-exp33.9%
log-pow16.4%
Applied egg-rr16.4%
Final simplification54.0%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(let* ((t_2 (/ (pow (sin k) 2.0) (pow l_m 2.0))))
(*
t_s
(if (<= k 6e-161)
(/ 2.0 (* (/ (* t_m (pow k 2.0)) (cos k)) 0.0))
(if (<= k 2e+147)
(/ 2.0 (* t_2 (* (pow k 2.0) (/ t_m (cos k)))))
(/ 2.0 (* (/ (exp (+ (log t_m) (* 2.0 (log k)))) (cos k)) t_2)))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double t_2 = pow(sin(k), 2.0) / pow(l_m, 2.0);
double tmp;
if (k <= 6e-161) {
tmp = 2.0 / (((t_m * pow(k, 2.0)) / cos(k)) * 0.0);
} else if (k <= 2e+147) {
tmp = 2.0 / (t_2 * (pow(k, 2.0) * (t_m / cos(k))));
} else {
tmp = 2.0 / ((exp((log(t_m) + (2.0 * log(k)))) / cos(k)) * t_2);
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = (sin(k) ** 2.0d0) / (l_m ** 2.0d0)
if (k <= 6d-161) then
tmp = 2.0d0 / (((t_m * (k ** 2.0d0)) / cos(k)) * 0.0d0)
else if (k <= 2d+147) then
tmp = 2.0d0 / (t_2 * ((k ** 2.0d0) * (t_m / cos(k))))
else
tmp = 2.0d0 / ((exp((log(t_m) + (2.0d0 * log(k)))) / cos(k)) * t_2)
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double t_2 = Math.pow(Math.sin(k), 2.0) / Math.pow(l_m, 2.0);
double tmp;
if (k <= 6e-161) {
tmp = 2.0 / (((t_m * Math.pow(k, 2.0)) / Math.cos(k)) * 0.0);
} else if (k <= 2e+147) {
tmp = 2.0 / (t_2 * (Math.pow(k, 2.0) * (t_m / Math.cos(k))));
} else {
tmp = 2.0 / ((Math.exp((Math.log(t_m) + (2.0 * Math.log(k)))) / Math.cos(k)) * t_2);
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): t_2 = math.pow(math.sin(k), 2.0) / math.pow(l_m, 2.0) tmp = 0 if k <= 6e-161: tmp = 2.0 / (((t_m * math.pow(k, 2.0)) / math.cos(k)) * 0.0) elif k <= 2e+147: tmp = 2.0 / (t_2 * (math.pow(k, 2.0) * (t_m / math.cos(k)))) else: tmp = 2.0 / ((math.exp((math.log(t_m) + (2.0 * math.log(k)))) / math.cos(k)) * t_2) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) t_2 = Float64((sin(k) ^ 2.0) / (l_m ^ 2.0)) tmp = 0.0 if (k <= 6e-161) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) / cos(k)) * 0.0)); elseif (k <= 2e+147) tmp = Float64(2.0 / Float64(t_2 * Float64((k ^ 2.0) * Float64(t_m / cos(k))))); else tmp = Float64(2.0 / Float64(Float64(exp(Float64(log(t_m) + Float64(2.0 * log(k)))) / cos(k)) * t_2)); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) t_2 = (sin(k) ^ 2.0) / (l_m ^ 2.0); tmp = 0.0; if (k <= 6e-161) tmp = 2.0 / (((t_m * (k ^ 2.0)) / cos(k)) * 0.0); elseif (k <= 2e+147) tmp = 2.0 / (t_2 * ((k ^ 2.0) * (t_m / cos(k)))); else tmp = 2.0 / ((exp((log(t_m) + (2.0 * log(k)))) / cos(k)) * t_2); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := Block[{t$95$2 = N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] / N[Power[l$95$m, 2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k, 6e-161], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * 0.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 2e+147], N[(2.0 / N[(t$95$2 * N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Exp[N[(N[Log[t$95$m], $MachinePrecision] + N[(2.0 * N[Log[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{{\sin k}^{2}}{{l\_m}^{2}}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 6 \cdot 10^{-161}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{2}}{\cos k} \cdot 0}\\
\mathbf{elif}\;k \leq 2 \cdot 10^{+147}:\\
\;\;\;\;\frac{2}{t\_2 \cdot \left({k}^{2} \cdot \frac{t\_m}{\cos k}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{e^{\log t\_m + 2 \cdot \log k}}{\cos k} \cdot t\_2}\\
\end{array}
\end{array}
\end{array}
if k < 5.99999999999999977e-161Initial program 29.8%
Simplified29.8%
Taylor expanded in t around 0 75.1%
associate-*r*75.1%
*-commutative75.1%
times-frac77.4%
Simplified77.4%
pow277.4%
add-log-exp62.3%
div-inv62.2%
exp-prod67.8%
pow267.8%
pow-flip67.8%
metadata-eval67.8%
Applied egg-rr67.8%
Taylor expanded in k around 0 47.6%
if 5.99999999999999977e-161 < k < 2e147Initial program 23.8%
Simplified23.8%
Taylor expanded in t around 0 75.3%
associate-*r*75.3%
*-commutative75.3%
times-frac81.0%
Simplified81.0%
add-log-exp34.6%
*-commutative34.6%
exp-prod24.9%
Applied egg-rr24.9%
Taylor expanded in t around 0 81.0%
associate-/l*81.1%
Simplified81.1%
if 2e147 < k Initial program 38.1%
Simplified38.1%
Taylor expanded in t around 0 67.5%
associate-*r*67.5%
*-commutative67.5%
times-frac67.5%
Simplified67.5%
add-log-exp67.5%
*-commutative67.5%
exp-prod38.7%
Applied egg-rr38.7%
pow-exp67.5%
rem-log-exp67.5%
add-exp-log43.4%
pow-to-exp43.4%
prod-exp47.4%
rem-log-exp43.4%
pow-to-exp43.4%
log-pow47.4%
Applied egg-rr47.4%
Final simplification56.5%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(*
t_s
(if (<= k 1.2e+147)
(/
2.0
(*
(/ (* t_m (pow k 2.0)) (cos k))
(exp (+ (log (pow (sin k) 2.0)) (* -2.0 (log l_m))))))
(/
2.0
(pow
(*
(/ t_m (pow (cbrt l_m) 2.0))
(cbrt (* (sin k) (* (tan k) (pow (/ k t_m) 2.0)))))
3.0)))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 1.2e+147) {
tmp = 2.0 / (((t_m * pow(k, 2.0)) / cos(k)) * exp((log(pow(sin(k), 2.0)) + (-2.0 * log(l_m)))));
} else {
tmp = 2.0 / pow(((t_m / pow(cbrt(l_m), 2.0)) * cbrt((sin(k) * (tan(k) * pow((k / t_m), 2.0))))), 3.0);
}
return t_s * tmp;
}
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 1.2e+147) {
tmp = 2.0 / (((t_m * Math.pow(k, 2.0)) / Math.cos(k)) * Math.exp((Math.log(Math.pow(Math.sin(k), 2.0)) + (-2.0 * Math.log(l_m)))));
} else {
tmp = 2.0 / Math.pow(((t_m / Math.pow(Math.cbrt(l_m), 2.0)) * Math.cbrt((Math.sin(k) * (Math.tan(k) * Math.pow((k / t_m), 2.0))))), 3.0);
}
return t_s * tmp;
}
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) tmp = 0.0 if (k <= 1.2e+147) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) / cos(k)) * exp(Float64(log((sin(k) ^ 2.0)) + Float64(-2.0 * log(l_m)))))); else tmp = Float64(2.0 / (Float64(Float64(t_m / (cbrt(l_m) ^ 2.0)) * cbrt(Float64(sin(k) * Float64(tan(k) * (Float64(k / t_m) ^ 2.0))))) ^ 3.0)); end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := N[(t$95$s * If[LessEqual[k, 1.2e+147], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * N[Exp[N[(N[Log[N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision] + N[(-2.0 * N[Log[l$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(t$95$m / N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[(N[Sin[k], $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 1.2 \cdot 10^{+147}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{2}}{\cos k} \cdot e^{\log \left({\sin k}^{2}\right) + -2 \cdot \log l\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{t\_m}{{\left(\sqrt[3]{l\_m}\right)}^{2}} \cdot \sqrt[3]{\sin k \cdot \left(\tan k \cdot {\left(\frac{k}{t\_m}\right)}^{2}\right)}\right)}^{3}}\\
\end{array}
\end{array}
if k < 1.20000000000000001e147Initial program 28.0%
Simplified28.0%
Taylor expanded in t around 0 75.1%
associate-*r*75.1%
*-commutative75.1%
times-frac78.4%
Simplified78.4%
pow278.4%
add-log-exp58.1%
div-inv58.1%
exp-prod61.1%
pow261.1%
pow-flip61.1%
metadata-eval61.1%
Applied egg-rr61.1%
pow-exp58.1%
rem-log-exp78.4%
pow-to-exp34.9%
pow-to-exp15.4%
prod-exp17.7%
rem-log-exp17.7%
pow-to-exp42.0%
rem-log-exp37.9%
pow-to-exp76.8%
log-pow42.0%
Applied egg-rr42.0%
if 1.20000000000000001e147 < k Initial program 38.1%
Simplified38.1%
add-cube-cbrt38.1%
pow338.1%
Applied egg-rr81.3%
Final simplification45.2%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(let* ((t_2 (pow (sin k) 2.0)))
(*
t_s
(if (<= k 2e+147)
(/
2.0
(*
(/ (* t_m (pow k 2.0)) (cos k))
(exp (+ (log t_2) (* -2.0 (log l_m))))))
(/
2.0
(*
(/ (exp (+ (log t_m) (* 2.0 (log k)))) (cos k))
(/ t_2 (pow l_m 2.0))))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double t_2 = pow(sin(k), 2.0);
double tmp;
if (k <= 2e+147) {
tmp = 2.0 / (((t_m * pow(k, 2.0)) / cos(k)) * exp((log(t_2) + (-2.0 * log(l_m)))));
} else {
tmp = 2.0 / ((exp((log(t_m) + (2.0 * log(k)))) / cos(k)) * (t_2 / pow(l_m, 2.0)));
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = sin(k) ** 2.0d0
if (k <= 2d+147) then
tmp = 2.0d0 / (((t_m * (k ** 2.0d0)) / cos(k)) * exp((log(t_2) + ((-2.0d0) * log(l_m)))))
else
tmp = 2.0d0 / ((exp((log(t_m) + (2.0d0 * log(k)))) / cos(k)) * (t_2 / (l_m ** 2.0d0)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double t_2 = Math.pow(Math.sin(k), 2.0);
double tmp;
if (k <= 2e+147) {
tmp = 2.0 / (((t_m * Math.pow(k, 2.0)) / Math.cos(k)) * Math.exp((Math.log(t_2) + (-2.0 * Math.log(l_m)))));
} else {
tmp = 2.0 / ((Math.exp((Math.log(t_m) + (2.0 * Math.log(k)))) / Math.cos(k)) * (t_2 / Math.pow(l_m, 2.0)));
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): t_2 = math.pow(math.sin(k), 2.0) tmp = 0 if k <= 2e+147: tmp = 2.0 / (((t_m * math.pow(k, 2.0)) / math.cos(k)) * math.exp((math.log(t_2) + (-2.0 * math.log(l_m))))) else: tmp = 2.0 / ((math.exp((math.log(t_m) + (2.0 * math.log(k)))) / math.cos(k)) * (t_2 / math.pow(l_m, 2.0))) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) t_2 = sin(k) ^ 2.0 tmp = 0.0 if (k <= 2e+147) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) / cos(k)) * exp(Float64(log(t_2) + Float64(-2.0 * log(l_m)))))); else tmp = Float64(2.0 / Float64(Float64(exp(Float64(log(t_m) + Float64(2.0 * log(k)))) / cos(k)) * Float64(t_2 / (l_m ^ 2.0)))); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) t_2 = sin(k) ^ 2.0; tmp = 0.0; if (k <= 2e+147) tmp = 2.0 / (((t_m * (k ^ 2.0)) / cos(k)) * exp((log(t_2) + (-2.0 * log(l_m))))); else tmp = 2.0 / ((exp((log(t_m) + (2.0 * log(k)))) / cos(k)) * (t_2 / (l_m ^ 2.0))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := Block[{t$95$2 = N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k, 2e+147], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * N[Exp[N[(N[Log[t$95$2], $MachinePrecision] + N[(-2.0 * N[Log[l$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Exp[N[(N[Log[t$95$m], $MachinePrecision] + N[(2.0 * N[Log[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * N[(t$95$2 / N[Power[l$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := {\sin k}^{2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 2 \cdot 10^{+147}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{2}}{\cos k} \cdot e^{\log t\_2 + -2 \cdot \log l\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{e^{\log t\_m + 2 \cdot \log k}}{\cos k} \cdot \frac{t\_2}{{l\_m}^{2}}}\\
\end{array}
\end{array}
\end{array}
if k < 2e147Initial program 28.0%
Simplified28.0%
Taylor expanded in t around 0 75.1%
associate-*r*75.1%
*-commutative75.1%
times-frac78.4%
Simplified78.4%
pow278.4%
add-log-exp58.1%
div-inv58.1%
exp-prod61.1%
pow261.1%
pow-flip61.1%
metadata-eval61.1%
Applied egg-rr61.1%
pow-exp58.1%
rem-log-exp78.4%
pow-to-exp34.9%
pow-to-exp15.4%
prod-exp17.7%
rem-log-exp17.7%
pow-to-exp42.0%
rem-log-exp37.9%
pow-to-exp76.8%
log-pow42.0%
Applied egg-rr42.0%
if 2e147 < k Initial program 38.1%
Simplified38.1%
Taylor expanded in t around 0 67.5%
associate-*r*67.5%
*-commutative67.5%
times-frac67.5%
Simplified67.5%
add-log-exp67.5%
*-commutative67.5%
exp-prod38.7%
Applied egg-rr38.7%
pow-exp67.5%
rem-log-exp67.5%
add-exp-log43.4%
pow-to-exp43.4%
prod-exp47.4%
rem-log-exp43.4%
pow-to-exp43.4%
log-pow47.4%
Applied egg-rr47.4%
Final simplification42.5%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(*
t_s
(if (<= k 8e-161)
(/ 2.0 (* (/ (* t_m (pow k 2.0)) (cos k)) 0.0))
(/
2.0
(*
(/ (pow (sin k) 2.0) (pow l_m 2.0))
(* (pow k 2.0) (/ t_m (cos k))))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 8e-161) {
tmp = 2.0 / (((t_m * pow(k, 2.0)) / cos(k)) * 0.0);
} else {
tmp = 2.0 / ((pow(sin(k), 2.0) / pow(l_m, 2.0)) * (pow(k, 2.0) * (t_m / cos(k))));
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 8d-161) then
tmp = 2.0d0 / (((t_m * (k ** 2.0d0)) / cos(k)) * 0.0d0)
else
tmp = 2.0d0 / (((sin(k) ** 2.0d0) / (l_m ** 2.0d0)) * ((k ** 2.0d0) * (t_m / cos(k))))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 8e-161) {
tmp = 2.0 / (((t_m * Math.pow(k, 2.0)) / Math.cos(k)) * 0.0);
} else {
tmp = 2.0 / ((Math.pow(Math.sin(k), 2.0) / Math.pow(l_m, 2.0)) * (Math.pow(k, 2.0) * (t_m / Math.cos(k))));
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): tmp = 0 if k <= 8e-161: tmp = 2.0 / (((t_m * math.pow(k, 2.0)) / math.cos(k)) * 0.0) else: tmp = 2.0 / ((math.pow(math.sin(k), 2.0) / math.pow(l_m, 2.0)) * (math.pow(k, 2.0) * (t_m / math.cos(k)))) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) tmp = 0.0 if (k <= 8e-161) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) / cos(k)) * 0.0)); else tmp = Float64(2.0 / Float64(Float64((sin(k) ^ 2.0) / (l_m ^ 2.0)) * Float64((k ^ 2.0) * Float64(t_m / cos(k))))); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) tmp = 0.0; if (k <= 8e-161) tmp = 2.0 / (((t_m * (k ^ 2.0)) / cos(k)) * 0.0); else tmp = 2.0 / (((sin(k) ^ 2.0) / (l_m ^ 2.0)) * ((k ^ 2.0) * (t_m / cos(k)))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := N[(t$95$s * If[LessEqual[k, 8e-161], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * 0.0), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] / N[Power[l$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 8 \cdot 10^{-161}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{2}}{\cos k} \cdot 0}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{{\sin k}^{2}}{{l\_m}^{2}} \cdot \left({k}^{2} \cdot \frac{t\_m}{\cos k}\right)}\\
\end{array}
\end{array}
if k < 8.00000000000000022e-161Initial program 29.8%
Simplified29.8%
Taylor expanded in t around 0 75.1%
associate-*r*75.1%
*-commutative75.1%
times-frac77.4%
Simplified77.4%
pow277.4%
add-log-exp62.3%
div-inv62.2%
exp-prod67.8%
pow267.8%
pow-flip67.8%
metadata-eval67.8%
Applied egg-rr67.8%
Taylor expanded in k around 0 47.6%
if 8.00000000000000022e-161 < k Initial program 27.2%
Simplified27.2%
Taylor expanded in t around 0 73.4%
associate-*r*73.5%
*-commutative73.5%
times-frac77.8%
Simplified77.8%
add-log-exp42.4%
*-commutative42.4%
exp-prod28.1%
Applied egg-rr28.1%
Taylor expanded in t around 0 77.8%
associate-/l*77.9%
Simplified77.9%
Final simplification58.1%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(let* ((t_2 (pow (sin k) 2.0)) (t_3 (* t_m (pow k 2.0))))
(*
t_s
(if (<= k 3.3e-161)
(/ 2.0 (* (/ t_3 (cos k)) 0.0))
(if (<= k 2e-24)
(/ 2.0 (* t_3 (/ t_2 (pow l_m 2.0))))
(* (/ (* 2.0 (cos k)) (* t_3 t_2)) (* l_m l_m)))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double t_2 = pow(sin(k), 2.0);
double t_3 = t_m * pow(k, 2.0);
double tmp;
if (k <= 3.3e-161) {
tmp = 2.0 / ((t_3 / cos(k)) * 0.0);
} else if (k <= 2e-24) {
tmp = 2.0 / (t_3 * (t_2 / pow(l_m, 2.0)));
} else {
tmp = ((2.0 * cos(k)) / (t_3 * t_2)) * (l_m * l_m);
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_2 = sin(k) ** 2.0d0
t_3 = t_m * (k ** 2.0d0)
if (k <= 3.3d-161) then
tmp = 2.0d0 / ((t_3 / cos(k)) * 0.0d0)
else if (k <= 2d-24) then
tmp = 2.0d0 / (t_3 * (t_2 / (l_m ** 2.0d0)))
else
tmp = ((2.0d0 * cos(k)) / (t_3 * t_2)) * (l_m * l_m)
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double t_2 = Math.pow(Math.sin(k), 2.0);
double t_3 = t_m * Math.pow(k, 2.0);
double tmp;
if (k <= 3.3e-161) {
tmp = 2.0 / ((t_3 / Math.cos(k)) * 0.0);
} else if (k <= 2e-24) {
tmp = 2.0 / (t_3 * (t_2 / Math.pow(l_m, 2.0)));
} else {
tmp = ((2.0 * Math.cos(k)) / (t_3 * t_2)) * (l_m * l_m);
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): t_2 = math.pow(math.sin(k), 2.0) t_3 = t_m * math.pow(k, 2.0) tmp = 0 if k <= 3.3e-161: tmp = 2.0 / ((t_3 / math.cos(k)) * 0.0) elif k <= 2e-24: tmp = 2.0 / (t_3 * (t_2 / math.pow(l_m, 2.0))) else: tmp = ((2.0 * math.cos(k)) / (t_3 * t_2)) * (l_m * l_m) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) t_2 = sin(k) ^ 2.0 t_3 = Float64(t_m * (k ^ 2.0)) tmp = 0.0 if (k <= 3.3e-161) tmp = Float64(2.0 / Float64(Float64(t_3 / cos(k)) * 0.0)); elseif (k <= 2e-24) tmp = Float64(2.0 / Float64(t_3 * Float64(t_2 / (l_m ^ 2.0)))); else tmp = Float64(Float64(Float64(2.0 * cos(k)) / Float64(t_3 * t_2)) * Float64(l_m * l_m)); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) t_2 = sin(k) ^ 2.0; t_3 = t_m * (k ^ 2.0); tmp = 0.0; if (k <= 3.3e-161) tmp = 2.0 / ((t_3 / cos(k)) * 0.0); elseif (k <= 2e-24) tmp = 2.0 / (t_3 * (t_2 / (l_m ^ 2.0))); else tmp = ((2.0 * cos(k)) / (t_3 * t_2)) * (l_m * l_m); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := Block[{t$95$2 = N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$3 = N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k, 3.3e-161], N[(2.0 / N[(N[(t$95$3 / N[Cos[k], $MachinePrecision]), $MachinePrecision] * 0.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 2e-24], N[(2.0 / N[(t$95$3 * N[(t$95$2 / N[Power[l$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 * N[Cos[k], $MachinePrecision]), $MachinePrecision] / N[(t$95$3 * t$95$2), $MachinePrecision]), $MachinePrecision] * N[(l$95$m * l$95$m), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := {\sin k}^{2}\\
t_3 := t\_m \cdot {k}^{2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 3.3 \cdot 10^{-161}:\\
\;\;\;\;\frac{2}{\frac{t\_3}{\cos k} \cdot 0}\\
\mathbf{elif}\;k \leq 2 \cdot 10^{-24}:\\
\;\;\;\;\frac{2}{t\_3 \cdot \frac{t\_2}{{l\_m}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \cos k}{t\_3 \cdot t\_2} \cdot \left(l\_m \cdot l\_m\right)\\
\end{array}
\end{array}
\end{array}
if k < 3.2999999999999998e-161Initial program 29.8%
Simplified29.8%
Taylor expanded in t around 0 75.1%
associate-*r*75.1%
*-commutative75.1%
times-frac77.4%
Simplified77.4%
pow277.4%
add-log-exp62.3%
div-inv62.2%
exp-prod67.8%
pow267.8%
pow-flip67.8%
metadata-eval67.8%
Applied egg-rr67.8%
Taylor expanded in k around 0 47.6%
if 3.2999999999999998e-161 < k < 1.99999999999999985e-24Initial program 25.2%
Simplified25.2%
Taylor expanded in t around 0 65.9%
associate-*r*65.9%
*-commutative65.9%
times-frac78.0%
Simplified78.0%
Taylor expanded in k around 0 78.0%
if 1.99999999999999985e-24 < k Initial program 28.3%
Simplified40.7%
Taylor expanded in t around 0 77.6%
associate-*r/77.6%
associate-*r*77.7%
Simplified77.7%
Final simplification58.1%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(let* ((t_2 (pow (sin k) 2.0)) (t_3 (* t_m (pow k 2.0))))
(*
t_s
(if (<= k 5.5e-161)
(/ 2.0 (* (/ t_3 (cos k)) 0.0))
(if (<= k 1e-24)
(/ 2.0 (* t_3 (/ t_2 (pow l_m 2.0))))
(* (* l_m l_m) (/ 2.0 (/ (* (* k k) (* t_m t_2)) (cos k)))))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double t_2 = pow(sin(k), 2.0);
double t_3 = t_m * pow(k, 2.0);
double tmp;
if (k <= 5.5e-161) {
tmp = 2.0 / ((t_3 / cos(k)) * 0.0);
} else if (k <= 1e-24) {
tmp = 2.0 / (t_3 * (t_2 / pow(l_m, 2.0)));
} else {
tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * t_2)) / cos(k)));
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_2 = sin(k) ** 2.0d0
t_3 = t_m * (k ** 2.0d0)
if (k <= 5.5d-161) then
tmp = 2.0d0 / ((t_3 / cos(k)) * 0.0d0)
else if (k <= 1d-24) then
tmp = 2.0d0 / (t_3 * (t_2 / (l_m ** 2.0d0)))
else
tmp = (l_m * l_m) * (2.0d0 / (((k * k) * (t_m * t_2)) / cos(k)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double t_2 = Math.pow(Math.sin(k), 2.0);
double t_3 = t_m * Math.pow(k, 2.0);
double tmp;
if (k <= 5.5e-161) {
tmp = 2.0 / ((t_3 / Math.cos(k)) * 0.0);
} else if (k <= 1e-24) {
tmp = 2.0 / (t_3 * (t_2 / Math.pow(l_m, 2.0)));
} else {
tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * t_2)) / Math.cos(k)));
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): t_2 = math.pow(math.sin(k), 2.0) t_3 = t_m * math.pow(k, 2.0) tmp = 0 if k <= 5.5e-161: tmp = 2.0 / ((t_3 / math.cos(k)) * 0.0) elif k <= 1e-24: tmp = 2.0 / (t_3 * (t_2 / math.pow(l_m, 2.0))) else: tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * t_2)) / math.cos(k))) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) t_2 = sin(k) ^ 2.0 t_3 = Float64(t_m * (k ^ 2.0)) tmp = 0.0 if (k <= 5.5e-161) tmp = Float64(2.0 / Float64(Float64(t_3 / cos(k)) * 0.0)); elseif (k <= 1e-24) tmp = Float64(2.0 / Float64(t_3 * Float64(t_2 / (l_m ^ 2.0)))); else tmp = Float64(Float64(l_m * l_m) * Float64(2.0 / Float64(Float64(Float64(k * k) * Float64(t_m * t_2)) / cos(k)))); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) t_2 = sin(k) ^ 2.0; t_3 = t_m * (k ^ 2.0); tmp = 0.0; if (k <= 5.5e-161) tmp = 2.0 / ((t_3 / cos(k)) * 0.0); elseif (k <= 1e-24) tmp = 2.0 / (t_3 * (t_2 / (l_m ^ 2.0))); else tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * t_2)) / cos(k))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := Block[{t$95$2 = N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$3 = N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k, 5.5e-161], N[(2.0 / N[(N[(t$95$3 / N[Cos[k], $MachinePrecision]), $MachinePrecision] * 0.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 1e-24], N[(2.0 / N[(t$95$3 * N[(t$95$2 / N[Power[l$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(2.0 / N[(N[(N[(k * k), $MachinePrecision] * N[(t$95$m * t$95$2), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := {\sin k}^{2}\\
t_3 := t\_m \cdot {k}^{2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 5.5 \cdot 10^{-161}:\\
\;\;\;\;\frac{2}{\frac{t\_3}{\cos k} \cdot 0}\\
\mathbf{elif}\;k \leq 10^{-24}:\\
\;\;\;\;\frac{2}{t\_3 \cdot \frac{t\_2}{{l\_m}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\left(l\_m \cdot l\_m\right) \cdot \frac{2}{\frac{\left(k \cdot k\right) \cdot \left(t\_m \cdot t\_2\right)}{\cos k}}\\
\end{array}
\end{array}
\end{array}
if k < 5.5e-161Initial program 29.8%
Simplified29.8%
Taylor expanded in t around 0 75.1%
associate-*r*75.1%
*-commutative75.1%
times-frac77.4%
Simplified77.4%
pow277.4%
add-log-exp62.3%
div-inv62.2%
exp-prod67.8%
pow267.8%
pow-flip67.8%
metadata-eval67.8%
Applied egg-rr67.8%
Taylor expanded in k around 0 47.6%
if 5.5e-161 < k < 9.99999999999999924e-25Initial program 25.2%
Simplified25.2%
Taylor expanded in t around 0 65.9%
associate-*r*65.9%
*-commutative65.9%
times-frac78.0%
Simplified78.0%
Taylor expanded in k around 0 78.0%
if 9.99999999999999924e-25 < k Initial program 28.3%
Simplified40.7%
Taylor expanded in t around 0 77.6%
unpow277.7%
Applied egg-rr77.6%
Final simplification58.1%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(*
t_s
(if (<= k 1.8e-161)
(/ 2.0 (* (/ (* t_m (pow k 2.0)) (cos k)) 0.0))
(/
2.0
(* (/ (pow (sin k) 2.0) (pow l_m 2.0)) (/ (* t_m (* k k)) (cos k)))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 1.8e-161) {
tmp = 2.0 / (((t_m * pow(k, 2.0)) / cos(k)) * 0.0);
} else {
tmp = 2.0 / ((pow(sin(k), 2.0) / pow(l_m, 2.0)) * ((t_m * (k * k)) / cos(k)));
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 1.8d-161) then
tmp = 2.0d0 / (((t_m * (k ** 2.0d0)) / cos(k)) * 0.0d0)
else
tmp = 2.0d0 / (((sin(k) ** 2.0d0) / (l_m ** 2.0d0)) * ((t_m * (k * k)) / cos(k)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 1.8e-161) {
tmp = 2.0 / (((t_m * Math.pow(k, 2.0)) / Math.cos(k)) * 0.0);
} else {
tmp = 2.0 / ((Math.pow(Math.sin(k), 2.0) / Math.pow(l_m, 2.0)) * ((t_m * (k * k)) / Math.cos(k)));
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): tmp = 0 if k <= 1.8e-161: tmp = 2.0 / (((t_m * math.pow(k, 2.0)) / math.cos(k)) * 0.0) else: tmp = 2.0 / ((math.pow(math.sin(k), 2.0) / math.pow(l_m, 2.0)) * ((t_m * (k * k)) / math.cos(k))) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) tmp = 0.0 if (k <= 1.8e-161) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) / cos(k)) * 0.0)); else tmp = Float64(2.0 / Float64(Float64((sin(k) ^ 2.0) / (l_m ^ 2.0)) * Float64(Float64(t_m * Float64(k * k)) / cos(k)))); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) tmp = 0.0; if (k <= 1.8e-161) tmp = 2.0 / (((t_m * (k ^ 2.0)) / cos(k)) * 0.0); else tmp = 2.0 / (((sin(k) ^ 2.0) / (l_m ^ 2.0)) * ((t_m * (k * k)) / cos(k))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := N[(t$95$s * If[LessEqual[k, 1.8e-161], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * 0.0), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] / N[Power[l$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(t$95$m * N[(k * k), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 1.8 \cdot 10^{-161}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{2}}{\cos k} \cdot 0}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{{\sin k}^{2}}{{l\_m}^{2}} \cdot \frac{t\_m \cdot \left(k \cdot k\right)}{\cos k}}\\
\end{array}
\end{array}
if k < 1.80000000000000009e-161Initial program 29.8%
Simplified29.8%
Taylor expanded in t around 0 75.1%
associate-*r*75.1%
*-commutative75.1%
times-frac77.4%
Simplified77.4%
pow277.4%
add-log-exp62.3%
div-inv62.2%
exp-prod67.8%
pow267.8%
pow-flip67.8%
metadata-eval67.8%
Applied egg-rr67.8%
Taylor expanded in k around 0 47.6%
if 1.80000000000000009e-161 < k Initial program 27.2%
Simplified27.2%
Taylor expanded in t around 0 73.4%
associate-*r*73.5%
*-commutative73.5%
times-frac77.8%
Simplified77.8%
unpow277.8%
Applied egg-rr77.8%
Final simplification58.1%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(*
t_s
(if (<= k 1e-160)
(/ 2.0 (* (/ (* t_m (pow k 2.0)) (cos k)) 0.0))
(*
(/ (pow l_m 2.0) (pow (sin k) 2.0))
(/ (* 2.0 (cos k)) (* t_m (* k k)))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 1e-160) {
tmp = 2.0 / (((t_m * pow(k, 2.0)) / cos(k)) * 0.0);
} else {
tmp = (pow(l_m, 2.0) / pow(sin(k), 2.0)) * ((2.0 * cos(k)) / (t_m * (k * k)));
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 1d-160) then
tmp = 2.0d0 / (((t_m * (k ** 2.0d0)) / cos(k)) * 0.0d0)
else
tmp = ((l_m ** 2.0d0) / (sin(k) ** 2.0d0)) * ((2.0d0 * cos(k)) / (t_m * (k * k)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 1e-160) {
tmp = 2.0 / (((t_m * Math.pow(k, 2.0)) / Math.cos(k)) * 0.0);
} else {
tmp = (Math.pow(l_m, 2.0) / Math.pow(Math.sin(k), 2.0)) * ((2.0 * Math.cos(k)) / (t_m * (k * k)));
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): tmp = 0 if k <= 1e-160: tmp = 2.0 / (((t_m * math.pow(k, 2.0)) / math.cos(k)) * 0.0) else: tmp = (math.pow(l_m, 2.0) / math.pow(math.sin(k), 2.0)) * ((2.0 * math.cos(k)) / (t_m * (k * k))) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) tmp = 0.0 if (k <= 1e-160) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) / cos(k)) * 0.0)); else tmp = Float64(Float64((l_m ^ 2.0) / (sin(k) ^ 2.0)) * Float64(Float64(2.0 * cos(k)) / Float64(t_m * Float64(k * k)))); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) tmp = 0.0; if (k <= 1e-160) tmp = 2.0 / (((t_m * (k ^ 2.0)) / cos(k)) * 0.0); else tmp = ((l_m ^ 2.0) / (sin(k) ^ 2.0)) * ((2.0 * cos(k)) / (t_m * (k * k))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := N[(t$95$s * If[LessEqual[k, 1e-160], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * 0.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[Power[l$95$m, 2.0], $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * N[Cos[k], $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 10^{-160}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{2}}{\cos k} \cdot 0}\\
\mathbf{else}:\\
\;\;\;\;\frac{{l\_m}^{2}}{{\sin k}^{2}} \cdot \frac{2 \cdot \cos k}{t\_m \cdot \left(k \cdot k\right)}\\
\end{array}
\end{array}
if k < 9.9999999999999999e-161Initial program 29.8%
Simplified29.8%
Taylor expanded in t around 0 75.1%
associate-*r*75.1%
*-commutative75.1%
times-frac77.4%
Simplified77.4%
pow277.4%
add-log-exp62.3%
div-inv62.2%
exp-prod67.8%
pow267.8%
pow-flip67.8%
metadata-eval67.8%
Applied egg-rr67.8%
Taylor expanded in k around 0 47.6%
if 9.9999999999999999e-161 < k Initial program 27.2%
Simplified27.2%
Applied egg-rr82.4%
associate-/r/82.4%
associate-/r*82.4%
associate-/r/83.5%
Simplified83.5%
Taylor expanded in k around inf 73.3%
associate-*r*73.3%
*-commutative73.3%
*-commutative73.3%
times-frac77.5%
*-commutative77.5%
unpow277.5%
rem-square-sqrt77.8%
Simplified77.8%
unpow277.8%
Applied egg-rr77.8%
Final simplification58.1%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(*
t_s
(if (<= k 1.26e-129)
(/ 2.0 (* (/ (* t_m (pow k 2.0)) (cos k)) 0.0))
(*
(* l_m l_m)
(/ 2.0 (/ (* (* k k) (* t_m (pow (sin k) 2.0))) (cos k)))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 1.26e-129) {
tmp = 2.0 / (((t_m * pow(k, 2.0)) / cos(k)) * 0.0);
} else {
tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * pow(sin(k), 2.0))) / cos(k)));
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 1.26d-129) then
tmp = 2.0d0 / (((t_m * (k ** 2.0d0)) / cos(k)) * 0.0d0)
else
tmp = (l_m * l_m) * (2.0d0 / (((k * k) * (t_m * (sin(k) ** 2.0d0))) / cos(k)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 1.26e-129) {
tmp = 2.0 / (((t_m * Math.pow(k, 2.0)) / Math.cos(k)) * 0.0);
} else {
tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * Math.pow(Math.sin(k), 2.0))) / Math.cos(k)));
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): tmp = 0 if k <= 1.26e-129: tmp = 2.0 / (((t_m * math.pow(k, 2.0)) / math.cos(k)) * 0.0) else: tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * math.pow(math.sin(k), 2.0))) / math.cos(k))) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) tmp = 0.0 if (k <= 1.26e-129) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) / cos(k)) * 0.0)); else tmp = Float64(Float64(l_m * l_m) * Float64(2.0 / Float64(Float64(Float64(k * k) * Float64(t_m * (sin(k) ^ 2.0))) / cos(k)))); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) tmp = 0.0; if (k <= 1.26e-129) tmp = 2.0 / (((t_m * (k ^ 2.0)) / cos(k)) * 0.0); else tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * (sin(k) ^ 2.0))) / cos(k))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := N[(t$95$s * If[LessEqual[k, 1.26e-129], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * 0.0), $MachinePrecision]), $MachinePrecision], N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(2.0 / N[(N[(N[(k * k), $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 1.26 \cdot 10^{-129}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{2}}{\cos k} \cdot 0}\\
\mathbf{else}:\\
\;\;\;\;\left(l\_m \cdot l\_m\right) \cdot \frac{2}{\frac{\left(k \cdot k\right) \cdot \left(t\_m \cdot {\sin k}^{2}\right)}{\cos k}}\\
\end{array}
\end{array}
if k < 1.2599999999999999e-129Initial program 29.9%
Simplified29.9%
Taylor expanded in t around 0 75.8%
associate-*r*75.8%
*-commutative75.8%
times-frac78.0%
Simplified78.0%
pow278.0%
add-log-exp63.7%
div-inv63.6%
exp-prod68.9%
pow268.9%
pow-flip68.9%
metadata-eval68.9%
Applied egg-rr68.9%
Taylor expanded in k around 0 49.7%
if 1.2599999999999999e-129 < k Initial program 26.5%
Simplified36.4%
Taylor expanded in t around 0 71.7%
unpow276.5%
Applied egg-rr71.7%
Final simplification56.6%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(let* ((t_2 (* t_m (pow k 2.0))))
(*
t_s
(if (<= k 1.35e-129)
(/ 2.0 (* (/ t_2 (cos k)) 0.0))
(if (<= k 2.05e-5)
(* (* l_m l_m) (/ 2.0 (/ (* t_2 (* k k)) (cos k))))
(*
(* l_m l_m)
(/
2.0
(/
(* (* k k) (* t_m (- 0.5 (/ (cos (* 2.0 k)) 2.0))))
(cos k)))))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double t_2 = t_m * pow(k, 2.0);
double tmp;
if (k <= 1.35e-129) {
tmp = 2.0 / ((t_2 / cos(k)) * 0.0);
} else if (k <= 2.05e-5) {
tmp = (l_m * l_m) * (2.0 / ((t_2 * (k * k)) / cos(k)));
} else {
tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * (0.5 - (cos((2.0 * k)) / 2.0)))) / cos(k)));
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = t_m * (k ** 2.0d0)
if (k <= 1.35d-129) then
tmp = 2.0d0 / ((t_2 / cos(k)) * 0.0d0)
else if (k <= 2.05d-5) then
tmp = (l_m * l_m) * (2.0d0 / ((t_2 * (k * k)) / cos(k)))
else
tmp = (l_m * l_m) * (2.0d0 / (((k * k) * (t_m * (0.5d0 - (cos((2.0d0 * k)) / 2.0d0)))) / cos(k)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double t_2 = t_m * Math.pow(k, 2.0);
double tmp;
if (k <= 1.35e-129) {
tmp = 2.0 / ((t_2 / Math.cos(k)) * 0.0);
} else if (k <= 2.05e-5) {
tmp = (l_m * l_m) * (2.0 / ((t_2 * (k * k)) / Math.cos(k)));
} else {
tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * (0.5 - (Math.cos((2.0 * k)) / 2.0)))) / Math.cos(k)));
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): t_2 = t_m * math.pow(k, 2.0) tmp = 0 if k <= 1.35e-129: tmp = 2.0 / ((t_2 / math.cos(k)) * 0.0) elif k <= 2.05e-5: tmp = (l_m * l_m) * (2.0 / ((t_2 * (k * k)) / math.cos(k))) else: tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * (0.5 - (math.cos((2.0 * k)) / 2.0)))) / math.cos(k))) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) t_2 = Float64(t_m * (k ^ 2.0)) tmp = 0.0 if (k <= 1.35e-129) tmp = Float64(2.0 / Float64(Float64(t_2 / cos(k)) * 0.0)); elseif (k <= 2.05e-5) tmp = Float64(Float64(l_m * l_m) * Float64(2.0 / Float64(Float64(t_2 * Float64(k * k)) / cos(k)))); else tmp = Float64(Float64(l_m * l_m) * Float64(2.0 / Float64(Float64(Float64(k * k) * Float64(t_m * Float64(0.5 - Float64(cos(Float64(2.0 * k)) / 2.0)))) / cos(k)))); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) t_2 = t_m * (k ^ 2.0); tmp = 0.0; if (k <= 1.35e-129) tmp = 2.0 / ((t_2 / cos(k)) * 0.0); elseif (k <= 2.05e-5) tmp = (l_m * l_m) * (2.0 / ((t_2 * (k * k)) / cos(k))); else tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * (0.5 - (cos((2.0 * k)) / 2.0)))) / cos(k))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := Block[{t$95$2 = N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k, 1.35e-129], N[(2.0 / N[(N[(t$95$2 / N[Cos[k], $MachinePrecision]), $MachinePrecision] * 0.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 2.05e-5], N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(2.0 / N[(N[(t$95$2 * N[(k * k), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(2.0 / N[(N[(N[(k * k), $MachinePrecision] * N[(t$95$m * N[(0.5 - N[(N[Cos[N[(2.0 * k), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := t\_m \cdot {k}^{2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 1.35 \cdot 10^{-129}:\\
\;\;\;\;\frac{2}{\frac{t\_2}{\cos k} \cdot 0}\\
\mathbf{elif}\;k \leq 2.05 \cdot 10^{-5}:\\
\;\;\;\;\left(l\_m \cdot l\_m\right) \cdot \frac{2}{\frac{t\_2 \cdot \left(k \cdot k\right)}{\cos k}}\\
\mathbf{else}:\\
\;\;\;\;\left(l\_m \cdot l\_m\right) \cdot \frac{2}{\frac{\left(k \cdot k\right) \cdot \left(t\_m \cdot \left(0.5 - \frac{\cos \left(2 \cdot k\right)}{2}\right)\right)}{\cos k}}\\
\end{array}
\end{array}
\end{array}
if k < 1.35e-129Initial program 29.9%
Simplified29.9%
Taylor expanded in t around 0 75.8%
associate-*r*75.8%
*-commutative75.8%
times-frac78.0%
Simplified78.0%
pow278.0%
add-log-exp63.7%
div-inv63.6%
exp-prod68.9%
pow268.9%
pow-flip68.9%
metadata-eval68.9%
Applied egg-rr68.9%
Taylor expanded in k around 0 49.7%
if 1.35e-129 < k < 2.05000000000000002e-5Initial program 21.8%
Simplified28.4%
Taylor expanded in t around 0 64.5%
unpow278.3%
Applied egg-rr64.5%
Taylor expanded in k around 0 64.4%
if 2.05000000000000002e-5 < k Initial program 29.0%
Simplified40.7%
Taylor expanded in t around 0 75.5%
unpow275.5%
sin-mult75.1%
Applied egg-rr75.1%
div-sub75.1%
+-inverses75.1%
cos-075.1%
metadata-eval75.1%
count-275.1%
Simplified75.1%
unpow275.6%
Applied egg-rr75.1%
Final simplification56.5%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(*
t_s
(if (<= k 2.05e-5)
(* (* l_m (/ l_m t_m)) (/ 2.0 (pow k 4.0)))
(*
(* l_m l_m)
(/
2.0
(/ (* (* k k) (* t_m (- 0.5 (/ (cos (* 2.0 k)) 2.0)))) (cos k)))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 2.05e-5) {
tmp = (l_m * (l_m / t_m)) * (2.0 / pow(k, 4.0));
} else {
tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * (0.5 - (cos((2.0 * k)) / 2.0)))) / cos(k)));
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 2.05d-5) then
tmp = (l_m * (l_m / t_m)) * (2.0d0 / (k ** 4.0d0))
else
tmp = (l_m * l_m) * (2.0d0 / (((k * k) * (t_m * (0.5d0 - (cos((2.0d0 * k)) / 2.0d0)))) / cos(k)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if (k <= 2.05e-5) {
tmp = (l_m * (l_m / t_m)) * (2.0 / Math.pow(k, 4.0));
} else {
tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * (0.5 - (Math.cos((2.0 * k)) / 2.0)))) / Math.cos(k)));
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): tmp = 0 if k <= 2.05e-5: tmp = (l_m * (l_m / t_m)) * (2.0 / math.pow(k, 4.0)) else: tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * (0.5 - (math.cos((2.0 * k)) / 2.0)))) / math.cos(k))) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) tmp = 0.0 if (k <= 2.05e-5) tmp = Float64(Float64(l_m * Float64(l_m / t_m)) * Float64(2.0 / (k ^ 4.0))); else tmp = Float64(Float64(l_m * l_m) * Float64(2.0 / Float64(Float64(Float64(k * k) * Float64(t_m * Float64(0.5 - Float64(cos(Float64(2.0 * k)) / 2.0)))) / cos(k)))); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) tmp = 0.0; if (k <= 2.05e-5) tmp = (l_m * (l_m / t_m)) * (2.0 / (k ^ 4.0)); else tmp = (l_m * l_m) * (2.0 / (((k * k) * (t_m * (0.5 - (cos((2.0 * k)) / 2.0)))) / cos(k))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := N[(t$95$s * If[LessEqual[k, 2.05e-5], N[(N[(l$95$m * N[(l$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] * N[(2.0 / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(2.0 / N[(N[(N[(k * k), $MachinePrecision] * N[(t$95$m * N[(0.5 - N[(N[Cos[N[(2.0 * k), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 2.05 \cdot 10^{-5}:\\
\;\;\;\;\left(l\_m \cdot \frac{l\_m}{t\_m}\right) \cdot \frac{2}{{k}^{4}}\\
\mathbf{else}:\\
\;\;\;\;\left(l\_m \cdot l\_m\right) \cdot \frac{2}{\frac{\left(k \cdot k\right) \cdot \left(t\_m \cdot \left(0.5 - \frac{\cos \left(2 \cdot k\right)}{2}\right)\right)}{\cos k}}\\
\end{array}
\end{array}
if k < 2.05000000000000002e-5Initial program 28.8%
Simplified36.8%
Taylor expanded in k around 0 64.2%
associate-*r/64.2%
*-commutative64.2%
*-commutative64.2%
times-frac63.4%
Simplified63.4%
pow263.4%
Applied egg-rr63.4%
associate-/l*68.1%
Applied egg-rr68.1%
if 2.05000000000000002e-5 < k Initial program 29.0%
Simplified40.7%
Taylor expanded in t around 0 75.5%
unpow275.5%
sin-mult75.1%
Applied egg-rr75.1%
div-sub75.1%
+-inverses75.1%
cos-075.1%
metadata-eval75.1%
count-275.1%
Simplified75.1%
unpow275.6%
Applied egg-rr75.1%
Final simplification69.6%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k)
:precision binary64
(*
t_s
(if (<= (* l_m l_m) 0.0)
(* (* l_m (/ l_m t_m)) (/ 2.0 (pow k 4.0)))
(* (* l_m l_m) (/ 2.0 (/ (* (* t_m (pow k 2.0)) (* k k)) (cos k)))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if ((l_m * l_m) <= 0.0) {
tmp = (l_m * (l_m / t_m)) * (2.0 / pow(k, 4.0));
} else {
tmp = (l_m * l_m) * (2.0 / (((t_m * pow(k, 2.0)) * (k * k)) / cos(k)));
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
real(8) :: tmp
if ((l_m * l_m) <= 0.0d0) then
tmp = (l_m * (l_m / t_m)) * (2.0d0 / (k ** 4.0d0))
else
tmp = (l_m * l_m) * (2.0d0 / (((t_m * (k ** 2.0d0)) * (k * k)) / cos(k)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
double tmp;
if ((l_m * l_m) <= 0.0) {
tmp = (l_m * (l_m / t_m)) * (2.0 / Math.pow(k, 4.0));
} else {
tmp = (l_m * l_m) * (2.0 / (((t_m * Math.pow(k, 2.0)) * (k * k)) / Math.cos(k)));
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): tmp = 0 if (l_m * l_m) <= 0.0: tmp = (l_m * (l_m / t_m)) * (2.0 / math.pow(k, 4.0)) else: tmp = (l_m * l_m) * (2.0 / (((t_m * math.pow(k, 2.0)) * (k * k)) / math.cos(k))) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) tmp = 0.0 if (Float64(l_m * l_m) <= 0.0) tmp = Float64(Float64(l_m * Float64(l_m / t_m)) * Float64(2.0 / (k ^ 4.0))); else tmp = Float64(Float64(l_m * l_m) * Float64(2.0 / Float64(Float64(Float64(t_m * (k ^ 2.0)) * Float64(k * k)) / cos(k)))); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k) tmp = 0.0; if ((l_m * l_m) <= 0.0) tmp = (l_m * (l_m / t_m)) * (2.0 / (k ^ 4.0)); else tmp = (l_m * l_m) * (2.0 / (((t_m * (k ^ 2.0)) * (k * k)) / cos(k))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := N[(t$95$s * If[LessEqual[N[(l$95$m * l$95$m), $MachinePrecision], 0.0], N[(N[(l$95$m * N[(l$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] * N[(2.0 / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(2.0 / N[(N[(N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;l\_m \cdot l\_m \leq 0:\\
\;\;\;\;\left(l\_m \cdot \frac{l\_m}{t\_m}\right) \cdot \frac{2}{{k}^{4}}\\
\mathbf{else}:\\
\;\;\;\;\left(l\_m \cdot l\_m\right) \cdot \frac{2}{\frac{\left(t\_m \cdot {k}^{2}\right) \cdot \left(k \cdot k\right)}{\cos k}}\\
\end{array}
\end{array}
if (*.f64 l l) < 0.0Initial program 25.9%
Simplified33.5%
Taylor expanded in k around 0 54.8%
associate-*r/54.8%
*-commutative54.8%
*-commutative54.8%
times-frac54.9%
Simplified54.9%
pow254.9%
Applied egg-rr54.9%
associate-/l*66.6%
Applied egg-rr66.6%
if 0.0 < (*.f64 l l) Initial program 29.9%
Simplified39.0%
Taylor expanded in t around 0 80.7%
unpow284.6%
Applied egg-rr80.7%
Taylor expanded in k around 0 65.9%
Final simplification66.1%
l_m = (fabs.f64 l) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l_m k) :precision binary64 (* t_s (* (* l_m (/ l_m t_m)) (/ 2.0 (pow k 4.0)))))
l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k) {
return t_s * ((l_m * (l_m / t_m)) * (2.0 / pow(k, 4.0)));
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k
code = t_s * ((l_m * (l_m / t_m)) * (2.0d0 / (k ** 4.0d0)))
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k) {
return t_s * ((l_m * (l_m / t_m)) * (2.0 / Math.pow(k, 4.0)));
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k): return t_s * ((l_m * (l_m / t_m)) * (2.0 / math.pow(k, 4.0)))
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k) return Float64(t_s * Float64(Float64(l_m * Float64(l_m / t_m)) * Float64(2.0 / (k ^ 4.0)))) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l_m, k) tmp = t_s * ((l_m * (l_m / t_m)) * (2.0 / (k ^ 4.0))); end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k_] := N[(t$95$s * N[(N[(l$95$m * N[(l$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] * N[(2.0 / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(l\_m \cdot \frac{l\_m}{t\_m}\right) \cdot \frac{2}{{k}^{4}}\right)
\end{array}
Initial program 28.9%
Simplified37.6%
Taylor expanded in k around 0 60.6%
associate-*r/60.6%
*-commutative60.6%
*-commutative60.6%
times-frac59.9%
Simplified59.9%
pow259.9%
Applied egg-rr59.9%
associate-/l*63.8%
Applied egg-rr63.8%
herbie shell --seed 2024137
(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))))