
(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}
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= (* l l) 1e-314)
(/ 2.0 (pow (* (/ (pow k 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+299)
(/
2.0
(*
(/ (* k (* (sqrt t_m) (sin k))) (pow l 2.0))
(/ (* (sqrt t_m) (* k (sin k))) (cos k))))
(*
(/ l k)
(* 2.0 (* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k)))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / pow(((pow(k, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+299) {
tmp = 2.0 / (((k * (sqrt(t_m) * sin(k))) / pow(l, 2.0)) * ((sqrt(t_m) * (k * sin(k))) / cos(k)));
} else {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if ((l * l) <= 1d-314) then
tmp = 2.0d0 / ((((k ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+299) then
tmp = 2.0d0 / (((k * (sqrt(t_m) * sin(k))) / (l ** 2.0d0)) * ((sqrt(t_m) * (k * sin(k))) / cos(k)))
else
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
end if
code = t_s * tmp
end function
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) {
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / Math.pow(((Math.pow(k, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+299) {
tmp = 2.0 / (((k * (Math.sqrt(t_m) * Math.sin(k))) / Math.pow(l, 2.0)) * ((Math.sqrt(t_m) * (k * Math.sin(k))) / Math.cos(k)));
} else {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if (l * l) <= 1e-314: tmp = 2.0 / math.pow(((math.pow(k, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+299: tmp = 2.0 / (((k * (math.sqrt(t_m) * math.sin(k))) / math.pow(l, 2.0)) * ((math.sqrt(t_m) * (k * math.sin(k))) / math.cos(k))) else: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-314) tmp = Float64(2.0 / (Float64(Float64((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+299) tmp = Float64(2.0 / Float64(Float64(Float64(k * Float64(sqrt(t_m) * sin(k))) / (l ^ 2.0)) * Float64(Float64(sqrt(t_m) * Float64(k * sin(k))) / cos(k)))); else tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if ((l * l) <= 1e-314) tmp = 2.0 / ((((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+299) tmp = 2.0 / (((k * (sqrt(t_m) * sin(k))) / (l ^ 2.0)) * ((sqrt(t_m) * (k * sin(k))) / cos(k))); else tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 1e-314], N[(2.0 / N[Power[N[(N[(N[Power[k, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+299], N[(2.0 / N[(N[(N[(k * N[(N[Sqrt[t$95$m], $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Sqrt[t$95$m], $MachinePrecision] * N[(k * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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 10^{-314}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+299}:\\
\;\;\;\;\frac{2}{\frac{k \cdot \left(\sqrt{t\_m} \cdot \sin k\right)}{{\ell}^{2}} \cdot \frac{\sqrt{t\_m} \cdot \left(k \cdot \sin k\right)}{\cos k}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 9.9999999996e-315Initial program 29.6%
Taylor expanded in t around 0 68.9%
associate-/l*68.9%
Simplified68.9%
add-sqr-sqrt30.1%
pow230.1%
Applied egg-rr32.7%
associate-/r*32.7%
Simplified32.7%
Taylor expanded in k around 0 52.9%
if 9.9999999996e-315 < (*.f64 l l) < 2.0000000000000001e299Initial program 39.5%
Taylor expanded in t around 0 84.1%
add-sqr-sqrt42.8%
times-frac42.8%
sqrt-prod42.8%
sqrt-pow136.9%
metadata-eval36.9%
pow136.9%
*-commutative36.9%
sqrt-prod29.3%
sqrt-pow137.7%
metadata-eval37.7%
pow137.7%
Applied egg-rr49.8%
Taylor expanded in k around inf 49.8%
if 2.0000000000000001e299 < (*.f64 l l) Initial program 44.5%
associate-/r*44.5%
associate-/l/44.5%
associate-*l/44.5%
associate-/r/44.5%
+-commutative44.5%
unpow244.5%
sqr-neg44.5%
distribute-frac-neg244.5%
distribute-frac-neg244.5%
unpow244.5%
+-rgt-identity44.5%
metadata-eval44.5%
associate--l+44.5%
+-commutative44.5%
associate--l+44.5%
Simplified44.5%
associate-*r*52.4%
unpow252.4%
times-frac67.4%
div-inv67.4%
frac-times67.4%
metadata-eval67.4%
Applied egg-rr67.4%
Taylor expanded in k around inf 63.8%
times-frac72.5%
*-commutative72.5%
Simplified72.5%
associate-*r/71.0%
associate-*r*71.0%
div-inv71.0%
pow-prod-down74.7%
pow-flip74.8%
metadata-eval74.8%
Applied egg-rr74.8%
associate-/l*76.3%
associate-/r/81.7%
*-commutative81.7%
associate-*l*83.4%
*-commutative83.4%
associate-*l*83.4%
*-commutative83.4%
Simplified83.4%
Final simplification57.7%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= (* l l) 1.5e-210)
(/ 2.0 (pow (* (/ (pow k 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+299)
(/
2.0
(*
(* k (* (sin k) (/ (sqrt t_m) (pow l 2.0))))
(* k (/ (* (sqrt t_m) (sin k)) (cos k)))))
(*
(/ l k)
(* 2.0 (* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k)))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if ((l * l) <= 1.5e-210) {
tmp = 2.0 / pow(((pow(k, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+299) {
tmp = 2.0 / ((k * (sin(k) * (sqrt(t_m) / pow(l, 2.0)))) * (k * ((sqrt(t_m) * sin(k)) / cos(k))));
} else {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if ((l * l) <= 1.5d-210) then
tmp = 2.0d0 / ((((k ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+299) then
tmp = 2.0d0 / ((k * (sin(k) * (sqrt(t_m) / (l ** 2.0d0)))) * (k * ((sqrt(t_m) * sin(k)) / cos(k))))
else
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
end if
code = t_s * tmp
end function
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) {
double tmp;
if ((l * l) <= 1.5e-210) {
tmp = 2.0 / Math.pow(((Math.pow(k, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+299) {
tmp = 2.0 / ((k * (Math.sin(k) * (Math.sqrt(t_m) / Math.pow(l, 2.0)))) * (k * ((Math.sqrt(t_m) * Math.sin(k)) / Math.cos(k))));
} else {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if (l * l) <= 1.5e-210: tmp = 2.0 / math.pow(((math.pow(k, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+299: tmp = 2.0 / ((k * (math.sin(k) * (math.sqrt(t_m) / math.pow(l, 2.0)))) * (k * ((math.sqrt(t_m) * math.sin(k)) / math.cos(k)))) else: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (Float64(l * l) <= 1.5e-210) tmp = Float64(2.0 / (Float64(Float64((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+299) tmp = Float64(2.0 / Float64(Float64(k * Float64(sin(k) * Float64(sqrt(t_m) / (l ^ 2.0)))) * Float64(k * Float64(Float64(sqrt(t_m) * sin(k)) / cos(k))))); else tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if ((l * l) <= 1.5e-210) tmp = 2.0 / ((((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+299) tmp = 2.0 / ((k * (sin(k) * (sqrt(t_m) / (l ^ 2.0)))) * (k * ((sqrt(t_m) * sin(k)) / cos(k)))); else tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 1.5e-210], N[(2.0 / N[Power[N[(N[(N[Power[k, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+299], N[(2.0 / N[(N[(k * N[(N[Sin[k], $MachinePrecision] * N[(N[Sqrt[t$95$m], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(k * N[(N[(N[Sqrt[t$95$m], $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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 1.5 \cdot 10^{-210}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+299}:\\
\;\;\;\;\frac{2}{\left(k \cdot \left(\sin k \cdot \frac{\sqrt{t\_m}}{{\ell}^{2}}\right)\right) \cdot \left(k \cdot \frac{\sqrt{t\_m} \cdot \sin k}{\cos k}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 1.5000000000000001e-210Initial program 31.9%
Taylor expanded in t around 0 74.9%
associate-/l*74.9%
Simplified74.9%
add-sqr-sqrt35.7%
pow235.7%
Applied egg-rr36.7%
associate-/r*36.7%
Simplified36.7%
Taylor expanded in k around 0 56.4%
if 1.5000000000000001e-210 < (*.f64 l l) < 2.0000000000000001e299Initial program 39.3%
Taylor expanded in t around 0 81.8%
add-sqr-sqrt37.8%
times-frac37.8%
sqrt-prod37.8%
sqrt-pow131.0%
metadata-eval31.0%
pow131.0%
*-commutative31.0%
sqrt-prod23.2%
sqrt-pow132.8%
metadata-eval32.8%
pow132.8%
Applied egg-rr45.8%
associate-/l*44.2%
associate-/l*44.2%
associate-/l*44.2%
Simplified44.2%
if 2.0000000000000001e299 < (*.f64 l l) Initial program 44.5%
associate-/r*44.5%
associate-/l/44.5%
associate-*l/44.5%
associate-/r/44.5%
+-commutative44.5%
unpow244.5%
sqr-neg44.5%
distribute-frac-neg244.5%
distribute-frac-neg244.5%
unpow244.5%
+-rgt-identity44.5%
metadata-eval44.5%
associate--l+44.5%
+-commutative44.5%
associate--l+44.5%
Simplified44.5%
associate-*r*52.4%
unpow252.4%
times-frac67.4%
div-inv67.4%
frac-times67.4%
metadata-eval67.4%
Applied egg-rr67.4%
Taylor expanded in k around inf 63.8%
times-frac72.5%
*-commutative72.5%
Simplified72.5%
associate-*r/71.0%
associate-*r*71.0%
div-inv71.0%
pow-prod-down74.7%
pow-flip74.8%
metadata-eval74.8%
Applied egg-rr74.8%
associate-/l*76.3%
associate-/r/81.7%
*-commutative81.7%
associate-*l*83.4%
*-commutative83.4%
associate-*l*83.4%
*-commutative83.4%
Simplified83.4%
Final simplification56.6%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= (* l l) 1e-314)
(/ 2.0 (pow (* (/ (pow k 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+299)
(/
2.0
(*
(/ (* (sqrt t_m) (* k (sin k))) (cos k))
(* (* k (* (sqrt t_m) (sin k))) (pow l -2.0))))
(*
(/ l k)
(* 2.0 (* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k)))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / pow(((pow(k, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+299) {
tmp = 2.0 / (((sqrt(t_m) * (k * sin(k))) / cos(k)) * ((k * (sqrt(t_m) * sin(k))) * pow(l, -2.0)));
} else {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if ((l * l) <= 1d-314) then
tmp = 2.0d0 / ((((k ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+299) then
tmp = 2.0d0 / (((sqrt(t_m) * (k * sin(k))) / cos(k)) * ((k * (sqrt(t_m) * sin(k))) * (l ** (-2.0d0))))
else
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
end if
code = t_s * tmp
end function
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) {
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / Math.pow(((Math.pow(k, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+299) {
tmp = 2.0 / (((Math.sqrt(t_m) * (k * Math.sin(k))) / Math.cos(k)) * ((k * (Math.sqrt(t_m) * Math.sin(k))) * Math.pow(l, -2.0)));
} else {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if (l * l) <= 1e-314: tmp = 2.0 / math.pow(((math.pow(k, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+299: tmp = 2.0 / (((math.sqrt(t_m) * (k * math.sin(k))) / math.cos(k)) * ((k * (math.sqrt(t_m) * math.sin(k))) * math.pow(l, -2.0))) else: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-314) tmp = Float64(2.0 / (Float64(Float64((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+299) tmp = Float64(2.0 / Float64(Float64(Float64(sqrt(t_m) * Float64(k * sin(k))) / cos(k)) * Float64(Float64(k * Float64(sqrt(t_m) * sin(k))) * (l ^ -2.0)))); else tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if ((l * l) <= 1e-314) tmp = 2.0 / ((((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+299) tmp = 2.0 / (((sqrt(t_m) * (k * sin(k))) / cos(k)) * ((k * (sqrt(t_m) * sin(k))) * (l ^ -2.0))); else tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 1e-314], N[(2.0 / N[Power[N[(N[(N[Power[k, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+299], N[(2.0 / N[(N[(N[(N[Sqrt[t$95$m], $MachinePrecision] * N[(k * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision] * N[(N[(k * N[(N[Sqrt[t$95$m], $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Power[l, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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 10^{-314}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+299}:\\
\;\;\;\;\frac{2}{\frac{\sqrt{t\_m} \cdot \left(k \cdot \sin k\right)}{\cos k} \cdot \left(\left(k \cdot \left(\sqrt{t\_m} \cdot \sin k\right)\right) \cdot {\ell}^{-2}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 9.9999999996e-315Initial program 29.6%
Taylor expanded in t around 0 68.9%
associate-/l*68.9%
Simplified68.9%
add-sqr-sqrt30.1%
pow230.1%
Applied egg-rr32.7%
associate-/r*32.7%
Simplified32.7%
Taylor expanded in k around 0 52.9%
if 9.9999999996e-315 < (*.f64 l l) < 2.0000000000000001e299Initial program 39.5%
Taylor expanded in t around 0 84.1%
add-sqr-sqrt42.8%
times-frac42.8%
sqrt-prod42.8%
sqrt-pow136.9%
metadata-eval36.9%
pow136.9%
*-commutative36.9%
sqrt-prod29.3%
sqrt-pow137.7%
metadata-eval37.7%
pow137.7%
Applied egg-rr49.8%
Taylor expanded in k around inf 49.8%
pow249.8%
div-inv49.8%
pow249.8%
pow-flip49.8%
metadata-eval49.8%
Applied egg-rr49.8%
if 2.0000000000000001e299 < (*.f64 l l) Initial program 44.5%
associate-/r*44.5%
associate-/l/44.5%
associate-*l/44.5%
associate-/r/44.5%
+-commutative44.5%
unpow244.5%
sqr-neg44.5%
distribute-frac-neg244.5%
distribute-frac-neg244.5%
unpow244.5%
+-rgt-identity44.5%
metadata-eval44.5%
associate--l+44.5%
+-commutative44.5%
associate--l+44.5%
Simplified44.5%
associate-*r*52.4%
unpow252.4%
times-frac67.4%
div-inv67.4%
frac-times67.4%
metadata-eval67.4%
Applied egg-rr67.4%
Taylor expanded in k around inf 63.8%
times-frac72.5%
*-commutative72.5%
Simplified72.5%
associate-*r/71.0%
associate-*r*71.0%
div-inv71.0%
pow-prod-down74.7%
pow-flip74.8%
metadata-eval74.8%
Applied egg-rr74.8%
associate-/l*76.3%
associate-/r/81.7%
*-commutative81.7%
associate-*l*83.4%
*-commutative83.4%
associate-*l*83.4%
*-commutative83.4%
Simplified83.4%
Final simplification57.7%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (* k (* (sqrt t_m) (sin k)))))
(*
t_s
(if (<= (* l l) 1e-314)
(/ 2.0 (pow (* (/ (pow k 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+299)
(/ 2.0 (* (/ t_2 (pow l 2.0)) (/ t_2 (cos k))))
(*
(/ l k)
(*
2.0
(* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k))))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = k * (sqrt(t_m) * sin(k));
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / pow(((pow(k, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+299) {
tmp = 2.0 / ((t_2 / pow(l, 2.0)) * (t_2 / cos(k)));
} else {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = k * (sqrt(t_m) * sin(k))
if ((l * l) <= 1d-314) then
tmp = 2.0d0 / ((((k ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+299) then
tmp = 2.0d0 / ((t_2 / (l ** 2.0d0)) * (t_2 / cos(k)))
else
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
end if
code = t_s * tmp
end function
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) {
double t_2 = k * (Math.sqrt(t_m) * Math.sin(k));
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / Math.pow(((Math.pow(k, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+299) {
tmp = 2.0 / ((t_2 / Math.pow(l, 2.0)) * (t_2 / Math.cos(k)));
} else {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): t_2 = k * (math.sqrt(t_m) * math.sin(k)) tmp = 0 if (l * l) <= 1e-314: tmp = 2.0 / math.pow(((math.pow(k, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+299: tmp = 2.0 / ((t_2 / math.pow(l, 2.0)) * (t_2 / math.cos(k))) else: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(k * Float64(sqrt(t_m) * sin(k))) tmp = 0.0 if (Float64(l * l) <= 1e-314) tmp = Float64(2.0 / (Float64(Float64((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+299) tmp = Float64(2.0 / Float64(Float64(t_2 / (l ^ 2.0)) * Float64(t_2 / cos(k)))); else tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) t_2 = k * (sqrt(t_m) * sin(k)); tmp = 0.0; if ((l * l) <= 1e-314) tmp = 2.0 / ((((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+299) tmp = 2.0 / ((t_2 / (l ^ 2.0)) * (t_2 / cos(k))); else tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); end tmp_2 = t_s * tmp; end
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_] := Block[{t$95$2 = N[(k * N[(N[Sqrt[t$95$m], $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 1e-314], N[(2.0 / N[Power[N[(N[(N[Power[k, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+299], N[(2.0 / N[(N[(t$95$2 / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$2 / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := k \cdot \left(\sqrt{t\_m} \cdot \sin k\right)\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-314}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+299}:\\
\;\;\;\;\frac{2}{\frac{t\_2}{{\ell}^{2}} \cdot \frac{t\_2}{\cos k}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\end{array}
\end{array}
\end{array}
if (*.f64 l l) < 9.9999999996e-315Initial program 29.6%
Taylor expanded in t around 0 68.9%
associate-/l*68.9%
Simplified68.9%
add-sqr-sqrt30.1%
pow230.1%
Applied egg-rr32.7%
associate-/r*32.7%
Simplified32.7%
Taylor expanded in k around 0 52.9%
if 9.9999999996e-315 < (*.f64 l l) < 2.0000000000000001e299Initial program 39.5%
Taylor expanded in t around 0 84.1%
add-sqr-sqrt42.8%
times-frac42.8%
sqrt-prod42.8%
sqrt-pow136.9%
metadata-eval36.9%
pow136.9%
*-commutative36.9%
sqrt-prod29.3%
sqrt-pow137.7%
metadata-eval37.7%
pow137.7%
Applied egg-rr49.8%
if 2.0000000000000001e299 < (*.f64 l l) Initial program 44.5%
associate-/r*44.5%
associate-/l/44.5%
associate-*l/44.5%
associate-/r/44.5%
+-commutative44.5%
unpow244.5%
sqr-neg44.5%
distribute-frac-neg244.5%
distribute-frac-neg244.5%
unpow244.5%
+-rgt-identity44.5%
metadata-eval44.5%
associate--l+44.5%
+-commutative44.5%
associate--l+44.5%
Simplified44.5%
associate-*r*52.4%
unpow252.4%
times-frac67.4%
div-inv67.4%
frac-times67.4%
metadata-eval67.4%
Applied egg-rr67.4%
Taylor expanded in k around inf 63.8%
times-frac72.5%
*-commutative72.5%
Simplified72.5%
associate-*r/71.0%
associate-*r*71.0%
div-inv71.0%
pow-prod-down74.7%
pow-flip74.8%
metadata-eval74.8%
Applied egg-rr74.8%
associate-/l*76.3%
associate-/r/81.7%
*-commutative81.7%
associate-*l*83.4%
*-commutative83.4%
associate-*l*83.4%
*-commutative83.4%
Simplified83.4%
Final simplification57.7%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= (* l l) 1e-314)
(/ 2.0 (pow (* (/ (pow k 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+105)
(/
2.0
(/ (* (pow l -2.0) (pow (* (sqrt t_m) (* k (sin k))) 2.0)) (cos k)))
(*
(/ l k)
(* 2.0 (* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k)))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / pow(((pow(k, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+105) {
tmp = 2.0 / ((pow(l, -2.0) * pow((sqrt(t_m) * (k * sin(k))), 2.0)) / cos(k));
} else {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if ((l * l) <= 1d-314) then
tmp = 2.0d0 / ((((k ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+105) then
tmp = 2.0d0 / (((l ** (-2.0d0)) * ((sqrt(t_m) * (k * sin(k))) ** 2.0d0)) / cos(k))
else
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
end if
code = t_s * tmp
end function
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) {
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / Math.pow(((Math.pow(k, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+105) {
tmp = 2.0 / ((Math.pow(l, -2.0) * Math.pow((Math.sqrt(t_m) * (k * Math.sin(k))), 2.0)) / Math.cos(k));
} else {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if (l * l) <= 1e-314: tmp = 2.0 / math.pow(((math.pow(k, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+105: tmp = 2.0 / ((math.pow(l, -2.0) * math.pow((math.sqrt(t_m) * (k * math.sin(k))), 2.0)) / math.cos(k)) else: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-314) tmp = Float64(2.0 / (Float64(Float64((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+105) tmp = Float64(2.0 / Float64(Float64((l ^ -2.0) * (Float64(sqrt(t_m) * Float64(k * sin(k))) ^ 2.0)) / cos(k))); else tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if ((l * l) <= 1e-314) tmp = 2.0 / ((((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+105) tmp = 2.0 / (((l ^ -2.0) * ((sqrt(t_m) * (k * sin(k))) ^ 2.0)) / cos(k)); else tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 1e-314], N[(2.0 / N[Power[N[(N[(N[Power[k, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+105], N[(2.0 / N[(N[(N[Power[l, -2.0], $MachinePrecision] * N[Power[N[(N[Sqrt[t$95$m], $MachinePrecision] * N[(k * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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 10^{-314}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+105}:\\
\;\;\;\;\frac{2}{\frac{{\ell}^{-2} \cdot {\left(\sqrt{t\_m} \cdot \left(k \cdot \sin k\right)\right)}^{2}}{\cos k}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 9.9999999996e-315Initial program 29.6%
Taylor expanded in t around 0 68.9%
associate-/l*68.9%
Simplified68.9%
add-sqr-sqrt30.1%
pow230.1%
Applied egg-rr32.7%
associate-/r*32.7%
Simplified32.7%
Taylor expanded in k around 0 52.9%
if 9.9999999996e-315 < (*.f64 l l) < 1.9999999999999999e105Initial program 36.0%
Taylor expanded in t around 0 91.5%
add-sqr-sqrt48.0%
times-frac48.0%
sqrt-prod47.9%
sqrt-pow143.8%
metadata-eval43.8%
pow143.8%
*-commutative43.8%
sqrt-prod36.1%
sqrt-pow142.8%
metadata-eval42.8%
pow142.8%
Applied egg-rr50.9%
Taylor expanded in k around inf 50.9%
associate-*r*50.9%
associate-*r/50.9%
div-inv50.9%
associate-*r*50.9%
*-commutative50.9%
pow-flip50.9%
metadata-eval50.9%
associate-*r*50.9%
*-commutative50.9%
Applied egg-rr50.9%
*-commutative50.9%
associate-*r*50.0%
unpow250.0%
Simplified50.0%
if 1.9999999999999999e105 < (*.f64 l l) Initial program 45.9%
associate-/r*45.8%
associate-/l/45.8%
associate-*l/45.8%
associate-/r/45.8%
+-commutative45.8%
unpow245.8%
sqr-neg45.8%
distribute-frac-neg245.8%
distribute-frac-neg245.8%
unpow245.8%
+-rgt-identity45.8%
metadata-eval45.8%
associate--l+45.8%
+-commutative45.8%
associate--l+45.8%
Simplified46.9%
associate-*r*51.5%
unpow251.5%
times-frac64.4%
div-inv64.4%
frac-times64.4%
metadata-eval64.4%
Applied egg-rr64.4%
Taylor expanded in k around inf 67.7%
times-frac72.7%
*-commutative72.7%
Simplified72.7%
associate-*r/71.9%
associate-*r*71.9%
div-inv72.0%
pow-prod-down75.1%
pow-flip75.9%
metadata-eval75.9%
Applied egg-rr75.9%
associate-/l*76.7%
associate-/r/81.9%
*-commutative81.9%
associate-*l*82.9%
*-commutative82.9%
associate-*l*82.9%
*-commutative82.9%
Simplified82.9%
Final simplification62.9%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= (* l l) 1e-314)
(/ 2.0 (pow (* (/ (pow k 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+105)
(/
2.0
(/ (/ (pow (* k (* (sqrt t_m) (sin k))) 2.0) (pow l 2.0)) (cos k)))
(*
(/ l k)
(* 2.0 (* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k)))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / pow(((pow(k, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+105) {
tmp = 2.0 / ((pow((k * (sqrt(t_m) * sin(k))), 2.0) / pow(l, 2.0)) / cos(k));
} else {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if ((l * l) <= 1d-314) then
tmp = 2.0d0 / ((((k ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+105) then
tmp = 2.0d0 / ((((k * (sqrt(t_m) * sin(k))) ** 2.0d0) / (l ** 2.0d0)) / cos(k))
else
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
end if
code = t_s * tmp
end function
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) {
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / Math.pow(((Math.pow(k, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+105) {
tmp = 2.0 / ((Math.pow((k * (Math.sqrt(t_m) * Math.sin(k))), 2.0) / Math.pow(l, 2.0)) / Math.cos(k));
} else {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if (l * l) <= 1e-314: tmp = 2.0 / math.pow(((math.pow(k, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+105: tmp = 2.0 / ((math.pow((k * (math.sqrt(t_m) * math.sin(k))), 2.0) / math.pow(l, 2.0)) / math.cos(k)) else: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-314) tmp = Float64(2.0 / (Float64(Float64((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+105) tmp = Float64(2.0 / Float64(Float64((Float64(k * Float64(sqrt(t_m) * sin(k))) ^ 2.0) / (l ^ 2.0)) / cos(k))); else tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if ((l * l) <= 1e-314) tmp = 2.0 / ((((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+105) tmp = 2.0 / ((((k * (sqrt(t_m) * sin(k))) ^ 2.0) / (l ^ 2.0)) / cos(k)); else tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 1e-314], N[(2.0 / N[Power[N[(N[(N[Power[k, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+105], N[(2.0 / N[(N[(N[Power[N[(k * N[(N[Sqrt[t$95$m], $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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 10^{-314}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+105}:\\
\;\;\;\;\frac{2}{\frac{\frac{{\left(k \cdot \left(\sqrt{t\_m} \cdot \sin k\right)\right)}^{2}}{{\ell}^{2}}}{\cos k}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 9.9999999996e-315Initial program 29.6%
Taylor expanded in t around 0 68.9%
associate-/l*68.9%
Simplified68.9%
add-sqr-sqrt30.1%
pow230.1%
Applied egg-rr32.7%
associate-/r*32.7%
Simplified32.7%
Taylor expanded in k around 0 52.9%
if 9.9999999996e-315 < (*.f64 l l) < 1.9999999999999999e105Initial program 36.0%
Taylor expanded in t around 0 91.5%
associate-/l*94.3%
Simplified94.3%
associate-*r/91.5%
associate-/r*91.6%
add-sqr-sqrt47.9%
pow247.9%
sqrt-prod47.9%
sqrt-pow150.0%
metadata-eval50.0%
pow150.0%
*-commutative50.0%
sqrt-prod50.0%
sqrt-pow150.0%
metadata-eval50.0%
pow150.0%
Applied egg-rr50.0%
if 1.9999999999999999e105 < (*.f64 l l) Initial program 45.9%
associate-/r*45.8%
associate-/l/45.8%
associate-*l/45.8%
associate-/r/45.8%
+-commutative45.8%
unpow245.8%
sqr-neg45.8%
distribute-frac-neg245.8%
distribute-frac-neg245.8%
unpow245.8%
+-rgt-identity45.8%
metadata-eval45.8%
associate--l+45.8%
+-commutative45.8%
associate--l+45.8%
Simplified46.9%
associate-*r*51.5%
unpow251.5%
times-frac64.4%
div-inv64.4%
frac-times64.4%
metadata-eval64.4%
Applied egg-rr64.4%
Taylor expanded in k around inf 67.7%
times-frac72.7%
*-commutative72.7%
Simplified72.7%
associate-*r/71.9%
associate-*r*71.9%
div-inv72.0%
pow-prod-down75.1%
pow-flip75.9%
metadata-eval75.9%
Applied egg-rr75.9%
associate-/l*76.7%
associate-/r/81.9%
*-commutative81.9%
associate-*l*82.9%
*-commutative82.9%
associate-*l*82.9%
*-commutative82.9%
Simplified82.9%
Final simplification62.9%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= (* l l) 1e-314)
(/ 2.0 (pow (* (/ (pow k 2.0) l) (sqrt t_m)) 2.0))
(if (<= (* l l) 2e+105)
(/ 2.0 (* (* (pow k 2.0) t_m) (* (tan k) (/ (sin k) (pow l 2.0)))))
(*
(/ l k)
(* 2.0 (* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k)))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / pow(((pow(k, 2.0) / l) * sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+105) {
tmp = 2.0 / ((pow(k, 2.0) * t_m) * (tan(k) * (sin(k) / pow(l, 2.0))));
} else {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if ((l * l) <= 1d-314) then
tmp = 2.0d0 / ((((k ** 2.0d0) / l) * sqrt(t_m)) ** 2.0d0)
else if ((l * l) <= 2d+105) then
tmp = 2.0d0 / (((k ** 2.0d0) * t_m) * (tan(k) * (sin(k) / (l ** 2.0d0))))
else
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
end if
code = t_s * tmp
end function
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) {
double tmp;
if ((l * l) <= 1e-314) {
tmp = 2.0 / Math.pow(((Math.pow(k, 2.0) / l) * Math.sqrt(t_m)), 2.0);
} else if ((l * l) <= 2e+105) {
tmp = 2.0 / ((Math.pow(k, 2.0) * t_m) * (Math.tan(k) * (Math.sin(k) / Math.pow(l, 2.0))));
} else {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if (l * l) <= 1e-314: tmp = 2.0 / math.pow(((math.pow(k, 2.0) / l) * math.sqrt(t_m)), 2.0) elif (l * l) <= 2e+105: tmp = 2.0 / ((math.pow(k, 2.0) * t_m) * (math.tan(k) * (math.sin(k) / math.pow(l, 2.0)))) else: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-314) tmp = Float64(2.0 / (Float64(Float64((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); elseif (Float64(l * l) <= 2e+105) tmp = Float64(2.0 / Float64(Float64((k ^ 2.0) * t_m) * Float64(tan(k) * Float64(sin(k) / (l ^ 2.0))))); else tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if ((l * l) <= 1e-314) tmp = 2.0 / ((((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0); elseif ((l * l) <= 2e+105) tmp = 2.0 / (((k ^ 2.0) * t_m) * (tan(k) * (sin(k) / (l ^ 2.0)))); else tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 1e-314], N[(2.0 / N[Power[N[(N[(N[Power[k, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+105], N[(2.0 / N[(N[(N[Power[k, 2.0], $MachinePrecision] * t$95$m), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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 10^{-314}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+105}:\\
\;\;\;\;\frac{2}{\left({k}^{2} \cdot t\_m\right) \cdot \left(\tan k \cdot \frac{\sin k}{{\ell}^{2}}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 9.9999999996e-315Initial program 29.6%
Taylor expanded in t around 0 68.9%
associate-/l*68.9%
Simplified68.9%
add-sqr-sqrt30.1%
pow230.1%
Applied egg-rr32.7%
associate-/r*32.7%
Simplified32.7%
Taylor expanded in k around 0 52.9%
if 9.9999999996e-315 < (*.f64 l l) < 1.9999999999999999e105Initial program 36.0%
Taylor expanded in t around 0 91.5%
associate-/l*94.3%
Simplified94.3%
pow194.3%
*-commutative94.3%
associate-/l*94.3%
unpow294.3%
*-commutative94.3%
times-frac94.2%
tan-quot94.3%
Applied egg-rr94.3%
unpow194.3%
*-commutative94.3%
associate-*r*91.6%
Simplified91.6%
if 1.9999999999999999e105 < (*.f64 l l) Initial program 45.9%
associate-/r*45.8%
associate-/l/45.8%
associate-*l/45.8%
associate-/r/45.8%
+-commutative45.8%
unpow245.8%
sqr-neg45.8%
distribute-frac-neg245.8%
distribute-frac-neg245.8%
unpow245.8%
+-rgt-identity45.8%
metadata-eval45.8%
associate--l+45.8%
+-commutative45.8%
associate--l+45.8%
Simplified46.9%
associate-*r*51.5%
unpow251.5%
times-frac64.4%
div-inv64.4%
frac-times64.4%
metadata-eval64.4%
Applied egg-rr64.4%
Taylor expanded in k around inf 67.7%
times-frac72.7%
*-commutative72.7%
Simplified72.7%
associate-*r/71.9%
associate-*r*71.9%
div-inv72.0%
pow-prod-down75.1%
pow-flip75.9%
metadata-eval75.9%
Applied egg-rr75.9%
associate-/l*76.7%
associate-/r/81.9%
*-commutative81.9%
associate-*l*82.9%
*-commutative82.9%
associate-*l*82.9%
*-commutative82.9%
Simplified82.9%
Final simplification77.9%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.5e-162)
(/ 2.0 (pow (* k (* (sqrt (/ t_m (cos k))) (/ (sin k) l))) 2.0))
(if (<= t_m 1.35e+219)
(*
(/ l k)
(* 2.0 (* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k)))))
(* 2.0 (pow (/ l (* (pow k 2.0) (sqrt t_m))) 2.0))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.5e-162) {
tmp = 2.0 / pow((k * (sqrt((t_m / cos(k))) * (sin(k) / l))), 2.0);
} else if (t_m <= 1.35e+219) {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
} else {
tmp = 2.0 * pow((l / (pow(k, 2.0) * sqrt(t_m))), 2.0);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.5d-162) then
tmp = 2.0d0 / ((k * (sqrt((t_m / cos(k))) * (sin(k) / l))) ** 2.0d0)
else if (t_m <= 1.35d+219) then
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
else
tmp = 2.0d0 * ((l / ((k ** 2.0d0) * sqrt(t_m))) ** 2.0d0)
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 1.5e-162) {
tmp = 2.0 / Math.pow((k * (Math.sqrt((t_m / Math.cos(k))) * (Math.sin(k) / l))), 2.0);
} else if (t_m <= 1.35e+219) {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
} else {
tmp = 2.0 * Math.pow((l / (Math.pow(k, 2.0) * Math.sqrt(t_m))), 2.0);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.5e-162: tmp = 2.0 / math.pow((k * (math.sqrt((t_m / math.cos(k))) * (math.sin(k) / l))), 2.0) elif t_m <= 1.35e+219: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) else: tmp = 2.0 * math.pow((l / (math.pow(k, 2.0) * math.sqrt(t_m))), 2.0) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.5e-162) tmp = Float64(2.0 / (Float64(k * Float64(sqrt(Float64(t_m / cos(k))) * Float64(sin(k) / l))) ^ 2.0)); elseif (t_m <= 1.35e+219) tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); else tmp = Float64(2.0 * (Float64(l / Float64((k ^ 2.0) * sqrt(t_m))) ^ 2.0)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.5e-162) tmp = 2.0 / ((k * (sqrt((t_m / cos(k))) * (sin(k) / l))) ^ 2.0); elseif (t_m <= 1.35e+219) tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); else tmp = 2.0 * ((l / ((k ^ 2.0) * sqrt(t_m))) ^ 2.0); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.5e-162], N[(2.0 / N[Power[N[(k * N[(N[Sqrt[N[(t$95$m / N[Cos[k], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 1.35e+219], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[Power[N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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.5 \cdot 10^{-162}:\\
\;\;\;\;\frac{2}{{\left(k \cdot \left(\sqrt{\frac{t\_m}{\cos k}} \cdot \frac{\sin k}{\ell}\right)\right)}^{2}}\\
\mathbf{elif}\;t\_m \leq 1.35 \cdot 10^{+219}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{{k}^{2} \cdot \sqrt{t\_m}}\right)}^{2}\\
\end{array}
\end{array}
if t < 1.49999999999999999e-162Initial program 35.2%
Taylor expanded in t around 0 71.2%
associate-/l*72.4%
Simplified72.4%
add-sqr-sqrt24.1%
pow224.1%
Applied egg-rr11.0%
associate-/r*11.0%
Simplified11.0%
Taylor expanded in k around inf 36.0%
if 1.49999999999999999e-162 < t < 1.3499999999999999e219Initial program 50.8%
associate-/r*50.8%
associate-/l/50.8%
associate-*l/53.3%
associate-/r/53.3%
+-commutative53.3%
unpow253.3%
sqr-neg53.3%
distribute-frac-neg253.3%
distribute-frac-neg253.3%
unpow253.3%
+-rgt-identity53.3%
metadata-eval53.3%
associate--l+53.3%
+-commutative53.3%
associate--l+53.3%
Simplified63.5%
associate-*r*71.8%
unpow271.8%
times-frac76.7%
div-inv76.7%
frac-times76.8%
metadata-eval76.8%
Applied egg-rr76.8%
Taylor expanded in k around inf 90.0%
times-frac91.1%
*-commutative91.1%
Simplified91.1%
associate-*r/91.0%
associate-*r*91.0%
div-inv91.0%
pow-prod-down92.1%
pow-flip92.8%
metadata-eval92.8%
Applied egg-rr92.8%
associate-/l*93.0%
associate-/r/94.2%
*-commutative94.2%
associate-*l*95.3%
*-commutative95.3%
associate-*l*95.3%
*-commutative95.3%
Simplified95.3%
if 1.3499999999999999e219 < t Initial program 0.0%
associate-/r*0.0%
associate-/l/0.0%
associate-*l/0.0%
associate-/r/0.0%
+-commutative0.0%
unpow20.0%
sqr-neg0.0%
distribute-frac-neg20.0%
distribute-frac-neg20.0%
unpow20.0%
+-rgt-identity0.0%
metadata-eval0.0%
associate--l+0.0%
+-commutative0.0%
associate--l+0.0%
Simplified11.1%
Taylor expanded in k around 0 73.3%
add-sqr-sqrt73.3%
sqrt-div73.3%
sqrt-pow162.3%
metadata-eval62.3%
pow162.3%
sqrt-prod62.3%
sqrt-pow162.3%
metadata-eval62.3%
sqrt-div62.3%
sqrt-pow179.0%
metadata-eval79.0%
pow179.0%
sqrt-prod79.0%
sqrt-pow184.4%
metadata-eval84.4%
Applied egg-rr84.4%
unpow284.4%
Simplified84.4%
Final simplification58.9%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.5e-162)
(/ 2.0 (pow (* (sqrt (/ t_m (cos k))) (* k (/ (sin k) l))) 2.0))
(if (<= t_m 1.5e+219)
(*
(/ l k)
(* 2.0 (* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k)))))
(* 2.0 (pow (/ l (* (pow k 2.0) (sqrt t_m))) 2.0))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.5e-162) {
tmp = 2.0 / pow((sqrt((t_m / cos(k))) * (k * (sin(k) / l))), 2.0);
} else if (t_m <= 1.5e+219) {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
} else {
tmp = 2.0 * pow((l / (pow(k, 2.0) * sqrt(t_m))), 2.0);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.5d-162) then
tmp = 2.0d0 / ((sqrt((t_m / cos(k))) * (k * (sin(k) / l))) ** 2.0d0)
else if (t_m <= 1.5d+219) then
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
else
tmp = 2.0d0 * ((l / ((k ** 2.0d0) * sqrt(t_m))) ** 2.0d0)
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 1.5e-162) {
tmp = 2.0 / Math.pow((Math.sqrt((t_m / Math.cos(k))) * (k * (Math.sin(k) / l))), 2.0);
} else if (t_m <= 1.5e+219) {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
} else {
tmp = 2.0 * Math.pow((l / (Math.pow(k, 2.0) * Math.sqrt(t_m))), 2.0);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.5e-162: tmp = 2.0 / math.pow((math.sqrt((t_m / math.cos(k))) * (k * (math.sin(k) / l))), 2.0) elif t_m <= 1.5e+219: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) else: tmp = 2.0 * math.pow((l / (math.pow(k, 2.0) * math.sqrt(t_m))), 2.0) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.5e-162) tmp = Float64(2.0 / (Float64(sqrt(Float64(t_m / cos(k))) * Float64(k * Float64(sin(k) / l))) ^ 2.0)); elseif (t_m <= 1.5e+219) tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); else tmp = Float64(2.0 * (Float64(l / Float64((k ^ 2.0) * sqrt(t_m))) ^ 2.0)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.5e-162) tmp = 2.0 / ((sqrt((t_m / cos(k))) * (k * (sin(k) / l))) ^ 2.0); elseif (t_m <= 1.5e+219) tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); else tmp = 2.0 * ((l / ((k ^ 2.0) * sqrt(t_m))) ^ 2.0); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.5e-162], N[(2.0 / N[Power[N[(N[Sqrt[N[(t$95$m / N[Cos[k], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(k * N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 1.5e+219], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[Power[N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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.5 \cdot 10^{-162}:\\
\;\;\;\;\frac{2}{{\left(\sqrt{\frac{t\_m}{\cos k}} \cdot \left(k \cdot \frac{\sin k}{\ell}\right)\right)}^{2}}\\
\mathbf{elif}\;t\_m \leq 1.5 \cdot 10^{+219}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{{k}^{2} \cdot \sqrt{t\_m}}\right)}^{2}\\
\end{array}
\end{array}
if t < 1.49999999999999999e-162Initial program 35.2%
Taylor expanded in t around 0 71.2%
associate-/l*72.4%
Simplified72.4%
add-sqr-sqrt24.1%
pow224.1%
Applied egg-rr11.0%
associate-/r*11.0%
Simplified11.0%
Taylor expanded in k around inf 36.7%
associate-/l*36.7%
Simplified36.7%
if 1.49999999999999999e-162 < t < 1.4999999999999999e219Initial program 50.8%
associate-/r*50.8%
associate-/l/50.8%
associate-*l/53.3%
associate-/r/53.3%
+-commutative53.3%
unpow253.3%
sqr-neg53.3%
distribute-frac-neg253.3%
distribute-frac-neg253.3%
unpow253.3%
+-rgt-identity53.3%
metadata-eval53.3%
associate--l+53.3%
+-commutative53.3%
associate--l+53.3%
Simplified63.5%
associate-*r*71.8%
unpow271.8%
times-frac76.7%
div-inv76.7%
frac-times76.8%
metadata-eval76.8%
Applied egg-rr76.8%
Taylor expanded in k around inf 90.0%
times-frac91.1%
*-commutative91.1%
Simplified91.1%
associate-*r/91.0%
associate-*r*91.0%
div-inv91.0%
pow-prod-down92.1%
pow-flip92.8%
metadata-eval92.8%
Applied egg-rr92.8%
associate-/l*93.0%
associate-/r/94.2%
*-commutative94.2%
associate-*l*95.3%
*-commutative95.3%
associate-*l*95.3%
*-commutative95.3%
Simplified95.3%
if 1.4999999999999999e219 < t Initial program 0.0%
associate-/r*0.0%
associate-/l/0.0%
associate-*l/0.0%
associate-/r/0.0%
+-commutative0.0%
unpow20.0%
sqr-neg0.0%
distribute-frac-neg20.0%
distribute-frac-neg20.0%
unpow20.0%
+-rgt-identity0.0%
metadata-eval0.0%
associate--l+0.0%
+-commutative0.0%
associate--l+0.0%
Simplified11.1%
Taylor expanded in k around 0 73.3%
add-sqr-sqrt73.3%
sqrt-div73.3%
sqrt-pow162.3%
metadata-eval62.3%
pow162.3%
sqrt-prod62.3%
sqrt-pow162.3%
metadata-eval62.3%
sqrt-div62.3%
sqrt-pow179.0%
metadata-eval79.0%
pow179.0%
sqrt-prod79.0%
sqrt-pow184.4%
metadata-eval84.4%
Applied egg-rr84.4%
unpow284.4%
Simplified84.4%
Final simplification59.2%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.5e-162)
(/ 2.0 (pow (* (/ (* k (sin k)) l) (sqrt (/ t_m (cos k)))) 2.0))
(if (<= t_m 2e+220)
(*
(/ l k)
(* 2.0 (* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k)))))
(* 2.0 (pow (/ l (* (pow k 2.0) (sqrt t_m))) 2.0))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.5e-162) {
tmp = 2.0 / pow((((k * sin(k)) / l) * sqrt((t_m / cos(k)))), 2.0);
} else if (t_m <= 2e+220) {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
} else {
tmp = 2.0 * pow((l / (pow(k, 2.0) * sqrt(t_m))), 2.0);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.5d-162) then
tmp = 2.0d0 / ((((k * sin(k)) / l) * sqrt((t_m / cos(k)))) ** 2.0d0)
else if (t_m <= 2d+220) then
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
else
tmp = 2.0d0 * ((l / ((k ** 2.0d0) * sqrt(t_m))) ** 2.0d0)
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 1.5e-162) {
tmp = 2.0 / Math.pow((((k * Math.sin(k)) / l) * Math.sqrt((t_m / Math.cos(k)))), 2.0);
} else if (t_m <= 2e+220) {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
} else {
tmp = 2.0 * Math.pow((l / (Math.pow(k, 2.0) * Math.sqrt(t_m))), 2.0);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.5e-162: tmp = 2.0 / math.pow((((k * math.sin(k)) / l) * math.sqrt((t_m / math.cos(k)))), 2.0) elif t_m <= 2e+220: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) else: tmp = 2.0 * math.pow((l / (math.pow(k, 2.0) * math.sqrt(t_m))), 2.0) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.5e-162) tmp = Float64(2.0 / (Float64(Float64(Float64(k * sin(k)) / l) * sqrt(Float64(t_m / cos(k)))) ^ 2.0)); elseif (t_m <= 2e+220) tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); else tmp = Float64(2.0 * (Float64(l / Float64((k ^ 2.0) * sqrt(t_m))) ^ 2.0)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.5e-162) tmp = 2.0 / ((((k * sin(k)) / l) * sqrt((t_m / cos(k)))) ^ 2.0); elseif (t_m <= 2e+220) tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); else tmp = 2.0 * ((l / ((k ^ 2.0) * sqrt(t_m))) ^ 2.0); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.5e-162], N[(2.0 / N[Power[N[(N[(N[(k * N[Sin[k], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 2e+220], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[Power[N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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.5 \cdot 10^{-162}:\\
\;\;\;\;\frac{2}{{\left(\frac{k \cdot \sin k}{\ell} \cdot \sqrt{\frac{t\_m}{\cos k}}\right)}^{2}}\\
\mathbf{elif}\;t\_m \leq 2 \cdot 10^{+220}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{{k}^{2} \cdot \sqrt{t\_m}}\right)}^{2}\\
\end{array}
\end{array}
if t < 1.49999999999999999e-162Initial program 35.2%
Taylor expanded in t around 0 71.2%
associate-/l*72.4%
Simplified72.4%
add-sqr-sqrt24.1%
pow224.1%
Applied egg-rr11.0%
associate-/r*11.0%
Simplified11.0%
Taylor expanded in k around inf 36.7%
if 1.49999999999999999e-162 < t < 2e220Initial program 50.8%
associate-/r*50.8%
associate-/l/50.8%
associate-*l/53.3%
associate-/r/53.3%
+-commutative53.3%
unpow253.3%
sqr-neg53.3%
distribute-frac-neg253.3%
distribute-frac-neg253.3%
unpow253.3%
+-rgt-identity53.3%
metadata-eval53.3%
associate--l+53.3%
+-commutative53.3%
associate--l+53.3%
Simplified63.5%
associate-*r*71.8%
unpow271.8%
times-frac76.7%
div-inv76.7%
frac-times76.8%
metadata-eval76.8%
Applied egg-rr76.8%
Taylor expanded in k around inf 90.0%
times-frac91.1%
*-commutative91.1%
Simplified91.1%
associate-*r/91.0%
associate-*r*91.0%
div-inv91.0%
pow-prod-down92.1%
pow-flip92.8%
metadata-eval92.8%
Applied egg-rr92.8%
associate-/l*93.0%
associate-/r/94.2%
*-commutative94.2%
associate-*l*95.3%
*-commutative95.3%
associate-*l*95.3%
*-commutative95.3%
Simplified95.3%
if 2e220 < t Initial program 0.0%
associate-/r*0.0%
associate-/l/0.0%
associate-*l/0.0%
associate-/r/0.0%
+-commutative0.0%
unpow20.0%
sqr-neg0.0%
distribute-frac-neg20.0%
distribute-frac-neg20.0%
unpow20.0%
+-rgt-identity0.0%
metadata-eval0.0%
associate--l+0.0%
+-commutative0.0%
associate--l+0.0%
Simplified11.1%
Taylor expanded in k around 0 73.3%
add-sqr-sqrt73.3%
sqrt-div73.3%
sqrt-pow162.3%
metadata-eval62.3%
pow162.3%
sqrt-prod62.3%
sqrt-pow162.3%
metadata-eval62.3%
sqrt-div62.3%
sqrt-pow179.0%
metadata-eval79.0%
pow179.0%
sqrt-prod79.0%
sqrt-pow184.4%
metadata-eval84.4%
Applied egg-rr84.4%
unpow284.4%
Simplified84.4%
Final simplification59.2%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.5e-162)
(* 2.0 (pow (/ (* l (sqrt (/ 1.0 t_m))) (pow k 2.0)) 2.0))
(if (<= t_m 2e+219)
(*
(* 2.0 (* (cos k) (* (/ l k) (pow (* t_m (sin k)) -2.0))))
(/ l (/ k t_m)))
(* 2.0 (pow (/ l (* (pow k 2.0) (sqrt t_m))) 2.0))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.5e-162) {
tmp = 2.0 * pow(((l * sqrt((1.0 / t_m))) / pow(k, 2.0)), 2.0);
} else if (t_m <= 2e+219) {
tmp = (2.0 * (cos(k) * ((l / k) * pow((t_m * sin(k)), -2.0)))) * (l / (k / t_m));
} else {
tmp = 2.0 * pow((l / (pow(k, 2.0) * sqrt(t_m))), 2.0);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.5d-162) then
tmp = 2.0d0 * (((l * sqrt((1.0d0 / t_m))) / (k ** 2.0d0)) ** 2.0d0)
else if (t_m <= 2d+219) then
tmp = (2.0d0 * (cos(k) * ((l / k) * ((t_m * sin(k)) ** (-2.0d0))))) * (l / (k / t_m))
else
tmp = 2.0d0 * ((l / ((k ** 2.0d0) * sqrt(t_m))) ** 2.0d0)
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 1.5e-162) {
tmp = 2.0 * Math.pow(((l * Math.sqrt((1.0 / t_m))) / Math.pow(k, 2.0)), 2.0);
} else if (t_m <= 2e+219) {
tmp = (2.0 * (Math.cos(k) * ((l / k) * Math.pow((t_m * Math.sin(k)), -2.0)))) * (l / (k / t_m));
} else {
tmp = 2.0 * Math.pow((l / (Math.pow(k, 2.0) * Math.sqrt(t_m))), 2.0);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.5e-162: tmp = 2.0 * math.pow(((l * math.sqrt((1.0 / t_m))) / math.pow(k, 2.0)), 2.0) elif t_m <= 2e+219: tmp = (2.0 * (math.cos(k) * ((l / k) * math.pow((t_m * math.sin(k)), -2.0)))) * (l / (k / t_m)) else: tmp = 2.0 * math.pow((l / (math.pow(k, 2.0) * math.sqrt(t_m))), 2.0) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.5e-162) tmp = Float64(2.0 * (Float64(Float64(l * sqrt(Float64(1.0 / t_m))) / (k ^ 2.0)) ^ 2.0)); elseif (t_m <= 2e+219) tmp = Float64(Float64(2.0 * Float64(cos(k) * Float64(Float64(l / k) * (Float64(t_m * sin(k)) ^ -2.0)))) * Float64(l / Float64(k / t_m))); else tmp = Float64(2.0 * (Float64(l / Float64((k ^ 2.0) * sqrt(t_m))) ^ 2.0)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.5e-162) tmp = 2.0 * (((l * sqrt((1.0 / t_m))) / (k ^ 2.0)) ^ 2.0); elseif (t_m <= 2e+219) tmp = (2.0 * (cos(k) * ((l / k) * ((t_m * sin(k)) ^ -2.0)))) * (l / (k / t_m)); else tmp = 2.0 * ((l / ((k ^ 2.0) * sqrt(t_m))) ^ 2.0); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.5e-162], N[(2.0 * N[Power[N[(N[(l * N[Sqrt[N[(1.0 / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 2e+219], N[(N[(2.0 * N[(N[Cos[k], $MachinePrecision] * N[(N[(l / k), $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[(k / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[Power[N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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.5 \cdot 10^{-162}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell \cdot \sqrt{\frac{1}{t\_m}}}{{k}^{2}}\right)}^{2}\\
\mathbf{elif}\;t\_m \leq 2 \cdot 10^{+219}:\\
\;\;\;\;\left(2 \cdot \left(\cos k \cdot \left(\frac{\ell}{k} \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right)\right)\right) \cdot \frac{\ell}{\frac{k}{t\_m}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{{k}^{2} \cdot \sqrt{t\_m}}\right)}^{2}\\
\end{array}
\end{array}
if t < 1.49999999999999999e-162Initial program 35.2%
associate-/r*35.2%
associate-/l/35.2%
associate-*l/35.8%
associate-/r/35.1%
+-commutative35.1%
unpow235.1%
sqr-neg35.1%
distribute-frac-neg235.1%
distribute-frac-neg235.1%
unpow235.1%
+-rgt-identity35.1%
metadata-eval35.1%
associate--l+35.1%
+-commutative35.1%
associate--l+35.1%
Simplified41.7%
Taylor expanded in k around 0 54.9%
add-sqr-sqrt32.7%
sqrt-div7.4%
sqrt-pow14.8%
metadata-eval4.8%
pow14.8%
sqrt-prod4.8%
sqrt-pow14.8%
metadata-eval4.8%
sqrt-div4.8%
sqrt-pow19.4%
metadata-eval9.4%
pow19.4%
sqrt-prod9.4%
sqrt-pow110.0%
metadata-eval10.0%
Applied egg-rr10.0%
unpow210.0%
Simplified10.0%
Taylor expanded in l around 0 10.0%
associate-*l/10.0%
Simplified10.0%
if 1.49999999999999999e-162 < t < 1.99999999999999993e219Initial program 50.8%
associate-/r*50.8%
associate-/l/50.8%
associate-*l/53.3%
associate-/r/53.3%
+-commutative53.3%
unpow253.3%
sqr-neg53.3%
distribute-frac-neg253.3%
distribute-frac-neg253.3%
unpow253.3%
+-rgt-identity53.3%
metadata-eval53.3%
associate--l+53.3%
+-commutative53.3%
associate--l+53.3%
Simplified63.5%
associate-*r*71.8%
unpow271.8%
times-frac76.7%
div-inv76.7%
frac-times76.8%
metadata-eval76.8%
Applied egg-rr76.8%
Taylor expanded in k around inf 90.0%
times-frac91.1%
*-commutative91.1%
Simplified91.1%
associate-*l/90.0%
div-inv90.0%
pow-prod-down91.2%
pow-flip91.9%
metadata-eval91.9%
Applied egg-rr91.9%
associate-*l/93.0%
associate-*r*93.1%
*-commutative93.1%
associate-*l*93.0%
*-commutative93.0%
Simplified93.0%
if 1.99999999999999993e219 < t Initial program 0.0%
associate-/r*0.0%
associate-/l/0.0%
associate-*l/0.0%
associate-/r/0.0%
+-commutative0.0%
unpow20.0%
sqr-neg0.0%
distribute-frac-neg20.0%
distribute-frac-neg20.0%
unpow20.0%
+-rgt-identity0.0%
metadata-eval0.0%
associate--l+0.0%
+-commutative0.0%
associate--l+0.0%
Simplified11.1%
Taylor expanded in k around 0 73.3%
add-sqr-sqrt73.3%
sqrt-div73.3%
sqrt-pow162.3%
metadata-eval62.3%
pow162.3%
sqrt-prod62.3%
sqrt-pow162.3%
metadata-eval62.3%
sqrt-div62.3%
sqrt-pow179.0%
metadata-eval79.0%
pow179.0%
sqrt-prod79.0%
sqrt-pow184.4%
metadata-eval84.4%
Applied egg-rr84.4%
unpow284.4%
Simplified84.4%
Final simplification42.5%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.5e-162)
(* 2.0 (pow (/ (* l (sqrt (/ 1.0 t_m))) (pow k 2.0)) 2.0))
(if (<= t_m 1.6e+219)
(*
(/ l k)
(* 2.0 (* (* (cos k) (pow (* t_m (sin k)) -2.0)) (* t_m (/ l k)))))
(* 2.0 (pow (/ l (* (pow k 2.0) (sqrt t_m))) 2.0))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.5e-162) {
tmp = 2.0 * pow(((l * sqrt((1.0 / t_m))) / pow(k, 2.0)), 2.0);
} else if (t_m <= 1.6e+219) {
tmp = (l / k) * (2.0 * ((cos(k) * pow((t_m * sin(k)), -2.0)) * (t_m * (l / k))));
} else {
tmp = 2.0 * pow((l / (pow(k, 2.0) * sqrt(t_m))), 2.0);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.5d-162) then
tmp = 2.0d0 * (((l * sqrt((1.0d0 / t_m))) / (k ** 2.0d0)) ** 2.0d0)
else if (t_m <= 1.6d+219) then
tmp = (l / k) * (2.0d0 * ((cos(k) * ((t_m * sin(k)) ** (-2.0d0))) * (t_m * (l / k))))
else
tmp = 2.0d0 * ((l / ((k ** 2.0d0) * sqrt(t_m))) ** 2.0d0)
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 1.5e-162) {
tmp = 2.0 * Math.pow(((l * Math.sqrt((1.0 / t_m))) / Math.pow(k, 2.0)), 2.0);
} else if (t_m <= 1.6e+219) {
tmp = (l / k) * (2.0 * ((Math.cos(k) * Math.pow((t_m * Math.sin(k)), -2.0)) * (t_m * (l / k))));
} else {
tmp = 2.0 * Math.pow((l / (Math.pow(k, 2.0) * Math.sqrt(t_m))), 2.0);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.5e-162: tmp = 2.0 * math.pow(((l * math.sqrt((1.0 / t_m))) / math.pow(k, 2.0)), 2.0) elif t_m <= 1.6e+219: tmp = (l / k) * (2.0 * ((math.cos(k) * math.pow((t_m * math.sin(k)), -2.0)) * (t_m * (l / k)))) else: tmp = 2.0 * math.pow((l / (math.pow(k, 2.0) * math.sqrt(t_m))), 2.0) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.5e-162) tmp = Float64(2.0 * (Float64(Float64(l * sqrt(Float64(1.0 / t_m))) / (k ^ 2.0)) ^ 2.0)); elseif (t_m <= 1.6e+219) tmp = Float64(Float64(l / k) * Float64(2.0 * Float64(Float64(cos(k) * (Float64(t_m * sin(k)) ^ -2.0)) * Float64(t_m * Float64(l / k))))); else tmp = Float64(2.0 * (Float64(l / Float64((k ^ 2.0) * sqrt(t_m))) ^ 2.0)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.5e-162) tmp = 2.0 * (((l * sqrt((1.0 / t_m))) / (k ^ 2.0)) ^ 2.0); elseif (t_m <= 1.6e+219) tmp = (l / k) * (2.0 * ((cos(k) * ((t_m * sin(k)) ^ -2.0)) * (t_m * (l / k)))); else tmp = 2.0 * ((l / ((k ^ 2.0) * sqrt(t_m))) ^ 2.0); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.5e-162], N[(2.0 * N[Power[N[(N[(l * N[Sqrt[N[(1.0 / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 1.6e+219], N[(N[(l / k), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[Power[N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
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.5 \cdot 10^{-162}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell \cdot \sqrt{\frac{1}{t\_m}}}{{k}^{2}}\right)}^{2}\\
\mathbf{elif}\;t\_m \leq 1.6 \cdot 10^{+219}:\\
\;\;\;\;\frac{\ell}{k} \cdot \left(2 \cdot \left(\left(\cos k \cdot {\left(t\_m \cdot \sin k\right)}^{-2}\right) \cdot \left(t\_m \cdot \frac{\ell}{k}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{{k}^{2} \cdot \sqrt{t\_m}}\right)}^{2}\\
\end{array}
\end{array}
if t < 1.49999999999999999e-162Initial program 35.2%
associate-/r*35.2%
associate-/l/35.2%
associate-*l/35.8%
associate-/r/35.1%
+-commutative35.1%
unpow235.1%
sqr-neg35.1%
distribute-frac-neg235.1%
distribute-frac-neg235.1%
unpow235.1%
+-rgt-identity35.1%
metadata-eval35.1%
associate--l+35.1%
+-commutative35.1%
associate--l+35.1%
Simplified41.7%
Taylor expanded in k around 0 54.9%
add-sqr-sqrt32.7%
sqrt-div7.4%
sqrt-pow14.8%
metadata-eval4.8%
pow14.8%
sqrt-prod4.8%
sqrt-pow14.8%
metadata-eval4.8%
sqrt-div4.8%
sqrt-pow19.4%
metadata-eval9.4%
pow19.4%
sqrt-prod9.4%
sqrt-pow110.0%
metadata-eval10.0%
Applied egg-rr10.0%
unpow210.0%
Simplified10.0%
Taylor expanded in l around 0 10.0%
associate-*l/10.0%
Simplified10.0%
if 1.49999999999999999e-162 < t < 1.60000000000000013e219Initial program 50.8%
associate-/r*50.8%
associate-/l/50.8%
associate-*l/53.3%
associate-/r/53.3%
+-commutative53.3%
unpow253.3%
sqr-neg53.3%
distribute-frac-neg253.3%
distribute-frac-neg253.3%
unpow253.3%
+-rgt-identity53.3%
metadata-eval53.3%
associate--l+53.3%
+-commutative53.3%
associate--l+53.3%
Simplified63.5%
associate-*r*71.8%
unpow271.8%
times-frac76.7%
div-inv76.7%
frac-times76.8%
metadata-eval76.8%
Applied egg-rr76.8%
Taylor expanded in k around inf 90.0%
times-frac91.1%
*-commutative91.1%
Simplified91.1%
associate-*r/91.0%
associate-*r*91.0%
div-inv91.0%
pow-prod-down92.1%
pow-flip92.8%
metadata-eval92.8%
Applied egg-rr92.8%
associate-/l*93.0%
associate-/r/94.2%
*-commutative94.2%
associate-*l*95.3%
*-commutative95.3%
associate-*l*95.3%
*-commutative95.3%
Simplified95.3%
if 1.60000000000000013e219 < t Initial program 0.0%
associate-/r*0.0%
associate-/l/0.0%
associate-*l/0.0%
associate-/r/0.0%
+-commutative0.0%
unpow20.0%
sqr-neg0.0%
distribute-frac-neg20.0%
distribute-frac-neg20.0%
unpow20.0%
+-rgt-identity0.0%
metadata-eval0.0%
associate--l+0.0%
+-commutative0.0%
associate--l+0.0%
Simplified11.1%
Taylor expanded in k around 0 73.3%
add-sqr-sqrt73.3%
sqrt-div73.3%
sqrt-pow162.3%
metadata-eval62.3%
pow162.3%
sqrt-prod62.3%
sqrt-pow162.3%
metadata-eval62.3%
sqrt-div62.3%
sqrt-pow179.0%
metadata-eval79.0%
pow179.0%
sqrt-prod79.0%
sqrt-pow184.4%
metadata-eval84.4%
Applied egg-rr84.4%
unpow284.4%
Simplified84.4%
Final simplification43.2%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 6.5e-23)
(* 2.0 (pow (/ l (* (pow k 2.0) (sqrt t_m))) 2.0))
(*
l
(/
(fma -0.3333333333333333 (/ l k) (* 2.0 (/ l (pow k 3.0))))
(* k t_m))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 6.5e-23) {
tmp = 2.0 * pow((l / (pow(k, 2.0) * sqrt(t_m))), 2.0);
} else {
tmp = l * (fma(-0.3333333333333333, (l / k), (2.0 * (l / pow(k, 3.0)))) / (k * t_m));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 6.5e-23) tmp = Float64(2.0 * (Float64(l / Float64((k ^ 2.0) * sqrt(t_m))) ^ 2.0)); else tmp = Float64(l * Float64(fma(-0.3333333333333333, Float64(l / k), Float64(2.0 * Float64(l / (k ^ 3.0)))) / Float64(k * t_m))); end return Float64(t_s * tmp) end
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_] := N[(t$95$s * If[LessEqual[k, 6.5e-23], N[(2.0 * N[Power[N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(l * N[(N[(-0.3333333333333333 * N[(l / k), $MachinePrecision] + N[(2.0 * N[(l / N[Power[k, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 6.5 \cdot 10^{-23}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{{k}^{2} \cdot \sqrt{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \frac{\mathsf{fma}\left(-0.3333333333333333, \frac{\ell}{k}, 2 \cdot \frac{\ell}{{k}^{3}}\right)}{k \cdot t\_m}\\
\end{array}
\end{array}
if k < 6.5e-23Initial program 40.8%
associate-/r*40.8%
associate-/l/40.8%
associate-*l/42.5%
associate-/r/42.5%
+-commutative42.5%
unpow242.5%
sqr-neg42.5%
distribute-frac-neg242.5%
distribute-frac-neg242.5%
unpow242.5%
+-rgt-identity42.5%
metadata-eval42.5%
associate--l+42.5%
+-commutative42.5%
associate--l+42.5%
Simplified51.3%
Taylor expanded in k around 0 67.3%
add-sqr-sqrt48.9%
sqrt-div39.0%
sqrt-pow125.2%
metadata-eval25.2%
pow125.2%
sqrt-prod25.2%
sqrt-pow125.2%
metadata-eval25.2%
sqrt-div25.2%
sqrt-pow145.3%
metadata-eval45.3%
pow145.3%
sqrt-prod45.2%
sqrt-pow146.3%
metadata-eval46.3%
Applied egg-rr46.3%
unpow246.3%
Simplified46.3%
if 6.5e-23 < k Initial program 31.4%
associate-/r*31.3%
associate-/l/31.3%
associate-*l/31.3%
associate-/r/30.1%
+-commutative30.1%
unpow230.1%
sqr-neg30.1%
distribute-frac-neg230.1%
distribute-frac-neg230.1%
unpow230.1%
+-rgt-identity30.1%
metadata-eval30.1%
associate--l+30.1%
+-commutative30.1%
associate--l+30.1%
Simplified36.5%
associate-*r*37.9%
unpow237.9%
times-frac51.8%
div-inv51.8%
frac-times51.8%
metadata-eval51.8%
Applied egg-rr51.8%
Taylor expanded in k around 0 44.5%
Taylor expanded in t around 0 54.9%
associate-/l*55.0%
fma-define55.0%
Simplified55.0%
Final simplification49.0%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 3.3e-19)
(/ 2.0 (pow (* (/ (pow k 2.0) l) (sqrt t_m)) 2.0))
(*
l
(/
(fma -0.3333333333333333 (/ l k) (* 2.0 (/ l (pow k 3.0))))
(* k t_m))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 3.3e-19) {
tmp = 2.0 / pow(((pow(k, 2.0) / l) * sqrt(t_m)), 2.0);
} else {
tmp = l * (fma(-0.3333333333333333, (l / k), (2.0 * (l / pow(k, 3.0)))) / (k * t_m));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 3.3e-19) tmp = Float64(2.0 / (Float64(Float64((k ^ 2.0) / l) * sqrt(t_m)) ^ 2.0)); else tmp = Float64(l * Float64(fma(-0.3333333333333333, Float64(l / k), Float64(2.0 * Float64(l / (k ^ 3.0)))) / Float64(k * t_m))); end return Float64(t_s * tmp) end
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_] := N[(t$95$s * If[LessEqual[k, 3.3e-19], N[(2.0 / N[Power[N[(N[(N[Power[k, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(l * N[(N[(-0.3333333333333333 * N[(l / k), $MachinePrecision] + N[(2.0 * N[(l / N[Power[k, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 3.3 \cdot 10^{-19}:\\
\;\;\;\;\frac{2}{{\left(\frac{{k}^{2}}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \frac{\mathsf{fma}\left(-0.3333333333333333, \frac{\ell}{k}, 2 \cdot \frac{\ell}{{k}^{3}}\right)}{k \cdot t\_m}\\
\end{array}
\end{array}
if k < 3.2999999999999998e-19Initial program 40.6%
Taylor expanded in t around 0 76.3%
associate-/l*77.7%
Simplified77.7%
add-sqr-sqrt40.3%
pow240.3%
Applied egg-rr42.8%
associate-/r*42.8%
Simplified42.8%
Taylor expanded in k around 0 46.6%
if 3.2999999999999998e-19 < k Initial program 31.8%
associate-/r*31.7%
associate-/l/31.7%
associate-*l/31.7%
associate-/r/30.5%
+-commutative30.5%
unpow230.5%
sqr-neg30.5%
distribute-frac-neg230.5%
distribute-frac-neg230.5%
unpow230.5%
+-rgt-identity30.5%
metadata-eval30.5%
associate--l+30.5%
+-commutative30.5%
associate--l+30.5%
Simplified36.9%
associate-*r*38.4%
unpow238.4%
times-frac52.4%
div-inv52.4%
frac-times52.4%
metadata-eval52.4%
Applied egg-rr52.4%
Taylor expanded in k around 0 43.8%
Taylor expanded in t around 0 54.4%
associate-/l*54.4%
fma-define54.4%
Simplified54.4%
Final simplification49.0%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 5e-24)
(* 2.0 (* l (/ l (* t_m (pow k 4.0)))))
(*
l
(/
(fma -0.3333333333333333 (/ l k) (* 2.0 (/ l (pow k 3.0))))
(* k t_m))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 5e-24) {
tmp = 2.0 * (l * (l / (t_m * pow(k, 4.0))));
} else {
tmp = l * (fma(-0.3333333333333333, (l / k), (2.0 * (l / pow(k, 3.0)))) / (k * t_m));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 5e-24) tmp = Float64(2.0 * Float64(l * Float64(l / Float64(t_m * (k ^ 4.0))))); else tmp = Float64(l * Float64(fma(-0.3333333333333333, Float64(l / k), Float64(2.0 * Float64(l / (k ^ 3.0)))) / Float64(k * t_m))); end return Float64(t_s * tmp) end
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_] := N[(t$95$s * If[LessEqual[k, 5e-24], N[(2.0 * N[(l * N[(l / N[(t$95$m * N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l * N[(N[(-0.3333333333333333 * N[(l / k), $MachinePrecision] + N[(2.0 * N[(l / N[Power[k, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 5 \cdot 10^{-24}:\\
\;\;\;\;2 \cdot \left(\ell \cdot \frac{\ell}{t\_m \cdot {k}^{4}}\right)\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \frac{\mathsf{fma}\left(-0.3333333333333333, \frac{\ell}{k}, 2 \cdot \frac{\ell}{{k}^{3}}\right)}{k \cdot t\_m}\\
\end{array}
\end{array}
if k < 4.9999999999999998e-24Initial program 40.8%
associate-/r*40.8%
associate-/l/40.8%
associate-*l/42.5%
associate-/r/42.5%
+-commutative42.5%
unpow242.5%
sqr-neg42.5%
distribute-frac-neg242.5%
distribute-frac-neg242.5%
unpow242.5%
+-rgt-identity42.5%
metadata-eval42.5%
associate--l+42.5%
+-commutative42.5%
associate--l+42.5%
Simplified51.3%
Taylor expanded in k around 0 67.3%
pow267.3%
*-un-lft-identity67.3%
times-frac77.7%
*-commutative77.7%
Applied egg-rr77.7%
if 4.9999999999999998e-24 < k Initial program 31.4%
associate-/r*31.3%
associate-/l/31.3%
associate-*l/31.3%
associate-/r/30.1%
+-commutative30.1%
unpow230.1%
sqr-neg30.1%
distribute-frac-neg230.1%
distribute-frac-neg230.1%
unpow230.1%
+-rgt-identity30.1%
metadata-eval30.1%
associate--l+30.1%
+-commutative30.1%
associate--l+30.1%
Simplified36.5%
associate-*r*37.9%
unpow237.9%
times-frac51.8%
div-inv51.8%
frac-times51.8%
metadata-eval51.8%
Applied egg-rr51.8%
Taylor expanded in k around 0 44.5%
Taylor expanded in t around 0 54.9%
associate-/l*55.0%
fma-define55.0%
Simplified55.0%
Final simplification70.6%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 2.55e-20)
(* 2.0 (* l (/ l (* t_m (pow k 4.0)))))
(/
(* l (+ (* 2.0 (/ l (pow k 3.0))) (* (/ l k) -0.3333333333333333)))
(* k t_m)))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 2.55e-20) {
tmp = 2.0 * (l * (l / (t_m * pow(k, 4.0))));
} else {
tmp = (l * ((2.0 * (l / pow(k, 3.0))) + ((l / k) * -0.3333333333333333))) / (k * t_m);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 2.55d-20) then
tmp = 2.0d0 * (l * (l / (t_m * (k ** 4.0d0))))
else
tmp = (l * ((2.0d0 * (l / (k ** 3.0d0))) + ((l / k) * (-0.3333333333333333d0)))) / (k * t_m)
end if
code = t_s * tmp
end function
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) {
double tmp;
if (k <= 2.55e-20) {
tmp = 2.0 * (l * (l / (t_m * Math.pow(k, 4.0))));
} else {
tmp = (l * ((2.0 * (l / Math.pow(k, 3.0))) + ((l / k) * -0.3333333333333333))) / (k * t_m);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if k <= 2.55e-20: tmp = 2.0 * (l * (l / (t_m * math.pow(k, 4.0)))) else: tmp = (l * ((2.0 * (l / math.pow(k, 3.0))) + ((l / k) * -0.3333333333333333))) / (k * t_m) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 2.55e-20) tmp = Float64(2.0 * Float64(l * Float64(l / Float64(t_m * (k ^ 4.0))))); else tmp = Float64(Float64(l * Float64(Float64(2.0 * Float64(l / (k ^ 3.0))) + Float64(Float64(l / k) * -0.3333333333333333))) / Float64(k * t_m)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (k <= 2.55e-20) tmp = 2.0 * (l * (l / (t_m * (k ^ 4.0)))); else tmp = (l * ((2.0 * (l / (k ^ 3.0))) + ((l / k) * -0.3333333333333333))) / (k * t_m); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[k, 2.55e-20], N[(2.0 * N[(l * N[(l / N[(t$95$m * N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * N[(N[(2.0 * N[(l / N[Power[k, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(l / k), $MachinePrecision] * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * t$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 2.55 \cdot 10^{-20}:\\
\;\;\;\;2 \cdot \left(\ell \cdot \frac{\ell}{t\_m \cdot {k}^{4}}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \left(2 \cdot \frac{\ell}{{k}^{3}} + \frac{\ell}{k} \cdot -0.3333333333333333\right)}{k \cdot t\_m}\\
\end{array}
\end{array}
if k < 2.55000000000000009e-20Initial program 40.8%
associate-/r*40.8%
associate-/l/40.8%
associate-*l/42.5%
associate-/r/42.5%
+-commutative42.5%
unpow242.5%
sqr-neg42.5%
distribute-frac-neg242.5%
distribute-frac-neg242.5%
unpow242.5%
+-rgt-identity42.5%
metadata-eval42.5%
associate--l+42.5%
+-commutative42.5%
associate--l+42.5%
Simplified51.3%
Taylor expanded in k around 0 67.3%
pow267.3%
*-un-lft-identity67.3%
times-frac77.7%
*-commutative77.7%
Applied egg-rr77.7%
if 2.55000000000000009e-20 < k Initial program 31.4%
associate-/r*31.3%
associate-/l/31.3%
associate-*l/31.3%
associate-/r/30.1%
+-commutative30.1%
unpow230.1%
sqr-neg30.1%
distribute-frac-neg230.1%
distribute-frac-neg230.1%
unpow230.1%
+-rgt-identity30.1%
metadata-eval30.1%
associate--l+30.1%
+-commutative30.1%
associate--l+30.1%
Simplified36.5%
associate-*r*37.9%
unpow237.9%
times-frac51.8%
div-inv51.8%
frac-times51.8%
metadata-eval51.8%
Applied egg-rr51.8%
Taylor expanded in k around 0 44.5%
Taylor expanded in t around 0 54.9%
Final simplification70.6%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 1 t) (FPCore (t_s t_m l k) :precision binary64 (* t_s (* 2.0 (* (* l (pow k -4.0)) (/ l t_m)))))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
return t_s * (2.0 * ((l * pow(k, -4.0)) * (l / t_m)));
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
code = t_s * (2.0d0 * ((l * (k ** (-4.0d0))) * (l / t_m)))
end function
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) {
return t_s * (2.0 * ((l * Math.pow(k, -4.0)) * (l / t_m)));
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): return t_s * (2.0 * ((l * math.pow(k, -4.0)) * (l / t_m)))
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) return Float64(t_s * Float64(2.0 * Float64(Float64(l * (k ^ -4.0)) * Float64(l / t_m)))) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k) tmp = t_s * (2.0 * ((l * (k ^ -4.0)) * (l / t_m))); end
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_] := N[(t$95$s * N[(2.0 * N[(N[(l * N[Power[k, -4.0], $MachinePrecision]), $MachinePrecision] * N[(l / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(2 \cdot \left(\left(\ell \cdot {k}^{-4}\right) \cdot \frac{\ell}{t\_m}\right)\right)
\end{array}
Initial program 37.9%
associate-/r*37.8%
associate-/l/37.8%
associate-*l/39.0%
associate-/r/38.6%
+-commutative38.6%
unpow238.6%
sqr-neg38.6%
distribute-frac-neg238.6%
distribute-frac-neg238.6%
unpow238.6%
+-rgt-identity38.6%
metadata-eval38.6%
associate--l+38.6%
+-commutative38.6%
associate--l+38.6%
Simplified46.7%
Taylor expanded in k around 0 62.7%
pow262.7%
times-frac69.5%
Applied egg-rr69.5%
associate-*r/70.0%
div-inv69.6%
pow-flip69.6%
metadata-eval69.6%
Applied egg-rr69.6%
associate-/l*69.1%
Simplified69.1%
Final simplification69.1%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 1 t) (FPCore (t_s t_m l k) :precision binary64 (* t_s (* 2.0 (* (/ l (pow k 4.0)) (/ l t_m)))))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
return t_s * (2.0 * ((l / pow(k, 4.0)) * (l / t_m)));
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
code = t_s * (2.0d0 * ((l / (k ** 4.0d0)) * (l / t_m)))
end function
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) {
return t_s * (2.0 * ((l / Math.pow(k, 4.0)) * (l / t_m)));
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): return t_s * (2.0 * ((l / math.pow(k, 4.0)) * (l / t_m)))
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) return Float64(t_s * Float64(2.0 * Float64(Float64(l / (k ^ 4.0)) * Float64(l / t_m)))) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k) tmp = t_s * (2.0 * ((l / (k ^ 4.0)) * (l / t_m))); end
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_] := N[(t$95$s * N[(2.0 * N[(N[(l / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision] * N[(l / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(2 \cdot \left(\frac{\ell}{{k}^{4}} \cdot \frac{\ell}{t\_m}\right)\right)
\end{array}
Initial program 37.9%
associate-/r*37.8%
associate-/l/37.8%
associate-*l/39.0%
associate-/r/38.6%
+-commutative38.6%
unpow238.6%
sqr-neg38.6%
distribute-frac-neg238.6%
distribute-frac-neg238.6%
unpow238.6%
+-rgt-identity38.6%
metadata-eval38.6%
associate--l+38.6%
+-commutative38.6%
associate--l+38.6%
Simplified46.7%
Taylor expanded in k around 0 62.7%
pow262.7%
times-frac69.5%
Applied egg-rr69.5%
Final simplification69.5%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 1 t) (FPCore (t_s t_m l k) :precision binary64 (* t_s (* 2.0 (/ (* l (* l (pow k -4.0))) t_m))))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
return t_s * (2.0 * ((l * (l * pow(k, -4.0))) / t_m));
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
code = t_s * (2.0d0 * ((l * (l * (k ** (-4.0d0)))) / t_m))
end function
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) {
return t_s * (2.0 * ((l * (l * Math.pow(k, -4.0))) / t_m));
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): return t_s * (2.0 * ((l * (l * math.pow(k, -4.0))) / t_m))
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) return Float64(t_s * Float64(2.0 * Float64(Float64(l * Float64(l * (k ^ -4.0))) / t_m))) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k) tmp = t_s * (2.0 * ((l * (l * (k ^ -4.0))) / t_m)); end
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_] := N[(t$95$s * N[(2.0 * N[(N[(l * N[(l * N[Power[k, -4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(2 \cdot \frac{\ell \cdot \left(\ell \cdot {k}^{-4}\right)}{t\_m}\right)
\end{array}
Initial program 37.9%
associate-/r*37.8%
associate-/l/37.8%
associate-*l/39.0%
associate-/r/38.6%
+-commutative38.6%
unpow238.6%
sqr-neg38.6%
distribute-frac-neg238.6%
distribute-frac-neg238.6%
unpow238.6%
+-rgt-identity38.6%
metadata-eval38.6%
associate--l+38.6%
+-commutative38.6%
associate--l+38.6%
Simplified46.7%
Taylor expanded in k around 0 62.7%
pow262.7%
times-frac69.5%
Applied egg-rr69.5%
associate-*r/70.0%
div-inv69.6%
pow-flip69.6%
metadata-eval69.6%
Applied egg-rr69.6%
Final simplification69.6%
herbie shell --seed 2024050
(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))))