
(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 15 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) - 1.0));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) - 1.0d0))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) - 1.0));
}
def code(t, l, k): return 2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t), 2.0)) - 1.0))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) - 1.0))) end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) - 1.0)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}
\end{array}
l_m = (fabs.f64 l)
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k_m)
:precision binary64
(let* ((t_2 (cbrt (* (tan k_m) (sin k_m)))))
(*
t_s
(if (<= k_m 0.000305)
(pow
(*
l_m
(*
(sqrt (fma -0.3333333333333333 (/ (pow k_m 2.0) t_m) (/ 2.0 t_m)))
(pow k_m -2.0)))
2.0)
(/
(pow (/ (/ (/ (sqrt 2.0) k_m) (pow (sqrt (cbrt l_m)) -4.0)) t_2) 2.0)
(* t_m (* (pow (cbrt l_m) -2.0) t_2)))))))l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
double t_2 = cbrt((tan(k_m) * sin(k_m)));
double tmp;
if (k_m <= 0.000305) {
tmp = pow((l_m * (sqrt(fma(-0.3333333333333333, (pow(k_m, 2.0) / t_m), (2.0 / t_m))) * pow(k_m, -2.0))), 2.0);
} else {
tmp = pow((((sqrt(2.0) / k_m) / pow(sqrt(cbrt(l_m)), -4.0)) / t_2), 2.0) / (t_m * (pow(cbrt(l_m), -2.0) * t_2));
}
return t_s * tmp;
}
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) t_2 = cbrt(Float64(tan(k_m) * sin(k_m))) tmp = 0.0 if (k_m <= 0.000305) tmp = Float64(l_m * Float64(sqrt(fma(-0.3333333333333333, Float64((k_m ^ 2.0) / t_m), Float64(2.0 / t_m))) * (k_m ^ -2.0))) ^ 2.0; else tmp = Float64((Float64(Float64(Float64(sqrt(2.0) / k_m) / (sqrt(cbrt(l_m)) ^ -4.0)) / t_2) ^ 2.0) / Float64(t_m * Float64((cbrt(l_m) ^ -2.0) * t_2))); end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := Block[{t$95$2 = N[Power[N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 0.000305], N[Power[N[(l$95$m * N[(N[Sqrt[N[(-0.3333333333333333 * N[(N[Power[k$95$m, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] + N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Power[k$95$m, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / N[Power[N[Sqrt[N[Power[l$95$m, 1/3], $MachinePrecision]], $MachinePrecision], -4.0], $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision], 2.0], $MachinePrecision] / N[(t$95$m * N[(N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], -2.0], $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sqrt[3]{\tan k\_m \cdot \sin k\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.000305:\\
\;\;\;\;{\left(l\_m \cdot \left(\sqrt{\mathsf{fma}\left(-0.3333333333333333, \frac{{k\_m}^{2}}{t\_m}, \frac{2}{t\_m}\right)} \cdot {k\_m}^{-2}\right)\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left(\frac{\frac{\frac{\sqrt{2}}{k\_m}}{{\left(\sqrt{\sqrt[3]{l\_m}}\right)}^{-4}}}{t\_2}\right)}^{2}}{t\_m \cdot \left({\left(\sqrt[3]{l\_m}\right)}^{-2} \cdot t\_2\right)}\\
\end{array}
\end{array}
\end{array}
if k < 3.04999999999999987e-4Initial program 35.4%
Simplified41.3%
Taylor expanded in k around 0 27.7%
fma-define27.7%
*-commutative27.7%
associate-/l*26.6%
associate-*r/26.6%
*-commutative26.6%
associate-/l*26.6%
Simplified26.6%
div-inv26.6%
pow226.6%
associate-*r/27.7%
pow227.7%
pow-prod-down47.5%
pow247.5%
*-commutative47.5%
pow247.5%
pow-flip47.5%
metadata-eval47.5%
Applied egg-rr47.5%
*-commutative47.5%
fma-define47.5%
associate-*r/47.5%
unpow247.5%
swap-sqr27.7%
unpow227.7%
unpow227.7%
*-commutative27.7%
associate-*r*27.7%
associate-*l/26.6%
associate-*r/26.6%
distribute-rgt-in43.8%
metadata-eval43.8%
associate-*r/43.8%
fma-define43.8%
Simplified43.8%
add-sqr-sqrt24.4%
pow224.4%
*-commutative24.4%
pow224.4%
sqrt-prod24.4%
sqrt-prod22.3%
sqrt-prod12.4%
add-sqr-sqrt28.8%
sqrt-pow134.4%
metadata-eval34.4%
Applied egg-rr34.4%
associate-*l*35.0%
Simplified35.0%
if 3.04999999999999987e-4 < k Initial program 32.2%
*-commutative32.2%
associate-/r*32.2%
Simplified38.1%
add-sqr-sqrt38.1%
add-cube-cbrt38.0%
times-frac38.0%
Applied egg-rr89.1%
associate-/r/89.1%
associate-/r*89.1%
associate-/r/89.1%
Simplified89.1%
frac-times81.9%
associate-/l*81.9%
div-inv81.8%
pow-flip81.8%
metadata-eval81.8%
Applied egg-rr81.8%
times-frac88.9%
associate-*r/89.0%
associate-/r*88.9%
associate-*l/81.8%
unpow281.8%
Simplified94.5%
add-sqr-sqrt44.7%
unpow-prod-down44.7%
Applied egg-rr44.7%
pow-sqr44.7%
metadata-eval44.7%
Simplified44.7%
Final simplification37.6%
l_m = (fabs.f64 l)
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k_m)
:precision binary64
(let* ((t_2 (* (cbrt (tan k_m)) (cbrt (sin k_m))))
(t_3 (pow (cbrt l_m) -2.0)))
(*
t_s
(/ (pow (/ (/ (/ (sqrt 2.0) k_m) t_3) t_2) 2.0) (* t_m (* t_3 t_2))))))l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
double t_2 = cbrt(tan(k_m)) * cbrt(sin(k_m));
double t_3 = pow(cbrt(l_m), -2.0);
return t_s * (pow((((sqrt(2.0) / k_m) / t_3) / t_2), 2.0) / (t_m * (t_3 * t_2)));
}
l_m = Math.abs(l);
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k_m) {
double t_2 = Math.cbrt(Math.tan(k_m)) * Math.cbrt(Math.sin(k_m));
double t_3 = Math.pow(Math.cbrt(l_m), -2.0);
return t_s * (Math.pow((((Math.sqrt(2.0) / k_m) / t_3) / t_2), 2.0) / (t_m * (t_3 * t_2)));
}
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) t_2 = Float64(cbrt(tan(k_m)) * cbrt(sin(k_m))) t_3 = cbrt(l_m) ^ -2.0 return Float64(t_s * Float64((Float64(Float64(Float64(sqrt(2.0) / k_m) / t_3) / t_2) ^ 2.0) / Float64(t_m * Float64(t_3 * t_2)))) end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := Block[{t$95$2 = N[(N[Power[N[Tan[k$95$m], $MachinePrecision], 1/3], $MachinePrecision] * N[Power[N[Sin[k$95$m], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], -2.0], $MachinePrecision]}, N[(t$95$s * N[(N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / t$95$3), $MachinePrecision] / t$95$2), $MachinePrecision], 2.0], $MachinePrecision] / N[(t$95$m * N[(t$95$3 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sqrt[3]{\tan k\_m} \cdot \sqrt[3]{\sin k\_m}\\
t_3 := {\left(\sqrt[3]{l\_m}\right)}^{-2}\\
t\_s \cdot \frac{{\left(\frac{\frac{\frac{\sqrt{2}}{k\_m}}{t\_3}}{t\_2}\right)}^{2}}{t\_m \cdot \left(t\_3 \cdot t\_2\right)}
\end{array}
\end{array}
Initial program 34.5%
*-commutative34.5%
associate-/r*34.6%
Simplified40.4%
add-sqr-sqrt40.4%
add-cube-cbrt40.4%
times-frac40.4%
Applied egg-rr82.8%
associate-/r/82.8%
associate-/r*82.9%
associate-/r/83.1%
Simplified83.1%
frac-times79.2%
associate-/l*79.2%
div-inv79.2%
pow-flip79.1%
metadata-eval79.1%
Applied egg-rr79.1%
times-frac83.0%
associate-*r/83.0%
associate-/r*83.1%
associate-*l/79.1%
unpow279.1%
Simplified92.5%
*-commutative92.5%
cbrt-prod92.5%
Applied egg-rr92.5%
*-commutative92.5%
cbrt-prod92.5%
Applied egg-rr93.2%
l_m = (fabs.f64 l)
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k_m)
:precision binary64
(let* ((t_2 (pow (cbrt l_m) -2.0)))
(*
t_s
(if (<= k_m 0.0011)
(pow
(*
l_m
(*
(sqrt (fma -0.3333333333333333 (/ (pow k_m 2.0) t_m) (/ 2.0 t_m)))
(pow k_m -2.0)))
2.0)
(/
(pow
(/ (/ (/ (sqrt 2.0) k_m) t_2) (* (cbrt (tan k_m)) (cbrt (sin k_m))))
2.0)
(* t_m (* t_2 (cbrt (* (tan k_m) (sin k_m))))))))))l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
double t_2 = pow(cbrt(l_m), -2.0);
double tmp;
if (k_m <= 0.0011) {
tmp = pow((l_m * (sqrt(fma(-0.3333333333333333, (pow(k_m, 2.0) / t_m), (2.0 / t_m))) * pow(k_m, -2.0))), 2.0);
} else {
tmp = pow((((sqrt(2.0) / k_m) / t_2) / (cbrt(tan(k_m)) * cbrt(sin(k_m)))), 2.0) / (t_m * (t_2 * cbrt((tan(k_m) * sin(k_m)))));
}
return t_s * tmp;
}
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) t_2 = cbrt(l_m) ^ -2.0 tmp = 0.0 if (k_m <= 0.0011) tmp = Float64(l_m * Float64(sqrt(fma(-0.3333333333333333, Float64((k_m ^ 2.0) / t_m), Float64(2.0 / t_m))) * (k_m ^ -2.0))) ^ 2.0; else tmp = Float64((Float64(Float64(Float64(sqrt(2.0) / k_m) / t_2) / Float64(cbrt(tan(k_m)) * cbrt(sin(k_m)))) ^ 2.0) / Float64(t_m * Float64(t_2 * cbrt(Float64(tan(k_m) * sin(k_m)))))); end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := Block[{t$95$2 = N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], -2.0], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 0.0011], N[Power[N[(l$95$m * N[(N[Sqrt[N[(-0.3333333333333333 * N[(N[Power[k$95$m, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] + N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Power[k$95$m, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / t$95$2), $MachinePrecision] / N[(N[Power[N[Tan[k$95$m], $MachinePrecision], 1/3], $MachinePrecision] * N[Power[N[Sin[k$95$m], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] / N[(t$95$m * N[(t$95$2 * N[Power[N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := {\left(\sqrt[3]{l\_m}\right)}^{-2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.0011:\\
\;\;\;\;{\left(l\_m \cdot \left(\sqrt{\mathsf{fma}\left(-0.3333333333333333, \frac{{k\_m}^{2}}{t\_m}, \frac{2}{t\_m}\right)} \cdot {k\_m}^{-2}\right)\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left(\frac{\frac{\frac{\sqrt{2}}{k\_m}}{t\_2}}{\sqrt[3]{\tan k\_m} \cdot \sqrt[3]{\sin k\_m}}\right)}^{2}}{t\_m \cdot \left(t\_2 \cdot \sqrt[3]{\tan k\_m \cdot \sin k\_m}\right)}\\
\end{array}
\end{array}
\end{array}
if k < 0.00110000000000000007Initial program 35.4%
Simplified41.3%
Taylor expanded in k around 0 27.7%
fma-define27.7%
*-commutative27.7%
associate-/l*26.6%
associate-*r/26.6%
*-commutative26.6%
associate-/l*26.6%
Simplified26.6%
div-inv26.6%
pow226.6%
associate-*r/27.7%
pow227.7%
pow-prod-down47.5%
pow247.5%
*-commutative47.5%
pow247.5%
pow-flip47.5%
metadata-eval47.5%
Applied egg-rr47.5%
*-commutative47.5%
fma-define47.5%
associate-*r/47.5%
unpow247.5%
swap-sqr27.7%
unpow227.7%
unpow227.7%
*-commutative27.7%
associate-*r*27.7%
associate-*l/26.6%
associate-*r/26.6%
distribute-rgt-in43.8%
metadata-eval43.8%
associate-*r/43.8%
fma-define43.8%
Simplified43.8%
add-sqr-sqrt24.4%
pow224.4%
*-commutative24.4%
pow224.4%
sqrt-prod24.4%
sqrt-prod22.3%
sqrt-prod12.4%
add-sqr-sqrt28.8%
sqrt-pow134.4%
metadata-eval34.4%
Applied egg-rr34.4%
associate-*l*35.0%
Simplified35.0%
if 0.00110000000000000007 < k Initial program 32.2%
*-commutative32.2%
associate-/r*32.2%
Simplified38.1%
add-sqr-sqrt38.1%
add-cube-cbrt38.0%
times-frac38.0%
Applied egg-rr89.1%
associate-/r/89.1%
associate-/r*89.1%
associate-/r/89.1%
Simplified89.1%
frac-times81.9%
associate-/l*81.9%
div-inv81.8%
pow-flip81.8%
metadata-eval81.8%
Applied egg-rr81.8%
times-frac88.9%
associate-*r/89.0%
associate-/r*88.9%
associate-*l/81.8%
unpow281.8%
Simplified94.5%
*-commutative94.5%
cbrt-prod94.5%
Applied egg-rr94.5%
Final simplification51.0%
l_m = (fabs.f64 l)
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k_m)
:precision binary64
(let* ((t_2 (cbrt (* (tan k_m) (sin k_m)))) (t_3 (pow (cbrt l_m) -2.0)))
(*
t_s
(if (<= k_m 0.00051)
(pow
(*
l_m
(*
(sqrt (fma -0.3333333333333333 (/ (pow k_m 2.0) t_m) (/ 2.0 t_m)))
(pow k_m -2.0)))
2.0)
(/ (pow (/ (/ (/ (sqrt 2.0) k_m) t_3) t_2) 2.0) (* t_m (* t_3 t_2)))))))l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
double t_2 = cbrt((tan(k_m) * sin(k_m)));
double t_3 = pow(cbrt(l_m), -2.0);
double tmp;
if (k_m <= 0.00051) {
tmp = pow((l_m * (sqrt(fma(-0.3333333333333333, (pow(k_m, 2.0) / t_m), (2.0 / t_m))) * pow(k_m, -2.0))), 2.0);
} else {
tmp = pow((((sqrt(2.0) / k_m) / t_3) / t_2), 2.0) / (t_m * (t_3 * t_2));
}
return t_s * tmp;
}
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) t_2 = cbrt(Float64(tan(k_m) * sin(k_m))) t_3 = cbrt(l_m) ^ -2.0 tmp = 0.0 if (k_m <= 0.00051) tmp = Float64(l_m * Float64(sqrt(fma(-0.3333333333333333, Float64((k_m ^ 2.0) / t_m), Float64(2.0 / t_m))) * (k_m ^ -2.0))) ^ 2.0; else tmp = Float64((Float64(Float64(Float64(sqrt(2.0) / k_m) / t_3) / t_2) ^ 2.0) / Float64(t_m * Float64(t_3 * t_2))); end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := Block[{t$95$2 = N[Power[N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]}, Block[{t$95$3 = N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], -2.0], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 0.00051], N[Power[N[(l$95$m * N[(N[Sqrt[N[(-0.3333333333333333 * N[(N[Power[k$95$m, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] + N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Power[k$95$m, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / t$95$3), $MachinePrecision] / t$95$2), $MachinePrecision], 2.0], $MachinePrecision] / N[(t$95$m * N[(t$95$3 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sqrt[3]{\tan k\_m \cdot \sin k\_m}\\
t_3 := {\left(\sqrt[3]{l\_m}\right)}^{-2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.00051:\\
\;\;\;\;{\left(l\_m \cdot \left(\sqrt{\mathsf{fma}\left(-0.3333333333333333, \frac{{k\_m}^{2}}{t\_m}, \frac{2}{t\_m}\right)} \cdot {k\_m}^{-2}\right)\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left(\frac{\frac{\frac{\sqrt{2}}{k\_m}}{t\_3}}{t\_2}\right)}^{2}}{t\_m \cdot \left(t\_3 \cdot t\_2\right)}\\
\end{array}
\end{array}
\end{array}
if k < 5.1e-4Initial program 35.4%
Simplified41.3%
Taylor expanded in k around 0 27.7%
fma-define27.7%
*-commutative27.7%
associate-/l*26.6%
associate-*r/26.6%
*-commutative26.6%
associate-/l*26.6%
Simplified26.6%
div-inv26.6%
pow226.6%
associate-*r/27.7%
pow227.7%
pow-prod-down47.5%
pow247.5%
*-commutative47.5%
pow247.5%
pow-flip47.5%
metadata-eval47.5%
Applied egg-rr47.5%
*-commutative47.5%
fma-define47.5%
associate-*r/47.5%
unpow247.5%
swap-sqr27.7%
unpow227.7%
unpow227.7%
*-commutative27.7%
associate-*r*27.7%
associate-*l/26.6%
associate-*r/26.6%
distribute-rgt-in43.8%
metadata-eval43.8%
associate-*r/43.8%
fma-define43.8%
Simplified43.8%
add-sqr-sqrt24.4%
pow224.4%
*-commutative24.4%
pow224.4%
sqrt-prod24.4%
sqrt-prod22.3%
sqrt-prod12.4%
add-sqr-sqrt28.8%
sqrt-pow134.4%
metadata-eval34.4%
Applied egg-rr34.4%
associate-*l*35.0%
Simplified35.0%
if 5.1e-4 < k Initial program 32.2%
*-commutative32.2%
associate-/r*32.2%
Simplified38.1%
add-sqr-sqrt38.1%
add-cube-cbrt38.0%
times-frac38.0%
Applied egg-rr89.1%
associate-/r/89.1%
associate-/r*89.1%
associate-/r/89.1%
Simplified89.1%
frac-times81.9%
associate-/l*81.9%
div-inv81.8%
pow-flip81.8%
metadata-eval81.8%
Applied egg-rr81.8%
times-frac88.9%
associate-*r/89.0%
associate-/r*88.9%
associate-*l/81.8%
unpow281.8%
Simplified94.5%
Final simplification51.0%
l_m = (fabs.f64 l)
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k_m)
:precision binary64
(let* ((t_2 (* (pow (cbrt l_m) -2.0) (cbrt (* (tan k_m) (sin k_m))))))
(*
t_s
(if (<= k_m 0.00122)
(pow
(*
l_m
(*
(sqrt (fma -0.3333333333333333 (/ (pow k_m 2.0) t_m) (/ 2.0 t_m)))
(pow k_m -2.0)))
2.0)
(/ (pow (/ (/ (sqrt 2.0) k_m) t_2) 2.0) (* t_m t_2))))))l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
double t_2 = pow(cbrt(l_m), -2.0) * cbrt((tan(k_m) * sin(k_m)));
double tmp;
if (k_m <= 0.00122) {
tmp = pow((l_m * (sqrt(fma(-0.3333333333333333, (pow(k_m, 2.0) / t_m), (2.0 / t_m))) * pow(k_m, -2.0))), 2.0);
} else {
tmp = pow(((sqrt(2.0) / k_m) / t_2), 2.0) / (t_m * t_2);
}
return t_s * tmp;
}
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) t_2 = Float64((cbrt(l_m) ^ -2.0) * cbrt(Float64(tan(k_m) * sin(k_m)))) tmp = 0.0 if (k_m <= 0.00122) tmp = Float64(l_m * Float64(sqrt(fma(-0.3333333333333333, Float64((k_m ^ 2.0) / t_m), Float64(2.0 / t_m))) * (k_m ^ -2.0))) ^ 2.0; else tmp = Float64((Float64(Float64(sqrt(2.0) / k_m) / t_2) ^ 2.0) / Float64(t_m * t_2)); end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := Block[{t$95$2 = N[(N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], -2.0], $MachinePrecision] * N[Power[N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 0.00122], N[Power[N[(l$95$m * N[(N[Sqrt[N[(-0.3333333333333333 * N[(N[Power[k$95$m, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] + N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Power[k$95$m, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[Power[N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / t$95$2), $MachinePrecision], 2.0], $MachinePrecision] / N[(t$95$m * t$95$2), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := {\left(\sqrt[3]{l\_m}\right)}^{-2} \cdot \sqrt[3]{\tan k\_m \cdot \sin k\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.00122:\\
\;\;\;\;{\left(l\_m \cdot \left(\sqrt{\mathsf{fma}\left(-0.3333333333333333, \frac{{k\_m}^{2}}{t\_m}, \frac{2}{t\_m}\right)} \cdot {k\_m}^{-2}\right)\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left(\frac{\frac{\sqrt{2}}{k\_m}}{t\_2}\right)}^{2}}{t\_m \cdot t\_2}\\
\end{array}
\end{array}
\end{array}
if k < 0.00121999999999999995Initial program 35.4%
Simplified41.3%
Taylor expanded in k around 0 27.7%
fma-define27.7%
*-commutative27.7%
associate-/l*26.6%
associate-*r/26.6%
*-commutative26.6%
associate-/l*26.6%
Simplified26.6%
div-inv26.6%
pow226.6%
associate-*r/27.7%
pow227.7%
pow-prod-down47.5%
pow247.5%
*-commutative47.5%
pow247.5%
pow-flip47.5%
metadata-eval47.5%
Applied egg-rr47.5%
*-commutative47.5%
fma-define47.5%
associate-*r/47.5%
unpow247.5%
swap-sqr27.7%
unpow227.7%
unpow227.7%
*-commutative27.7%
associate-*r*27.7%
associate-*l/26.6%
associate-*r/26.6%
distribute-rgt-in43.8%
metadata-eval43.8%
associate-*r/43.8%
fma-define43.8%
Simplified43.8%
add-sqr-sqrt24.4%
pow224.4%
*-commutative24.4%
pow224.4%
sqrt-prod24.4%
sqrt-prod22.3%
sqrt-prod12.4%
add-sqr-sqrt28.8%
sqrt-pow134.4%
metadata-eval34.4%
Applied egg-rr34.4%
associate-*l*35.0%
Simplified35.0%
if 0.00121999999999999995 < k Initial program 32.2%
*-commutative32.2%
associate-/r*32.2%
Simplified38.1%
add-sqr-sqrt38.1%
add-cube-cbrt38.0%
times-frac38.0%
Applied egg-rr89.1%
associate-/r/89.1%
associate-/r*89.1%
associate-/r/89.1%
Simplified89.1%
frac-times81.9%
associate-/l*81.9%
div-inv81.8%
pow-flip81.8%
metadata-eval81.8%
Applied egg-rr81.8%
times-frac88.9%
associate-*r/89.0%
associate-/r*88.9%
associate-*l/81.8%
unpow281.8%
Simplified94.5%
unpow294.5%
associate-/l/94.5%
*-commutative94.5%
associate-/l/94.5%
*-commutative94.5%
Applied egg-rr94.5%
unpow294.5%
Simplified94.5%
Final simplification51.0%
l_m = (fabs.f64 l)
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k_m)
:precision binary64
(*
t_s
(if (<= k_m 0.0048)
(pow
(*
l_m
(*
(sqrt (fma -0.3333333333333333 (/ (pow k_m 2.0) t_m) (/ 2.0 t_m)))
(pow k_m -2.0)))
2.0)
(if (<= k_m 2.35e+151)
(*
(/ 2.0 (* (pow k_m 2.0) (* t_m (/ (pow (sin k_m) 2.0) (cos k_m)))))
(* l_m l_m))
(pow
(/
(cbrt (* 2.0 (pow (/ k_m t_m) -2.0)))
(* (cbrt (* (tan k_m) (sin k_m))) (* (pow (cbrt l_m) -2.0) t_m)))
3.0)))))l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
double tmp;
if (k_m <= 0.0048) {
tmp = pow((l_m * (sqrt(fma(-0.3333333333333333, (pow(k_m, 2.0) / t_m), (2.0 / t_m))) * pow(k_m, -2.0))), 2.0);
} else if (k_m <= 2.35e+151) {
tmp = (2.0 / (pow(k_m, 2.0) * (t_m * (pow(sin(k_m), 2.0) / cos(k_m))))) * (l_m * l_m);
} else {
tmp = pow((cbrt((2.0 * pow((k_m / t_m), -2.0))) / (cbrt((tan(k_m) * sin(k_m))) * (pow(cbrt(l_m), -2.0) * t_m))), 3.0);
}
return t_s * tmp;
}
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) tmp = 0.0 if (k_m <= 0.0048) tmp = Float64(l_m * Float64(sqrt(fma(-0.3333333333333333, Float64((k_m ^ 2.0) / t_m), Float64(2.0 / t_m))) * (k_m ^ -2.0))) ^ 2.0; elseif (k_m <= 2.35e+151) tmp = Float64(Float64(2.0 / Float64((k_m ^ 2.0) * Float64(t_m * Float64((sin(k_m) ^ 2.0) / cos(k_m))))) * Float64(l_m * l_m)); else tmp = Float64(cbrt(Float64(2.0 * (Float64(k_m / t_m) ^ -2.0))) / Float64(cbrt(Float64(tan(k_m) * sin(k_m))) * Float64((cbrt(l_m) ^ -2.0) * t_m))) ^ 3.0; end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 0.0048], N[Power[N[(l$95$m * N[(N[Sqrt[N[(-0.3333333333333333 * N[(N[Power[k$95$m, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] + N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Power[k$95$m, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 2.35e+151], N[(N[(2.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l$95$m * l$95$m), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[Power[N[(2.0 * N[Power[N[(k$95$m / t$95$m), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] / N[(N[Power[N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * N[(N[Power[N[Power[l$95$m, 1/3], $MachinePrecision], -2.0], $MachinePrecision] * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.0048:\\
\;\;\;\;{\left(l\_m \cdot \left(\sqrt{\mathsf{fma}\left(-0.3333333333333333, \frac{{k\_m}^{2}}{t\_m}, \frac{2}{t\_m}\right)} \cdot {k\_m}^{-2}\right)\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 2.35 \cdot 10^{+151}:\\
\;\;\;\;\frac{2}{{k\_m}^{2} \cdot \left(t\_m \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}\right)} \cdot \left(l\_m \cdot l\_m\right)\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{\sqrt[3]{2 \cdot {\left(\frac{k\_m}{t\_m}\right)}^{-2}}}{\sqrt[3]{\tan k\_m \cdot \sin k\_m} \cdot \left({\left(\sqrt[3]{l\_m}\right)}^{-2} \cdot t\_m\right)}\right)}^{3}\\
\end{array}
\end{array}
if k < 0.00479999999999999958Initial program 35.4%
Simplified41.3%
Taylor expanded in k around 0 27.7%
fma-define27.7%
*-commutative27.7%
associate-/l*26.6%
associate-*r/26.6%
*-commutative26.6%
associate-/l*26.6%
Simplified26.6%
div-inv26.6%
pow226.6%
associate-*r/27.7%
pow227.7%
pow-prod-down47.5%
pow247.5%
*-commutative47.5%
pow247.5%
pow-flip47.5%
metadata-eval47.5%
Applied egg-rr47.5%
*-commutative47.5%
fma-define47.5%
associate-*r/47.5%
unpow247.5%
swap-sqr27.7%
unpow227.7%
unpow227.7%
*-commutative27.7%
associate-*r*27.7%
associate-*l/26.6%
associate-*r/26.6%
distribute-rgt-in43.8%
metadata-eval43.8%
associate-*r/43.8%
fma-define43.8%
Simplified43.8%
add-sqr-sqrt24.4%
pow224.4%
*-commutative24.4%
pow224.4%
sqrt-prod24.4%
sqrt-prod22.3%
sqrt-prod12.4%
add-sqr-sqrt28.8%
sqrt-pow134.4%
metadata-eval34.4%
Applied egg-rr34.4%
associate-*l*35.0%
Simplified35.0%
if 0.00479999999999999958 < k < 2.34999999999999995e151Initial program 25.6%
Simplified38.4%
Taylor expanded in t around 0 91.1%
associate-/l*91.0%
associate-/l*91.1%
Simplified91.1%
if 2.34999999999999995e151 < k Initial program 37.8%
*-commutative37.8%
associate-/r*37.8%
Simplified40.5%
add-cube-cbrt40.4%
pow340.4%
cbrt-prod40.5%
cbrt-div40.4%
rem-cbrt-cube57.2%
cbrt-prod67.4%
pow267.4%
Applied egg-rr67.4%
add-cube-cbrt67.3%
Applied egg-rr80.9%
unpow280.9%
unpow380.9%
Simplified80.9%
Final simplification48.6%
l_m = (fabs.f64 l)
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k_m)
:precision binary64
(*
t_s
(if (<= k_m 0.00185)
(pow
(*
l_m
(*
(sqrt (fma -0.3333333333333333 (/ (pow k_m 2.0) t_m) (/ 2.0 t_m)))
(pow k_m -2.0)))
2.0)
(*
(/ 2.0 (* (pow k_m 2.0) (* t_m (/ (pow (sin k_m) 2.0) (cos k_m)))))
(* l_m l_m)))))l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
double tmp;
if (k_m <= 0.00185) {
tmp = pow((l_m * (sqrt(fma(-0.3333333333333333, (pow(k_m, 2.0) / t_m), (2.0 / t_m))) * pow(k_m, -2.0))), 2.0);
} else {
tmp = (2.0 / (pow(k_m, 2.0) * (t_m * (pow(sin(k_m), 2.0) / cos(k_m))))) * (l_m * l_m);
}
return t_s * tmp;
}
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) tmp = 0.0 if (k_m <= 0.00185) tmp = Float64(l_m * Float64(sqrt(fma(-0.3333333333333333, Float64((k_m ^ 2.0) / t_m), Float64(2.0 / t_m))) * (k_m ^ -2.0))) ^ 2.0; else tmp = Float64(Float64(2.0 / Float64((k_m ^ 2.0) * Float64(t_m * Float64((sin(k_m) ^ 2.0) / cos(k_m))))) * Float64(l_m * l_m)); end return Float64(t_s * tmp) end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 0.00185], N[Power[N[(l$95$m * N[(N[Sqrt[N[(-0.3333333333333333 * N[(N[Power[k$95$m, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] + N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Power[k$95$m, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[(2.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l$95$m * l$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.00185:\\
\;\;\;\;{\left(l\_m \cdot \left(\sqrt{\mathsf{fma}\left(-0.3333333333333333, \frac{{k\_m}^{2}}{t\_m}, \frac{2}{t\_m}\right)} \cdot {k\_m}^{-2}\right)\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{k\_m}^{2} \cdot \left(t\_m \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}\right)} \cdot \left(l\_m \cdot l\_m\right)\\
\end{array}
\end{array}
if k < 0.0018500000000000001Initial program 35.4%
Simplified41.3%
Taylor expanded in k around 0 27.7%
fma-define27.7%
*-commutative27.7%
associate-/l*26.6%
associate-*r/26.6%
*-commutative26.6%
associate-/l*26.6%
Simplified26.6%
div-inv26.6%
pow226.6%
associate-*r/27.7%
pow227.7%
pow-prod-down47.5%
pow247.5%
*-commutative47.5%
pow247.5%
pow-flip47.5%
metadata-eval47.5%
Applied egg-rr47.5%
*-commutative47.5%
fma-define47.5%
associate-*r/47.5%
unpow247.5%
swap-sqr27.7%
unpow227.7%
unpow227.7%
*-commutative27.7%
associate-*r*27.7%
associate-*l/26.6%
associate-*r/26.6%
distribute-rgt-in43.8%
metadata-eval43.8%
associate-*r/43.8%
fma-define43.8%
Simplified43.8%
add-sqr-sqrt24.4%
pow224.4%
*-commutative24.4%
pow224.4%
sqrt-prod24.4%
sqrt-prod22.3%
sqrt-prod12.4%
add-sqr-sqrt28.8%
sqrt-pow134.4%
metadata-eval34.4%
Applied egg-rr34.4%
associate-*l*35.0%
Simplified35.0%
if 0.0018500000000000001 < k Initial program 32.2%
Simplified38.2%
Taylor expanded in t around 0 77.2%
associate-/l*77.2%
associate-/l*77.2%
Simplified77.2%
l_m = (fabs.f64 l) k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l_m k_m) :precision binary64 (* t_s (* (/ 2.0 (* (pow k_m 2.0) (* t_m (/ (pow (sin k_m) 2.0) (cos k_m))))) (* l_m l_m))))
l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
return t_s * ((2.0 / (pow(k_m, 2.0) * (t_m * (pow(sin(k_m), 2.0) / cos(k_m))))) * (l_m * l_m));
}
l_m = abs(l)
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k_m
code = t_s * ((2.0d0 / ((k_m ** 2.0d0) * (t_m * ((sin(k_m) ** 2.0d0) / cos(k_m))))) * (l_m * l_m))
end function
l_m = Math.abs(l);
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k_m) {
return t_s * ((2.0 / (Math.pow(k_m, 2.0) * (t_m * (Math.pow(Math.sin(k_m), 2.0) / Math.cos(k_m))))) * (l_m * l_m));
}
l_m = math.fabs(l) k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k_m): return t_s * ((2.0 / (math.pow(k_m, 2.0) * (t_m * (math.pow(math.sin(k_m), 2.0) / math.cos(k_m))))) * (l_m * l_m))
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) return Float64(t_s * Float64(Float64(2.0 / Float64((k_m ^ 2.0) * Float64(t_m * Float64((sin(k_m) ^ 2.0) / cos(k_m))))) * Float64(l_m * l_m))) end
l_m = abs(l); k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l_m, k_m) tmp = t_s * ((2.0 / ((k_m ^ 2.0) * (t_m * ((sin(k_m) ^ 2.0) / cos(k_m))))) * (l_m * l_m)); end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := N[(t$95$s * N[(N[(2.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l$95$m * l$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\frac{2}{{k\_m}^{2} \cdot \left(t\_m \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}\right)} \cdot \left(l\_m \cdot l\_m\right)\right)
\end{array}
Initial program 34.5%
Simplified40.5%
Taylor expanded in t around 0 74.9%
associate-/l*74.9%
associate-/l*74.9%
Simplified74.9%
l_m = (fabs.f64 l) k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l_m k_m) :precision binary64 (* t_s (* (* l_m l_m) (* 2.0 (/ (cos k_m) (* (pow k_m 2.0) (* t_m (pow (sin k_m) 2.0))))))))
l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
return t_s * ((l_m * l_m) * (2.0 * (cos(k_m) / (pow(k_m, 2.0) * (t_m * pow(sin(k_m), 2.0))))));
}
l_m = abs(l)
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k_m
code = t_s * ((l_m * l_m) * (2.0d0 * (cos(k_m) / ((k_m ** 2.0d0) * (t_m * (sin(k_m) ** 2.0d0))))))
end function
l_m = Math.abs(l);
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k_m) {
return t_s * ((l_m * l_m) * (2.0 * (Math.cos(k_m) / (Math.pow(k_m, 2.0) * (t_m * Math.pow(Math.sin(k_m), 2.0))))));
}
l_m = math.fabs(l) k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k_m): return t_s * ((l_m * l_m) * (2.0 * (math.cos(k_m) / (math.pow(k_m, 2.0) * (t_m * math.pow(math.sin(k_m), 2.0))))))
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) return Float64(t_s * Float64(Float64(l_m * l_m) * Float64(2.0 * Float64(cos(k_m) / Float64((k_m ^ 2.0) * Float64(t_m * (sin(k_m) ^ 2.0))))))) end
l_m = abs(l); k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l_m, k_m) tmp = t_s * ((l_m * l_m) * (2.0 * (cos(k_m) / ((k_m ^ 2.0) * (t_m * (sin(k_m) ^ 2.0)))))); end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := N[(t$95$s * N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(2.0 * N[(N[Cos[k$95$m], $MachinePrecision] / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(l\_m \cdot l\_m\right) \cdot \left(2 \cdot \frac{\cos k\_m}{{k\_m}^{2} \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}\right)\right)
\end{array}
Initial program 34.5%
Simplified40.5%
Taylor expanded in t around 0 74.9%
Final simplification74.9%
l_m = (fabs.f64 l) k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l_m k_m) :precision binary64 (* t_s (* (* l_m l_m) (* 2.0 (/ 1.0 (* (pow k_m 2.0) (* t_m (pow (sin k_m) 2.0))))))))
l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
return t_s * ((l_m * l_m) * (2.0 * (1.0 / (pow(k_m, 2.0) * (t_m * pow(sin(k_m), 2.0))))));
}
l_m = abs(l)
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k_m
code = t_s * ((l_m * l_m) * (2.0d0 * (1.0d0 / ((k_m ** 2.0d0) * (t_m * (sin(k_m) ** 2.0d0))))))
end function
l_m = Math.abs(l);
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k_m) {
return t_s * ((l_m * l_m) * (2.0 * (1.0 / (Math.pow(k_m, 2.0) * (t_m * Math.pow(Math.sin(k_m), 2.0))))));
}
l_m = math.fabs(l) k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k_m): return t_s * ((l_m * l_m) * (2.0 * (1.0 / (math.pow(k_m, 2.0) * (t_m * math.pow(math.sin(k_m), 2.0))))))
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) return Float64(t_s * Float64(Float64(l_m * l_m) * Float64(2.0 * Float64(1.0 / Float64((k_m ^ 2.0) * Float64(t_m * (sin(k_m) ^ 2.0))))))) end
l_m = abs(l); k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l_m, k_m) tmp = t_s * ((l_m * l_m) * (2.0 * (1.0 / ((k_m ^ 2.0) * (t_m * (sin(k_m) ^ 2.0)))))); end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := N[(t$95$s * N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(2.0 * N[(1.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(l\_m \cdot l\_m\right) \cdot \left(2 \cdot \frac{1}{{k\_m}^{2} \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}\right)\right)
\end{array}
Initial program 34.5%
Simplified40.5%
Taylor expanded in t around 0 74.9%
Taylor expanded in k around 0 66.4%
Final simplification66.4%
l_m = (fabs.f64 l)
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k_m)
:precision binary64
(*
t_s
(if (<= k_m 4.6e+51)
(* 2.0 (/ (* (pow l_m 2.0) (pow k_m -4.0)) t_m))
(* -0.3333333333333333 (/ (pow l_m 2.0) (* t_m (pow k_m 2.0)))))))l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
double tmp;
if (k_m <= 4.6e+51) {
tmp = 2.0 * ((pow(l_m, 2.0) * pow(k_m, -4.0)) / t_m);
} else {
tmp = -0.3333333333333333 * (pow(l_m, 2.0) / (t_m * pow(k_m, 2.0)));
}
return t_s * tmp;
}
l_m = abs(l)
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 4.6d+51) then
tmp = 2.0d0 * (((l_m ** 2.0d0) * (k_m ** (-4.0d0))) / t_m)
else
tmp = (-0.3333333333333333d0) * ((l_m ** 2.0d0) / (t_m * (k_m ** 2.0d0)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k_m) {
double tmp;
if (k_m <= 4.6e+51) {
tmp = 2.0 * ((Math.pow(l_m, 2.0) * Math.pow(k_m, -4.0)) / t_m);
} else {
tmp = -0.3333333333333333 * (Math.pow(l_m, 2.0) / (t_m * Math.pow(k_m, 2.0)));
}
return t_s * tmp;
}
l_m = math.fabs(l) k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k_m): tmp = 0 if k_m <= 4.6e+51: tmp = 2.0 * ((math.pow(l_m, 2.0) * math.pow(k_m, -4.0)) / t_m) else: tmp = -0.3333333333333333 * (math.pow(l_m, 2.0) / (t_m * math.pow(k_m, 2.0))) return t_s * tmp
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) tmp = 0.0 if (k_m <= 4.6e+51) tmp = Float64(2.0 * Float64(Float64((l_m ^ 2.0) * (k_m ^ -4.0)) / t_m)); else tmp = Float64(-0.3333333333333333 * Float64((l_m ^ 2.0) / Float64(t_m * (k_m ^ 2.0)))); end return Float64(t_s * tmp) end
l_m = abs(l); k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k_m) tmp = 0.0; if (k_m <= 4.6e+51) tmp = 2.0 * (((l_m ^ 2.0) * (k_m ^ -4.0)) / t_m); else tmp = -0.3333333333333333 * ((l_m ^ 2.0) / (t_m * (k_m ^ 2.0))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 4.6e+51], N[(2.0 * N[(N[(N[Power[l$95$m, 2.0], $MachinePrecision] * N[Power[k$95$m, -4.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision], N[(-0.3333333333333333 * N[(N[Power[l$95$m, 2.0], $MachinePrecision] / N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 4.6 \cdot 10^{+51}:\\
\;\;\;\;2 \cdot \frac{{l\_m}^{2} \cdot {k\_m}^{-4}}{t\_m}\\
\mathbf{else}:\\
\;\;\;\;-0.3333333333333333 \cdot \frac{{l\_m}^{2}}{t\_m \cdot {k\_m}^{2}}\\
\end{array}
\end{array}
if k < 4.6000000000000001e51Initial program 33.6%
Simplified39.7%
Taylor expanded in k around 0 63.4%
*-commutative63.4%
associate-/r*63.4%
Simplified63.4%
div-inv63.4%
pow-flip63.4%
metadata-eval63.4%
Applied egg-rr63.4%
pow163.4%
associate-*l*63.5%
pow263.5%
Applied egg-rr63.5%
unpow163.5%
associate-*l/63.4%
associate-/l*63.4%
*-commutative63.4%
Simplified63.4%
if 4.6000000000000001e51 < k Initial program 37.9%
Simplified43.3%
Taylor expanded in k around 0 15.0%
fma-define15.0%
*-commutative15.0%
associate-/l*10.9%
associate-*r/10.9%
*-commutative10.9%
associate-/l*10.9%
Simplified10.9%
Taylor expanded in k around inf 62.3%
Final simplification63.2%
l_m = (fabs.f64 l)
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k_m)
:precision binary64
(*
t_s
(if (<= k_m 4.6e+51)
(* (* l_m l_m) (* (/ 2.0 t_m) (pow k_m -4.0)))
(* -0.3333333333333333 (/ (pow l_m 2.0) (* t_m (pow k_m 2.0)))))))l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
double tmp;
if (k_m <= 4.6e+51) {
tmp = (l_m * l_m) * ((2.0 / t_m) * pow(k_m, -4.0));
} else {
tmp = -0.3333333333333333 * (pow(l_m, 2.0) / (t_m * pow(k_m, 2.0)));
}
return t_s * tmp;
}
l_m = abs(l)
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 4.6d+51) then
tmp = (l_m * l_m) * ((2.0d0 / t_m) * (k_m ** (-4.0d0)))
else
tmp = (-0.3333333333333333d0) * ((l_m ** 2.0d0) / (t_m * (k_m ** 2.0d0)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k_m) {
double tmp;
if (k_m <= 4.6e+51) {
tmp = (l_m * l_m) * ((2.0 / t_m) * Math.pow(k_m, -4.0));
} else {
tmp = -0.3333333333333333 * (Math.pow(l_m, 2.0) / (t_m * Math.pow(k_m, 2.0)));
}
return t_s * tmp;
}
l_m = math.fabs(l) k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k_m): tmp = 0 if k_m <= 4.6e+51: tmp = (l_m * l_m) * ((2.0 / t_m) * math.pow(k_m, -4.0)) else: tmp = -0.3333333333333333 * (math.pow(l_m, 2.0) / (t_m * math.pow(k_m, 2.0))) return t_s * tmp
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) tmp = 0.0 if (k_m <= 4.6e+51) tmp = Float64(Float64(l_m * l_m) * Float64(Float64(2.0 / t_m) * (k_m ^ -4.0))); else tmp = Float64(-0.3333333333333333 * Float64((l_m ^ 2.0) / Float64(t_m * (k_m ^ 2.0)))); end return Float64(t_s * tmp) end
l_m = abs(l); k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k_m) tmp = 0.0; if (k_m <= 4.6e+51) tmp = (l_m * l_m) * ((2.0 / t_m) * (k_m ^ -4.0)); else tmp = -0.3333333333333333 * ((l_m ^ 2.0) / (t_m * (k_m ^ 2.0))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 4.6e+51], N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(N[(2.0 / t$95$m), $MachinePrecision] * N[Power[k$95$m, -4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-0.3333333333333333 * N[(N[Power[l$95$m, 2.0], $MachinePrecision] / N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 4.6 \cdot 10^{+51}:\\
\;\;\;\;\left(l\_m \cdot l\_m\right) \cdot \left(\frac{2}{t\_m} \cdot {k\_m}^{-4}\right)\\
\mathbf{else}:\\
\;\;\;\;-0.3333333333333333 \cdot \frac{{l\_m}^{2}}{t\_m \cdot {k\_m}^{2}}\\
\end{array}
\end{array}
if k < 4.6000000000000001e51Initial program 33.6%
Simplified39.7%
Taylor expanded in k around 0 63.4%
*-commutative63.4%
associate-/r*63.4%
Simplified63.4%
div-inv63.4%
pow-flip63.4%
metadata-eval63.4%
Applied egg-rr63.4%
if 4.6000000000000001e51 < k Initial program 37.9%
Simplified43.3%
Taylor expanded in k around 0 15.0%
fma-define15.0%
*-commutative15.0%
associate-/l*10.9%
associate-*r/10.9%
*-commutative10.9%
associate-/l*10.9%
Simplified10.9%
Taylor expanded in k around inf 62.3%
Final simplification63.1%
l_m = (fabs.f64 l)
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l_m k_m)
:precision binary64
(*
t_s
(if (<= k_m 4.6e+51)
(* (* l_m l_m) (* (/ 2.0 t_m) (pow k_m -4.0)))
(* -0.3333333333333333 (/ (/ (pow l_m 2.0) (pow k_m 2.0)) t_m)))))l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
double tmp;
if (k_m <= 4.6e+51) {
tmp = (l_m * l_m) * ((2.0 / t_m) * pow(k_m, -4.0));
} else {
tmp = -0.3333333333333333 * ((pow(l_m, 2.0) / pow(k_m, 2.0)) / t_m);
}
return t_s * tmp;
}
l_m = abs(l)
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 4.6d+51) then
tmp = (l_m * l_m) * ((2.0d0 / t_m) * (k_m ** (-4.0d0)))
else
tmp = (-0.3333333333333333d0) * (((l_m ** 2.0d0) / (k_m ** 2.0d0)) / t_m)
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k_m) {
double tmp;
if (k_m <= 4.6e+51) {
tmp = (l_m * l_m) * ((2.0 / t_m) * Math.pow(k_m, -4.0));
} else {
tmp = -0.3333333333333333 * ((Math.pow(l_m, 2.0) / Math.pow(k_m, 2.0)) / t_m);
}
return t_s * tmp;
}
l_m = math.fabs(l) k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k_m): tmp = 0 if k_m <= 4.6e+51: tmp = (l_m * l_m) * ((2.0 / t_m) * math.pow(k_m, -4.0)) else: tmp = -0.3333333333333333 * ((math.pow(l_m, 2.0) / math.pow(k_m, 2.0)) / t_m) return t_s * tmp
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) tmp = 0.0 if (k_m <= 4.6e+51) tmp = Float64(Float64(l_m * l_m) * Float64(Float64(2.0 / t_m) * (k_m ^ -4.0))); else tmp = Float64(-0.3333333333333333 * Float64(Float64((l_m ^ 2.0) / (k_m ^ 2.0)) / t_m)); end return Float64(t_s * tmp) end
l_m = abs(l); k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l_m, k_m) tmp = 0.0; if (k_m <= 4.6e+51) tmp = (l_m * l_m) * ((2.0 / t_m) * (k_m ^ -4.0)); else tmp = -0.3333333333333333 * (((l_m ^ 2.0) / (k_m ^ 2.0)) / t_m); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 4.6e+51], N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(N[(2.0 / t$95$m), $MachinePrecision] * N[Power[k$95$m, -4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-0.3333333333333333 * N[(N[(N[Power[l$95$m, 2.0], $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 4.6 \cdot 10^{+51}:\\
\;\;\;\;\left(l\_m \cdot l\_m\right) \cdot \left(\frac{2}{t\_m} \cdot {k\_m}^{-4}\right)\\
\mathbf{else}:\\
\;\;\;\;-0.3333333333333333 \cdot \frac{\frac{{l\_m}^{2}}{{k\_m}^{2}}}{t\_m}\\
\end{array}
\end{array}
if k < 4.6000000000000001e51Initial program 33.6%
Simplified39.7%
Taylor expanded in k around 0 63.4%
*-commutative63.4%
associate-/r*63.4%
Simplified63.4%
div-inv63.4%
pow-flip63.4%
metadata-eval63.4%
Applied egg-rr63.4%
if 4.6000000000000001e51 < k Initial program 37.9%
Simplified43.3%
Taylor expanded in k around 0 15.0%
fma-define15.0%
*-commutative15.0%
associate-/l*10.9%
associate-*r/10.9%
*-commutative10.9%
associate-/l*10.9%
Simplified10.9%
div-inv10.9%
pow210.9%
associate-*r/15.0%
pow215.0%
pow-prod-down40.0%
pow240.0%
*-commutative40.0%
pow240.0%
pow-flip40.0%
metadata-eval40.0%
Applied egg-rr40.0%
*-commutative40.0%
fma-define40.0%
associate-*r/40.0%
unpow240.0%
swap-sqr15.0%
unpow215.0%
unpow215.0%
*-commutative15.0%
associate-*r*15.0%
associate-*l/10.9%
associate-*r/10.9%
distribute-rgt-in12.7%
metadata-eval12.7%
associate-*r/12.7%
fma-define12.7%
Simplified12.7%
Taylor expanded in k around inf 62.3%
associate-/r*62.4%
Simplified62.4%
Final simplification63.2%
l_m = (fabs.f64 l) k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l_m k_m) :precision binary64 (* t_s (* (* l_m l_m) (* (/ 2.0 t_m) (pow k_m -4.0)))))
l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
return t_s * ((l_m * l_m) * ((2.0 / t_m) * pow(k_m, -4.0)));
}
l_m = abs(l)
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k_m
code = t_s * ((l_m * l_m) * ((2.0d0 / t_m) * (k_m ** (-4.0d0))))
end function
l_m = Math.abs(l);
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k_m) {
return t_s * ((l_m * l_m) * ((2.0 / t_m) * Math.pow(k_m, -4.0)));
}
l_m = math.fabs(l) k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k_m): return t_s * ((l_m * l_m) * ((2.0 / t_m) * math.pow(k_m, -4.0)))
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) return Float64(t_s * Float64(Float64(l_m * l_m) * Float64(Float64(2.0 / t_m) * (k_m ^ -4.0)))) end
l_m = abs(l); k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l_m, k_m) tmp = t_s * ((l_m * l_m) * ((2.0 / t_m) * (k_m ^ -4.0))); end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := N[(t$95$s * N[(N[(l$95$m * l$95$m), $MachinePrecision] * N[(N[(2.0 / t$95$m), $MachinePrecision] * N[Power[k$95$m, -4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(l\_m \cdot l\_m\right) \cdot \left(\frac{2}{t\_m} \cdot {k\_m}^{-4}\right)\right)
\end{array}
Initial program 34.5%
Simplified40.5%
Taylor expanded in k around 0 62.2%
*-commutative62.2%
associate-/r*62.2%
Simplified62.2%
div-inv62.2%
pow-flip62.2%
metadata-eval62.2%
Applied egg-rr62.2%
Final simplification62.2%
l_m = (fabs.f64 l) k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l_m k_m) :precision binary64 (* t_s 0.0))
l_m = fabs(l);
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l_m, double k_m) {
return t_s * 0.0;
}
l_m = abs(l)
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l_m, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l_m
real(8), intent (in) :: k_m
code = t_s * 0.0d0
end function
l_m = Math.abs(l);
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l_m, double k_m) {
return t_s * 0.0;
}
l_m = math.fabs(l) k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l_m, k_m): return t_s * 0.0
l_m = abs(l) k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l_m, k_m) return Float64(t_s * 0.0) end
l_m = abs(l); k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l_m, k_m) tmp = t_s * 0.0; end
l_m = N[Abs[l], $MachinePrecision]
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l$95$m_, k$95$m_] := N[(t$95$s * 0.0), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot 0
\end{array}
Initial program 34.5%
Simplified40.5%
Taylor expanded in k around 0 62.2%
*-commutative62.2%
associate-/r*62.2%
Simplified62.2%
add-log-exp58.4%
div-inv58.4%
exp-prod46.8%
pow-flip46.8%
metadata-eval46.8%
Applied egg-rr46.8%
Taylor expanded in t around inf 33.3%
Taylor expanded in l around 0 33.9%
herbie shell --seed 2024129
(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))))