
(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 17 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 2.5e-43)
(/ 2.0 (pow (* (/ k_m l) (* k_m (sqrt t_m))) 2.0))
(/ (pow (* (sqrt (/ 2.0 t_m)) (/ l k_m)) 2.0) (* (tan k_m) (sin k_m))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.5e-43) {
tmp = 2.0 / pow(((k_m / l) * (k_m * sqrt(t_m))), 2.0);
} else {
tmp = pow((sqrt((2.0 / t_m)) * (l / k_m)), 2.0) / (tan(k_m) * sin(k_m));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 2.5d-43) then
tmp = 2.0d0 / (((k_m / l) * (k_m * sqrt(t_m))) ** 2.0d0)
else
tmp = ((sqrt((2.0d0 / t_m)) * (l / k_m)) ** 2.0d0) / (tan(k_m) * sin(k_m))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.5e-43) {
tmp = 2.0 / Math.pow(((k_m / l) * (k_m * Math.sqrt(t_m))), 2.0);
} else {
tmp = Math.pow((Math.sqrt((2.0 / t_m)) * (l / k_m)), 2.0) / (Math.tan(k_m) * Math.sin(k_m));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 2.5e-43: tmp = 2.0 / math.pow(((k_m / l) * (k_m * math.sqrt(t_m))), 2.0) else: tmp = math.pow((math.sqrt((2.0 / t_m)) * (l / k_m)), 2.0) / (math.tan(k_m) * math.sin(k_m)) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 2.5e-43) tmp = Float64(2.0 / (Float64(Float64(k_m / l) * Float64(k_m * sqrt(t_m))) ^ 2.0)); else tmp = Float64((Float64(sqrt(Float64(2.0 / t_m)) * Float64(l / k_m)) ^ 2.0) / Float64(tan(k_m) * sin(k_m))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 2.5e-43) tmp = 2.0 / (((k_m / l) * (k_m * sqrt(t_m))) ^ 2.0); else tmp = ((sqrt((2.0 / t_m)) * (l / k_m)) ^ 2.0) / (tan(k_m) * sin(k_m)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 2.5e-43], N[(2.0 / N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[Power[N[(N[Sqrt[N[(2.0 / t$95$m), $MachinePrecision]], $MachinePrecision] * N[(l / k$95$m), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] / N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$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 2.5 \cdot 10^{-43}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{\ell} \cdot \left(k\_m \cdot \sqrt{t\_m}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left(\sqrt{\frac{2}{t\_m}} \cdot \frac{\ell}{k\_m}\right)}^{2}}{\tan k\_m \cdot \sin k\_m}\\
\end{array}
\end{array}
if k < 2.50000000000000009e-43Initial program 38.4%
Simplified44.6%
Taylor expanded in t around 0 73.2%
times-frac72.4%
Simplified72.4%
associate-*l/73.2%
associate-/l*73.2%
unpow273.2%
*-un-lft-identity73.2%
times-frac73.2%
tan-quot73.2%
Applied egg-rr73.2%
*-commutative73.2%
associate-/l*72.4%
unpow272.4%
associate-/r*84.1%
unpow284.1%
associate-*r/87.3%
associate-*l/89.3%
unpow289.3%
*-commutative89.3%
/-rgt-identity89.3%
Simplified89.3%
Taylor expanded in k around 0 76.8%
add-sqr-sqrt45.6%
pow245.6%
sqrt-prod45.6%
unpow245.6%
sqrt-prod26.2%
add-sqr-sqrt45.6%
*-commutative45.6%
sqrt-prod45.7%
unpow245.7%
sqrt-prod18.0%
add-sqr-sqrt47.7%
Applied egg-rr47.7%
if 2.50000000000000009e-43 < k Initial program 30.8%
Simplified47.3%
Taylor expanded in t around 0 72.3%
times-frac72.4%
Simplified72.4%
div-inv72.4%
add-sqr-sqrt72.3%
pow272.3%
sqrt-div72.4%
unpow272.4%
sqrt-prod76.4%
add-sqr-sqrt76.6%
unpow276.6%
sqrt-prod51.9%
add-sqr-sqrt93.6%
associate-/l*93.6%
unpow293.6%
*-un-lft-identity93.6%
times-frac93.6%
tan-quot93.7%
Applied egg-rr93.7%
associate-*r/93.7%
metadata-eval93.7%
associate-/r*93.6%
/-rgt-identity93.6%
associate-/r*93.6%
Simplified93.6%
add-sqr-sqrt73.2%
sqrt-div44.4%
sqrt-div44.4%
unpow244.4%
sqrt-prod28.7%
add-sqr-sqrt40.4%
sqrt-div40.4%
sqrt-div40.4%
unpow240.4%
sqrt-prod31.3%
add-sqr-sqrt49.1%
Applied egg-rr49.1%
unpow249.1%
associate-/l/48.9%
Simplified48.9%
*-un-lft-identity48.9%
associate-/r*49.0%
sqrt-undiv49.0%
Applied egg-rr49.0%
associate-*r/49.0%
associate-*l/49.0%
*-commutative49.0%
associate-/r/49.0%
associate-*l/48.9%
*-lft-identity48.9%
Simplified48.9%
Final simplification48.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 (let* ((t_2 (/ (sqrt (/ 2.0 t_m)) (/ k_m l)))) (* t_s (* (/ t_2 (tan k_m)) (/ t_2 (sin k_m))))))
k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = sqrt((2.0 / t_m)) / (k_m / l);
return t_s * ((t_2 / tan(k_m)) * (t_2 / sin(k_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
real(8) :: t_2
t_2 = sqrt((2.0d0 / t_m)) / (k_m / l)
code = t_s * ((t_2 / tan(k_m)) * (t_2 / sin(k_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) {
double t_2 = Math.sqrt((2.0 / t_m)) / (k_m / l);
return t_s * ((t_2 / Math.tan(k_m)) * (t_2 / Math.sin(k_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): t_2 = math.sqrt((2.0 / t_m)) / (k_m / l) return t_s * ((t_2 / math.tan(k_m)) * (t_2 / math.sin(k_m)))
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(sqrt(Float64(2.0 / t_m)) / Float64(k_m / l)) return Float64(t_s * Float64(Float64(t_2 / tan(k_m)) * Float64(t_2 / sin(k_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) t_2 = sqrt((2.0 / t_m)) / (k_m / l); tmp = t_s * ((t_2 / tan(k_m)) * (t_2 / sin(k_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_] := Block[{t$95$2 = N[(N[Sqrt[N[(2.0 / t$95$m), $MachinePrecision]], $MachinePrecision] / N[(k$95$m / l), $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * N[(N[(t$95$2 / N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(t$95$2 / N[Sin[k$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)
\\
\begin{array}{l}
t_2 := \frac{\sqrt{\frac{2}{t\_m}}}{\frac{k\_m}{\ell}}\\
t\_s \cdot \left(\frac{t\_2}{\tan k\_m} \cdot \frac{t\_2}{\sin k\_m}\right)
\end{array}
\end{array}
Initial program 36.2%
Simplified45.3%
Taylor expanded in t around 0 73.0%
times-frac72.4%
Simplified72.4%
div-inv72.4%
add-sqr-sqrt72.3%
pow272.3%
sqrt-div72.4%
unpow272.4%
sqrt-prod41.7%
add-sqr-sqrt75.2%
unpow275.2%
sqrt-prod50.8%
add-sqr-sqrt90.6%
associate-/l*90.6%
unpow290.6%
*-un-lft-identity90.6%
times-frac90.5%
tan-quot90.6%
Applied egg-rr90.6%
associate-*r/90.6%
metadata-eval90.6%
associate-/r*90.7%
/-rgt-identity90.7%
associate-/r*90.8%
Simplified90.8%
add-sqr-sqrt63.4%
sqrt-div51.2%
sqrt-div51.2%
unpow251.2%
sqrt-prod28.7%
add-sqr-sqrt31.4%
sqrt-div31.4%
sqrt-div31.4%
unpow231.4%
sqrt-prod26.6%
add-sqr-sqrt54.3%
Applied egg-rr54.3%
unpow254.3%
associate-/l/54.3%
Simplified54.3%
unpow254.3%
*-commutative54.3%
times-frac55.6%
associate-/r*55.6%
sqrt-undiv55.6%
associate-/r*55.6%
sqrt-undiv55.6%
Applied egg-rr55.6%
Final simplification55.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 6.8e-48)
(/ 2.0 (pow (* (/ k_m l) (* k_m (sqrt t_m))) 2.0))
(/ 2.0 (* (* (tan k_m) (sin k_m)) (pow (* (/ k_m l) (sqrt t_m)) 2.0))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 6.8e-48) {
tmp = 2.0 / pow(((k_m / l) * (k_m * sqrt(t_m))), 2.0);
} else {
tmp = 2.0 / ((tan(k_m) * sin(k_m)) * pow(((k_m / l) * sqrt(t_m)), 2.0));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 6.8d-48) then
tmp = 2.0d0 / (((k_m / l) * (k_m * sqrt(t_m))) ** 2.0d0)
else
tmp = 2.0d0 / ((tan(k_m) * sin(k_m)) * (((k_m / l) * sqrt(t_m)) ** 2.0d0))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 6.8e-48) {
tmp = 2.0 / Math.pow(((k_m / l) * (k_m * Math.sqrt(t_m))), 2.0);
} else {
tmp = 2.0 / ((Math.tan(k_m) * Math.sin(k_m)) * Math.pow(((k_m / l) * Math.sqrt(t_m)), 2.0));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 6.8e-48: tmp = 2.0 / math.pow(((k_m / l) * (k_m * math.sqrt(t_m))), 2.0) else: tmp = 2.0 / ((math.tan(k_m) * math.sin(k_m)) * math.pow(((k_m / l) * math.sqrt(t_m)), 2.0)) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 6.8e-48) tmp = Float64(2.0 / (Float64(Float64(k_m / l) * Float64(k_m * sqrt(t_m))) ^ 2.0)); else tmp = Float64(2.0 / Float64(Float64(tan(k_m) * sin(k_m)) * (Float64(Float64(k_m / l) * sqrt(t_m)) ^ 2.0))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 6.8e-48) tmp = 2.0 / (((k_m / l) * (k_m * sqrt(t_m))) ^ 2.0); else tmp = 2.0 / ((tan(k_m) * sin(k_m)) * (((k_m / l) * sqrt(t_m)) ^ 2.0)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 6.8e-48], N[(2.0 / N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] * N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 6.8 \cdot 10^{-48}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{\ell} \cdot \left(k\_m \cdot \sqrt{t\_m}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\tan k\_m \cdot \sin k\_m\right) \cdot {\left(\frac{k\_m}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}\\
\end{array}
\end{array}
if k < 6.80000000000000056e-48Initial program 38.4%
Simplified44.6%
Taylor expanded in t around 0 73.2%
times-frac72.4%
Simplified72.4%
associate-*l/73.2%
associate-/l*73.2%
unpow273.2%
*-un-lft-identity73.2%
times-frac73.2%
tan-quot73.2%
Applied egg-rr73.2%
*-commutative73.2%
associate-/l*72.4%
unpow272.4%
associate-/r*84.1%
unpow284.1%
associate-*r/87.3%
associate-*l/89.3%
unpow289.3%
*-commutative89.3%
/-rgt-identity89.3%
Simplified89.3%
Taylor expanded in k around 0 76.8%
add-sqr-sqrt45.6%
pow245.6%
sqrt-prod45.6%
unpow245.6%
sqrt-prod26.2%
add-sqr-sqrt45.6%
*-commutative45.6%
sqrt-prod45.7%
unpow245.7%
sqrt-prod18.0%
add-sqr-sqrt47.7%
Applied egg-rr47.7%
if 6.80000000000000056e-48 < k Initial program 30.8%
Simplified47.3%
Taylor expanded in t around 0 72.3%
times-frac72.4%
Simplified72.4%
div-inv72.4%
add-sqr-sqrt72.3%
pow272.3%
sqrt-div72.4%
unpow272.4%
sqrt-prod76.4%
add-sqr-sqrt76.6%
unpow276.6%
sqrt-prod51.9%
add-sqr-sqrt93.6%
associate-/l*93.6%
unpow293.6%
*-un-lft-identity93.6%
times-frac93.6%
tan-quot93.7%
Applied egg-rr93.7%
associate-*r/93.7%
metadata-eval93.7%
associate-/r*93.6%
/-rgt-identity93.6%
associate-/r*93.6%
Simplified93.6%
add-sqr-sqrt73.2%
sqrt-div44.4%
sqrt-div44.4%
unpow244.4%
sqrt-prod28.7%
add-sqr-sqrt40.4%
sqrt-div40.4%
sqrt-div40.4%
unpow240.4%
sqrt-prod31.3%
add-sqr-sqrt49.1%
Applied egg-rr49.1%
unpow249.1%
associate-/l/48.9%
Simplified48.9%
div-inv48.9%
unpow248.9%
frac-times48.1%
rem-square-sqrt48.2%
pow248.2%
Applied egg-rr48.2%
associate-*r/48.2%
*-rgt-identity48.2%
associate-/l/48.2%
*-commutative48.2%
Simplified48.2%
Final simplification47.8%
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 (* (/ 1.0 (tan k_m)) (/ (/ 2.0 (pow (* (/ k_m l) (sqrt t_m)) 2.0)) (sin k_m)))))
k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((1.0 / tan(k_m)) * ((2.0 / pow(((k_m / l) * sqrt(t_m)), 2.0)) / sin(k_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 * ((1.0d0 / tan(k_m)) * ((2.0d0 / (((k_m / l) * sqrt(t_m)) ** 2.0d0)) / sin(k_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 * ((1.0 / Math.tan(k_m)) * ((2.0 / Math.pow(((k_m / l) * Math.sqrt(t_m)), 2.0)) / Math.sin(k_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 * ((1.0 / math.tan(k_m)) * ((2.0 / math.pow(((k_m / l) * math.sqrt(t_m)), 2.0)) / math.sin(k_m)))
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64(1.0 / tan(k_m)) * Float64(Float64(2.0 / (Float64(Float64(k_m / l) * sqrt(t_m)) ^ 2.0)) / sin(k_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 * ((1.0 / tan(k_m)) * ((2.0 / (((k_m / l) * sqrt(t_m)) ^ 2.0)) / sin(k_m))); end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(N[(1.0 / N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 / N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[Sin[k$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(\frac{1}{\tan k\_m} \cdot \frac{\frac{2}{{\left(\frac{k\_m}{\ell} \cdot \sqrt{t\_m}\right)}^{2}}}{\sin k\_m}\right)
\end{array}
Initial program 36.2%
Simplified45.3%
Taylor expanded in t around 0 73.0%
times-frac72.4%
Simplified72.4%
div-inv72.4%
add-sqr-sqrt72.3%
pow272.3%
sqrt-div72.4%
unpow272.4%
sqrt-prod41.7%
add-sqr-sqrt75.2%
unpow275.2%
sqrt-prod50.8%
add-sqr-sqrt90.6%
associate-/l*90.6%
unpow290.6%
*-un-lft-identity90.6%
times-frac90.5%
tan-quot90.6%
Applied egg-rr90.6%
associate-*r/90.6%
metadata-eval90.6%
associate-/r*90.7%
/-rgt-identity90.7%
associate-/r*90.8%
Simplified90.8%
add-sqr-sqrt63.4%
sqrt-div51.2%
sqrt-div51.2%
unpow251.2%
sqrt-prod28.7%
add-sqr-sqrt31.4%
sqrt-div31.4%
sqrt-div31.4%
unpow231.4%
sqrt-prod26.6%
add-sqr-sqrt54.3%
Applied egg-rr54.3%
unpow254.3%
associate-/l/54.3%
Simplified54.3%
*-un-lft-identity54.3%
*-commutative54.3%
times-frac55.6%
unpow255.6%
frac-times55.1%
rem-square-sqrt55.2%
pow255.2%
Applied egg-rr55.2%
Final simplification55.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 2.5e-43)
(/ 2.0 (pow (* (/ k_m l) (* k_m (sqrt t_m))) 2.0))
(if (<= k_m 1.5e+83)
(* l (/ 2.0 (* (/ t_m l) (/ (pow (* k_m (sin k_m)) 2.0) (cos k_m)))))
(/ (/ 2.0 (pow (/ k_m l) 2.0)) (* (tan k_m) (* t_m (sin k_m))))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.5e-43) {
tmp = 2.0 / pow(((k_m / l) * (k_m * sqrt(t_m))), 2.0);
} else if (k_m <= 1.5e+83) {
tmp = l * (2.0 / ((t_m / l) * (pow((k_m * sin(k_m)), 2.0) / cos(k_m))));
} else {
tmp = (2.0 / pow((k_m / l), 2.0)) / (tan(k_m) * (t_m * sin(k_m)));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 2.5d-43) then
tmp = 2.0d0 / (((k_m / l) * (k_m * sqrt(t_m))) ** 2.0d0)
else if (k_m <= 1.5d+83) then
tmp = l * (2.0d0 / ((t_m / l) * (((k_m * sin(k_m)) ** 2.0d0) / cos(k_m))))
else
tmp = (2.0d0 / ((k_m / l) ** 2.0d0)) / (tan(k_m) * (t_m * sin(k_m)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.5e-43) {
tmp = 2.0 / Math.pow(((k_m / l) * (k_m * Math.sqrt(t_m))), 2.0);
} else if (k_m <= 1.5e+83) {
tmp = l * (2.0 / ((t_m / l) * (Math.pow((k_m * Math.sin(k_m)), 2.0) / Math.cos(k_m))));
} else {
tmp = (2.0 / Math.pow((k_m / l), 2.0)) / (Math.tan(k_m) * (t_m * Math.sin(k_m)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 2.5e-43: tmp = 2.0 / math.pow(((k_m / l) * (k_m * math.sqrt(t_m))), 2.0) elif k_m <= 1.5e+83: tmp = l * (2.0 / ((t_m / l) * (math.pow((k_m * math.sin(k_m)), 2.0) / math.cos(k_m)))) else: tmp = (2.0 / math.pow((k_m / l), 2.0)) / (math.tan(k_m) * (t_m * math.sin(k_m))) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 2.5e-43) tmp = Float64(2.0 / (Float64(Float64(k_m / l) * Float64(k_m * sqrt(t_m))) ^ 2.0)); elseif (k_m <= 1.5e+83) tmp = Float64(l * Float64(2.0 / Float64(Float64(t_m / l) * Float64((Float64(k_m * sin(k_m)) ^ 2.0) / cos(k_m))))); else tmp = Float64(Float64(2.0 / (Float64(k_m / l) ^ 2.0)) / Float64(tan(k_m) * Float64(t_m * sin(k_m)))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 2.5e-43) tmp = 2.0 / (((k_m / l) * (k_m * sqrt(t_m))) ^ 2.0); elseif (k_m <= 1.5e+83) tmp = l * (2.0 / ((t_m / l) * (((k_m * sin(k_m)) ^ 2.0) / cos(k_m)))); else tmp = (2.0 / ((k_m / l) ^ 2.0)) / (tan(k_m) * (t_m * sin(k_m))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 2.5e-43], N[(2.0 / N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.5e+83], N[(l * N[(2.0 / N[(N[(t$95$m / l), $MachinePrecision] * N[(N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[Power[N[(k$95$m / l), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[(N[Tan[k$95$m], $MachinePrecision] * N[(t$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 2.5 \cdot 10^{-43}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{\ell} \cdot \left(k\_m \cdot \sqrt{t\_m}\right)\right)}^{2}}\\
\mathbf{elif}\;k\_m \leq 1.5 \cdot 10^{+83}:\\
\;\;\;\;\ell \cdot \frac{2}{\frac{t\_m}{\ell} \cdot \frac{{\left(k\_m \cdot \sin k\_m\right)}^{2}}{\cos k\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{{\left(\frac{k\_m}{\ell}\right)}^{2}}}{\tan k\_m \cdot \left(t\_m \cdot \sin k\_m\right)}\\
\end{array}
\end{array}
if k < 2.50000000000000009e-43Initial program 38.4%
Simplified44.6%
Taylor expanded in t around 0 73.2%
times-frac72.4%
Simplified72.4%
associate-*l/73.2%
associate-/l*73.2%
unpow273.2%
*-un-lft-identity73.2%
times-frac73.2%
tan-quot73.2%
Applied egg-rr73.2%
*-commutative73.2%
associate-/l*72.4%
unpow272.4%
associate-/r*84.1%
unpow284.1%
associate-*r/87.3%
associate-*l/89.3%
unpow289.3%
*-commutative89.3%
/-rgt-identity89.3%
Simplified89.3%
Taylor expanded in k around 0 76.8%
add-sqr-sqrt45.6%
pow245.6%
sqrt-prod45.6%
unpow245.6%
sqrt-prod26.2%
add-sqr-sqrt45.6%
*-commutative45.6%
sqrt-prod45.7%
unpow245.7%
sqrt-prod18.0%
add-sqr-sqrt47.7%
Applied egg-rr47.7%
if 2.50000000000000009e-43 < k < 1.5e83Initial program 30.2%
Simplified51.0%
associate-/r*54.8%
associate-*l/54.8%
+-rgt-identity54.8%
associate-*l*54.8%
Applied egg-rr54.8%
Taylor expanded in t around 0 91.8%
associate-*r*91.7%
*-commutative91.7%
*-commutative91.7%
Simplified91.7%
associate-/r/91.8%
associate-*l*91.9%
pow-prod-down91.7%
*-commutative91.7%
Applied egg-rr91.7%
*-commutative91.7%
times-frac99.6%
Simplified99.6%
if 1.5e83 < k Initial program 31.1%
Simplified45.4%
Taylor expanded in t around 0 64.8%
times-frac66.7%
Simplified66.7%
div-inv66.7%
add-sqr-sqrt66.7%
pow266.7%
sqrt-div66.7%
unpow266.7%
sqrt-prod72.8%
add-sqr-sqrt72.9%
unpow272.9%
sqrt-prod56.8%
add-sqr-sqrt96.4%
associate-/l*96.4%
unpow296.4%
*-un-lft-identity96.4%
times-frac96.3%
tan-quot96.4%
Applied egg-rr96.4%
associate-*r/96.4%
metadata-eval96.4%
associate-/r*96.5%
/-rgt-identity96.5%
associate-/r*96.5%
Simplified96.5%
associate-/l/96.5%
*-commutative96.5%
associate-/r*96.4%
div-inv96.4%
Applied egg-rr96.4%
associate-*r/96.4%
metadata-eval96.4%
associate-/r*96.5%
associate-*r*96.6%
Simplified96.6%
Final simplification61.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 2.4e-35)
(/ 2.0 (pow (* (/ k_m l) (* k_m (sqrt t_m))) 2.0))
(/ (/ 2.0 (pow (/ k_m l) 2.0)) (* (tan k_m) (* t_m (sin k_m)))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.4e-35) {
tmp = 2.0 / pow(((k_m / l) * (k_m * sqrt(t_m))), 2.0);
} else {
tmp = (2.0 / pow((k_m / l), 2.0)) / (tan(k_m) * (t_m * sin(k_m)));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 2.4d-35) then
tmp = 2.0d0 / (((k_m / l) * (k_m * sqrt(t_m))) ** 2.0d0)
else
tmp = (2.0d0 / ((k_m / l) ** 2.0d0)) / (tan(k_m) * (t_m * sin(k_m)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.4e-35) {
tmp = 2.0 / Math.pow(((k_m / l) * (k_m * Math.sqrt(t_m))), 2.0);
} else {
tmp = (2.0 / Math.pow((k_m / l), 2.0)) / (Math.tan(k_m) * (t_m * Math.sin(k_m)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 2.4e-35: tmp = 2.0 / math.pow(((k_m / l) * (k_m * math.sqrt(t_m))), 2.0) else: tmp = (2.0 / math.pow((k_m / l), 2.0)) / (math.tan(k_m) * (t_m * math.sin(k_m))) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 2.4e-35) tmp = Float64(2.0 / (Float64(Float64(k_m / l) * Float64(k_m * sqrt(t_m))) ^ 2.0)); else tmp = Float64(Float64(2.0 / (Float64(k_m / l) ^ 2.0)) / Float64(tan(k_m) * Float64(t_m * sin(k_m)))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 2.4e-35) tmp = 2.0 / (((k_m / l) * (k_m * sqrt(t_m))) ^ 2.0); else tmp = (2.0 / ((k_m / l) ^ 2.0)) / (tan(k_m) * (t_m * sin(k_m))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 2.4e-35], N[(2.0 / N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[Power[N[(k$95$m / l), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[(N[Tan[k$95$m], $MachinePrecision] * N[(t$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 2.4 \cdot 10^{-35}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{\ell} \cdot \left(k\_m \cdot \sqrt{t\_m}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{{\left(\frac{k\_m}{\ell}\right)}^{2}}}{\tan k\_m \cdot \left(t\_m \cdot \sin k\_m\right)}\\
\end{array}
\end{array}
if k < 2.4000000000000001e-35Initial program 38.2%
Simplified44.4%
Taylor expanded in t around 0 72.9%
times-frac72.5%
Simplified72.5%
associate-*l/72.8%
associate-/l*72.9%
unpow272.9%
*-un-lft-identity72.9%
times-frac72.9%
tan-quot72.8%
Applied egg-rr72.8%
*-commutative72.8%
associate-/l*72.5%
unpow272.5%
associate-/r*84.2%
unpow284.2%
associate-*r/87.4%
associate-*l/89.4%
unpow289.4%
*-commutative89.4%
/-rgt-identity89.4%
Simplified89.4%
Taylor expanded in k around 0 76.9%
add-sqr-sqrt45.4%
pow245.4%
sqrt-prod45.4%
unpow245.4%
sqrt-prod26.1%
add-sqr-sqrt45.4%
*-commutative45.4%
sqrt-prod45.5%
unpow245.5%
sqrt-prod17.9%
add-sqr-sqrt47.4%
Applied egg-rr47.4%
if 2.4000000000000001e-35 < k Initial program 31.1%
Simplified47.8%
Taylor expanded in t around 0 73.2%
times-frac72.0%
Simplified72.0%
div-inv72.0%
add-sqr-sqrt72.0%
pow272.0%
sqrt-div72.0%
unpow272.0%
sqrt-prod76.2%
add-sqr-sqrt76.3%
unpow276.3%
sqrt-prod51.2%
add-sqr-sqrt93.5%
associate-/l*93.6%
unpow293.6%
*-un-lft-identity93.6%
times-frac93.5%
tan-quot93.6%
Applied egg-rr93.6%
associate-*r/93.6%
metadata-eval93.6%
associate-/r*93.6%
/-rgt-identity93.6%
associate-/r*93.6%
Simplified93.6%
associate-/l/93.6%
*-commutative93.6%
associate-/r*93.6%
div-inv93.6%
Applied egg-rr93.6%
associate-*r/93.6%
metadata-eval93.6%
associate-/r*93.6%
associate-*r*93.7%
Simplified93.7%
Final simplification60.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 2.4e-35)
(/ 2.0 (pow (* (/ k_m l) (* k_m (sqrt t_m))) 2.0))
(/ (/ (* 2.0 (pow (/ k_m l) -2.0)) t_m) (* (tan k_m) (sin k_m))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.4e-35) {
tmp = 2.0 / pow(((k_m / l) * (k_m * sqrt(t_m))), 2.0);
} else {
tmp = ((2.0 * pow((k_m / l), -2.0)) / t_m) / (tan(k_m) * sin(k_m));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 2.4d-35) then
tmp = 2.0d0 / (((k_m / l) * (k_m * sqrt(t_m))) ** 2.0d0)
else
tmp = ((2.0d0 * ((k_m / l) ** (-2.0d0))) / t_m) / (tan(k_m) * sin(k_m))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.4e-35) {
tmp = 2.0 / Math.pow(((k_m / l) * (k_m * Math.sqrt(t_m))), 2.0);
} else {
tmp = ((2.0 * Math.pow((k_m / l), -2.0)) / t_m) / (Math.tan(k_m) * Math.sin(k_m));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 2.4e-35: tmp = 2.0 / math.pow(((k_m / l) * (k_m * math.sqrt(t_m))), 2.0) else: tmp = ((2.0 * math.pow((k_m / l), -2.0)) / t_m) / (math.tan(k_m) * math.sin(k_m)) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 2.4e-35) tmp = Float64(2.0 / (Float64(Float64(k_m / l) * Float64(k_m * sqrt(t_m))) ^ 2.0)); else tmp = Float64(Float64(Float64(2.0 * (Float64(k_m / l) ^ -2.0)) / t_m) / Float64(tan(k_m) * sin(k_m))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 2.4e-35) tmp = 2.0 / (((k_m / l) * (k_m * sqrt(t_m))) ^ 2.0); else tmp = ((2.0 * ((k_m / l) ^ -2.0)) / t_m) / (tan(k_m) * sin(k_m)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 2.4e-35], N[(2.0 / N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 * N[Power[N[(k$95$m / l), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$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 2.4 \cdot 10^{-35}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{\ell} \cdot \left(k\_m \cdot \sqrt{t\_m}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2 \cdot {\left(\frac{k\_m}{\ell}\right)}^{-2}}{t\_m}}{\tan k\_m \cdot \sin k\_m}\\
\end{array}
\end{array}
if k < 2.4000000000000001e-35Initial program 38.2%
Simplified44.4%
Taylor expanded in t around 0 72.9%
times-frac72.5%
Simplified72.5%
associate-*l/72.8%
associate-/l*72.9%
unpow272.9%
*-un-lft-identity72.9%
times-frac72.9%
tan-quot72.8%
Applied egg-rr72.8%
*-commutative72.8%
associate-/l*72.5%
unpow272.5%
associate-/r*84.2%
unpow284.2%
associate-*r/87.4%
associate-*l/89.4%
unpow289.4%
*-commutative89.4%
/-rgt-identity89.4%
Simplified89.4%
Taylor expanded in k around 0 76.9%
add-sqr-sqrt45.4%
pow245.4%
sqrt-prod45.4%
unpow245.4%
sqrt-prod26.1%
add-sqr-sqrt45.4%
*-commutative45.4%
sqrt-prod45.5%
unpow245.5%
sqrt-prod17.9%
add-sqr-sqrt47.4%
Applied egg-rr47.4%
if 2.4000000000000001e-35 < k Initial program 31.1%
Simplified47.8%
Taylor expanded in t around 0 73.2%
times-frac72.0%
Simplified72.0%
div-inv72.0%
add-sqr-sqrt72.0%
pow272.0%
sqrt-div72.0%
unpow272.0%
sqrt-prod76.2%
add-sqr-sqrt76.3%
unpow276.3%
sqrt-prod51.2%
add-sqr-sqrt93.5%
associate-/l*93.6%
unpow293.6%
*-un-lft-identity93.6%
times-frac93.5%
tan-quot93.6%
Applied egg-rr93.6%
associate-*r/93.6%
metadata-eval93.6%
associate-/r*93.6%
/-rgt-identity93.6%
associate-/r*93.6%
Simplified93.6%
div-inv93.6%
div-inv93.6%
pow-flip93.5%
metadata-eval93.5%
Applied egg-rr93.5%
associate-*r/93.6%
*-rgt-identity93.6%
Simplified93.6%
Final simplification60.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 2.8e-43)
(/ 2.0 (pow (* (/ k_m l) (* k_m (sqrt t_m))) 2.0))
(/
(/ (/ 2.0 (/ 1.0 (* (/ l k_m) (/ l k_m)))) t_m)
(* (tan k_m) (sin k_m))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.8e-43) {
tmp = 2.0 / pow(((k_m / l) * (k_m * sqrt(t_m))), 2.0);
} else {
tmp = ((2.0 / (1.0 / ((l / k_m) * (l / k_m)))) / t_m) / (tan(k_m) * sin(k_m));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 2.8d-43) then
tmp = 2.0d0 / (((k_m / l) * (k_m * sqrt(t_m))) ** 2.0d0)
else
tmp = ((2.0d0 / (1.0d0 / ((l / k_m) * (l / k_m)))) / t_m) / (tan(k_m) * sin(k_m))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.8e-43) {
tmp = 2.0 / Math.pow(((k_m / l) * (k_m * Math.sqrt(t_m))), 2.0);
} else {
tmp = ((2.0 / (1.0 / ((l / k_m) * (l / k_m)))) / t_m) / (Math.tan(k_m) * Math.sin(k_m));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 2.8e-43: tmp = 2.0 / math.pow(((k_m / l) * (k_m * math.sqrt(t_m))), 2.0) else: tmp = ((2.0 / (1.0 / ((l / k_m) * (l / k_m)))) / t_m) / (math.tan(k_m) * math.sin(k_m)) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 2.8e-43) tmp = Float64(2.0 / (Float64(Float64(k_m / l) * Float64(k_m * sqrt(t_m))) ^ 2.0)); else tmp = Float64(Float64(Float64(2.0 / Float64(1.0 / Float64(Float64(l / k_m) * Float64(l / k_m)))) / t_m) / Float64(tan(k_m) * sin(k_m))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 2.8e-43) tmp = 2.0 / (((k_m / l) * (k_m * sqrt(t_m))) ^ 2.0); else tmp = ((2.0 / (1.0 / ((l / k_m) * (l / k_m)))) / t_m) / (tan(k_m) * sin(k_m)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 2.8e-43], N[(2.0 / N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / N[(1.0 / N[(N[(l / k$95$m), $MachinePrecision] * N[(l / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$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 2.8 \cdot 10^{-43}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{\ell} \cdot \left(k\_m \cdot \sqrt{t\_m}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{2}{\frac{1}{\frac{\ell}{k\_m} \cdot \frac{\ell}{k\_m}}}}{t\_m}}{\tan k\_m \cdot \sin k\_m}\\
\end{array}
\end{array}
if k < 2.7999999999999998e-43Initial program 38.4%
Simplified44.6%
Taylor expanded in t around 0 73.2%
times-frac72.4%
Simplified72.4%
associate-*l/73.2%
associate-/l*73.2%
unpow273.2%
*-un-lft-identity73.2%
times-frac73.2%
tan-quot73.2%
Applied egg-rr73.2%
*-commutative73.2%
associate-/l*72.4%
unpow272.4%
associate-/r*84.1%
unpow284.1%
associate-*r/87.3%
associate-*l/89.3%
unpow289.3%
*-commutative89.3%
/-rgt-identity89.3%
Simplified89.3%
Taylor expanded in k around 0 76.8%
add-sqr-sqrt45.6%
pow245.6%
sqrt-prod45.6%
unpow245.6%
sqrt-prod26.2%
add-sqr-sqrt45.6%
*-commutative45.6%
sqrt-prod45.7%
unpow245.7%
sqrt-prod18.0%
add-sqr-sqrt47.7%
Applied egg-rr47.7%
if 2.7999999999999998e-43 < k Initial program 30.8%
Simplified47.3%
Taylor expanded in t around 0 72.3%
times-frac72.4%
Simplified72.4%
div-inv72.4%
add-sqr-sqrt72.3%
pow272.3%
sqrt-div72.4%
unpow272.4%
sqrt-prod76.4%
add-sqr-sqrt76.6%
unpow276.6%
sqrt-prod51.9%
add-sqr-sqrt93.6%
associate-/l*93.6%
unpow293.6%
*-un-lft-identity93.6%
times-frac93.6%
tan-quot93.7%
Applied egg-rr93.7%
associate-*r/93.7%
metadata-eval93.7%
associate-/r*93.6%
/-rgt-identity93.6%
associate-/r*93.6%
Simplified93.6%
unpow293.6%
clear-num93.7%
clear-num93.6%
frac-times93.7%
metadata-eval93.7%
Applied egg-rr93.7%
Final simplification60.8%
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 2.4e-35)
(/ 2.0 (pow (* (/ k_m l) (* k_m (sqrt t_m))) 2.0))
(/ 2.0 (* (* (/ k_m l) (/ k_m l)) (* t_m (* (tan k_m) (sin k_m))))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.4e-35) {
tmp = 2.0 / pow(((k_m / l) * (k_m * sqrt(t_m))), 2.0);
} else {
tmp = 2.0 / (((k_m / l) * (k_m / l)) * (t_m * (tan(k_m) * sin(k_m))));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 2.4d-35) then
tmp = 2.0d0 / (((k_m / l) * (k_m * sqrt(t_m))) ** 2.0d0)
else
tmp = 2.0d0 / (((k_m / l) * (k_m / l)) * (t_m * (tan(k_m) * sin(k_m))))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2.4e-35) {
tmp = 2.0 / Math.pow(((k_m / l) * (k_m * Math.sqrt(t_m))), 2.0);
} else {
tmp = 2.0 / (((k_m / l) * (k_m / l)) * (t_m * (Math.tan(k_m) * Math.sin(k_m))));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 2.4e-35: tmp = 2.0 / math.pow(((k_m / l) * (k_m * math.sqrt(t_m))), 2.0) else: tmp = 2.0 / (((k_m / l) * (k_m / l)) * (t_m * (math.tan(k_m) * math.sin(k_m)))) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 2.4e-35) tmp = Float64(2.0 / (Float64(Float64(k_m / l) * Float64(k_m * sqrt(t_m))) ^ 2.0)); else tmp = Float64(2.0 / Float64(Float64(Float64(k_m / l) * Float64(k_m / l)) * Float64(t_m * Float64(tan(k_m) * sin(k_m))))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 2.4e-35) tmp = 2.0 / (((k_m / l) * (k_m * sqrt(t_m))) ^ 2.0); else tmp = 2.0 / (((k_m / l) * (k_m / l)) * (t_m * (tan(k_m) * sin(k_m)))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 2.4e-35], N[(2.0 / N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$95$m], $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 2.4 \cdot 10^{-35}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{\ell} \cdot \left(k\_m \cdot \sqrt{t\_m}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\frac{k\_m}{\ell} \cdot \frac{k\_m}{\ell}\right) \cdot \left(t\_m \cdot \left(\tan k\_m \cdot \sin k\_m\right)\right)}\\
\end{array}
\end{array}
if k < 2.4000000000000001e-35Initial program 38.2%
Simplified44.4%
Taylor expanded in t around 0 72.9%
times-frac72.5%
Simplified72.5%
associate-*l/72.8%
associate-/l*72.9%
unpow272.9%
*-un-lft-identity72.9%
times-frac72.9%
tan-quot72.8%
Applied egg-rr72.8%
*-commutative72.8%
associate-/l*72.5%
unpow272.5%
associate-/r*84.2%
unpow284.2%
associate-*r/87.4%
associate-*l/89.4%
unpow289.4%
*-commutative89.4%
/-rgt-identity89.4%
Simplified89.4%
Taylor expanded in k around 0 76.9%
add-sqr-sqrt45.4%
pow245.4%
sqrt-prod45.4%
unpow245.4%
sqrt-prod26.1%
add-sqr-sqrt45.4%
*-commutative45.4%
sqrt-prod45.5%
unpow245.5%
sqrt-prod17.9%
add-sqr-sqrt47.4%
Applied egg-rr47.4%
if 2.4000000000000001e-35 < k Initial program 31.1%
Simplified47.8%
Taylor expanded in t around 0 73.2%
times-frac72.0%
Simplified72.0%
associate-*l/73.2%
associate-/l*73.2%
unpow273.2%
*-un-lft-identity73.2%
times-frac73.1%
tan-quot73.2%
Applied egg-rr73.2%
*-commutative73.2%
associate-/l*72.1%
unpow272.1%
associate-/r*79.1%
unpow279.1%
associate-*r/91.1%
associate-*l/93.6%
unpow293.6%
*-commutative93.6%
/-rgt-identity93.6%
Simplified93.6%
unpow293.6%
Applied egg-rr93.6%
Final simplification60.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.8e-35)
(/ 2.0 (pow (* (/ k_m l) (* k_m (sqrt t_m))) 2.0))
(/ (/ (/ 2.0 (* (/ k_m l) (/ k_m l))) t_m) (* (tan k_m) (sin k_m))))))k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1.8e-35) {
tmp = 2.0 / pow(((k_m / l) * (k_m * sqrt(t_m))), 2.0);
} else {
tmp = ((2.0 / ((k_m / l) * (k_m / l))) / t_m) / (tan(k_m) * sin(k_m));
}
return t_s * tmp;
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 1.8d-35) then
tmp = 2.0d0 / (((k_m / l) * (k_m * sqrt(t_m))) ** 2.0d0)
else
tmp = ((2.0d0 / ((k_m / l) * (k_m / l))) / t_m) / (tan(k_m) * sin(k_m))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1.8e-35) {
tmp = 2.0 / Math.pow(((k_m / l) * (k_m * Math.sqrt(t_m))), 2.0);
} else {
tmp = ((2.0 / ((k_m / l) * (k_m / l))) / t_m) / (Math.tan(k_m) * Math.sin(k_m));
}
return t_s * tmp;
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 1.8e-35: tmp = 2.0 / math.pow(((k_m / l) * (k_m * math.sqrt(t_m))), 2.0) else: tmp = ((2.0 / ((k_m / l) * (k_m / l))) / t_m) / (math.tan(k_m) * math.sin(k_m)) return t_s * tmp
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 1.8e-35) tmp = Float64(2.0 / (Float64(Float64(k_m / l) * Float64(k_m * sqrt(t_m))) ^ 2.0)); else tmp = Float64(Float64(Float64(2.0 / Float64(Float64(k_m / l) * Float64(k_m / l))) / t_m) / Float64(tan(k_m) * sin(k_m))); end return Float64(t_s * tmp) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 1.8e-35) tmp = 2.0 / (((k_m / l) * (k_m * sqrt(t_m))) ^ 2.0); else tmp = ((2.0 / ((k_m / l) * (k_m / l))) / t_m) / (tan(k_m) * sin(k_m)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 1.8e-35], N[(2.0 / N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / N[(N[Tan[k$95$m], $MachinePrecision] * N[Sin[k$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.8 \cdot 10^{-35}:\\
\;\;\;\;\frac{2}{{\left(\frac{k\_m}{\ell} \cdot \left(k\_m \cdot \sqrt{t\_m}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{2}{\frac{k\_m}{\ell} \cdot \frac{k\_m}{\ell}}}{t\_m}}{\tan k\_m \cdot \sin k\_m}\\
\end{array}
\end{array}
if k < 1.80000000000000009e-35Initial program 38.2%
Simplified44.4%
Taylor expanded in t around 0 72.9%
times-frac72.5%
Simplified72.5%
associate-*l/72.8%
associate-/l*72.9%
unpow272.9%
*-un-lft-identity72.9%
times-frac72.9%
tan-quot72.8%
Applied egg-rr72.8%
*-commutative72.8%
associate-/l*72.5%
unpow272.5%
associate-/r*84.2%
unpow284.2%
associate-*r/87.4%
associate-*l/89.4%
unpow289.4%
*-commutative89.4%
/-rgt-identity89.4%
Simplified89.4%
Taylor expanded in k around 0 76.9%
add-sqr-sqrt45.4%
pow245.4%
sqrt-prod45.4%
unpow245.4%
sqrt-prod26.1%
add-sqr-sqrt45.4%
*-commutative45.4%
sqrt-prod45.5%
unpow245.5%
sqrt-prod17.9%
add-sqr-sqrt47.4%
Applied egg-rr47.4%
if 1.80000000000000009e-35 < k Initial program 31.1%
Simplified47.8%
Taylor expanded in t around 0 73.2%
times-frac72.0%
Simplified72.0%
div-inv72.0%
add-sqr-sqrt72.0%
pow272.0%
sqrt-div72.0%
unpow272.0%
sqrt-prod76.2%
add-sqr-sqrt76.3%
unpow276.3%
sqrt-prod51.2%
add-sqr-sqrt93.5%
associate-/l*93.6%
unpow293.6%
*-un-lft-identity93.6%
times-frac93.5%
tan-quot93.6%
Applied egg-rr93.6%
associate-*r/93.6%
metadata-eval93.6%
associate-/r*93.6%
/-rgt-identity93.6%
associate-/r*93.6%
Simplified93.6%
unpow293.6%
Applied egg-rr93.6%
Final simplification60.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 (* k_m (* k_m (/ (sqrt t_m) l))) 2.0))))
k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 / pow((k_m * (k_m * (sqrt(t_m) / l))), 2.0));
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 / ((k_m * (k_m * (sqrt(t_m) / l))) ** 2.0d0))
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 / Math.pow((k_m * (k_m * (Math.sqrt(t_m) / l))), 2.0));
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 / math.pow((k_m * (k_m * (math.sqrt(t_m) / l))), 2.0))
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 / (Float64(k_m * Float64(k_m * Float64(sqrt(t_m) / l))) ^ 2.0))) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 / ((k_m * (k_m * (sqrt(t_m) / l))) ^ 2.0)); end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(2.0 / N[Power[N[(k$95$m * N[(k$95$m * N[(N[Sqrt[t$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \frac{2}{{\left(k\_m \cdot \left(k\_m \cdot \frac{\sqrt{t\_m}}{\ell}\right)\right)}^{2}}
\end{array}
Initial program 36.2%
Simplified45.3%
Taylor expanded in t around 0 73.0%
times-frac72.4%
Simplified72.4%
associate-*l/72.9%
associate-/l*72.9%
unpow272.9%
*-un-lft-identity72.9%
times-frac72.9%
tan-quot72.9%
Applied egg-rr72.9%
*-commutative72.9%
associate-/l*72.4%
unpow272.4%
associate-/r*82.8%
unpow282.8%
associate-*r/88.4%
associate-*l/90.6%
unpow290.6%
*-commutative90.6%
/-rgt-identity90.6%
Simplified90.6%
Taylor expanded in k around 0 72.7%
add-sqr-sqrt41.2%
pow241.2%
sqrt-prod41.2%
unpow241.2%
sqrt-prod23.6%
add-sqr-sqrt41.2%
*-commutative41.2%
sqrt-prod41.3%
unpow241.3%
sqrt-prod21.4%
add-sqr-sqrt42.6%
Applied egg-rr42.6%
associate-*l/41.6%
associate-*r/42.6%
associate-*r/42.6%
*-commutative42.6%
associate-*l/42.6%
associate-/l*43.4%
Simplified43.4%
Final simplification43.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 (* (/ k_m l) (* k_m (sqrt t_m))) 2.0))))
k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 / pow(((k_m / l) * (k_m * sqrt(t_m))), 2.0));
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 / (((k_m / l) * (k_m * sqrt(t_m))) ** 2.0d0))
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 / Math.pow(((k_m / l) * (k_m * Math.sqrt(t_m))), 2.0));
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 / math.pow(((k_m / l) * (k_m * math.sqrt(t_m))), 2.0))
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 / (Float64(Float64(k_m / l) * Float64(k_m * sqrt(t_m))) ^ 2.0))) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 / (((k_m / l) * (k_m * sqrt(t_m))) ^ 2.0)); end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(2.0 / N[Power[N[(N[(k$95$m / l), $MachinePrecision] * N[(k$95$m * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \frac{2}{{\left(\frac{k\_m}{\ell} \cdot \left(k\_m \cdot \sqrt{t\_m}\right)\right)}^{2}}
\end{array}
Initial program 36.2%
Simplified45.3%
Taylor expanded in t around 0 73.0%
times-frac72.4%
Simplified72.4%
associate-*l/72.9%
associate-/l*72.9%
unpow272.9%
*-un-lft-identity72.9%
times-frac72.9%
tan-quot72.9%
Applied egg-rr72.9%
*-commutative72.9%
associate-/l*72.4%
unpow272.4%
associate-/r*82.8%
unpow282.8%
associate-*r/88.4%
associate-*l/90.6%
unpow290.6%
*-commutative90.6%
/-rgt-identity90.6%
Simplified90.6%
Taylor expanded in k around 0 72.7%
add-sqr-sqrt41.2%
pow241.2%
sqrt-prod41.2%
unpow241.2%
sqrt-prod23.6%
add-sqr-sqrt41.2%
*-commutative41.2%
sqrt-prod41.3%
unpow241.3%
sqrt-prod21.4%
add-sqr-sqrt42.6%
Applied egg-rr42.6%
Final simplification42.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 (/ 1.0 (* t_m (pow (* k_m (/ k_m l)) 2.0))))))
k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (1.0 / (t_m * pow((k_m * (k_m / l)), 2.0))));
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 * (1.0d0 / (t_m * ((k_m * (k_m / l)) ** 2.0d0))))
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (1.0 / (t_m * Math.pow((k_m * (k_m / l)), 2.0))));
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 * (1.0 / (t_m * math.pow((k_m * (k_m / l)), 2.0))))
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 * Float64(1.0 / Float64(t_m * (Float64(k_m * Float64(k_m / l)) ^ 2.0))))) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 * (1.0 / (t_m * ((k_m * (k_m / l)) ^ 2.0)))); end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(2.0 * N[(1.0 / N[(t$95$m * N[Power[N[(k$95$m * N[(k$95$m / l), $MachinePrecision]), $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 \left(2 \cdot \frac{1}{t\_m \cdot {\left(k\_m \cdot \frac{k\_m}{\ell}\right)}^{2}}\right)
\end{array}
Initial program 36.2%
Simplified45.3%
Taylor expanded in t around 0 73.0%
times-frac72.4%
Simplified72.4%
associate-*l/72.9%
associate-/l*72.9%
unpow272.9%
*-un-lft-identity72.9%
times-frac72.9%
tan-quot72.9%
Applied egg-rr72.9%
*-commutative72.9%
associate-/l*72.4%
unpow272.4%
associate-/r*82.8%
unpow282.8%
associate-*r/88.4%
associate-*l/90.6%
unpow290.6%
*-commutative90.6%
/-rgt-identity90.6%
Simplified90.6%
Taylor expanded in k around 0 72.7%
div-inv72.7%
associate-*r*72.0%
pow-prod-down73.7%
Applied egg-rr73.7%
Final simplification73.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 (* 2.0 (* l (* (/ l t_m) (pow k_m -4.0))))))
k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (l * ((l / t_m) * pow(k_m, -4.0))));
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 * (l * ((l / t_m) * (k_m ** (-4.0d0)))))
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (l * ((l / t_m) * Math.pow(k_m, -4.0))));
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 * (l * ((l / t_m) * math.pow(k_m, -4.0))))
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 * Float64(l * Float64(Float64(l / t_m) * (k_m ^ -4.0))))) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 * (l * ((l / t_m) * (k_m ^ -4.0)))); end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(2.0 * N[(l * N[(N[(l / t$95$m), $MachinePrecision] * N[Power[k$95$m, -4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(2 \cdot \left(\ell \cdot \left(\frac{\ell}{t\_m} \cdot {k\_m}^{-4}\right)\right)\right)
\end{array}
Initial program 36.2%
Simplified45.3%
Taylor expanded in k around 0 64.5%
unpow264.5%
times-frac70.5%
Applied egg-rr70.5%
associate-*r/70.8%
div-inv70.8%
pow-flip70.8%
metadata-eval70.8%
Applied egg-rr70.8%
associate-/l*70.5%
Simplified70.5%
associate-*r/70.8%
Applied egg-rr70.8%
associate-/l*70.5%
associate-*r*70.5%
Simplified70.5%
Final simplification70.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 (* (/ l t_m) (* l (pow k_m -4.0))))))
k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * ((l / t_m) * (l * pow(k_m, -4.0))));
}
k_m = abs(k)
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 * ((l / t_m) * (l * (k_m ** (-4.0d0)))))
end function
k_m = Math.abs(k);
t_m = Math.abs(t);
t_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * ((l / t_m) * (l * Math.pow(k_m, -4.0))));
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 * ((l / t_m) * (l * math.pow(k_m, -4.0))))
k_m = abs(k) t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 * Float64(Float64(l / t_m) * Float64(l * (k_m ^ -4.0))))) end
k_m = abs(k); t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 * ((l / t_m) * (l * (k_m ^ -4.0)))); end
k_m = N[Abs[k], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
t_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(2.0 * N[(N[(l / t$95$m), $MachinePrecision] * N[(l * N[Power[k$95$m, -4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(2 \cdot \left(\frac{\ell}{t\_m} \cdot \left(\ell \cdot {k\_m}^{-4}\right)\right)\right)
\end{array}
Initial program 36.2%
Simplified45.3%
Taylor expanded in k around 0 64.5%
unpow264.5%
times-frac70.5%
Applied egg-rr70.5%
associate-*r/70.8%
div-inv70.8%
pow-flip70.8%
metadata-eval70.8%
Applied egg-rr70.8%
associate-/l*70.5%
Simplified70.5%
Final simplification70.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 (* (/ l (pow k_m 4.0)) (/ l 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 * ((l / pow(k_m, 4.0)) * (l / 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 ** 4.0d0)) * (l / 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 * ((l / Math.pow(k_m, 4.0)) * (l / 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 * ((l / math.pow(k_m, 4.0)) * (l / 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 ^ 4.0)) * Float64(l / 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 ^ 4.0)) * (l / 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[(l / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision] * N[(l / 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(\frac{\ell}{{k\_m}^{4}} \cdot \frac{\ell}{t\_m}\right)\right)
\end{array}
Initial program 36.2%
Simplified45.3%
Taylor expanded in k around 0 64.5%
unpow264.5%
times-frac70.5%
Applied egg-rr70.5%
Final simplification70.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 (/ l (* t_m (/ (pow k_m 4.0) l))))))
k_m = fabs(k);
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (l / (t_m * (pow(k_m, 4.0) / l))));
}
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 / (t_m * ((k_m ** 4.0d0) / l))))
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 * (l / (t_m * (Math.pow(k_m, 4.0) / l))));
}
k_m = math.fabs(k) t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 * (l / (t_m * (math.pow(k_m, 4.0) / l))))
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(l / Float64(t_m * Float64((k_m ^ 4.0) / l))))) 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 / (t_m * ((k_m ^ 4.0) / l)))); 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[(l / N[(t$95$m * N[(N[Power[k$95$m, 4.0], $MachinePrecision] / l), $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 \frac{\ell}{t\_m \cdot \frac{{k\_m}^{4}}{\ell}}\right)
\end{array}
Initial program 36.2%
Simplified45.3%
Taylor expanded in k around 0 64.5%
unpow264.5%
times-frac70.5%
Applied egg-rr70.5%
clear-num70.5%
frac-times70.8%
*-un-lft-identity70.8%
Applied egg-rr70.8%
Final simplification70.8%
herbie shell --seed 2024040
(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))))