
(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 13 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) - 1.0));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) - 1.0d0))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) - 1.0));
}
def code(t, l, k): return 2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t), 2.0)) - 1.0))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) - 1.0))) end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) - 1.0)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}
\end{array}
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 1e-33)
(*
2.0
(pow (* (/ l (* k_m (sqrt t_m))) (/ (sqrt (cos k_m)) (sin k_m))) 2.0))
(*
2.0
(*
(pow (/ (/ l k_m) (sqrt t_m)) 2.0)
(/ (cos k_m) (pow (sin k_m) 2.0)))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1e-33) {
tmp = 2.0 * pow(((l / (k_m * sqrt(t_m))) * (sqrt(cos(k_m)) / sin(k_m))), 2.0);
} else {
tmp = 2.0 * (pow(((l / k_m) / sqrt(t_m)), 2.0) * (cos(k_m) / pow(sin(k_m), 2.0)));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 1d-33) then
tmp = 2.0d0 * (((l / (k_m * sqrt(t_m))) * (sqrt(cos(k_m)) / sin(k_m))) ** 2.0d0)
else
tmp = 2.0d0 * ((((l / k_m) / sqrt(t_m)) ** 2.0d0) * (cos(k_m) / (sin(k_m) ** 2.0d0)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1e-33) {
tmp = 2.0 * Math.pow(((l / (k_m * Math.sqrt(t_m))) * (Math.sqrt(Math.cos(k_m)) / Math.sin(k_m))), 2.0);
} else {
tmp = 2.0 * (Math.pow(((l / k_m) / Math.sqrt(t_m)), 2.0) * (Math.cos(k_m) / Math.pow(Math.sin(k_m), 2.0)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 1e-33: tmp = 2.0 * math.pow(((l / (k_m * math.sqrt(t_m))) * (math.sqrt(math.cos(k_m)) / math.sin(k_m))), 2.0) else: tmp = 2.0 * (math.pow(((l / k_m) / math.sqrt(t_m)), 2.0) * (math.cos(k_m) / math.pow(math.sin(k_m), 2.0))) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 1e-33) tmp = Float64(2.0 * (Float64(Float64(l / Float64(k_m * sqrt(t_m))) * Float64(sqrt(cos(k_m)) / sin(k_m))) ^ 2.0)); else tmp = Float64(2.0 * Float64((Float64(Float64(l / k_m) / sqrt(t_m)) ^ 2.0) * Float64(cos(k_m) / (sin(k_m) ^ 2.0)))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 1e-33) tmp = 2.0 * (((l / (k_m * sqrt(t_m))) * (sqrt(cos(k_m)) / sin(k_m))) ^ 2.0); else tmp = 2.0 * ((((l / k_m) / sqrt(t_m)) ^ 2.0) * (cos(k_m) / (sin(k_m) ^ 2.0))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 1e-33], N[(2.0 * N[Power[N[(N[(l / N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sqrt[N[Cos[k$95$m], $MachinePrecision]], $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[Power[N[(N[(l / k$95$m), $MachinePrecision] / N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 10^{-33}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{k\_m \cdot \sqrt{t\_m}} \cdot \frac{\sqrt{\cos k\_m}}{\sin k\_m}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left({\left(\frac{\frac{\ell}{k\_m}}{\sqrt{t\_m}}\right)}^{2} \cdot \frac{\cos k\_m}{{\sin k\_m}^{2}}\right)\\
\end{array}
\end{array}
if k < 1.0000000000000001e-33Initial program 35.6%
associate-*l*35.6%
associate-/r*35.6%
sub-neg35.6%
distribute-rgt-in27.9%
unpow227.9%
times-frac21.9%
sqr-neg21.9%
times-frac27.9%
unpow227.9%
distribute-rgt-in35.6%
+-commutative35.6%
associate-+l+42.5%
Simplified42.5%
Taylor expanded in t around 0 73.2%
times-frac75.7%
Simplified75.7%
Taylor expanded in l around 0 73.2%
associate-*r*73.2%
times-frac76.1%
*-commutative76.1%
associate-/r*73.8%
Simplified73.8%
add-sqr-sqrt47.2%
pow247.2%
Applied egg-rr43.9%
associate-/l/44.9%
Simplified44.9%
if 1.0000000000000001e-33 < k Initial program 33.8%
associate-*l*33.8%
associate-/r*33.8%
sub-neg33.8%
distribute-rgt-in33.8%
unpow233.8%
times-frac26.7%
sqr-neg26.7%
times-frac33.8%
unpow233.8%
distribute-rgt-in33.8%
+-commutative33.8%
associate-+l+41.3%
Simplified41.3%
Taylor expanded in t around 0 73.3%
times-frac75.6%
Simplified75.6%
Taylor expanded in l around 0 73.3%
associate-*r*73.3%
times-frac73.3%
*-commutative73.3%
associate-/r*75.4%
Simplified75.4%
*-un-lft-identity75.4%
add-sqr-sqrt58.5%
pow258.5%
sqrt-div51.4%
sqrt-div37.4%
unpow237.4%
sqrt-prod12.4%
add-sqr-sqrt38.8%
unpow238.8%
sqrt-prod44.1%
add-sqr-sqrt44.2%
Applied egg-rr44.2%
*-lft-identity44.2%
associate-/l/44.2%
associate-/r*44.2%
Simplified44.2%
Final simplification44.7%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (/ l (* k_m (sqrt t_m)))))
(*
t_s
(if (<= k_m 1.8e-64)
(* 2.0 (pow (* t_2 (/ (sqrt (cos k_m)) (sin k_m))) 2.0))
(* 2.0 (* (/ (cos k_m) (pow (sin k_m) 2.0)) (pow t_2 2.0)))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = l / (k_m * sqrt(t_m));
double tmp;
if (k_m <= 1.8e-64) {
tmp = 2.0 * pow((t_2 * (sqrt(cos(k_m)) / sin(k_m))), 2.0);
} else {
tmp = 2.0 * ((cos(k_m) / pow(sin(k_m), 2.0)) * pow(t_2, 2.0));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: t_2
real(8) :: tmp
t_2 = l / (k_m * sqrt(t_m))
if (k_m <= 1.8d-64) then
tmp = 2.0d0 * ((t_2 * (sqrt(cos(k_m)) / sin(k_m))) ** 2.0d0)
else
tmp = 2.0d0 * ((cos(k_m) / (sin(k_m) ** 2.0d0)) * (t_2 ** 2.0d0))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = l / (k_m * Math.sqrt(t_m));
double tmp;
if (k_m <= 1.8e-64) {
tmp = 2.0 * Math.pow((t_2 * (Math.sqrt(Math.cos(k_m)) / Math.sin(k_m))), 2.0);
} else {
tmp = 2.0 * ((Math.cos(k_m) / Math.pow(Math.sin(k_m), 2.0)) * Math.pow(t_2, 2.0));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): t_2 = l / (k_m * math.sqrt(t_m)) tmp = 0 if k_m <= 1.8e-64: tmp = 2.0 * math.pow((t_2 * (math.sqrt(math.cos(k_m)) / math.sin(k_m))), 2.0) else: tmp = 2.0 * ((math.cos(k_m) / math.pow(math.sin(k_m), 2.0)) * math.pow(t_2, 2.0)) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = Float64(l / Float64(k_m * sqrt(t_m))) tmp = 0.0 if (k_m <= 1.8e-64) tmp = Float64(2.0 * (Float64(t_2 * Float64(sqrt(cos(k_m)) / sin(k_m))) ^ 2.0)); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) / (sin(k_m) ^ 2.0)) * (t_2 ^ 2.0))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) t_2 = l / (k_m * sqrt(t_m)); tmp = 0.0; if (k_m <= 1.8e-64) tmp = 2.0 * ((t_2 * (sqrt(cos(k_m)) / sin(k_m))) ^ 2.0); else tmp = 2.0 * ((cos(k_m) / (sin(k_m) ^ 2.0)) * (t_2 ^ 2.0)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[(l / N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 1.8e-64], N[(2.0 * N[Power[N[(t$95$2 * N[(N[Sqrt[N[Cos[k$95$m], $MachinePrecision]], $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[t$95$2, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{\ell}{k\_m \cdot \sqrt{t\_m}}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 1.8 \cdot 10^{-64}:\\
\;\;\;\;2 \cdot {\left(t\_2 \cdot \frac{\sqrt{\cos k\_m}}{\sin k\_m}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{\cos k\_m}{{\sin k\_m}^{2}} \cdot {t\_2}^{2}\right)\\
\end{array}
\end{array}
\end{array}
if k < 1.7999999999999999e-64Initial program 36.4%
associate-*l*36.4%
associate-/r*36.4%
sub-neg36.4%
distribute-rgt-in28.2%
unpow228.2%
times-frac21.9%
sqr-neg21.9%
times-frac28.2%
unpow228.2%
distribute-rgt-in36.4%
+-commutative36.4%
associate-+l+43.2%
Simplified43.2%
Taylor expanded in t around 0 72.8%
times-frac74.9%
Simplified74.9%
Taylor expanded in l around 0 72.8%
associate-*r*72.8%
times-frac74.7%
*-commutative74.7%
associate-/r*72.3%
Simplified72.3%
add-sqr-sqrt48.2%
pow248.2%
Applied egg-rr44.7%
associate-/l/45.8%
Simplified45.8%
if 1.7999999999999999e-64 < k Initial program 32.3%
associate-*l*32.3%
associate-/r*32.3%
sub-neg32.3%
distribute-rgt-in32.3%
unpow232.3%
times-frac26.0%
sqr-neg26.0%
times-frac32.3%
unpow232.3%
distribute-rgt-in32.3%
+-commutative32.3%
associate-+l+40.1%
Simplified40.1%
Taylor expanded in t around 0 74.2%
times-frac77.4%
Simplified77.4%
Taylor expanded in l around 0 74.2%
associate-*r*74.1%
times-frac76.4%
*-commutative76.4%
associate-/r*78.3%
Simplified78.3%
*-un-lft-identity78.3%
add-sqr-sqrt55.0%
pow255.0%
sqrt-div48.8%
sqrt-div36.5%
unpow236.5%
sqrt-prod10.9%
add-sqr-sqrt37.7%
unpow237.7%
sqrt-prod42.4%
add-sqr-sqrt42.4%
Applied egg-rr42.4%
*-lft-identity42.4%
associate-/l/42.4%
Simplified42.4%
Final simplification44.7%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 0.008)
(*
2.0
(pow (* (/ l (* k_m (sqrt t_m))) (/ (sqrt (cos k_m)) (sin k_m))) 2.0))
(*
2.0
(* (/ (cos k_m) (pow (sin k_m) 2.0)) (/ (pow (/ l k_m) 2.0) t_m))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 0.008) {
tmp = 2.0 * pow(((l / (k_m * sqrt(t_m))) * (sqrt(cos(k_m)) / sin(k_m))), 2.0);
} else {
tmp = 2.0 * ((cos(k_m) / pow(sin(k_m), 2.0)) * (pow((l / k_m), 2.0) / t_m));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 0.008d0) then
tmp = 2.0d0 * (((l / (k_m * sqrt(t_m))) * (sqrt(cos(k_m)) / sin(k_m))) ** 2.0d0)
else
tmp = 2.0d0 * ((cos(k_m) / (sin(k_m) ** 2.0d0)) * (((l / k_m) ** 2.0d0) / t_m))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 0.008) {
tmp = 2.0 * Math.pow(((l / (k_m * Math.sqrt(t_m))) * (Math.sqrt(Math.cos(k_m)) / Math.sin(k_m))), 2.0);
} else {
tmp = 2.0 * ((Math.cos(k_m) / Math.pow(Math.sin(k_m), 2.0)) * (Math.pow((l / k_m), 2.0) / t_m));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 0.008: tmp = 2.0 * math.pow(((l / (k_m * math.sqrt(t_m))) * (math.sqrt(math.cos(k_m)) / math.sin(k_m))), 2.0) else: tmp = 2.0 * ((math.cos(k_m) / math.pow(math.sin(k_m), 2.0)) * (math.pow((l / k_m), 2.0) / t_m)) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 0.008) tmp = Float64(2.0 * (Float64(Float64(l / Float64(k_m * sqrt(t_m))) * Float64(sqrt(cos(k_m)) / sin(k_m))) ^ 2.0)); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) / (sin(k_m) ^ 2.0)) * Float64((Float64(l / k_m) ^ 2.0) / t_m))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 0.008) tmp = 2.0 * (((l / (k_m * sqrt(t_m))) * (sqrt(cos(k_m)) / sin(k_m))) ^ 2.0); else tmp = 2.0 * ((cos(k_m) / (sin(k_m) ^ 2.0)) * (((l / k_m) ^ 2.0) / t_m)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 0.008], N[(2.0 * N[Power[N[(N[(l / N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sqrt[N[Cos[k$95$m], $MachinePrecision]], $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.008:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{k\_m \cdot \sqrt{t\_m}} \cdot \frac{\sqrt{\cos k\_m}}{\sin k\_m}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{\cos k\_m}{{\sin k\_m}^{2}} \cdot \frac{{\left(\frac{\ell}{k\_m}\right)}^{2}}{t\_m}\right)\\
\end{array}
\end{array}
if k < 0.0080000000000000002Initial program 35.4%
associate-*l*35.4%
associate-/r*35.4%
sub-neg35.4%
distribute-rgt-in27.9%
unpow227.9%
times-frac21.6%
sqr-neg21.6%
times-frac27.9%
unpow227.9%
distribute-rgt-in35.4%
+-commutative35.4%
associate-+l+42.0%
Simplified42.0%
Taylor expanded in t around 0 73.6%
times-frac76.0%
Simplified76.0%
Taylor expanded in l around 0 73.6%
associate-*r*73.6%
times-frac76.4%
*-commutative76.4%
associate-/r*74.2%
Simplified74.2%
add-sqr-sqrt46.5%
pow246.5%
Applied egg-rr43.3%
associate-/l/44.3%
Simplified44.3%
if 0.0080000000000000002 < k Initial program 34.4%
associate-*l*34.3%
associate-/r*34.4%
sub-neg34.4%
distribute-rgt-in34.4%
unpow234.4%
times-frac28.1%
sqr-neg28.1%
times-frac34.4%
unpow234.4%
distribute-rgt-in34.4%
+-commutative34.4%
associate-+l+42.7%
Simplified42.7%
Taylor expanded in t around 0 72.0%
times-frac74.6%
Simplified74.6%
Taylor expanded in l around 0 72.0%
associate-*r*71.9%
times-frac72.0%
*-commutative72.0%
associate-/r*74.4%
Simplified74.4%
associate-*r/74.3%
add-sqr-sqrt61.8%
pow261.8%
sqrt-div53.9%
sqrt-div38.4%
unpow238.4%
sqrt-prod13.8%
add-sqr-sqrt39.9%
unpow239.9%
sqrt-prod45.8%
add-sqr-sqrt45.9%
Applied egg-rr45.9%
*-commutative45.9%
associate-/l*45.9%
associate-/l/45.9%
associate-/r*45.9%
Simplified45.9%
Taylor expanded in k around inf 72.0%
associate-*r*71.9%
times-frac72.0%
associate-/r*74.6%
unpow274.6%
unpow274.6%
times-frac95.2%
unpow295.2%
Simplified95.2%
Final simplification57.3%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 14.0)
(*
2.0
(pow (/ (* (sqrt (cos k_m)) (/ (/ l k_m) (sqrt t_m))) (sin k_m)) 2.0))
(*
2.0
(* (/ (cos k_m) (pow (sin k_m) 2.0)) (/ (pow (/ l k_m) 2.0) t_m))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 14.0) {
tmp = 2.0 * pow(((sqrt(cos(k_m)) * ((l / k_m) / sqrt(t_m))) / sin(k_m)), 2.0);
} else {
tmp = 2.0 * ((cos(k_m) / pow(sin(k_m), 2.0)) * (pow((l / k_m), 2.0) / t_m));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 14.0d0) then
tmp = 2.0d0 * (((sqrt(cos(k_m)) * ((l / k_m) / sqrt(t_m))) / sin(k_m)) ** 2.0d0)
else
tmp = 2.0d0 * ((cos(k_m) / (sin(k_m) ** 2.0d0)) * (((l / k_m) ** 2.0d0) / t_m))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 14.0) {
tmp = 2.0 * Math.pow(((Math.sqrt(Math.cos(k_m)) * ((l / k_m) / Math.sqrt(t_m))) / Math.sin(k_m)), 2.0);
} else {
tmp = 2.0 * ((Math.cos(k_m) / Math.pow(Math.sin(k_m), 2.0)) * (Math.pow((l / k_m), 2.0) / t_m));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 14.0: tmp = 2.0 * math.pow(((math.sqrt(math.cos(k_m)) * ((l / k_m) / math.sqrt(t_m))) / math.sin(k_m)), 2.0) else: tmp = 2.0 * ((math.cos(k_m) / math.pow(math.sin(k_m), 2.0)) * (math.pow((l / k_m), 2.0) / t_m)) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 14.0) tmp = Float64(2.0 * (Float64(Float64(sqrt(cos(k_m)) * Float64(Float64(l / k_m) / sqrt(t_m))) / sin(k_m)) ^ 2.0)); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) / (sin(k_m) ^ 2.0)) * Float64((Float64(l / k_m) ^ 2.0) / t_m))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 14.0) tmp = 2.0 * (((sqrt(cos(k_m)) * ((l / k_m) / sqrt(t_m))) / sin(k_m)) ^ 2.0); else tmp = 2.0 * ((cos(k_m) / (sin(k_m) ^ 2.0)) * (((l / k_m) ^ 2.0) / t_m)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 14.0], N[(2.0 * N[Power[N[(N[(N[Sqrt[N[Cos[k$95$m], $MachinePrecision]], $MachinePrecision] * N[(N[(l / k$95$m), $MachinePrecision] / N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 14:\\
\;\;\;\;2 \cdot {\left(\frac{\sqrt{\cos k\_m} \cdot \frac{\frac{\ell}{k\_m}}{\sqrt{t\_m}}}{\sin k\_m}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{\cos k\_m}{{\sin k\_m}^{2}} \cdot \frac{{\left(\frac{\ell}{k\_m}\right)}^{2}}{t\_m}\right)\\
\end{array}
\end{array}
if k < 14Initial program 35.2%
associate-*l*35.2%
associate-/r*35.2%
sub-neg35.2%
distribute-rgt-in27.8%
unpow227.8%
times-frac21.5%
sqr-neg21.5%
times-frac27.8%
unpow227.8%
distribute-rgt-in35.2%
+-commutative35.2%
associate-+l+42.3%
Simplified42.3%
Taylor expanded in t around 0 73.8%
times-frac76.2%
Simplified76.2%
Taylor expanded in l around 0 73.8%
associate-*r*73.8%
times-frac76.5%
*-commutative76.5%
associate-/r*74.4%
Simplified74.4%
add-sqr-sqrt46.2%
pow246.2%
Applied egg-rr43.1%
associate-*r/43.1%
associate-/l/44.1%
associate-/r*44.1%
Simplified44.1%
if 14 < k Initial program 34.9%
associate-*l*34.8%
associate-/r*34.9%
sub-neg34.9%
distribute-rgt-in34.9%
unpow234.9%
times-frac28.5%
sqr-neg28.5%
times-frac34.9%
unpow234.9%
distribute-rgt-in34.9%
+-commutative34.9%
associate-+l+41.8%
Simplified41.8%
Taylor expanded in t around 0 71.6%
times-frac74.2%
Simplified74.2%
Taylor expanded in l around 0 71.6%
associate-*r*71.5%
times-frac71.6%
*-commutative71.6%
associate-/r*74.0%
Simplified74.0%
associate-*r/73.9%
add-sqr-sqrt62.7%
pow262.7%
sqrt-div54.7%
sqrt-div39.0%
unpow239.0%
sqrt-prod14.0%
add-sqr-sqrt40.6%
unpow240.6%
sqrt-prod46.5%
add-sqr-sqrt46.6%
Applied egg-rr46.6%
*-commutative46.6%
associate-/l*46.6%
associate-/l/46.6%
associate-/r*46.6%
Simplified46.6%
Taylor expanded in k around inf 71.6%
associate-*r*71.5%
times-frac71.6%
associate-/r*74.2%
unpow274.2%
unpow274.2%
times-frac95.2%
unpow295.2%
Simplified95.2%
Final simplification56.9%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 0.00048)
(* 2.0 (pow (* (/ l (pow k_m 2.0)) (sqrt (/ 1.0 t_m))) 2.0))
(*
2.0
(* (/ (cos k_m) (pow (sin k_m) 2.0)) (/ (pow (/ l k_m) 2.0) t_m))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 0.00048) {
tmp = 2.0 * pow(((l / pow(k_m, 2.0)) * sqrt((1.0 / t_m))), 2.0);
} else {
tmp = 2.0 * ((cos(k_m) / pow(sin(k_m), 2.0)) * (pow((l / k_m), 2.0) / t_m));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 0.00048d0) then
tmp = 2.0d0 * (((l / (k_m ** 2.0d0)) * sqrt((1.0d0 / t_m))) ** 2.0d0)
else
tmp = 2.0d0 * ((cos(k_m) / (sin(k_m) ** 2.0d0)) * (((l / k_m) ** 2.0d0) / t_m))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 0.00048) {
tmp = 2.0 * Math.pow(((l / Math.pow(k_m, 2.0)) * Math.sqrt((1.0 / t_m))), 2.0);
} else {
tmp = 2.0 * ((Math.cos(k_m) / Math.pow(Math.sin(k_m), 2.0)) * (Math.pow((l / k_m), 2.0) / t_m));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 0.00048: tmp = 2.0 * math.pow(((l / math.pow(k_m, 2.0)) * math.sqrt((1.0 / t_m))), 2.0) else: tmp = 2.0 * ((math.cos(k_m) / math.pow(math.sin(k_m), 2.0)) * (math.pow((l / k_m), 2.0) / t_m)) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 0.00048) tmp = Float64(2.0 * (Float64(Float64(l / (k_m ^ 2.0)) * sqrt(Float64(1.0 / t_m))) ^ 2.0)); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) / (sin(k_m) ^ 2.0)) * Float64((Float64(l / k_m) ^ 2.0) / t_m))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 0.00048) tmp = 2.0 * (((l / (k_m ^ 2.0)) * sqrt((1.0 / t_m))) ^ 2.0); else tmp = 2.0 * ((cos(k_m) / (sin(k_m) ^ 2.0)) * (((l / k_m) ^ 2.0) / t_m)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 0.00048], N[(2.0 * N[Power[N[(N[(l / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.00048:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{{k\_m}^{2}} \cdot \sqrt{\frac{1}{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{\cos k\_m}{{\sin k\_m}^{2}} \cdot \frac{{\left(\frac{\ell}{k\_m}\right)}^{2}}{t\_m}\right)\\
\end{array}
\end{array}
if k < 4.80000000000000012e-4Initial program 35.4%
associate-*l*35.4%
associate-/r*35.4%
sub-neg35.4%
distribute-rgt-in27.9%
unpow227.9%
times-frac21.6%
sqr-neg21.6%
times-frac27.9%
unpow227.9%
distribute-rgt-in35.4%
+-commutative35.4%
associate-+l+42.0%
Simplified42.0%
Taylor expanded in k around 0 65.3%
*-commutative65.3%
associate-/r*65.8%
Simplified65.8%
add-sqr-sqrt41.7%
pow241.7%
div-inv41.7%
pow-flip41.7%
metadata-eval41.7%
Applied egg-rr41.7%
Taylor expanded in l around 0 42.7%
if 4.80000000000000012e-4 < k Initial program 34.4%
associate-*l*34.3%
associate-/r*34.4%
sub-neg34.4%
distribute-rgt-in34.4%
unpow234.4%
times-frac28.1%
sqr-neg28.1%
times-frac34.4%
unpow234.4%
distribute-rgt-in34.4%
+-commutative34.4%
associate-+l+42.7%
Simplified42.7%
Taylor expanded in t around 0 72.0%
times-frac74.6%
Simplified74.6%
Taylor expanded in l around 0 72.0%
associate-*r*71.9%
times-frac72.0%
*-commutative72.0%
associate-/r*74.4%
Simplified74.4%
associate-*r/74.3%
add-sqr-sqrt61.8%
pow261.8%
sqrt-div53.9%
sqrt-div38.4%
unpow238.4%
sqrt-prod13.8%
add-sqr-sqrt39.9%
unpow239.9%
sqrt-prod45.8%
add-sqr-sqrt45.9%
Applied egg-rr45.9%
*-commutative45.9%
associate-/l*45.9%
associate-/l/45.9%
associate-/r*45.9%
Simplified45.9%
Taylor expanded in k around inf 72.0%
associate-*r*71.9%
times-frac72.0%
associate-/r*74.6%
unpow274.6%
unpow274.6%
times-frac95.2%
unpow295.2%
Simplified95.2%
Final simplification56.0%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 27000000.0)
(* 2.0 (pow (* (/ l (pow k_m 2.0)) (sqrt (/ 1.0 t_m))) 2.0))
(*
2.0
(*
(/ (pow (/ l k_m) 2.0) t_m)
(+ -0.16666666666666666 (/ 1.0 (pow k_m 2.0))))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 27000000.0) {
tmp = 2.0 * pow(((l / pow(k_m, 2.0)) * sqrt((1.0 / t_m))), 2.0);
} else {
tmp = 2.0 * ((pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / pow(k_m, 2.0))));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 27000000.0d0) then
tmp = 2.0d0 * (((l / (k_m ** 2.0d0)) * sqrt((1.0d0 / t_m))) ** 2.0d0)
else
tmp = 2.0d0 * ((((l / k_m) ** 2.0d0) / t_m) * ((-0.16666666666666666d0) + (1.0d0 / (k_m ** 2.0d0))))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 27000000.0) {
tmp = 2.0 * Math.pow(((l / Math.pow(k_m, 2.0)) * Math.sqrt((1.0 / t_m))), 2.0);
} else {
tmp = 2.0 * ((Math.pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / Math.pow(k_m, 2.0))));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 27000000.0: tmp = 2.0 * math.pow(((l / math.pow(k_m, 2.0)) * math.sqrt((1.0 / t_m))), 2.0) else: tmp = 2.0 * ((math.pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / math.pow(k_m, 2.0)))) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 27000000.0) tmp = Float64(2.0 * (Float64(Float64(l / (k_m ^ 2.0)) * sqrt(Float64(1.0 / t_m))) ^ 2.0)); else tmp = Float64(2.0 * Float64(Float64((Float64(l / k_m) ^ 2.0) / t_m) * Float64(-0.16666666666666666 + Float64(1.0 / (k_m ^ 2.0))))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 27000000.0) tmp = 2.0 * (((l / (k_m ^ 2.0)) * sqrt((1.0 / t_m))) ^ 2.0); else tmp = 2.0 * ((((l / k_m) ^ 2.0) / t_m) * (-0.16666666666666666 + (1.0 / (k_m ^ 2.0)))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 27000000.0], N[(2.0 * N[Power[N[(N[(l / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(1.0 / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 27000000:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{{k\_m}^{2}} \cdot \sqrt{\frac{1}{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{{\left(\frac{\ell}{k\_m}\right)}^{2}}{t\_m} \cdot \left(-0.16666666666666666 + \frac{1}{{k\_m}^{2}}\right)\right)\\
\end{array}
\end{array}
if k < 2.7e7Initial program 35.3%
associate-*l*35.3%
associate-/r*35.3%
sub-neg35.3%
distribute-rgt-in28.0%
unpow228.0%
times-frac21.8%
sqr-neg21.8%
times-frac28.0%
unpow228.0%
distribute-rgt-in35.3%
+-commutative35.3%
associate-+l+42.4%
Simplified42.4%
Taylor expanded in k around 0 65.1%
*-commutative65.1%
associate-/r*65.5%
Simplified65.5%
add-sqr-sqrt41.5%
pow241.5%
div-inv41.5%
pow-flip41.5%
metadata-eval41.5%
Applied egg-rr41.5%
Taylor expanded in l around 0 42.0%
if 2.7e7 < k Initial program 34.4%
associate-*l*34.4%
associate-/r*34.4%
sub-neg34.4%
distribute-rgt-in34.4%
unpow234.4%
times-frac27.8%
sqr-neg27.8%
times-frac34.4%
unpow234.4%
distribute-rgt-in34.4%
+-commutative34.4%
associate-+l+41.5%
Simplified41.5%
Taylor expanded in t around 0 70.6%
times-frac73.4%
Simplified73.4%
Taylor expanded in k around 0 57.4%
sub-neg57.4%
*-commutative57.4%
associate-/r*57.4%
associate-*r/57.4%
metadata-eval57.4%
distribute-neg-frac57.4%
metadata-eval57.4%
Simplified57.4%
Taylor expanded in l around 0 57.2%
associate-/l*57.0%
associate-/r*57.0%
associate-*r/57.0%
metadata-eval57.0%
div-sub57.0%
associate-/l*57.0%
associate-/r/57.0%
associate-/r*57.4%
unpow257.4%
unpow257.4%
times-frac61.5%
unpow261.5%
sub-neg61.5%
metadata-eval61.5%
+-commutative61.5%
Simplified61.5%
Final simplification46.7%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 27000000.0)
(* 2.0 (pow (* (pow k_m -2.0) (/ l (sqrt t_m))) 2.0))
(*
2.0
(*
(/ (pow (/ l k_m) 2.0) t_m)
(+ -0.16666666666666666 (/ 1.0 (pow k_m 2.0))))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 27000000.0) {
tmp = 2.0 * pow((pow(k_m, -2.0) * (l / sqrt(t_m))), 2.0);
} else {
tmp = 2.0 * ((pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / pow(k_m, 2.0))));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 27000000.0d0) then
tmp = 2.0d0 * (((k_m ** (-2.0d0)) * (l / sqrt(t_m))) ** 2.0d0)
else
tmp = 2.0d0 * ((((l / k_m) ** 2.0d0) / t_m) * ((-0.16666666666666666d0) + (1.0d0 / (k_m ** 2.0d0))))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 27000000.0) {
tmp = 2.0 * Math.pow((Math.pow(k_m, -2.0) * (l / Math.sqrt(t_m))), 2.0);
} else {
tmp = 2.0 * ((Math.pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / Math.pow(k_m, 2.0))));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 27000000.0: tmp = 2.0 * math.pow((math.pow(k_m, -2.0) * (l / math.sqrt(t_m))), 2.0) else: tmp = 2.0 * ((math.pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / math.pow(k_m, 2.0)))) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 27000000.0) tmp = Float64(2.0 * (Float64((k_m ^ -2.0) * Float64(l / sqrt(t_m))) ^ 2.0)); else tmp = Float64(2.0 * Float64(Float64((Float64(l / k_m) ^ 2.0) / t_m) * Float64(-0.16666666666666666 + Float64(1.0 / (k_m ^ 2.0))))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 27000000.0) tmp = 2.0 * (((k_m ^ -2.0) * (l / sqrt(t_m))) ^ 2.0); else tmp = 2.0 * ((((l / k_m) ^ 2.0) / t_m) * (-0.16666666666666666 + (1.0 / (k_m ^ 2.0)))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 27000000.0], N[(2.0 * N[Power[N[(N[Power[k$95$m, -2.0], $MachinePrecision] * N[(l / N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(1.0 / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 27000000:\\
\;\;\;\;2 \cdot {\left({k\_m}^{-2} \cdot \frac{\ell}{\sqrt{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{{\left(\frac{\ell}{k\_m}\right)}^{2}}{t\_m} \cdot \left(-0.16666666666666666 + \frac{1}{{k\_m}^{2}}\right)\right)\\
\end{array}
\end{array}
if k < 2.7e7Initial program 35.3%
associate-*l*35.3%
associate-/r*35.3%
sub-neg35.3%
distribute-rgt-in28.0%
unpow228.0%
times-frac21.8%
sqr-neg21.8%
times-frac28.0%
unpow228.0%
distribute-rgt-in35.3%
+-commutative35.3%
associate-+l+42.4%
Simplified42.4%
Taylor expanded in k around 0 65.1%
*-commutative65.1%
associate-/r*65.5%
Simplified65.5%
add-sqr-sqrt41.5%
pow241.5%
div-inv41.5%
pow-flip41.5%
metadata-eval41.5%
Applied egg-rr41.5%
*-commutative41.5%
sqrt-prod38.8%
sqrt-pow140.5%
metadata-eval40.5%
sqrt-div36.5%
unpow236.5%
sqrt-prod16.6%
add-sqr-sqrt42.0%
Applied egg-rr42.0%
if 2.7e7 < k Initial program 34.4%
associate-*l*34.4%
associate-/r*34.4%
sub-neg34.4%
distribute-rgt-in34.4%
unpow234.4%
times-frac27.8%
sqr-neg27.8%
times-frac34.4%
unpow234.4%
distribute-rgt-in34.4%
+-commutative34.4%
associate-+l+41.5%
Simplified41.5%
Taylor expanded in t around 0 70.6%
times-frac73.4%
Simplified73.4%
Taylor expanded in k around 0 57.4%
sub-neg57.4%
*-commutative57.4%
associate-/r*57.4%
associate-*r/57.4%
metadata-eval57.4%
distribute-neg-frac57.4%
metadata-eval57.4%
Simplified57.4%
Taylor expanded in l around 0 57.2%
associate-/l*57.0%
associate-/r*57.0%
associate-*r/57.0%
metadata-eval57.0%
div-sub57.0%
associate-/l*57.0%
associate-/r/57.0%
associate-/r*57.4%
unpow257.4%
unpow257.4%
times-frac61.5%
unpow261.5%
sub-neg61.5%
metadata-eval61.5%
+-commutative61.5%
Simplified61.5%
Final simplification46.7%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 27000000.0)
(* 2.0 (pow (/ (* l (pow k_m -2.0)) (sqrt t_m)) 2.0))
(*
2.0
(*
(/ (pow (/ l k_m) 2.0) t_m)
(+ -0.16666666666666666 (/ 1.0 (pow k_m 2.0))))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 27000000.0) {
tmp = 2.0 * pow(((l * pow(k_m, -2.0)) / sqrt(t_m)), 2.0);
} else {
tmp = 2.0 * ((pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / pow(k_m, 2.0))));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 27000000.0d0) then
tmp = 2.0d0 * (((l * (k_m ** (-2.0d0))) / sqrt(t_m)) ** 2.0d0)
else
tmp = 2.0d0 * ((((l / k_m) ** 2.0d0) / t_m) * ((-0.16666666666666666d0) + (1.0d0 / (k_m ** 2.0d0))))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 27000000.0) {
tmp = 2.0 * Math.pow(((l * Math.pow(k_m, -2.0)) / Math.sqrt(t_m)), 2.0);
} else {
tmp = 2.0 * ((Math.pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / Math.pow(k_m, 2.0))));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 27000000.0: tmp = 2.0 * math.pow(((l * math.pow(k_m, -2.0)) / math.sqrt(t_m)), 2.0) else: tmp = 2.0 * ((math.pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / math.pow(k_m, 2.0)))) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 27000000.0) tmp = Float64(2.0 * (Float64(Float64(l * (k_m ^ -2.0)) / sqrt(t_m)) ^ 2.0)); else tmp = Float64(2.0 * Float64(Float64((Float64(l / k_m) ^ 2.0) / t_m) * Float64(-0.16666666666666666 + Float64(1.0 / (k_m ^ 2.0))))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 27000000.0) tmp = 2.0 * (((l * (k_m ^ -2.0)) / sqrt(t_m)) ^ 2.0); else tmp = 2.0 * ((((l / k_m) ^ 2.0) / t_m) * (-0.16666666666666666 + (1.0 / (k_m ^ 2.0)))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 27000000.0], N[(2.0 * N[Power[N[(N[(l * N[Power[k$95$m, -2.0], $MachinePrecision]), $MachinePrecision] / N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(1.0 / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 27000000:\\
\;\;\;\;2 \cdot {\left(\frac{\ell \cdot {k\_m}^{-2}}{\sqrt{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{{\left(\frac{\ell}{k\_m}\right)}^{2}}{t\_m} \cdot \left(-0.16666666666666666 + \frac{1}{{k\_m}^{2}}\right)\right)\\
\end{array}
\end{array}
if k < 2.7e7Initial program 35.3%
associate-*l*35.3%
associate-/r*35.3%
sub-neg35.3%
distribute-rgt-in28.0%
unpow228.0%
times-frac21.8%
sqr-neg21.8%
times-frac28.0%
unpow228.0%
distribute-rgt-in35.3%
+-commutative35.3%
associate-+l+42.4%
Simplified42.4%
Taylor expanded in k around 0 65.1%
*-commutative65.1%
associate-/r*65.5%
Simplified65.5%
add-sqr-sqrt41.5%
pow241.5%
div-inv41.5%
pow-flip41.5%
metadata-eval41.5%
Applied egg-rr41.5%
add-cbrt-cube40.1%
pow1/340.0%
pow340.0%
sqrt-unprod37.3%
sqrt-div32.8%
unpow232.8%
sqrt-prod14.6%
add-sqr-sqrt35.1%
sqrt-pow137.2%
metadata-eval37.2%
Applied egg-rr37.2%
pow-pow42.0%
metadata-eval42.0%
pow142.0%
associate-*l/41.6%
Applied egg-rr41.6%
if 2.7e7 < k Initial program 34.4%
associate-*l*34.4%
associate-/r*34.4%
sub-neg34.4%
distribute-rgt-in34.4%
unpow234.4%
times-frac27.8%
sqr-neg27.8%
times-frac34.4%
unpow234.4%
distribute-rgt-in34.4%
+-commutative34.4%
associate-+l+41.5%
Simplified41.5%
Taylor expanded in t around 0 70.6%
times-frac73.4%
Simplified73.4%
Taylor expanded in k around 0 57.4%
sub-neg57.4%
*-commutative57.4%
associate-/r*57.4%
associate-*r/57.4%
metadata-eval57.4%
distribute-neg-frac57.4%
metadata-eval57.4%
Simplified57.4%
Taylor expanded in l around 0 57.2%
associate-/l*57.0%
associate-/r*57.0%
associate-*r/57.0%
metadata-eval57.0%
div-sub57.0%
associate-/l*57.0%
associate-/r/57.0%
associate-/r*57.4%
unpow257.4%
unpow257.4%
times-frac61.5%
unpow261.5%
sub-neg61.5%
metadata-eval61.5%
+-commutative61.5%
Simplified61.5%
Final simplification46.4%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(*
2.0
(*
(pow (/ l k_m) 2.0)
(/ (- (- 0.16666666666666666) (/ -1.0 (pow k_m 2.0))) t_m)))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (pow((l / k_m), 2.0) * ((-0.16666666666666666 - (-1.0 / pow(k_m, 2.0))) / t_m)));
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 * (((l / k_m) ** 2.0d0) * ((-0.16666666666666666d0 - ((-1.0d0) / (k_m ** 2.0d0))) / t_m)))
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (Math.pow((l / k_m), 2.0) * ((-0.16666666666666666 - (-1.0 / Math.pow(k_m, 2.0))) / t_m)));
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 * (math.pow((l / k_m), 2.0) * ((-0.16666666666666666 - (-1.0 / math.pow(k_m, 2.0))) / t_m)))
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 * Float64((Float64(l / k_m) ^ 2.0) * Float64(Float64(Float64(-0.16666666666666666) - Float64(-1.0 / (k_m ^ 2.0))) / t_m)))) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 * (((l / k_m) ^ 2.0) * ((-0.16666666666666666 - (-1.0 / (k_m ^ 2.0))) / t_m))); end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(2.0 * N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] * N[(N[((-0.16666666666666666) - N[(-1.0 / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(2 \cdot \left({\left(\frac{\ell}{k\_m}\right)}^{2} \cdot \frac{\left(-0.16666666666666666\right) - \frac{-1}{{k\_m}^{2}}}{t\_m}\right)\right)
\end{array}
Initial program 35.1%
associate-*l*35.1%
associate-/r*35.1%
sub-neg35.1%
distribute-rgt-in29.5%
unpow229.5%
times-frac23.2%
sqr-neg23.2%
times-frac29.5%
unpow229.5%
distribute-rgt-in35.1%
+-commutative35.1%
associate-+l+42.2%
Simplified42.2%
Taylor expanded in t around 0 73.2%
times-frac75.7%
Simplified75.7%
Taylor expanded in k around 0 68.0%
sub-neg68.0%
*-commutative68.0%
associate-/r*68.0%
associate-*r/68.0%
metadata-eval68.0%
distribute-neg-frac68.0%
metadata-eval68.0%
Simplified68.0%
Taylor expanded in t around -inf 68.2%
mul-1-neg68.2%
times-frac68.0%
distribute-rgt-neg-in68.0%
unpow268.0%
unpow268.0%
times-frac74.6%
unpow274.6%
sub-neg74.6%
distribute-neg-frac74.6%
metadata-eval74.6%
Simplified74.6%
Final simplification74.6%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(*
2.0
(*
(/ (pow (/ l k_m) 2.0) t_m)
(+ -0.16666666666666666 (/ 1.0 (pow k_m 2.0)))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * ((pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / pow(k_m, 2.0)))));
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 * ((((l / k_m) ** 2.0d0) / t_m) * ((-0.16666666666666666d0) + (1.0d0 / (k_m ** 2.0d0)))))
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * ((Math.pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / Math.pow(k_m, 2.0)))));
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 * ((math.pow((l / k_m), 2.0) / t_m) * (-0.16666666666666666 + (1.0 / math.pow(k_m, 2.0)))))
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 * Float64(Float64((Float64(l / k_m) ^ 2.0) / t_m) * Float64(-0.16666666666666666 + Float64(1.0 / (k_m ^ 2.0)))))) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 * ((((l / k_m) ^ 2.0) / t_m) * (-0.16666666666666666 + (1.0 / (k_m ^ 2.0))))); end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(2.0 * N[(N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(1.0 / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(2 \cdot \left(\frac{{\left(\frac{\ell}{k\_m}\right)}^{2}}{t\_m} \cdot \left(-0.16666666666666666 + \frac{1}{{k\_m}^{2}}\right)\right)\right)
\end{array}
Initial program 35.1%
associate-*l*35.1%
associate-/r*35.1%
sub-neg35.1%
distribute-rgt-in29.5%
unpow229.5%
times-frac23.2%
sqr-neg23.2%
times-frac29.5%
unpow229.5%
distribute-rgt-in35.1%
+-commutative35.1%
associate-+l+42.2%
Simplified42.2%
Taylor expanded in t around 0 73.2%
times-frac75.7%
Simplified75.7%
Taylor expanded in k around 0 68.0%
sub-neg68.0%
*-commutative68.0%
associate-/r*68.0%
associate-*r/68.0%
metadata-eval68.0%
distribute-neg-frac68.0%
metadata-eval68.0%
Simplified68.0%
Taylor expanded in l around 0 67.9%
associate-/l*66.1%
associate-/r*66.1%
associate-*r/66.1%
metadata-eval66.1%
div-sub66.1%
associate-/l*66.1%
associate-/r/68.2%
associate-/r*68.7%
unpow268.7%
unpow268.7%
times-frac74.5%
unpow274.5%
sub-neg74.5%
metadata-eval74.5%
+-commutative74.5%
Simplified74.5%
Final simplification74.5%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(*
2.0
(*
(pow (/ l k_m) 2.0)
(* (/ 1.0 t_m) (+ -0.16666666666666666 (pow k_m -2.0)))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (pow((l / k_m), 2.0) * ((1.0 / t_m) * (-0.16666666666666666 + pow(k_m, -2.0)))));
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 * (((l / k_m) ** 2.0d0) * ((1.0d0 / t_m) * ((-0.16666666666666666d0) + (k_m ** (-2.0d0))))))
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (Math.pow((l / k_m), 2.0) * ((1.0 / t_m) * (-0.16666666666666666 + Math.pow(k_m, -2.0)))));
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 * (math.pow((l / k_m), 2.0) * ((1.0 / t_m) * (-0.16666666666666666 + math.pow(k_m, -2.0)))))
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 * Float64((Float64(l / k_m) ^ 2.0) * Float64(Float64(1.0 / t_m) * Float64(-0.16666666666666666 + (k_m ^ -2.0)))))) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 * (((l / k_m) ^ 2.0) * ((1.0 / t_m) * (-0.16666666666666666 + (k_m ^ -2.0))))); end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(2.0 * N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] * N[(N[(1.0 / t$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[Power[k$95$m, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(2 \cdot \left({\left(\frac{\ell}{k\_m}\right)}^{2} \cdot \left(\frac{1}{t\_m} \cdot \left(-0.16666666666666666 + {k\_m}^{-2}\right)\right)\right)\right)
\end{array}
Initial program 35.1%
associate-*l*35.1%
associate-/r*35.1%
sub-neg35.1%
distribute-rgt-in29.5%
unpow229.5%
times-frac23.2%
sqr-neg23.2%
times-frac29.5%
unpow229.5%
distribute-rgt-in35.1%
+-commutative35.1%
associate-+l+42.2%
Simplified42.2%
Taylor expanded in t around 0 73.2%
times-frac75.7%
Simplified75.7%
Taylor expanded in k around 0 68.0%
sub-neg68.0%
*-commutative68.0%
associate-/r*68.0%
associate-*r/68.0%
metadata-eval68.0%
distribute-neg-frac68.0%
metadata-eval68.0%
Simplified68.0%
associate-*l/67.9%
div-inv67.9%
fma-define67.9%
pow-flip67.9%
metadata-eval67.9%
Applied egg-rr67.9%
associate-/l*66.1%
associate-/r/68.0%
unpow268.0%
unpow268.0%
times-frac74.6%
unpow274.6%
fma-undefine74.6%
metadata-eval74.6%
associate-*l/74.6%
distribute-lft-out74.6%
Simplified74.6%
Final simplification74.6%
k_m = (fabs.f64 k)
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 5e+56)
(* 2.0 (* (/ (pow l 2.0) t_m) (pow k_m -4.0)))
(* 2.0 (* (/ (pow (/ l k_m) 2.0) t_m) -0.16666666666666666)))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 5e+56) {
tmp = 2.0 * ((pow(l, 2.0) / t_m) * pow(k_m, -4.0));
} else {
tmp = 2.0 * ((pow((l / k_m), 2.0) / t_m) * -0.16666666666666666);
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 5d+56) then
tmp = 2.0d0 * (((l ** 2.0d0) / t_m) * (k_m ** (-4.0d0)))
else
tmp = 2.0d0 * ((((l / k_m) ** 2.0d0) / t_m) * (-0.16666666666666666d0))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 5e+56) {
tmp = 2.0 * ((Math.pow(l, 2.0) / t_m) * Math.pow(k_m, -4.0));
} else {
tmp = 2.0 * ((Math.pow((l / k_m), 2.0) / t_m) * -0.16666666666666666);
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 5e+56: tmp = 2.0 * ((math.pow(l, 2.0) / t_m) * math.pow(k_m, -4.0)) else: tmp = 2.0 * ((math.pow((l / k_m), 2.0) / t_m) * -0.16666666666666666) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 5e+56) tmp = Float64(2.0 * Float64(Float64((l ^ 2.0) / t_m) * (k_m ^ -4.0))); else tmp = Float64(2.0 * Float64(Float64((Float64(l / k_m) ^ 2.0) / t_m) * -0.16666666666666666)); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 5e+56) tmp = 2.0 * (((l ^ 2.0) / t_m) * (k_m ^ -4.0)); else tmp = 2.0 * ((((l / k_m) ^ 2.0) / t_m) * -0.16666666666666666); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 5e+56], N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] * N[Power[k$95$m, -4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 5 \cdot 10^{+56}:\\
\;\;\;\;2 \cdot \left(\frac{{\ell}^{2}}{t\_m} \cdot {k\_m}^{-4}\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{{\left(\frac{\ell}{k\_m}\right)}^{2}}{t\_m} \cdot -0.16666666666666666\right)\\
\end{array}
\end{array}
if k < 5.00000000000000024e56Initial program 35.6%
associate-*l*35.6%
associate-/r*35.6%
sub-neg35.6%
distribute-rgt-in28.5%
unpow228.5%
times-frac22.5%
sqr-neg22.5%
times-frac28.5%
unpow228.5%
distribute-rgt-in35.6%
+-commutative35.6%
associate-+l+42.4%
Simplified42.4%
Taylor expanded in k around 0 64.4%
*-commutative64.4%
associate-/r*65.3%
Simplified65.3%
div-inv65.3%
pow-flip65.3%
metadata-eval65.3%
Applied egg-rr65.3%
if 5.00000000000000024e56 < k Initial program 33.3%
associate-*l*33.3%
associate-/r*33.3%
sub-neg33.3%
distribute-rgt-in33.3%
unpow233.3%
times-frac25.9%
sqr-neg25.9%
times-frac33.3%
unpow233.3%
distribute-rgt-in33.3%
+-commutative33.3%
associate-+l+41.4%
Simplified41.4%
Taylor expanded in t around 0 66.9%
times-frac70.0%
Simplified70.0%
Taylor expanded in k around 0 58.9%
sub-neg58.9%
*-commutative58.9%
associate-/r*58.9%
associate-*r/58.9%
metadata-eval58.9%
distribute-neg-frac58.9%
metadata-eval58.9%
Simplified58.9%
Taylor expanded in k around inf 58.4%
*-commutative58.4%
associate-/r*58.9%
unpow258.9%
unpow258.9%
times-frac63.5%
unpow263.5%
Simplified63.5%
Final simplification64.9%
k_m = (fabs.f64 k) t_m = (fabs.f64 t) t_s = (copysign.f64 1 t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* 2.0 (* (/ (pow (/ l k_m) 2.0) t_m) -0.16666666666666666))))
k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * ((pow((l / k_m), 2.0) / t_m) * -0.16666666666666666));
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 * ((((l / k_m) ** 2.0d0) / t_m) * (-0.16666666666666666d0)))
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * ((Math.pow((l / k_m), 2.0) / t_m) * -0.16666666666666666));
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 * ((math.pow((l / k_m), 2.0) / t_m) * -0.16666666666666666))
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 * Float64(Float64((Float64(l / k_m) ^ 2.0) / t_m) * -0.16666666666666666))) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 * ((((l / k_m) ^ 2.0) / t_m) * -0.16666666666666666)); end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(2.0 * N[(N[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(2 \cdot \left(\frac{{\left(\frac{\ell}{k\_m}\right)}^{2}}{t\_m} \cdot -0.16666666666666666\right)\right)
\end{array}
Initial program 35.1%
associate-*l*35.1%
associate-/r*35.1%
sub-neg35.1%
distribute-rgt-in29.5%
unpow229.5%
times-frac23.2%
sqr-neg23.2%
times-frac29.5%
unpow229.5%
distribute-rgt-in35.1%
+-commutative35.1%
associate-+l+42.2%
Simplified42.2%
Taylor expanded in t around 0 73.2%
times-frac75.7%
Simplified75.7%
Taylor expanded in k around 0 68.0%
sub-neg68.0%
*-commutative68.0%
associate-/r*68.0%
associate-*r/68.0%
metadata-eval68.0%
distribute-neg-frac68.0%
metadata-eval68.0%
Simplified68.0%
Taylor expanded in k around inf 31.7%
*-commutative31.7%
associate-/r*31.8%
unpow231.8%
unpow231.8%
times-frac33.4%
unpow233.4%
Simplified33.4%
Final simplification33.4%
herbie shell --seed 2024039
(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))))