
(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 19 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)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 1.9e-14)
(pow
(/ (sqrt (/ 2.0 t_m)) (* k_m (/ (sin k_m) (* l (sqrt (cos k_m))))))
2.0)
(if (<= k_m 2.5e+152)
(/
2.0
(*
t_m
(* (/ (pow (sin k_m) 2.0) (* (cos k_m) (pow l 2.0))) (pow k_m 2.0))))
(/
2.0
(*
(pow (* (/ k_m t_m) (/ (pow t_m 1.5) l)) 2.0)
(* (sin k_m) (tan k_m))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1.9e-14) {
tmp = pow((sqrt((2.0 / t_m)) / (k_m * (sin(k_m) / (l * sqrt(cos(k_m)))))), 2.0);
} else if (k_m <= 2.5e+152) {
tmp = 2.0 / (t_m * ((pow(sin(k_m), 2.0) / (cos(k_m) * pow(l, 2.0))) * pow(k_m, 2.0)));
} else {
tmp = 2.0 / (pow(((k_m / t_m) * (pow(t_m, 1.5) / l)), 2.0) * (sin(k_m) * tan(k_m)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 1.9d-14) then
tmp = (sqrt((2.0d0 / t_m)) / (k_m * (sin(k_m) / (l * sqrt(cos(k_m)))))) ** 2.0d0
else if (k_m <= 2.5d+152) then
tmp = 2.0d0 / (t_m * (((sin(k_m) ** 2.0d0) / (cos(k_m) * (l ** 2.0d0))) * (k_m ** 2.0d0)))
else
tmp = 2.0d0 / ((((k_m / t_m) * ((t_m ** 1.5d0) / l)) ** 2.0d0) * (sin(k_m) * tan(k_m)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1.9e-14) {
tmp = Math.pow((Math.sqrt((2.0 / t_m)) / (k_m * (Math.sin(k_m) / (l * Math.sqrt(Math.cos(k_m)))))), 2.0);
} else if (k_m <= 2.5e+152) {
tmp = 2.0 / (t_m * ((Math.pow(Math.sin(k_m), 2.0) / (Math.cos(k_m) * Math.pow(l, 2.0))) * Math.pow(k_m, 2.0)));
} else {
tmp = 2.0 / (Math.pow(((k_m / t_m) * (Math.pow(t_m, 1.5) / l)), 2.0) * (Math.sin(k_m) * Math.tan(k_m)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 1.9e-14: tmp = math.pow((math.sqrt((2.0 / t_m)) / (k_m * (math.sin(k_m) / (l * math.sqrt(math.cos(k_m)))))), 2.0) elif k_m <= 2.5e+152: tmp = 2.0 / (t_m * ((math.pow(math.sin(k_m), 2.0) / (math.cos(k_m) * math.pow(l, 2.0))) * math.pow(k_m, 2.0))) else: tmp = 2.0 / (math.pow(((k_m / t_m) * (math.pow(t_m, 1.5) / l)), 2.0) * (math.sin(k_m) * math.tan(k_m))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 1.9e-14) tmp = Float64(sqrt(Float64(2.0 / t_m)) / Float64(k_m * Float64(sin(k_m) / Float64(l * sqrt(cos(k_m)))))) ^ 2.0; elseif (k_m <= 2.5e+152) tmp = Float64(2.0 / Float64(t_m * Float64(Float64((sin(k_m) ^ 2.0) / Float64(cos(k_m) * (l ^ 2.0))) * (k_m ^ 2.0)))); else tmp = Float64(2.0 / Float64((Float64(Float64(k_m / t_m) * Float64((t_m ^ 1.5) / l)) ^ 2.0) * Float64(sin(k_m) * tan(k_m)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 1.9e-14) tmp = (sqrt((2.0 / t_m)) / (k_m * (sin(k_m) / (l * sqrt(cos(k_m)))))) ^ 2.0; elseif (k_m <= 2.5e+152) tmp = 2.0 / (t_m * (((sin(k_m) ^ 2.0) / (cos(k_m) * (l ^ 2.0))) * (k_m ^ 2.0))); else tmp = 2.0 / ((((k_m / t_m) * ((t_m ^ 1.5) / l)) ^ 2.0) * (sin(k_m) * tan(k_m))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 1.9e-14], N[Power[N[(N[Sqrt[N[(2.0 / t$95$m), $MachinePrecision]], $MachinePrecision] / N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / N[(l * N[Sqrt[N[Cos[k$95$m], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 2.5e+152], N[(2.0 / N[(t$95$m * N[(N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Power[N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 1.9 \cdot 10^{-14}:\\
\;\;\;\;{\left(\frac{\sqrt{\frac{2}{t\_m}}}{k\_m \cdot \frac{\sin k\_m}{\ell \cdot \sqrt{\cos k\_m}}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 2.5 \cdot 10^{+152}:\\
\;\;\;\;\frac{2}{t\_m \cdot \left(\frac{{\sin k\_m}^{2}}{\cos k\_m \cdot {\ell}^{2}} \cdot {k\_m}^{2}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{t\_m} \cdot \frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(\sin k\_m \cdot \tan k\_m\right)}\\
\end{array}
\end{array}
if k < 1.9000000000000001e-14Initial program 41.4%
associate-/r*46.8%
div-inv46.8%
Applied egg-rr46.8%
add-exp-log46.5%
expm1-define46.5%
log1p-define51.3%
unpow251.3%
expm1-log1p-u51.6%
clear-num51.6%
un-div-inv51.6%
Applied egg-rr51.6%
Taylor expanded in t around 0 76.7%
associate-/l*77.4%
*-commutative77.4%
associate-*r/78.0%
*-commutative78.0%
associate-*l*77.3%
*-commutative77.3%
Simplified77.3%
add-sqr-sqrt48.7%
pow248.7%
Applied egg-rr38.5%
if 1.9000000000000001e-14 < k < 2.5e152Initial program 16.2%
associate-/r*22.2%
div-inv22.2%
Applied egg-rr22.2%
add-exp-log21.6%
expm1-define21.6%
log1p-define31.8%
unpow231.8%
expm1-log1p-u32.5%
clear-num32.5%
un-div-inv32.5%
Applied egg-rr32.5%
Taylor expanded in t around 0 78.3%
associate-/l*78.5%
*-commutative78.5%
associate-*r/75.9%
*-commutative75.9%
associate-*l*81.9%
*-commutative81.9%
Simplified81.9%
if 2.5e152 < k Initial program 40.5%
add-sqr-sqrt17.3%
pow217.3%
Applied egg-rr14.2%
associate-*r*14.2%
unpow-prod-down14.2%
pow214.2%
add-sqr-sqrt34.3%
Applied egg-rr34.3%
Final simplification43.3%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 1.9e-14)
(pow
(/ (sqrt (/ 2.0 t_m)) (* k_m (/ (sin k_m) (* l (sqrt (cos k_m))))))
2.0)
(if (<= k_m 4.05e+152)
(/ 2.0 (* t_m (/ (pow (* k_m (sin k_m)) 2.0) (* (cos k_m) (pow l 2.0)))))
(/
2.0
(*
(pow (* (/ k_m t_m) (/ (pow t_m 1.5) l)) 2.0)
(* (sin k_m) (tan k_m))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1.9e-14) {
tmp = pow((sqrt((2.0 / t_m)) / (k_m * (sin(k_m) / (l * sqrt(cos(k_m)))))), 2.0);
} else if (k_m <= 4.05e+152) {
tmp = 2.0 / (t_m * (pow((k_m * sin(k_m)), 2.0) / (cos(k_m) * pow(l, 2.0))));
} else {
tmp = 2.0 / (pow(((k_m / t_m) * (pow(t_m, 1.5) / l)), 2.0) * (sin(k_m) * tan(k_m)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 1.9d-14) then
tmp = (sqrt((2.0d0 / t_m)) / (k_m * (sin(k_m) / (l * sqrt(cos(k_m)))))) ** 2.0d0
else if (k_m <= 4.05d+152) then
tmp = 2.0d0 / (t_m * (((k_m * sin(k_m)) ** 2.0d0) / (cos(k_m) * (l ** 2.0d0))))
else
tmp = 2.0d0 / ((((k_m / t_m) * ((t_m ** 1.5d0) / l)) ** 2.0d0) * (sin(k_m) * tan(k_m)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1.9e-14) {
tmp = Math.pow((Math.sqrt((2.0 / t_m)) / (k_m * (Math.sin(k_m) / (l * Math.sqrt(Math.cos(k_m)))))), 2.0);
} else if (k_m <= 4.05e+152) {
tmp = 2.0 / (t_m * (Math.pow((k_m * Math.sin(k_m)), 2.0) / (Math.cos(k_m) * Math.pow(l, 2.0))));
} else {
tmp = 2.0 / (Math.pow(((k_m / t_m) * (Math.pow(t_m, 1.5) / l)), 2.0) * (Math.sin(k_m) * Math.tan(k_m)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 1.9e-14: tmp = math.pow((math.sqrt((2.0 / t_m)) / (k_m * (math.sin(k_m) / (l * math.sqrt(math.cos(k_m)))))), 2.0) elif k_m <= 4.05e+152: tmp = 2.0 / (t_m * (math.pow((k_m * math.sin(k_m)), 2.0) / (math.cos(k_m) * math.pow(l, 2.0)))) else: tmp = 2.0 / (math.pow(((k_m / t_m) * (math.pow(t_m, 1.5) / l)), 2.0) * (math.sin(k_m) * math.tan(k_m))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 1.9e-14) tmp = Float64(sqrt(Float64(2.0 / t_m)) / Float64(k_m * Float64(sin(k_m) / Float64(l * sqrt(cos(k_m)))))) ^ 2.0; elseif (k_m <= 4.05e+152) tmp = Float64(2.0 / Float64(t_m * Float64((Float64(k_m * sin(k_m)) ^ 2.0) / Float64(cos(k_m) * (l ^ 2.0))))); else tmp = Float64(2.0 / Float64((Float64(Float64(k_m / t_m) * Float64((t_m ^ 1.5) / l)) ^ 2.0) * Float64(sin(k_m) * tan(k_m)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 1.9e-14) tmp = (sqrt((2.0 / t_m)) / (k_m * (sin(k_m) / (l * sqrt(cos(k_m)))))) ^ 2.0; elseif (k_m <= 4.05e+152) tmp = 2.0 / (t_m * (((k_m * sin(k_m)) ^ 2.0) / (cos(k_m) * (l ^ 2.0)))); else tmp = 2.0 / ((((k_m / t_m) * ((t_m ^ 1.5) / l)) ^ 2.0) * (sin(k_m) * tan(k_m))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 1.9e-14], N[Power[N[(N[Sqrt[N[(2.0 / t$95$m), $MachinePrecision]], $MachinePrecision] / N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / N[(l * N[Sqrt[N[Cos[k$95$m], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 4.05e+152], N[(2.0 / N[(t$95$m * N[(N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Power[N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 1.9 \cdot 10^{-14}:\\
\;\;\;\;{\left(\frac{\sqrt{\frac{2}{t\_m}}}{k\_m \cdot \frac{\sin k\_m}{\ell \cdot \sqrt{\cos k\_m}}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 4.05 \cdot 10^{+152}:\\
\;\;\;\;\frac{2}{t\_m \cdot \frac{{\left(k\_m \cdot \sin k\_m\right)}^{2}}{\cos k\_m \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{t\_m} \cdot \frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(\sin k\_m \cdot \tan k\_m\right)}\\
\end{array}
\end{array}
if k < 1.9000000000000001e-14Initial program 41.4%
associate-/r*46.8%
div-inv46.8%
Applied egg-rr46.8%
add-exp-log46.5%
expm1-define46.5%
log1p-define51.3%
unpow251.3%
expm1-log1p-u51.6%
clear-num51.6%
un-div-inv51.6%
Applied egg-rr51.6%
Taylor expanded in t around 0 76.7%
associate-/l*77.4%
*-commutative77.4%
associate-*r/78.0%
*-commutative78.0%
associate-*l*77.3%
*-commutative77.3%
Simplified77.3%
add-sqr-sqrt48.7%
pow248.7%
Applied egg-rr38.5%
if 1.9000000000000001e-14 < k < 4.04999999999999999e152Initial program 16.2%
associate-/r*22.2%
div-inv22.2%
Applied egg-rr22.2%
add-exp-log21.6%
expm1-define21.6%
log1p-define31.8%
unpow231.8%
expm1-log1p-u32.5%
clear-num32.5%
un-div-inv32.5%
Applied egg-rr32.5%
Taylor expanded in t around 0 78.3%
associate-/l*78.5%
*-commutative78.5%
associate-*r/75.9%
*-commutative75.9%
associate-*l*81.9%
*-commutative81.9%
Simplified81.9%
associate-*l/81.9%
pow-prod-down81.8%
Applied egg-rr81.8%
if 4.04999999999999999e152 < k Initial program 40.5%
add-sqr-sqrt17.3%
pow217.3%
Applied egg-rr14.2%
associate-*r*14.2%
unpow-prod-down14.2%
pow214.2%
add-sqr-sqrt34.3%
Applied egg-rr34.3%
Final simplification43.3%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 2e-219)
(/ 2.0 (pow (* (/ (pow k_m 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+300)
(/
2.0
(*
(* k_m k_m)
(/ (* t_m (pow (sin k_m) 2.0)) (* (cos k_m) (pow l 2.0)))))
(/
2.0
(pow (* (sqrt (/ t_m (cos k_m))) (/ (* k_m (sin k_m)) l)) 2.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 2e-219) {
tmp = 2.0 / pow(((pow(k_m, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+300) {
tmp = 2.0 / ((k_m * k_m) * ((t_m * pow(sin(k_m), 2.0)) / (cos(k_m) * pow(l, 2.0))));
} else {
tmp = 2.0 / pow((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 2d-219) then
tmp = 2.0d0 / ((((k_m ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+300) then
tmp = 2.0d0 / ((k_m * k_m) * ((t_m * (sin(k_m) ** 2.0d0)) / (cos(k_m) * (l ** 2.0d0))))
else
tmp = 2.0d0 / ((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)) ** 2.0d0)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 2e-219) {
tmp = 2.0 / Math.pow(((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+300) {
tmp = 2.0 / ((k_m * k_m) * ((t_m * Math.pow(Math.sin(k_m), 2.0)) / (Math.cos(k_m) * Math.pow(l, 2.0))));
} else {
tmp = 2.0 / Math.pow((Math.sqrt((t_m / Math.cos(k_m))) * ((k_m * Math.sin(k_m)) / l)), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 2e-219: tmp = 2.0 / math.pow(((math.pow(k_m, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+300: tmp = 2.0 / ((k_m * k_m) * ((t_m * math.pow(math.sin(k_m), 2.0)) / (math.cos(k_m) * math.pow(l, 2.0)))) else: tmp = 2.0 / math.pow((math.sqrt((t_m / math.cos(k_m))) * ((k_m * math.sin(k_m)) / l)), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 2e-219) tmp = Float64(2.0 / (Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+300) tmp = Float64(2.0 / Float64(Float64(k_m * k_m) * Float64(Float64(t_m * (sin(k_m) ^ 2.0)) / Float64(cos(k_m) * (l ^ 2.0))))); else tmp = Float64(2.0 / (Float64(sqrt(Float64(t_m / cos(k_m))) * Float64(Float64(k_m * sin(k_m)) / l)) ^ 2.0)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 2e-219) tmp = 2.0 / ((((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+300) tmp = 2.0 / ((k_m * k_m) * ((t_m * (sin(k_m) ^ 2.0)) / (cos(k_m) * (l ^ 2.0)))); else tmp = 2.0 / ((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)) ^ 2.0); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 2e-219], N[(2.0 / N[Power[N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+300], N[(2.0 / N[(N[(k$95$m * k$95$m), $MachinePrecision] * N[(N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 2 \cdot 10^{-219}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+300}:\\
\;\;\;\;\frac{2}{\left(k\_m \cdot k\_m\right) \cdot \frac{t\_m \cdot {\sin k\_m}^{2}}{\cos k\_m \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\sqrt{\frac{t\_m}{\cos k\_m}} \cdot \frac{k\_m \cdot \sin k\_m}{\ell}\right)}^{2}}\\
\end{array}
\end{array}
if (*.f64 l l) < 2.0000000000000001e-219Initial program 28.8%
add-sqr-sqrt11.0%
pow211.0%
Applied egg-rr13.9%
Taylor expanded in k around 0 38.2%
if 2.0000000000000001e-219 < (*.f64 l l) < 2.0000000000000001e300Initial program 42.8%
Taylor expanded in t around 0 87.5%
associate-/l*88.6%
*-commutative88.6%
Simplified88.6%
unpow288.6%
Applied egg-rr88.6%
if 2.0000000000000001e300 < (*.f64 l l) Initial program 40.2%
add-sqr-sqrt18.7%
pow218.7%
Applied egg-rr28.4%
Taylor expanded in k around inf 56.8%
Final simplification65.5%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 2e-219)
(/ 2.0 (pow (* (/ (pow k_m 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+300)
(*
(* l l)
(/ 2.0 (/ (* (pow (sin k_m) 2.0) (* t_m (pow k_m 2.0))) (cos k_m))))
(/
2.0
(pow (* (sqrt (/ t_m (cos k_m))) (/ (* k_m (sin k_m)) l)) 2.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 2e-219) {
tmp = 2.0 / pow(((pow(k_m, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+300) {
tmp = (l * l) * (2.0 / ((pow(sin(k_m), 2.0) * (t_m * pow(k_m, 2.0))) / cos(k_m)));
} else {
tmp = 2.0 / pow((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 2d-219) then
tmp = 2.0d0 / ((((k_m ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+300) then
tmp = (l * l) * (2.0d0 / (((sin(k_m) ** 2.0d0) * (t_m * (k_m ** 2.0d0))) / cos(k_m)))
else
tmp = 2.0d0 / ((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)) ** 2.0d0)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 2e-219) {
tmp = 2.0 / Math.pow(((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+300) {
tmp = (l * l) * (2.0 / ((Math.pow(Math.sin(k_m), 2.0) * (t_m * Math.pow(k_m, 2.0))) / Math.cos(k_m)));
} else {
tmp = 2.0 / Math.pow((Math.sqrt((t_m / Math.cos(k_m))) * ((k_m * Math.sin(k_m)) / l)), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 2e-219: tmp = 2.0 / math.pow(((math.pow(k_m, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+300: tmp = (l * l) * (2.0 / ((math.pow(math.sin(k_m), 2.0) * (t_m * math.pow(k_m, 2.0))) / math.cos(k_m))) else: tmp = 2.0 / math.pow((math.sqrt((t_m / math.cos(k_m))) * ((k_m * math.sin(k_m)) / l)), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 2e-219) tmp = Float64(2.0 / (Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+300) tmp = Float64(Float64(l * l) * Float64(2.0 / Float64(Float64((sin(k_m) ^ 2.0) * Float64(t_m * (k_m ^ 2.0))) / cos(k_m)))); else tmp = Float64(2.0 / (Float64(sqrt(Float64(t_m / cos(k_m))) * Float64(Float64(k_m * sin(k_m)) / l)) ^ 2.0)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 2e-219) tmp = 2.0 / ((((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+300) tmp = (l * l) * (2.0 / (((sin(k_m) ^ 2.0) * (t_m * (k_m ^ 2.0))) / cos(k_m))); else tmp = 2.0 / ((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)) ^ 2.0); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 2e-219], N[(2.0 / N[Power[N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+300], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 2 \cdot 10^{-219}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+300}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{\frac{{\sin k\_m}^{2} \cdot \left(t\_m \cdot {k\_m}^{2}\right)}{\cos k\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\sqrt{\frac{t\_m}{\cos k\_m}} \cdot \frac{k\_m \cdot \sin k\_m}{\ell}\right)}^{2}}\\
\end{array}
\end{array}
if (*.f64 l l) < 2.0000000000000001e-219Initial program 28.8%
add-sqr-sqrt11.0%
pow211.0%
Applied egg-rr13.9%
Taylor expanded in k around 0 38.2%
if 2.0000000000000001e-219 < (*.f64 l l) < 2.0000000000000001e300Initial program 42.8%
Simplified47.5%
Taylor expanded in t around 0 87.5%
associate-*r*87.5%
Simplified87.5%
if 2.0000000000000001e300 < (*.f64 l l) Initial program 40.2%
add-sqr-sqrt18.7%
pow218.7%
Applied egg-rr28.4%
Taylor expanded in k around inf 56.8%
Final simplification65.1%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 2e-219)
(/ 2.0 (pow (* (/ (pow k_m 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+300)
(*
(* l l)
(/ 2.0 (* (pow k_m 2.0) (/ (* t_m (pow (sin k_m) 2.0)) (cos k_m)))))
(/
2.0
(pow (* (sqrt (/ t_m (cos k_m))) (/ (* k_m (sin k_m)) l)) 2.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 2e-219) {
tmp = 2.0 / pow(((pow(k_m, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+300) {
tmp = (l * l) * (2.0 / (pow(k_m, 2.0) * ((t_m * pow(sin(k_m), 2.0)) / cos(k_m))));
} else {
tmp = 2.0 / pow((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 2d-219) then
tmp = 2.0d0 / ((((k_m ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+300) then
tmp = (l * l) * (2.0d0 / ((k_m ** 2.0d0) * ((t_m * (sin(k_m) ** 2.0d0)) / cos(k_m))))
else
tmp = 2.0d0 / ((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)) ** 2.0d0)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 2e-219) {
tmp = 2.0 / Math.pow(((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+300) {
tmp = (l * l) * (2.0 / (Math.pow(k_m, 2.0) * ((t_m * Math.pow(Math.sin(k_m), 2.0)) / Math.cos(k_m))));
} else {
tmp = 2.0 / Math.pow((Math.sqrt((t_m / Math.cos(k_m))) * ((k_m * Math.sin(k_m)) / l)), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 2e-219: tmp = 2.0 / math.pow(((math.pow(k_m, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+300: tmp = (l * l) * (2.0 / (math.pow(k_m, 2.0) * ((t_m * math.pow(math.sin(k_m), 2.0)) / math.cos(k_m)))) else: tmp = 2.0 / math.pow((math.sqrt((t_m / math.cos(k_m))) * ((k_m * math.sin(k_m)) / l)), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 2e-219) tmp = Float64(2.0 / (Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+300) tmp = Float64(Float64(l * l) * Float64(2.0 / Float64((k_m ^ 2.0) * Float64(Float64(t_m * (sin(k_m) ^ 2.0)) / cos(k_m))))); else tmp = Float64(2.0 / (Float64(sqrt(Float64(t_m / cos(k_m))) * Float64(Float64(k_m * sin(k_m)) / l)) ^ 2.0)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 2e-219) tmp = 2.0 / ((((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+300) tmp = (l * l) * (2.0 / ((k_m ^ 2.0) * ((t_m * (sin(k_m) ^ 2.0)) / cos(k_m)))); else tmp = 2.0 / ((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)) ^ 2.0); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 2e-219], N[(2.0 / N[Power[N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+300], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 2 \cdot 10^{-219}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+300}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{{k\_m}^{2} \cdot \frac{t\_m \cdot {\sin k\_m}^{2}}{\cos k\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\sqrt{\frac{t\_m}{\cos k\_m}} \cdot \frac{k\_m \cdot \sin k\_m}{\ell}\right)}^{2}}\\
\end{array}
\end{array}
if (*.f64 l l) < 2.0000000000000001e-219Initial program 28.8%
add-sqr-sqrt11.0%
pow211.0%
Applied egg-rr13.9%
Taylor expanded in k around 0 38.2%
if 2.0000000000000001e-219 < (*.f64 l l) < 2.0000000000000001e300Initial program 42.8%
Simplified47.5%
Taylor expanded in t around 0 87.5%
associate-/l*87.5%
Simplified87.5%
if 2.0000000000000001e300 < (*.f64 l l) Initial program 40.2%
add-sqr-sqrt18.7%
pow218.7%
Applied egg-rr28.4%
Taylor expanded in k around inf 56.8%
Final simplification65.1%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (* k_m (sin k_m))))
(*
t_s
(if (<= (* l l) 5e-324)
(/ 2.0 (pow (* (/ (pow k_m 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+300)
(/ 2.0 (* t_m (/ (pow t_2 2.0) (* (cos k_m) (pow l 2.0)))))
(/ 2.0 (pow (* (sqrt (/ t_m (cos k_m))) (/ t_2 l)) 2.0)))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = k_m * sin(k_m);
double tmp;
if ((l * l) <= 5e-324) {
tmp = 2.0 / pow(((pow(k_m, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+300) {
tmp = 2.0 / (t_m * (pow(t_2, 2.0) / (cos(k_m) * pow(l, 2.0))));
} else {
tmp = 2.0 / pow((sqrt((t_m / cos(k_m))) * (t_2 / l)), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: t_2
real(8) :: tmp
t_2 = k_m * sin(k_m)
if ((l * l) <= 5d-324) then
tmp = 2.0d0 / ((((k_m ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+300) then
tmp = 2.0d0 / (t_m * ((t_2 ** 2.0d0) / (cos(k_m) * (l ** 2.0d0))))
else
tmp = 2.0d0 / ((sqrt((t_m / cos(k_m))) * (t_2 / l)) ** 2.0d0)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = k_m * Math.sin(k_m);
double tmp;
if ((l * l) <= 5e-324) {
tmp = 2.0 / Math.pow(((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+300) {
tmp = 2.0 / (t_m * (Math.pow(t_2, 2.0) / (Math.cos(k_m) * Math.pow(l, 2.0))));
} else {
tmp = 2.0 / Math.pow((Math.sqrt((t_m / Math.cos(k_m))) * (t_2 / l)), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): t_2 = k_m * math.sin(k_m) tmp = 0 if (l * l) <= 5e-324: tmp = 2.0 / math.pow(((math.pow(k_m, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+300: tmp = 2.0 / (t_m * (math.pow(t_2, 2.0) / (math.cos(k_m) * math.pow(l, 2.0)))) else: tmp = 2.0 / math.pow((math.sqrt((t_m / math.cos(k_m))) * (t_2 / l)), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = Float64(k_m * sin(k_m)) tmp = 0.0 if (Float64(l * l) <= 5e-324) tmp = Float64(2.0 / (Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+300) tmp = Float64(2.0 / Float64(t_m * Float64((t_2 ^ 2.0) / Float64(cos(k_m) * (l ^ 2.0))))); else tmp = Float64(2.0 / (Float64(sqrt(Float64(t_m / cos(k_m))) * Float64(t_2 / l)) ^ 2.0)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) t_2 = k_m * sin(k_m); tmp = 0.0; if ((l * l) <= 5e-324) tmp = 2.0 / ((((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+300) tmp = 2.0 / (t_m * ((t_2 ^ 2.0) / (cos(k_m) * (l ^ 2.0)))); else tmp = 2.0 / ((sqrt((t_m / cos(k_m))) * (t_2 / l)) ^ 2.0); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := Block[{t$95$2 = N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 5e-324], N[(2.0 / N[Power[N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+300], N[(2.0 / N[(t$95$m * N[(N[Power[t$95$2, 2.0], $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(t$95$2 / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := k\_m \cdot \sin k\_m\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 5 \cdot 10^{-324}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+300}:\\
\;\;\;\;\frac{2}{t\_m \cdot \frac{{t\_2}^{2}}{\cos k\_m \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\sqrt{\frac{t\_m}{\cos k\_m}} \cdot \frac{t\_2}{\ell}\right)}^{2}}\\
\end{array}
\end{array}
\end{array}
if (*.f64 l l) < 4.94066e-324Initial program 25.5%
add-sqr-sqrt9.1%
pow29.1%
Applied egg-rr13.0%
Taylor expanded in k around 0 36.3%
if 4.94066e-324 < (*.f64 l l) < 2.0000000000000001e300Initial program 42.3%
associate-/r*42.3%
div-inv42.3%
Applied egg-rr42.3%
add-exp-log41.9%
expm1-define41.9%
log1p-define48.1%
unpow248.1%
expm1-log1p-u48.5%
clear-num48.5%
un-div-inv48.5%
Applied egg-rr48.5%
Taylor expanded in t around 0 88.4%
associate-/l*89.3%
*-commutative89.3%
associate-*r/89.4%
*-commutative89.4%
associate-*l*89.9%
*-commutative89.9%
Simplified89.9%
associate-*l/88.5%
pow-prod-down88.4%
Applied egg-rr88.4%
if 2.0000000000000001e300 < (*.f64 l l) Initial program 40.2%
add-sqr-sqrt18.7%
pow218.7%
Applied egg-rr28.4%
Taylor expanded in k around inf 56.8%
Final simplification68.6%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 1.9e-14)
(/ 2.0 (pow (* k_m (* (/ (sin k_m) l) (sqrt (/ t_m (cos k_m))))) 2.0))
(if (<= k_m 4.1e+148)
(/ 2.0 (* t_m (/ (pow (* k_m (sin k_m)) 2.0) (* (cos k_m) (pow l 2.0)))))
(/
2.0
(*
(pow (* (/ k_m t_m) (/ (pow t_m 1.5) l)) 2.0)
(* (sin k_m) (tan k_m))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1.9e-14) {
tmp = 2.0 / pow((k_m * ((sin(k_m) / l) * sqrt((t_m / cos(k_m))))), 2.0);
} else if (k_m <= 4.1e+148) {
tmp = 2.0 / (t_m * (pow((k_m * sin(k_m)), 2.0) / (cos(k_m) * pow(l, 2.0))));
} else {
tmp = 2.0 / (pow(((k_m / t_m) * (pow(t_m, 1.5) / l)), 2.0) * (sin(k_m) * tan(k_m)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 1.9d-14) then
tmp = 2.0d0 / ((k_m * ((sin(k_m) / l) * sqrt((t_m / cos(k_m))))) ** 2.0d0)
else if (k_m <= 4.1d+148) then
tmp = 2.0d0 / (t_m * (((k_m * sin(k_m)) ** 2.0d0) / (cos(k_m) * (l ** 2.0d0))))
else
tmp = 2.0d0 / ((((k_m / t_m) * ((t_m ** 1.5d0) / l)) ** 2.0d0) * (sin(k_m) * tan(k_m)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1.9e-14) {
tmp = 2.0 / Math.pow((k_m * ((Math.sin(k_m) / l) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else if (k_m <= 4.1e+148) {
tmp = 2.0 / (t_m * (Math.pow((k_m * Math.sin(k_m)), 2.0) / (Math.cos(k_m) * Math.pow(l, 2.0))));
} else {
tmp = 2.0 / (Math.pow(((k_m / t_m) * (Math.pow(t_m, 1.5) / l)), 2.0) * (Math.sin(k_m) * Math.tan(k_m)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 1.9e-14: tmp = 2.0 / math.pow((k_m * ((math.sin(k_m) / l) * math.sqrt((t_m / math.cos(k_m))))), 2.0) elif k_m <= 4.1e+148: tmp = 2.0 / (t_m * (math.pow((k_m * math.sin(k_m)), 2.0) / (math.cos(k_m) * math.pow(l, 2.0)))) else: tmp = 2.0 / (math.pow(((k_m / t_m) * (math.pow(t_m, 1.5) / l)), 2.0) * (math.sin(k_m) * math.tan(k_m))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 1.9e-14) tmp = Float64(2.0 / (Float64(k_m * Float64(Float64(sin(k_m) / l) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); elseif (k_m <= 4.1e+148) tmp = Float64(2.0 / Float64(t_m * Float64((Float64(k_m * sin(k_m)) ^ 2.0) / Float64(cos(k_m) * (l ^ 2.0))))); else tmp = Float64(2.0 / Float64((Float64(Float64(k_m / t_m) * Float64((t_m ^ 1.5) / l)) ^ 2.0) * Float64(sin(k_m) * tan(k_m)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 1.9e-14) tmp = 2.0 / ((k_m * ((sin(k_m) / l) * sqrt((t_m / cos(k_m))))) ^ 2.0); elseif (k_m <= 4.1e+148) tmp = 2.0 / (t_m * (((k_m * sin(k_m)) ^ 2.0) / (cos(k_m) * (l ^ 2.0)))); else tmp = 2.0 / ((((k_m / t_m) * ((t_m ^ 1.5) / l)) ^ 2.0) * (sin(k_m) * tan(k_m))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 1.9e-14], N[(2.0 / N[Power[N[(k$95$m * N[(N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 4.1e+148], N[(2.0 / N[(t$95$m * N[(N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Power[N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 1.9 \cdot 10^{-14}:\\
\;\;\;\;\frac{2}{{\left(k\_m \cdot \left(\frac{\sin k\_m}{\ell} \cdot \sqrt{\frac{t\_m}{\cos k\_m}}\right)\right)}^{2}}\\
\mathbf{elif}\;k\_m \leq 4.1 \cdot 10^{+148}:\\
\;\;\;\;\frac{2}{t\_m \cdot \frac{{\left(k\_m \cdot \sin k\_m\right)}^{2}}{\cos k\_m \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{t\_m} \cdot \frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(\sin k\_m \cdot \tan k\_m\right)}\\
\end{array}
\end{array}
if k < 1.9000000000000001e-14Initial program 41.4%
add-sqr-sqrt20.7%
pow220.7%
Applied egg-rr25.5%
Taylor expanded in k around inf 49.5%
associate-/l*50.0%
associate-*l*49.2%
Simplified49.2%
if 1.9000000000000001e-14 < k < 4.0999999999999998e148Initial program 16.2%
associate-/r*22.2%
div-inv22.2%
Applied egg-rr22.2%
add-exp-log21.6%
expm1-define21.6%
log1p-define31.8%
unpow231.8%
expm1-log1p-u32.5%
clear-num32.5%
un-div-inv32.5%
Applied egg-rr32.5%
Taylor expanded in t around 0 78.3%
associate-/l*78.5%
*-commutative78.5%
associate-*r/75.9%
*-commutative75.9%
associate-*l*81.9%
*-commutative81.9%
Simplified81.9%
associate-*l/81.9%
pow-prod-down81.8%
Applied egg-rr81.8%
if 4.0999999999999998e148 < k Initial program 40.5%
add-sqr-sqrt17.3%
pow217.3%
Applied egg-rr14.2%
associate-*r*14.2%
unpow-prod-down14.2%
pow214.2%
add-sqr-sqrt34.3%
Applied egg-rr34.3%
Final simplification51.2%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 5e+211)
(/ 2.0 (pow (* (/ (pow k_m 2.0) l) (sqrt t_m)) 2.0))
(/ 2.0 (pow (* (sqrt (/ t_m (cos k_m))) (/ (* k_m (sin k_m)) l)) 2.0)))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 5e+211) {
tmp = 2.0 / pow(((pow(k_m, 2.0) / l) * sqrt(t_m)), 2.0);
} else {
tmp = 2.0 / pow((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 5d+211) then
tmp = 2.0d0 / ((((k_m ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else
tmp = 2.0d0 / ((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)) ** 2.0d0)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 5e+211) {
tmp = 2.0 / Math.pow(((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else {
tmp = 2.0 / Math.pow((Math.sqrt((t_m / Math.cos(k_m))) * ((k_m * Math.sin(k_m)) / l)), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 5e+211: tmp = 2.0 / math.pow(((math.pow(k_m, 2.0) / l) * math.sqrt(t_m)), 2.0) else: tmp = 2.0 / math.pow((math.sqrt((t_m / math.cos(k_m))) * ((k_m * math.sin(k_m)) / l)), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 5e+211) tmp = Float64(2.0 / (Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); else tmp = Float64(2.0 / (Float64(sqrt(Float64(t_m / cos(k_m))) * Float64(Float64(k_m * sin(k_m)) / l)) ^ 2.0)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 5e+211) tmp = 2.0 / ((((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); else tmp = 2.0 / ((sqrt((t_m / cos(k_m))) * ((k_m * sin(k_m)) / l)) ^ 2.0); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 5e+211], N[(2.0 / N[Power[N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 5 \cdot 10^{+211}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\sqrt{\frac{t\_m}{\cos k\_m}} \cdot \frac{k\_m \cdot \sin k\_m}{\ell}\right)}^{2}}\\
\end{array}
\end{array}
if (*.f64 l l) < 4.9999999999999995e211Initial program 36.7%
add-sqr-sqrt17.8%
pow217.8%
Applied egg-rr20.2%
Taylor expanded in k around 0 37.4%
if 4.9999999999999995e211 < (*.f64 l l) Initial program 40.5%
add-sqr-sqrt20.8%
pow220.8%
Applied egg-rr28.2%
Taylor expanded in k around inf 56.2%
Final simplification44.1%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 5e+211)
(/ 2.0 (pow (* (/ (pow k_m 2.0) l) (sqrt t_m)) 2.0))
(/ 2.0 (* t_m (pow (* k_m (/ (/ (sin k_m) l) (sqrt (cos k_m)))) 2.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 5e+211) {
tmp = 2.0 / pow(((pow(k_m, 2.0) / l) * sqrt(t_m)), 2.0);
} else {
tmp = 2.0 / (t_m * pow((k_m * ((sin(k_m) / l) / sqrt(cos(k_m)))), 2.0));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 5d+211) then
tmp = 2.0d0 / ((((k_m ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else
tmp = 2.0d0 / (t_m * ((k_m * ((sin(k_m) / l) / sqrt(cos(k_m)))) ** 2.0d0))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 5e+211) {
tmp = 2.0 / Math.pow(((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else {
tmp = 2.0 / (t_m * Math.pow((k_m * ((Math.sin(k_m) / l) / Math.sqrt(Math.cos(k_m)))), 2.0));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 5e+211: tmp = 2.0 / math.pow(((math.pow(k_m, 2.0) / l) * math.sqrt(t_m)), 2.0) else: tmp = 2.0 / (t_m * math.pow((k_m * ((math.sin(k_m) / l) / math.sqrt(math.cos(k_m)))), 2.0)) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 5e+211) tmp = Float64(2.0 / (Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); else tmp = Float64(2.0 / Float64(t_m * (Float64(k_m * Float64(Float64(sin(k_m) / l) / sqrt(cos(k_m)))) ^ 2.0))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 5e+211) tmp = 2.0 / ((((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); else tmp = 2.0 / (t_m * ((k_m * ((sin(k_m) / l) / sqrt(cos(k_m)))) ^ 2.0)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 5e+211], N[(2.0 / N[Power[N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(t$95$m * N[Power[N[(k$95$m * N[(N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision] / N[Sqrt[N[Cos[k$95$m], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 5 \cdot 10^{+211}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t\_m \cdot {\left(k\_m \cdot \frac{\frac{\sin k\_m}{\ell}}{\sqrt{\cos k\_m}}\right)}^{2}}\\
\end{array}
\end{array}
if (*.f64 l l) < 4.9999999999999995e211Initial program 36.7%
add-sqr-sqrt17.8%
pow217.8%
Applied egg-rr20.2%
Taylor expanded in k around 0 37.4%
if 4.9999999999999995e211 < (*.f64 l l) Initial program 40.5%
associate-/r*42.8%
div-inv42.9%
Applied egg-rr42.9%
add-exp-log42.6%
expm1-define42.6%
log1p-define42.8%
unpow242.8%
expm1-log1p-u43.1%
clear-num43.1%
un-div-inv43.1%
Applied egg-rr43.1%
Taylor expanded in t around 0 64.6%
associate-/l*64.1%
*-commutative64.1%
associate-*r/64.1%
*-commutative64.1%
associate-*l*65.8%
*-commutative65.8%
Simplified65.8%
pow165.8%
Applied egg-rr71.4%
unpow171.4%
associate-/r*71.5%
Simplified71.5%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= t_m 1.85e-118)
(pow (/ (/ (sqrt 2.0) k_m) (* k_m (/ (sqrt t_m) l))) 2.0)
(if (<= t_m 7e+76)
(/
2.0
(/
(* (/ k_m t_m) (* (* (sin k_m) (tan k_m)) (/ (/ (pow t_m 3.0) l) l)))
(/ t_m k_m)))
(/ 2.0 (pow (* (/ (pow k_m 2.0) l) (sqrt t_m)) 2.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (t_m <= 1.85e-118) {
tmp = pow(((sqrt(2.0) / k_m) / (k_m * (sqrt(t_m) / l))), 2.0);
} else if (t_m <= 7e+76) {
tmp = 2.0 / (((k_m / t_m) * ((sin(k_m) * tan(k_m)) * ((pow(t_m, 3.0) / l) / l))) / (t_m / k_m));
} else {
tmp = 2.0 / pow(((pow(k_m, 2.0) / l) * sqrt(t_m)), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (t_m <= 1.85d-118) then
tmp = ((sqrt(2.0d0) / k_m) / (k_m * (sqrt(t_m) / l))) ** 2.0d0
else if (t_m <= 7d+76) then
tmp = 2.0d0 / (((k_m / t_m) * ((sin(k_m) * tan(k_m)) * (((t_m ** 3.0d0) / l) / l))) / (t_m / k_m))
else
tmp = 2.0d0 / ((((k_m ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (t_m <= 1.85e-118) {
tmp = Math.pow(((Math.sqrt(2.0) / k_m) / (k_m * (Math.sqrt(t_m) / l))), 2.0);
} else if (t_m <= 7e+76) {
tmp = 2.0 / (((k_m / t_m) * ((Math.sin(k_m) * Math.tan(k_m)) * ((Math.pow(t_m, 3.0) / l) / l))) / (t_m / k_m));
} else {
tmp = 2.0 / Math.pow(((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m)), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if t_m <= 1.85e-118: tmp = math.pow(((math.sqrt(2.0) / k_m) / (k_m * (math.sqrt(t_m) / l))), 2.0) elif t_m <= 7e+76: tmp = 2.0 / (((k_m / t_m) * ((math.sin(k_m) * math.tan(k_m)) * ((math.pow(t_m, 3.0) / l) / l))) / (t_m / k_m)) else: tmp = 2.0 / math.pow(((math.pow(k_m, 2.0) / l) * math.sqrt(t_m)), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (t_m <= 1.85e-118) tmp = Float64(Float64(sqrt(2.0) / k_m) / Float64(k_m * Float64(sqrt(t_m) / l))) ^ 2.0; elseif (t_m <= 7e+76) tmp = Float64(2.0 / Float64(Float64(Float64(k_m / t_m) * Float64(Float64(sin(k_m) * tan(k_m)) * Float64(Float64((t_m ^ 3.0) / l) / l))) / Float64(t_m / k_m))); else tmp = Float64(2.0 / (Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (t_m <= 1.85e-118) tmp = ((sqrt(2.0) / k_m) / (k_m * (sqrt(t_m) / l))) ^ 2.0; elseif (t_m <= 7e+76) tmp = 2.0 / (((k_m / t_m) * ((sin(k_m) * tan(k_m)) * (((t_m ^ 3.0) / l) / l))) / (t_m / k_m)); else tmp = 2.0 / ((((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 1.85e-118], N[Power[N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / N[(k$95$m * N[(N[Sqrt[t$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[t$95$m, 7e+76], N[(2.0 / N[(N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[t$95$m, 3.0], $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.85 \cdot 10^{-118}:\\
\;\;\;\;{\left(\frac{\frac{\sqrt{2}}{k\_m}}{k\_m \cdot \frac{\sqrt{t\_m}}{\ell}}\right)}^{2}\\
\mathbf{elif}\;t\_m \leq 7 \cdot 10^{+76}:\\
\;\;\;\;\frac{2}{\frac{\frac{k\_m}{t\_m} \cdot \left(\left(\sin k\_m \cdot \tan k\_m\right) \cdot \frac{\frac{{t\_m}^{3}}{\ell}}{\ell}\right)}{\frac{t\_m}{k\_m}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\end{array}
\end{array}
if t < 1.85000000000000007e-118Initial program 34.0%
Taylor expanded in t around 0 72.5%
associate-/l*73.0%
*-commutative73.0%
Simplified73.0%
*-un-lft-identity73.0%
associate-/r*73.0%
associate-/l*73.1%
Applied egg-rr73.1%
Taylor expanded in k around 0 62.9%
associate-/l*61.9%
Simplified61.9%
add-sqr-sqrt34.3%
pow234.3%
sqrt-div13.5%
sqrt-div13.5%
sqrt-pow113.5%
metadata-eval13.5%
pow113.5%
sqrt-prod13.5%
sqrt-pow114.9%
metadata-eval14.9%
pow114.9%
sqrt-div13.0%
sqrt-pow115.8%
metadata-eval15.8%
pow115.8%
Applied egg-rr15.8%
if 1.85000000000000007e-118 < t < 7.00000000000000001e76Initial program 74.4%
associate-/r*74.4%
div-inv74.4%
Applied egg-rr74.4%
add-exp-log73.8%
expm1-define73.8%
log1p-define76.1%
unpow276.1%
expm1-log1p-u76.7%
clear-num76.7%
un-div-inv76.7%
Applied egg-rr76.7%
associate-*r/77.1%
associate-*l*77.2%
un-div-inv77.2%
Applied egg-rr77.2%
if 7.00000000000000001e76 < t Initial program 11.1%
add-sqr-sqrt5.4%
pow25.4%
Applied egg-rr40.5%
Taylor expanded in k around 0 85.6%
Final simplification37.2%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= t_m 1.85e-118)
(pow (/ (/ (sqrt 2.0) k_m) (* k_m (/ (sqrt t_m) l))) 2.0)
(if (<= t_m 1.2e+65)
(/
2.0
(*
(* (tan k_m) (* (sin k_m) (/ (/ (pow t_m 3.0) l) l)))
(/ (/ k_m t_m) (/ t_m k_m))))
(/ 2.0 (pow (* (/ (pow k_m 2.0) l) (sqrt t_m)) 2.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (t_m <= 1.85e-118) {
tmp = pow(((sqrt(2.0) / k_m) / (k_m * (sqrt(t_m) / l))), 2.0);
} else if (t_m <= 1.2e+65) {
tmp = 2.0 / ((tan(k_m) * (sin(k_m) * ((pow(t_m, 3.0) / l) / l))) * ((k_m / t_m) / (t_m / k_m)));
} else {
tmp = 2.0 / pow(((pow(k_m, 2.0) / l) * sqrt(t_m)), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (t_m <= 1.85d-118) then
tmp = ((sqrt(2.0d0) / k_m) / (k_m * (sqrt(t_m) / l))) ** 2.0d0
else if (t_m <= 1.2d+65) then
tmp = 2.0d0 / ((tan(k_m) * (sin(k_m) * (((t_m ** 3.0d0) / l) / l))) * ((k_m / t_m) / (t_m / k_m)))
else
tmp = 2.0d0 / ((((k_m ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (t_m <= 1.85e-118) {
tmp = Math.pow(((Math.sqrt(2.0) / k_m) / (k_m * (Math.sqrt(t_m) / l))), 2.0);
} else if (t_m <= 1.2e+65) {
tmp = 2.0 / ((Math.tan(k_m) * (Math.sin(k_m) * ((Math.pow(t_m, 3.0) / l) / l))) * ((k_m / t_m) / (t_m / k_m)));
} else {
tmp = 2.0 / Math.pow(((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m)), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if t_m <= 1.85e-118: tmp = math.pow(((math.sqrt(2.0) / k_m) / (k_m * (math.sqrt(t_m) / l))), 2.0) elif t_m <= 1.2e+65: tmp = 2.0 / ((math.tan(k_m) * (math.sin(k_m) * ((math.pow(t_m, 3.0) / l) / l))) * ((k_m / t_m) / (t_m / k_m))) else: tmp = 2.0 / math.pow(((math.pow(k_m, 2.0) / l) * math.sqrt(t_m)), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (t_m <= 1.85e-118) tmp = Float64(Float64(sqrt(2.0) / k_m) / Float64(k_m * Float64(sqrt(t_m) / l))) ^ 2.0; elseif (t_m <= 1.2e+65) tmp = Float64(2.0 / Float64(Float64(tan(k_m) * Float64(sin(k_m) * Float64(Float64((t_m ^ 3.0) / l) / l))) * Float64(Float64(k_m / t_m) / Float64(t_m / k_m)))); else tmp = Float64(2.0 / (Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (t_m <= 1.85e-118) tmp = ((sqrt(2.0) / k_m) / (k_m * (sqrt(t_m) / l))) ^ 2.0; elseif (t_m <= 1.2e+65) tmp = 2.0 / ((tan(k_m) * (sin(k_m) * (((t_m ^ 3.0) / l) / l))) * ((k_m / t_m) / (t_m / k_m))); else tmp = 2.0 / ((((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 1.85e-118], N[Power[N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / N[(k$95$m * N[(N[Sqrt[t$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[t$95$m, 1.2e+65], N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(N[Power[t$95$m, 3.0], $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(k$95$m / t$95$m), $MachinePrecision] / N[(t$95$m / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.85 \cdot 10^{-118}:\\
\;\;\;\;{\left(\frac{\frac{\sqrt{2}}{k\_m}}{k\_m \cdot \frac{\sqrt{t\_m}}{\ell}}\right)}^{2}\\
\mathbf{elif}\;t\_m \leq 1.2 \cdot 10^{+65}:\\
\;\;\;\;\frac{2}{\left(\tan k\_m \cdot \left(\sin k\_m \cdot \frac{\frac{{t\_m}^{3}}{\ell}}{\ell}\right)\right) \cdot \frac{\frac{k\_m}{t\_m}}{\frac{t\_m}{k\_m}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\end{array}
\end{array}
if t < 1.85000000000000007e-118Initial program 34.0%
Taylor expanded in t around 0 72.5%
associate-/l*73.0%
*-commutative73.0%
Simplified73.0%
*-un-lft-identity73.0%
associate-/r*73.0%
associate-/l*73.1%
Applied egg-rr73.1%
Taylor expanded in k around 0 62.9%
associate-/l*61.9%
Simplified61.9%
add-sqr-sqrt34.3%
pow234.3%
sqrt-div13.5%
sqrt-div13.5%
sqrt-pow113.5%
metadata-eval13.5%
pow113.5%
sqrt-prod13.5%
sqrt-pow114.9%
metadata-eval14.9%
pow114.9%
sqrt-div13.0%
sqrt-pow115.8%
metadata-eval15.8%
pow115.8%
Applied egg-rr15.8%
if 1.85000000000000007e-118 < t < 1.2000000000000001e65Initial program 76.1%
associate-/r*76.1%
div-inv76.1%
Applied egg-rr76.1%
add-exp-log75.5%
expm1-define75.5%
log1p-define77.8%
unpow277.8%
expm1-log1p-u78.4%
clear-num78.4%
un-div-inv78.4%
Applied egg-rr78.4%
un-div-inv78.3%
Applied egg-rr78.3%
if 1.2000000000000001e65 < t Initial program 16.9%
add-sqr-sqrt9.5%
pow29.5%
Applied egg-rr40.6%
Taylor expanded in k around 0 85.1%
Final simplification37.4%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 6.2e+14)
(/ 2.0 (pow (* (/ (pow k_m 2.0) l) (sqrt t_m)) 2.0))
(if (<= k_m 1.45e+71)
(*
(* l l)
(/
(+ (* -0.3333333333333333 (/ (pow k_m 2.0) t_m)) (* 2.0 (/ 1.0 t_m)))
(pow k_m 4.0)))
(pow (/ (/ (sqrt 2.0) k_m) (* k_m (/ (sqrt t_m) l))) 2.0)))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 6.2e+14) {
tmp = 2.0 / pow(((pow(k_m, 2.0) / l) * sqrt(t_m)), 2.0);
} else if (k_m <= 1.45e+71) {
tmp = (l * l) * (((-0.3333333333333333 * (pow(k_m, 2.0) / t_m)) + (2.0 * (1.0 / t_m))) / pow(k_m, 4.0));
} else {
tmp = pow(((sqrt(2.0) / k_m) / (k_m * (sqrt(t_m) / l))), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 6.2d+14) then
tmp = 2.0d0 / ((((k_m ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if (k_m <= 1.45d+71) then
tmp = (l * l) * ((((-0.3333333333333333d0) * ((k_m ** 2.0d0) / t_m)) + (2.0d0 * (1.0d0 / t_m))) / (k_m ** 4.0d0))
else
tmp = ((sqrt(2.0d0) / k_m) / (k_m * (sqrt(t_m) / l))) ** 2.0d0
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 6.2e+14) {
tmp = 2.0 / Math.pow(((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if (k_m <= 1.45e+71) {
tmp = (l * l) * (((-0.3333333333333333 * (Math.pow(k_m, 2.0) / t_m)) + (2.0 * (1.0 / t_m))) / Math.pow(k_m, 4.0));
} else {
tmp = Math.pow(((Math.sqrt(2.0) / k_m) / (k_m * (Math.sqrt(t_m) / l))), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 6.2e+14: tmp = 2.0 / math.pow(((math.pow(k_m, 2.0) / l) * math.sqrt(t_m)), 2.0) elif k_m <= 1.45e+71: tmp = (l * l) * (((-0.3333333333333333 * (math.pow(k_m, 2.0) / t_m)) + (2.0 * (1.0 / t_m))) / math.pow(k_m, 4.0)) else: tmp = math.pow(((math.sqrt(2.0) / k_m) / (k_m * (math.sqrt(t_m) / l))), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 6.2e+14) tmp = Float64(2.0 / (Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (k_m <= 1.45e+71) tmp = Float64(Float64(l * l) * Float64(Float64(Float64(-0.3333333333333333 * Float64((k_m ^ 2.0) / t_m)) + Float64(2.0 * Float64(1.0 / t_m))) / (k_m ^ 4.0))); else tmp = Float64(Float64(sqrt(2.0) / k_m) / Float64(k_m * Float64(sqrt(t_m) / l))) ^ 2.0; end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 6.2e+14) tmp = 2.0 / ((((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif (k_m <= 1.45e+71) tmp = (l * l) * (((-0.3333333333333333 * ((k_m ^ 2.0) / t_m)) + (2.0 * (1.0 / t_m))) / (k_m ^ 4.0)); else tmp = ((sqrt(2.0) / k_m) / (k_m * (sqrt(t_m) / l))) ^ 2.0; end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 6.2e+14], N[(2.0 / N[Power[N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.45e+71], N[(N[(l * l), $MachinePrecision] * N[(N[(N[(-0.3333333333333333 * N[(N[Power[k$95$m, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision] + N[(2.0 * N[(1.0 / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / N[(k$95$m * N[(N[Sqrt[t$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 6.2 \cdot 10^{+14}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;k\_m \leq 1.45 \cdot 10^{+71}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{-0.3333333333333333 \cdot \frac{{k\_m}^{2}}{t\_m} + 2 \cdot \frac{1}{t\_m}}{{k\_m}^{4}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{\frac{\sqrt{2}}{k\_m}}{k\_m \cdot \frac{\sqrt{t\_m}}{\ell}}\right)}^{2}\\
\end{array}
\end{array}
if k < 6.2e14Initial program 40.3%
add-sqr-sqrt20.2%
pow220.2%
Applied egg-rr26.4%
Taylor expanded in k around 0 39.4%
if 6.2e14 < k < 1.45000000000000004e71Initial program 23.3%
Simplified38.5%
Taylor expanded in k around 0 53.1%
if 1.45000000000000004e71 < k Initial program 33.3%
Taylor expanded in t around 0 56.2%
associate-/l*54.4%
*-commutative54.4%
Simplified54.4%
*-un-lft-identity54.4%
associate-/r*54.4%
associate-/l*54.4%
Applied egg-rr54.4%
Taylor expanded in k around 0 50.3%
associate-/l*50.4%
Simplified50.4%
add-sqr-sqrt50.2%
pow250.2%
sqrt-div25.2%
sqrt-div25.2%
sqrt-pow125.2%
metadata-eval25.2%
pow125.2%
sqrt-prod25.2%
sqrt-pow125.7%
metadata-eval25.7%
pow125.7%
sqrt-div25.3%
sqrt-pow123.6%
metadata-eval23.6%
pow123.6%
Applied egg-rr23.6%
Final simplification37.1%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (/ 2.0 (pow (* (/ (pow k_m 2.0) l) (sqrt t_m)) 2.0))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 / pow(((pow(k_m, 2.0) / l) * sqrt(t_m)), 2.0));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 / ((((k_m ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 / Math.pow(((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m)), 2.0));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 / math.pow(((math.pow(k_m, 2.0) / l) * math.sqrt(t_m)), 2.0))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 / (Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 / ((((k_m ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * N[(2.0 / N[Power[N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \frac{2}{{\left(\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}
\end{array}
Initial program 38.1%
add-sqr-sqrt18.9%
pow218.9%
Applied egg-rr23.1%
Taylor expanded in k around 0 34.5%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 0.0)
(/ 2.0 (pow (* (/ k_m t_m) (* k_m (/ (pow t_m 1.5) l))) 2.0))
(/ (/ 2.0 (pow k_m 2.0)) (* (* k_m k_m) (/ t_m (pow l 2.0)))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / pow(((k_m / t_m) * (k_m * (pow(t_m, 1.5) / l))), 2.0);
} else {
tmp = (2.0 / pow(k_m, 2.0)) / ((k_m * k_m) * (t_m / pow(l, 2.0)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 0.0d0) then
tmp = 2.0d0 / (((k_m / t_m) * (k_m * ((t_m ** 1.5d0) / l))) ** 2.0d0)
else
tmp = (2.0d0 / (k_m ** 2.0d0)) / ((k_m * k_m) * (t_m / (l ** 2.0d0)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / Math.pow(((k_m / t_m) * (k_m * (Math.pow(t_m, 1.5) / l))), 2.0);
} else {
tmp = (2.0 / Math.pow(k_m, 2.0)) / ((k_m * k_m) * (t_m / Math.pow(l, 2.0)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 0.0: tmp = 2.0 / math.pow(((k_m / t_m) * (k_m * (math.pow(t_m, 1.5) / l))), 2.0) else: tmp = (2.0 / math.pow(k_m, 2.0)) / ((k_m * k_m) * (t_m / math.pow(l, 2.0))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 0.0) tmp = Float64(2.0 / (Float64(Float64(k_m / t_m) * Float64(k_m * Float64((t_m ^ 1.5) / l))) ^ 2.0)); else tmp = Float64(Float64(2.0 / (k_m ^ 2.0)) / Float64(Float64(k_m * k_m) * Float64(t_m / (l ^ 2.0)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 0.0) tmp = 2.0 / (((k_m / t_m) * (k_m * ((t_m ^ 1.5) / l))) ^ 2.0); else tmp = (2.0 / (k_m ^ 2.0)) / ((k_m * k_m) * (t_m / (l ^ 2.0))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 0.0], N[(2.0 / N[Power[N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(k$95$m * N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] * N[(t$95$m / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 0:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{t\_m} \cdot \left(k\_m \cdot \frac{{t\_m}^{1.5}}{\ell}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{{k\_m}^{2}}}{\left(k\_m \cdot k\_m\right) \cdot \frac{t\_m}{{\ell}^{2}}}\\
\end{array}
\end{array}
if (*.f64 l l) < 0.0Initial program 25.9%
add-sqr-sqrt9.3%
pow29.3%
Applied egg-rr13.2%
Taylor expanded in k around 0 20.6%
if 0.0 < (*.f64 l l) Initial program 41.3%
Taylor expanded in t around 0 78.8%
associate-/l*79.4%
*-commutative79.4%
Simplified79.4%
*-un-lft-identity79.4%
associate-/r*79.6%
associate-/l*79.7%
Applied egg-rr79.7%
Taylor expanded in k around 0 67.3%
associate-/l*65.8%
Simplified65.8%
unpow279.4%
Applied egg-rr65.8%
Final simplification56.3%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 0.0)
(/ 2.0 (pow (* (/ k_m t_m) (* k_m (/ (pow t_m 1.5) l))) 2.0))
(* (* l l) (* (/ 2.0 t_m) (pow k_m -4.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / pow(((k_m / t_m) * (k_m * (pow(t_m, 1.5) / l))), 2.0);
} else {
tmp = (l * l) * ((2.0 / t_m) * pow(k_m, -4.0));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 0.0d0) then
tmp = 2.0d0 / (((k_m / t_m) * (k_m * ((t_m ** 1.5d0) / l))) ** 2.0d0)
else
tmp = (l * l) * ((2.0d0 / t_m) * (k_m ** (-4.0d0)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / Math.pow(((k_m / t_m) * (k_m * (Math.pow(t_m, 1.5) / l))), 2.0);
} else {
tmp = (l * l) * ((2.0 / t_m) * Math.pow(k_m, -4.0));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 0.0: tmp = 2.0 / math.pow(((k_m / t_m) * (k_m * (math.pow(t_m, 1.5) / l))), 2.0) else: tmp = (l * l) * ((2.0 / t_m) * math.pow(k_m, -4.0)) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 0.0) tmp = Float64(2.0 / (Float64(Float64(k_m / t_m) * Float64(k_m * Float64((t_m ^ 1.5) / l))) ^ 2.0)); else tmp = Float64(Float64(l * l) * Float64(Float64(2.0 / t_m) * (k_m ^ -4.0))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 0.0) tmp = 2.0 / (((k_m / t_m) * (k_m * ((t_m ^ 1.5) / l))) ^ 2.0); else tmp = (l * l) * ((2.0 / t_m) * (k_m ^ -4.0)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 0.0], N[(2.0 / N[Power[N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(k$95$m * N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(N[(2.0 / t$95$m), $MachinePrecision] * N[Power[k$95$m, -4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 0:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{t\_m} \cdot \left(k\_m \cdot \frac{{t\_m}^{1.5}}{\ell}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \left(\frac{2}{t\_m} \cdot {k\_m}^{-4}\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 0.0Initial program 25.9%
add-sqr-sqrt9.3%
pow29.3%
Applied egg-rr13.2%
Taylor expanded in k around 0 20.6%
if 0.0 < (*.f64 l l) Initial program 41.3%
Simplified44.9%
Taylor expanded in k around 0 64.5%
*-un-lft-identity64.5%
associate-/r*65.0%
Applied egg-rr65.0%
Taylor expanded in k around 0 64.5%
*-commutative64.5%
associate-/r*65.0%
Simplified65.0%
associate-/l/64.5%
*-un-lft-identity64.5%
associate-/l/65.0%
div-inv65.0%
pow-flip65.0%
metadata-eval65.0%
Applied egg-rr65.0%
*-lft-identity65.0%
Simplified65.0%
Final simplification55.7%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (/ 2.0 (* t_m (/ (pow k_m 4.0) (pow l 2.0))))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 / (t_m * (pow(k_m, 4.0) / pow(l, 2.0))));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 / (t_m * ((k_m ** 4.0d0) / (l ** 2.0d0))))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 / (t_m * (Math.pow(k_m, 4.0) / Math.pow(l, 2.0))));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 / (t_m * (math.pow(k_m, 4.0) / math.pow(l, 2.0))))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 / Float64(t_m * Float64((k_m ^ 4.0) / (l ^ 2.0))))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 / (t_m * ((k_m ^ 4.0) / (l ^ 2.0)))); end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * N[(2.0 / N[(t$95$m * N[(N[Power[k$95$m, 4.0], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \frac{2}{t\_m \cdot \frac{{k\_m}^{4}}{{\ell}^{2}}}
\end{array}
Initial program 38.1%
associate-/r*42.9%
div-inv42.9%
Applied egg-rr42.9%
add-exp-log42.6%
expm1-define42.6%
log1p-define47.9%
unpow247.9%
expm1-log1p-u48.1%
clear-num48.1%
un-div-inv48.1%
Applied egg-rr48.1%
Taylor expanded in t around 0 74.0%
associate-/l*74.5%
*-commutative74.5%
associate-*r/74.6%
*-commutative74.6%
associate-*l*74.9%
*-commutative74.9%
Simplified74.9%
Taylor expanded in k around 0 63.1%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* (* l l) (* (/ 2.0 t_m) (pow k_m -4.0)))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * ((2.0 / t_m) * pow(k_m, -4.0)));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * l) * ((2.0d0 / t_m) * (k_m ** (-4.0d0))))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * ((2.0 / t_m) * Math.pow(k_m, -4.0)));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * ((l * l) * ((2.0 / t_m) * math.pow(k_m, -4.0)))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64(l * l) * Float64(Float64(2.0 / t_m) * (k_m ^ -4.0)))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * l) * ((2.0 / t_m) * (k_m ^ -4.0))); end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * N[(N[(l * l), $MachinePrecision] * N[(N[(2.0 / t$95$m), $MachinePrecision] * N[Power[k$95$m, -4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(\ell \cdot \ell\right) \cdot \left(\frac{2}{t\_m} \cdot {k\_m}^{-4}\right)\right)
\end{array}
Initial program 38.1%
Simplified42.9%
Taylor expanded in k around 0 62.7%
*-un-lft-identity62.7%
associate-/r*63.1%
Applied egg-rr63.1%
Taylor expanded in k around 0 62.7%
*-commutative62.7%
associate-/r*63.0%
Simplified63.0%
associate-/l/62.7%
*-un-lft-identity62.7%
associate-/l/63.0%
div-inv63.1%
pow-flip63.1%
metadata-eval63.1%
Applied egg-rr63.1%
*-lft-identity63.1%
Simplified63.1%
Final simplification63.1%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* (* l l) (/ (/ 2.0 (pow k_m 4.0)) t_m))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * ((2.0 / pow(k_m, 4.0)) / t_m));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * l) * ((2.0d0 / (k_m ** 4.0d0)) / t_m))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * ((2.0 / Math.pow(k_m, 4.0)) / t_m));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * ((l * l) * ((2.0 / math.pow(k_m, 4.0)) / t_m))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64(l * l) * Float64(Float64(2.0 / (k_m ^ 4.0)) / t_m))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * l) * ((2.0 / (k_m ^ 4.0)) / t_m)); end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * N[(N[(l * l), $MachinePrecision] * N[(N[(2.0 / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(\ell \cdot \ell\right) \cdot \frac{\frac{2}{{k\_m}^{4}}}{t\_m}\right)
\end{array}
Initial program 38.1%
Simplified42.9%
Taylor expanded in k around 0 62.7%
*-un-lft-identity62.7%
associate-/r*63.1%
Applied egg-rr63.1%
Taylor expanded in k around 0 62.7%
*-commutative62.7%
associate-/r*63.0%
Simplified63.0%
*-un-lft-identity63.0%
associate-/l/62.7%
associate-/r*63.1%
Applied egg-rr63.1%
Final simplification63.1%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* (* l l) (/ 2.0 (* t_m (pow k_m 4.0))))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (2.0 / (t_m * pow(k_m, 4.0))));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * l) * (2.0d0 / (t_m * (k_m ** 4.0d0))))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (2.0 / (t_m * Math.pow(k_m, 4.0))));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * ((l * l) * (2.0 / (t_m * math.pow(k_m, 4.0))))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64(l * l) * Float64(2.0 / Float64(t_m * (k_m ^ 4.0))))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * l) * (2.0 / (t_m * (k_m ^ 4.0)))); end
k_m = N[Abs[k], $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_, k$95$m_] := N[(t$95$s * N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(t$95$m * N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(\ell \cdot \ell\right) \cdot \frac{2}{t\_m \cdot {k\_m}^{4}}\right)
\end{array}
Initial program 38.1%
Simplified42.9%
Taylor expanded in k around 0 62.7%
Final simplification62.7%
herbie shell --seed 2024112
(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))))