
(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 9 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 5.2e-38)
(* 2.0 (pow (* (/ l k_m) (/ (sqrt (/ (cos k_m) t_m)) (sin k_m))) 2.0))
(*
2.0
(/ (* (cos k_m) (* (/ l k_m) (/ l k_m))) (* t_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 <= 5.2e-38) {
tmp = 2.0 * pow(((l / k_m) * (sqrt((cos(k_m) / t_m)) / sin(k_m))), 2.0);
} else {
tmp = 2.0 * ((cos(k_m) * ((l / k_m) * (l / k_m))) / (t_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 <= 5.2d-38) then
tmp = 2.0d0 * (((l / k_m) * (sqrt((cos(k_m) / t_m)) / sin(k_m))) ** 2.0d0)
else
tmp = 2.0d0 * ((cos(k_m) * ((l / k_m) * (l / k_m))) / (t_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 <= 5.2e-38) {
tmp = 2.0 * Math.pow(((l / k_m) * (Math.sqrt((Math.cos(k_m) / t_m)) / Math.sin(k_m))), 2.0);
} else {
tmp = 2.0 * ((Math.cos(k_m) * ((l / k_m) * (l / k_m))) / (t_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 <= 5.2e-38: tmp = 2.0 * math.pow(((l / k_m) * (math.sqrt((math.cos(k_m) / t_m)) / math.sin(k_m))), 2.0) else: tmp = 2.0 * ((math.cos(k_m) * ((l / k_m) * (l / k_m))) / (t_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 <= 5.2e-38) tmp = Float64(2.0 * (Float64(Float64(l / k_m) * Float64(sqrt(Float64(cos(k_m) / t_m)) / sin(k_m))) ^ 2.0)); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) * Float64(Float64(l / k_m) * Float64(l / k_m))) / Float64(t_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 <= 5.2e-38) tmp = 2.0 * (((l / k_m) * (sqrt((cos(k_m) / t_m)) / sin(k_m))) ^ 2.0); else tmp = 2.0 * ((cos(k_m) * ((l / k_m) * (l / k_m))) / (t_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, 5.2e-38], N[(2.0 * N[Power[N[(N[(l / k$95$m), $MachinePrecision] * N[(N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$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[(N[(l / k$95$m), $MachinePrecision] * N[(l / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * 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 5.2 \cdot 10^{-38}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell}{k_m} \cdot \frac{\sqrt{\frac{\cos k_m}{t_m}}}{\sin k_m}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{\cos k_m \cdot \left(\frac{\ell}{k_m} \cdot \frac{\ell}{k_m}\right)}{t_m \cdot {\sin k_m}^{2}}\\
\end{array}
\end{array}
if k < 5.20000000000000022e-38Initial program 39.4%
associate-*l*39.4%
associate-/r*39.4%
sub-neg39.4%
distribute-rgt-in33.6%
unpow233.6%
times-frac23.2%
sqr-neg23.2%
times-frac33.6%
unpow233.6%
distribute-rgt-in39.4%
+-commutative39.4%
associate-+l+44.5%
Simplified44.5%
Taylor expanded in t around 0 76.9%
times-frac78.9%
Simplified78.9%
associate-*l/77.3%
Applied egg-rr77.3%
expm1-log1p-u44.1%
expm1-udef39.2%
Applied egg-rr37.0%
expm1-def42.4%
expm1-log1p43.1%
associate-/l*40.5%
associate-/r/43.2%
Simplified43.2%
if 5.20000000000000022e-38 < k Initial program 29.2%
associate-*l*29.2%
associate-/r*29.2%
sub-neg29.2%
distribute-rgt-in29.2%
unpow229.2%
times-frac25.5%
sqr-neg25.5%
times-frac29.2%
unpow229.2%
distribute-rgt-in29.2%
+-commutative29.2%
associate-+l+46.0%
Simplified46.0%
Taylor expanded in t around 0 73.0%
times-frac73.8%
Simplified73.8%
associate-*r/73.7%
div-inv73.7%
pow-flip73.8%
metadata-eval73.8%
Applied egg-rr73.8%
Taylor expanded in l around 0 73.7%
unpow273.7%
unpow273.7%
times-frac95.1%
*-rgt-identity95.1%
associate-*r/95.1%
*-rgt-identity95.1%
associate-*r/95.0%
unpow295.0%
associate-*r/95.1%
*-rgt-identity95.1%
Simplified95.1%
unpow295.1%
Applied egg-rr95.1%
Final simplification60.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 (<= (* l l) 5e-315)
(* 2.0 (pow (/ (* l (pow k_m -2.0)) (sqrt t_m)) 2.0))
(*
2.0
(/ (* (cos k_m) (* (/ l k_m) (/ l k_m))) (* t_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 ((l * l) <= 5e-315) {
tmp = 2.0 * pow(((l * pow(k_m, -2.0)) / sqrt(t_m)), 2.0);
} else {
tmp = 2.0 * ((cos(k_m) * ((l / k_m) * (l / k_m))) / (t_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 ((l * l) <= 5d-315) then
tmp = 2.0d0 * (((l * (k_m ** (-2.0d0))) / sqrt(t_m)) ** 2.0d0)
else
tmp = 2.0d0 * ((cos(k_m) * ((l / k_m) * (l / k_m))) / (t_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 ((l * l) <= 5e-315) {
tmp = 2.0 * Math.pow(((l * Math.pow(k_m, -2.0)) / Math.sqrt(t_m)), 2.0);
} else {
tmp = 2.0 * ((Math.cos(k_m) * ((l / k_m) * (l / k_m))) / (t_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 (l * l) <= 5e-315: tmp = 2.0 * math.pow(((l * math.pow(k_m, -2.0)) / math.sqrt(t_m)), 2.0) else: tmp = 2.0 * ((math.cos(k_m) * ((l / k_m) * (l / k_m))) / (t_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 (Float64(l * l) <= 5e-315) tmp = Float64(2.0 * (Float64(Float64(l * (k_m ^ -2.0)) / sqrt(t_m)) ^ 2.0)); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) * Float64(Float64(l / k_m) * Float64(l / k_m))) / Float64(t_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 ((l * l) <= 5e-315) tmp = 2.0 * (((l * (k_m ^ -2.0)) / sqrt(t_m)) ^ 2.0); else tmp = 2.0 * ((cos(k_m) * ((l / k_m) * (l / k_m))) / (t_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[N[(l * l), $MachinePrecision], 5e-315], 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[Cos[k$95$m], $MachinePrecision] * N[(N[(l / k$95$m), $MachinePrecision] * N[(l / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * 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}\;\ell \cdot \ell \leq 5 \cdot 10^{-315}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell \cdot {k_m}^{-2}}{\sqrt{t_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{\cos k_m \cdot \left(\frac{\ell}{k_m} \cdot \frac{\ell}{k_m}\right)}{t_m \cdot {\sin k_m}^{2}}\\
\end{array}
\end{array}
if (*.f64 l l) < 5.0000000023e-315Initial program 23.1%
associate-*l*23.1%
associate-/r*23.1%
sub-neg23.1%
distribute-rgt-in23.1%
unpow223.1%
times-frac21.5%
sqr-neg21.5%
times-frac23.1%
unpow223.1%
distribute-rgt-in23.1%
+-commutative23.1%
associate-+l+43.7%
Simplified43.7%
Taylor expanded in k around 0 59.4%
*-commutative59.4%
associate-/r*60.5%
Simplified60.5%
expm1-log1p-u60.5%
expm1-udef60.4%
div-inv60.4%
pow-flip60.4%
metadata-eval60.4%
Applied egg-rr60.4%
expm1-def60.4%
expm1-log1p60.4%
Simplified60.4%
rem-cbrt-cube59.4%
unpow1/359.4%
sqr-pow59.4%
pow259.4%
sqrt-pow159.4%
pow-pow60.4%
metadata-eval60.4%
pow160.4%
*-commutative60.4%
sqrt-prod58.8%
sqrt-pow159.9%
metadata-eval59.9%
sqrt-div33.4%
pow233.4%
sqrt-prod14.9%
add-sqr-sqrt44.5%
Applied egg-rr44.5%
associate-*r/44.5%
Simplified44.5%
if 5.0000000023e-315 < (*.f64 l l) Initial program 40.2%
associate-*l*40.2%
associate-/r*40.2%
sub-neg40.2%
distribute-rgt-in35.0%
unpow235.0%
times-frac24.7%
sqr-neg24.7%
times-frac35.0%
unpow235.0%
distribute-rgt-in40.2%
+-commutative40.2%
associate-+l+45.4%
Simplified45.4%
Taylor expanded in t around 0 80.5%
times-frac82.5%
Simplified82.5%
associate-*r/82.5%
div-inv82.5%
pow-flip82.5%
metadata-eval82.5%
Applied egg-rr82.5%
Taylor expanded in l around 0 82.5%
unpow282.5%
unpow282.5%
times-frac98.1%
*-rgt-identity98.1%
associate-*r/98.2%
*-rgt-identity98.2%
associate-*r/98.1%
unpow298.1%
associate-*r/98.1%
*-rgt-identity98.1%
Simplified98.1%
unpow298.1%
Applied egg-rr98.1%
Final simplification85.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
(if (<= k_m 1.1e+74)
(* 2.0 (pow (* (pow k_m -2.0) (/ l (sqrt t_m))) 2.0))
(* 2.0 (* -0.16666666666666666 (/ (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 <= 1.1e+74) {
tmp = 2.0 * pow((pow(k_m, -2.0) * (l / sqrt(t_m))), 2.0);
} else {
tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1d+74) then
tmp = 2.0d0 * (((k_m ** (-2.0d0)) * (l / sqrt(t_m))) ** 2.0d0)
else
tmp = 2.0d0 * ((-0.16666666666666666d0) * (((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 <= 1.1e+74) {
tmp = 2.0 * Math.pow((Math.pow(k_m, -2.0) * (l / Math.sqrt(t_m))), 2.0);
} else {
tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1e+74: tmp = 2.0 * math.pow((math.pow(k_m, -2.0) * (l / math.sqrt(t_m))), 2.0) else: tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1e+74) tmp = Float64(2.0 * (Float64((k_m ^ -2.0) * Float64(l / sqrt(t_m))) ^ 2.0)); else tmp = Float64(2.0 * Float64(-0.16666666666666666 * 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 <= 1.1e+74) tmp = 2.0 * (((k_m ^ -2.0) * (l / sqrt(t_m))) ^ 2.0); else tmp = 2.0 * (-0.16666666666666666 * (((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, 1.1e+74], 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[(-0.16666666666666666 * 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 1.1 \cdot 10^{+74}:\\
\;\;\;\;2 \cdot {\left({k_m}^{-2} \cdot \frac{\ell}{\sqrt{t_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(-0.16666666666666666 \cdot \frac{{\left(\frac{\ell}{k_m}\right)}^{2}}{t_m}\right)\\
\end{array}
\end{array}
if k < 1.1000000000000001e74Initial program 37.4%
associate-*l*37.4%
associate-/r*37.4%
sub-neg37.4%
distribute-rgt-in32.4%
unpow232.4%
times-frac23.5%
sqr-neg23.5%
times-frac32.4%
unpow232.4%
distribute-rgt-in37.4%
+-commutative37.4%
associate-+l+44.8%
Simplified44.8%
Taylor expanded in k around 0 67.4%
*-commutative67.4%
associate-/r*67.5%
Simplified67.5%
expm1-log1p-u40.1%
expm1-udef38.6%
div-inv38.5%
pow-flip38.5%
metadata-eval38.5%
Applied egg-rr38.5%
expm1-def40.0%
expm1-log1p67.5%
Simplified67.5%
rem-cbrt-cube64.8%
unpow1/339.1%
sqr-pow39.1%
pow239.1%
sqrt-pow139.1%
pow-pow38.8%
metadata-eval38.8%
pow138.8%
*-commutative38.8%
sqrt-prod34.7%
sqrt-pow135.6%
metadata-eval35.6%
sqrt-div29.5%
pow229.5%
sqrt-prod15.0%
add-sqr-sqrt32.9%
Applied egg-rr32.9%
if 1.1000000000000001e74 < k Initial program 31.3%
associate-*l*31.3%
associate-/r*31.3%
sub-neg31.3%
distribute-rgt-in31.3%
unpow231.3%
times-frac25.7%
sqr-neg25.7%
times-frac31.3%
unpow231.3%
distribute-rgt-in31.3%
+-commutative31.3%
associate-+l+45.8%
Simplified45.8%
Taylor expanded in t around 0 64.0%
times-frac62.4%
Simplified62.4%
associate-*l/62.3%
Applied egg-rr62.3%
Taylor expanded in k around 0 55.6%
+-commutative55.6%
associate--l+55.6%
associate-/r*55.6%
unpow255.6%
unpow255.6%
times-frac55.6%
*-rgt-identity55.6%
associate-*r/55.6%
*-rgt-identity55.6%
associate-*r/55.6%
unpow255.6%
associate-*r/55.6%
*-rgt-identity55.6%
distribute-rgt-out--56.0%
metadata-eval56.0%
associate-*l/56.0%
Simplified56.0%
Taylor expanded in k around inf 56.3%
associate-/r*56.5%
unpow256.5%
unpow256.5%
times-frac66.1%
unpow266.1%
Simplified66.1%
Final simplification40.1%
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 1.1e+74)
(* 2.0 (pow (/ (* l (pow k_m -2.0)) (sqrt t_m)) 2.0))
(* 2.0 (* -0.16666666666666666 (/ (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 <= 1.1e+74) {
tmp = 2.0 * pow(((l * pow(k_m, -2.0)) / sqrt(t_m)), 2.0);
} else {
tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1d+74) then
tmp = 2.0d0 * (((l * (k_m ** (-2.0d0))) / sqrt(t_m)) ** 2.0d0)
else
tmp = 2.0d0 * ((-0.16666666666666666d0) * (((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 <= 1.1e+74) {
tmp = 2.0 * Math.pow(((l * Math.pow(k_m, -2.0)) / Math.sqrt(t_m)), 2.0);
} else {
tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1e+74: tmp = 2.0 * math.pow(((l * math.pow(k_m, -2.0)) / math.sqrt(t_m)), 2.0) else: tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1e+74) tmp = Float64(2.0 * (Float64(Float64(l * (k_m ^ -2.0)) / sqrt(t_m)) ^ 2.0)); else tmp = Float64(2.0 * Float64(-0.16666666666666666 * 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 <= 1.1e+74) tmp = 2.0 * (((l * (k_m ^ -2.0)) / sqrt(t_m)) ^ 2.0); else tmp = 2.0 * (-0.16666666666666666 * (((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, 1.1e+74], 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[(-0.16666666666666666 * 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 1.1 \cdot 10^{+74}:\\
\;\;\;\;2 \cdot {\left(\frac{\ell \cdot {k_m}^{-2}}{\sqrt{t_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(-0.16666666666666666 \cdot \frac{{\left(\frac{\ell}{k_m}\right)}^{2}}{t_m}\right)\\
\end{array}
\end{array}
if k < 1.1000000000000001e74Initial program 37.4%
associate-*l*37.4%
associate-/r*37.4%
sub-neg37.4%
distribute-rgt-in32.4%
unpow232.4%
times-frac23.5%
sqr-neg23.5%
times-frac32.4%
unpow232.4%
distribute-rgt-in37.4%
+-commutative37.4%
associate-+l+44.8%
Simplified44.8%
Taylor expanded in k around 0 67.4%
*-commutative67.4%
associate-/r*67.5%
Simplified67.5%
expm1-log1p-u40.1%
expm1-udef38.6%
div-inv38.5%
pow-flip38.5%
metadata-eval38.5%
Applied egg-rr38.5%
expm1-def40.0%
expm1-log1p67.5%
Simplified67.5%
rem-cbrt-cube64.8%
unpow1/339.1%
sqr-pow39.1%
pow239.1%
sqrt-pow139.1%
pow-pow38.8%
metadata-eval38.8%
pow138.8%
*-commutative38.8%
sqrt-prod34.7%
sqrt-pow135.6%
metadata-eval35.6%
sqrt-div29.5%
pow229.5%
sqrt-prod15.0%
add-sqr-sqrt32.9%
Applied egg-rr32.9%
associate-*r/32.9%
Simplified32.9%
if 1.1000000000000001e74 < k Initial program 31.3%
associate-*l*31.3%
associate-/r*31.3%
sub-neg31.3%
distribute-rgt-in31.3%
unpow231.3%
times-frac25.7%
sqr-neg25.7%
times-frac31.3%
unpow231.3%
distribute-rgt-in31.3%
+-commutative31.3%
associate-+l+45.8%
Simplified45.8%
Taylor expanded in t around 0 64.0%
times-frac62.4%
Simplified62.4%
associate-*l/62.3%
Applied egg-rr62.3%
Taylor expanded in k around 0 55.6%
+-commutative55.6%
associate--l+55.6%
associate-/r*55.6%
unpow255.6%
unpow255.6%
times-frac55.6%
*-rgt-identity55.6%
associate-*r/55.6%
*-rgt-identity55.6%
associate-*r/55.6%
unpow255.6%
associate-*r/55.6%
*-rgt-identity55.6%
distribute-rgt-out--56.0%
metadata-eval56.0%
associate-*l/56.0%
Simplified56.0%
Taylor expanded in k around inf 56.3%
associate-/r*56.5%
unpow256.5%
unpow256.5%
times-frac66.1%
unpow266.1%
Simplified66.1%
Final simplification40.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 1.1e+74)
(* 2.0 (/ (* (pow l 2.0) (/ 1.0 t_m)) (pow k_m 4.0)))
(* 2.0 (* -0.16666666666666666 (/ (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 <= 1.1e+74) {
tmp = 2.0 * ((pow(l, 2.0) * (1.0 / t_m)) / pow(k_m, 4.0));
} else {
tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1d+74) then
tmp = 2.0d0 * (((l ** 2.0d0) * (1.0d0 / t_m)) / (k_m ** 4.0d0))
else
tmp = 2.0d0 * ((-0.16666666666666666d0) * (((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 <= 1.1e+74) {
tmp = 2.0 * ((Math.pow(l, 2.0) * (1.0 / t_m)) / Math.pow(k_m, 4.0));
} else {
tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1e+74: tmp = 2.0 * ((math.pow(l, 2.0) * (1.0 / t_m)) / math.pow(k_m, 4.0)) else: tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1e+74) tmp = Float64(2.0 * Float64(Float64((l ^ 2.0) * Float64(1.0 / t_m)) / (k_m ^ 4.0))); else tmp = Float64(2.0 * Float64(-0.16666666666666666 * 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 <= 1.1e+74) tmp = 2.0 * (((l ^ 2.0) * (1.0 / t_m)) / (k_m ^ 4.0)); else tmp = 2.0 * (-0.16666666666666666 * (((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, 1.1e+74], N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] * N[(1.0 / t$95$m), $MachinePrecision]), $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(-0.16666666666666666 * 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 1.1 \cdot 10^{+74}:\\
\;\;\;\;2 \cdot \frac{{\ell}^{2} \cdot \frac{1}{t_m}}{{k_m}^{4}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(-0.16666666666666666 \cdot \frac{{\left(\frac{\ell}{k_m}\right)}^{2}}{t_m}\right)\\
\end{array}
\end{array}
if k < 1.1000000000000001e74Initial program 37.4%
associate-*l*37.4%
associate-/r*37.4%
sub-neg37.4%
distribute-rgt-in32.4%
unpow232.4%
times-frac23.5%
sqr-neg23.5%
times-frac32.4%
unpow232.4%
distribute-rgt-in37.4%
+-commutative37.4%
associate-+l+44.8%
Simplified44.8%
Taylor expanded in k around 0 67.4%
*-commutative67.4%
associate-/r*67.5%
Simplified67.5%
div-inv67.5%
Applied egg-rr67.5%
if 1.1000000000000001e74 < k Initial program 31.3%
associate-*l*31.3%
associate-/r*31.3%
sub-neg31.3%
distribute-rgt-in31.3%
unpow231.3%
times-frac25.7%
sqr-neg25.7%
times-frac31.3%
unpow231.3%
distribute-rgt-in31.3%
+-commutative31.3%
associate-+l+45.8%
Simplified45.8%
Taylor expanded in t around 0 64.0%
times-frac62.4%
Simplified62.4%
associate-*l/62.3%
Applied egg-rr62.3%
Taylor expanded in k around 0 55.6%
+-commutative55.6%
associate--l+55.6%
associate-/r*55.6%
unpow255.6%
unpow255.6%
times-frac55.6%
*-rgt-identity55.6%
associate-*r/55.6%
*-rgt-identity55.6%
associate-*r/55.6%
unpow255.6%
associate-*r/55.6%
*-rgt-identity55.6%
distribute-rgt-out--56.0%
metadata-eval56.0%
associate-*l/56.0%
Simplified56.0%
Taylor expanded in k around inf 56.3%
associate-/r*56.5%
unpow256.5%
unpow256.5%
times-frac66.1%
unpow266.1%
Simplified66.1%
Final simplification67.2%
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 (/ (pow (/ l k_m) 2.0) t_m)))
(*
t_s
(if (<= k_m 1.1e+74)
(* 2.0 (/ t_2 (pow k_m 2.0)))
(* 2.0 (* -0.16666666666666666 t_2))))))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 = pow((l / k_m), 2.0) / t_m;
double tmp;
if (k_m <= 1.1e+74) {
tmp = 2.0 * (t_2 / pow(k_m, 2.0));
} else {
tmp = 2.0 * (-0.16666666666666666 * t_2);
}
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) ** 2.0d0) / t_m
if (k_m <= 1.1d+74) then
tmp = 2.0d0 * (t_2 / (k_m ** 2.0d0))
else
tmp = 2.0d0 * ((-0.16666666666666666d0) * t_2)
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 = Math.pow((l / k_m), 2.0) / t_m;
double tmp;
if (k_m <= 1.1e+74) {
tmp = 2.0 * (t_2 / Math.pow(k_m, 2.0));
} else {
tmp = 2.0 * (-0.16666666666666666 * t_2);
}
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 = math.pow((l / k_m), 2.0) / t_m tmp = 0 if k_m <= 1.1e+74: tmp = 2.0 * (t_2 / math.pow(k_m, 2.0)) else: tmp = 2.0 * (-0.16666666666666666 * t_2) 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((Float64(l / k_m) ^ 2.0) / t_m) tmp = 0.0 if (k_m <= 1.1e+74) tmp = Float64(2.0 * Float64(t_2 / (k_m ^ 2.0))); else tmp = Float64(2.0 * Float64(-0.16666666666666666 * t_2)); 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) ^ 2.0) / t_m; tmp = 0.0; if (k_m <= 1.1e+74) tmp = 2.0 * (t_2 / (k_m ^ 2.0)); else tmp = 2.0 * (-0.16666666666666666 * t_2); 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[(N[Power[N[(l / k$95$m), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 1.1e+74], N[(2.0 * N[(t$95$2 / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(-0.16666666666666666 * t$95$2), $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{{\left(\frac{\ell}{k_m}\right)}^{2}}{t_m}\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;k_m \leq 1.1 \cdot 10^{+74}:\\
\;\;\;\;2 \cdot \frac{t_2}{{k_m}^{2}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(-0.16666666666666666 \cdot t_2\right)\\
\end{array}
\end{array}
\end{array}
if k < 1.1000000000000001e74Initial program 37.4%
associate-*l*37.4%
associate-/r*37.4%
sub-neg37.4%
distribute-rgt-in32.4%
unpow232.4%
times-frac23.5%
sqr-neg23.5%
times-frac32.4%
unpow232.4%
distribute-rgt-in37.4%
+-commutative37.4%
associate-+l+44.8%
Simplified44.8%
Taylor expanded in t around 0 78.9%
times-frac81.3%
Simplified81.3%
associate-*l/80.1%
Applied egg-rr80.1%
Taylor expanded in k around 0 69.0%
associate-/r*69.4%
unpow269.4%
unpow269.4%
times-frac76.4%
*-rgt-identity76.4%
associate-*r/76.4%
*-rgt-identity76.4%
associate-*r/76.4%
unpow276.4%
associate-*r/76.4%
*-rgt-identity76.4%
Simplified76.4%
if 1.1000000000000001e74 < k Initial program 31.3%
associate-*l*31.3%
associate-/r*31.3%
sub-neg31.3%
distribute-rgt-in31.3%
unpow231.3%
times-frac25.7%
sqr-neg25.7%
times-frac31.3%
unpow231.3%
distribute-rgt-in31.3%
+-commutative31.3%
associate-+l+45.8%
Simplified45.8%
Taylor expanded in t around 0 64.0%
times-frac62.4%
Simplified62.4%
associate-*l/62.3%
Applied egg-rr62.3%
Taylor expanded in k around 0 55.6%
+-commutative55.6%
associate--l+55.6%
associate-/r*55.6%
unpow255.6%
unpow255.6%
times-frac55.6%
*-rgt-identity55.6%
associate-*r/55.6%
*-rgt-identity55.6%
associate-*r/55.6%
unpow255.6%
associate-*r/55.6%
*-rgt-identity55.6%
distribute-rgt-out--56.0%
metadata-eval56.0%
associate-*l/56.0%
Simplified56.0%
Taylor expanded in k around inf 56.3%
associate-/r*56.5%
unpow256.5%
unpow256.5%
times-frac66.1%
unpow266.1%
Simplified66.1%
Final simplification74.2%
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 1.1e+74)
(* 2.0 (* (/ (pow l 2.0) t_m) (pow k_m -4.0)))
(* 2.0 (* -0.16666666666666666 (/ (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 <= 1.1e+74) {
tmp = 2.0 * ((pow(l, 2.0) / t_m) * pow(k_m, -4.0));
} else {
tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1d+74) then
tmp = 2.0d0 * (((l ** 2.0d0) / t_m) * (k_m ** (-4.0d0)))
else
tmp = 2.0d0 * ((-0.16666666666666666d0) * (((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 <= 1.1e+74) {
tmp = 2.0 * ((Math.pow(l, 2.0) / t_m) * Math.pow(k_m, -4.0));
} else {
tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1e+74: tmp = 2.0 * ((math.pow(l, 2.0) / t_m) * math.pow(k_m, -4.0)) else: tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1e+74) tmp = Float64(2.0 * Float64(Float64((l ^ 2.0) / t_m) * (k_m ^ -4.0))); else tmp = Float64(2.0 * Float64(-0.16666666666666666 * 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 <= 1.1e+74) tmp = 2.0 * (((l ^ 2.0) / t_m) * (k_m ^ -4.0)); else tmp = 2.0 * (-0.16666666666666666 * (((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, 1.1e+74], 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[(-0.16666666666666666 * 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 1.1 \cdot 10^{+74}:\\
\;\;\;\;2 \cdot \left(\frac{{\ell}^{2}}{t_m} \cdot {k_m}^{-4}\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(-0.16666666666666666 \cdot \frac{{\left(\frac{\ell}{k_m}\right)}^{2}}{t_m}\right)\\
\end{array}
\end{array}
if k < 1.1000000000000001e74Initial program 37.4%
associate-*l*37.4%
associate-/r*37.4%
sub-neg37.4%
distribute-rgt-in32.4%
unpow232.4%
times-frac23.5%
sqr-neg23.5%
times-frac32.4%
unpow232.4%
distribute-rgt-in37.4%
+-commutative37.4%
associate-+l+44.8%
Simplified44.8%
Taylor expanded in k around 0 67.4%
*-commutative67.4%
associate-/r*67.5%
Simplified67.5%
expm1-log1p-u40.1%
expm1-udef38.6%
div-inv38.5%
pow-flip38.5%
metadata-eval38.5%
Applied egg-rr38.5%
expm1-def40.0%
expm1-log1p67.5%
Simplified67.5%
if 1.1000000000000001e74 < k Initial program 31.3%
associate-*l*31.3%
associate-/r*31.3%
sub-neg31.3%
distribute-rgt-in31.3%
unpow231.3%
times-frac25.7%
sqr-neg25.7%
times-frac31.3%
unpow231.3%
distribute-rgt-in31.3%
+-commutative31.3%
associate-+l+45.8%
Simplified45.8%
Taylor expanded in t around 0 64.0%
times-frac62.4%
Simplified62.4%
associate-*l/62.3%
Applied egg-rr62.3%
Taylor expanded in k around 0 55.6%
+-commutative55.6%
associate--l+55.6%
associate-/r*55.6%
unpow255.6%
unpow255.6%
times-frac55.6%
*-rgt-identity55.6%
associate-*r/55.6%
*-rgt-identity55.6%
associate-*r/55.6%
unpow255.6%
associate-*r/55.6%
*-rgt-identity55.6%
distribute-rgt-out--56.0%
metadata-eval56.0%
associate-*l/56.0%
Simplified56.0%
Taylor expanded in k around inf 56.3%
associate-/r*56.5%
unpow256.5%
unpow256.5%
times-frac66.1%
unpow266.1%
Simplified66.1%
Final simplification67.2%
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 1.1e+74)
(* 2.0 (/ (/ (pow l 2.0) t_m) (pow k_m 4.0)))
(* 2.0 (* -0.16666666666666666 (/ (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 <= 1.1e+74) {
tmp = 2.0 * ((pow(l, 2.0) / t_m) / pow(k_m, 4.0));
} else {
tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1d+74) then
tmp = 2.0d0 * (((l ** 2.0d0) / t_m) / (k_m ** 4.0d0))
else
tmp = 2.0d0 * ((-0.16666666666666666d0) * (((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 <= 1.1e+74) {
tmp = 2.0 * ((Math.pow(l, 2.0) / t_m) / Math.pow(k_m, 4.0));
} else {
tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1e+74: tmp = 2.0 * ((math.pow(l, 2.0) / t_m) / math.pow(k_m, 4.0)) else: tmp = 2.0 * (-0.16666666666666666 * (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 <= 1.1e+74) tmp = Float64(2.0 * Float64(Float64((l ^ 2.0) / t_m) / (k_m ^ 4.0))); else tmp = Float64(2.0 * Float64(-0.16666666666666666 * 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 <= 1.1e+74) tmp = 2.0 * (((l ^ 2.0) / t_m) / (k_m ^ 4.0)); else tmp = 2.0 * (-0.16666666666666666 * (((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, 1.1e+74], 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[(-0.16666666666666666 * 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 1.1 \cdot 10^{+74}:\\
\;\;\;\;2 \cdot \frac{\frac{{\ell}^{2}}{t_m}}{{k_m}^{4}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(-0.16666666666666666 \cdot \frac{{\left(\frac{\ell}{k_m}\right)}^{2}}{t_m}\right)\\
\end{array}
\end{array}
if k < 1.1000000000000001e74Initial program 37.4%
associate-*l*37.4%
associate-/r*37.4%
sub-neg37.4%
distribute-rgt-in32.4%
unpow232.4%
times-frac23.5%
sqr-neg23.5%
times-frac32.4%
unpow232.4%
distribute-rgt-in37.4%
+-commutative37.4%
associate-+l+44.8%
Simplified44.8%
Taylor expanded in k around 0 67.4%
*-commutative67.4%
associate-/r*67.5%
Simplified67.5%
if 1.1000000000000001e74 < k Initial program 31.3%
associate-*l*31.3%
associate-/r*31.3%
sub-neg31.3%
distribute-rgt-in31.3%
unpow231.3%
times-frac25.7%
sqr-neg25.7%
times-frac31.3%
unpow231.3%
distribute-rgt-in31.3%
+-commutative31.3%
associate-+l+45.8%
Simplified45.8%
Taylor expanded in t around 0 64.0%
times-frac62.4%
Simplified62.4%
associate-*l/62.3%
Applied egg-rr62.3%
Taylor expanded in k around 0 55.6%
+-commutative55.6%
associate--l+55.6%
associate-/r*55.6%
unpow255.6%
unpow255.6%
times-frac55.6%
*-rgt-identity55.6%
associate-*r/55.6%
*-rgt-identity55.6%
associate-*r/55.6%
unpow255.6%
associate-*r/55.6%
*-rgt-identity55.6%
distribute-rgt-out--56.0%
metadata-eval56.0%
associate-*l/56.0%
Simplified56.0%
Taylor expanded in k around inf 56.3%
associate-/r*56.5%
unpow256.5%
unpow256.5%
times-frac66.1%
unpow266.1%
Simplified66.1%
Final simplification67.2%
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 (* -0.16666666666666666 (/ (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) {
return t_s * (2.0 * (-0.16666666666666666 * (pow((l / 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 * ((-0.16666666666666666d0) * (((l / 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 * (-0.16666666666666666 * (Math.pow((l / 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 * (-0.16666666666666666 * (math.pow((l / 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(-0.16666666666666666 * Float64((Float64(l / 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 * (-0.16666666666666666 * (((l / 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[(-0.16666666666666666 * 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 \left(2 \cdot \left(-0.16666666666666666 \cdot \frac{{\left(\frac{\ell}{k_m}\right)}^{2}}{t_m}\right)\right)
\end{array}
Initial program 36.1%
associate-*l*36.1%
associate-/r*36.1%
sub-neg36.1%
distribute-rgt-in32.2%
unpow232.2%
times-frac24.0%
sqr-neg24.0%
times-frac32.2%
unpow232.2%
distribute-rgt-in36.1%
+-commutative36.1%
associate-+l+45.0%
Simplified45.0%
Taylor expanded in t around 0 75.7%
times-frac77.2%
Simplified77.2%
associate-*l/76.3%
Applied egg-rr76.3%
Taylor expanded in k around 0 47.7%
+-commutative47.7%
associate--l+47.7%
associate-/r*48.0%
unpow248.0%
unpow248.0%
times-frac53.4%
*-rgt-identity53.4%
associate-*r/53.4%
*-rgt-identity53.4%
associate-*r/53.4%
unpow253.4%
associate-*r/53.4%
*-rgt-identity53.4%
distribute-rgt-out--53.5%
metadata-eval53.5%
associate-*l/53.5%
Simplified53.5%
Taylor expanded in k around inf 31.2%
associate-/r*31.3%
unpow231.3%
unpow231.3%
times-frac34.0%
unpow234.0%
Simplified34.0%
Final simplification34.0%
herbie shell --seed 2024019
(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))))