
(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)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 2.35e-147)
(/ 2.0 (* (log (pow (exp (pow k_m 4.0)) (pow l -2.0))) t))
(/
(* 2.0 (* (/ (pow l 2.0) k_m) (/ (cos k_m) (* t (pow (sin k_m) 2.0)))))
k_m)))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 2.35e-147) {
tmp = 2.0 / (log(pow(exp(pow(k_m, 4.0)), pow(l, -2.0))) * t);
} else {
tmp = (2.0 * ((pow(l, 2.0) / k_m) * (cos(k_m) / (t * pow(sin(k_m), 2.0))))) / k_m;
}
return tmp;
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 2.35d-147) then
tmp = 2.0d0 / (log((exp((k_m ** 4.0d0)) ** (l ** (-2.0d0)))) * t)
else
tmp = (2.0d0 * (((l ** 2.0d0) / k_m) * (cos(k_m) / (t * (sin(k_m) ** 2.0d0))))) / k_m
end if
code = tmp
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 2.35e-147) {
tmp = 2.0 / (Math.log(Math.pow(Math.exp(Math.pow(k_m, 4.0)), Math.pow(l, -2.0))) * t);
} else {
tmp = (2.0 * ((Math.pow(l, 2.0) / k_m) * (Math.cos(k_m) / (t * Math.pow(Math.sin(k_m), 2.0))))) / k_m;
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 2.35e-147: tmp = 2.0 / (math.log(math.pow(math.exp(math.pow(k_m, 4.0)), math.pow(l, -2.0))) * t) else: tmp = (2.0 * ((math.pow(l, 2.0) / k_m) * (math.cos(k_m) / (t * math.pow(math.sin(k_m), 2.0))))) / k_m return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 2.35e-147) tmp = Float64(2.0 / Float64(log((exp((k_m ^ 4.0)) ^ (l ^ -2.0))) * t)); else tmp = Float64(Float64(2.0 * Float64(Float64((l ^ 2.0) / k_m) * Float64(cos(k_m) / Float64(t * (sin(k_m) ^ 2.0))))) / k_m); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 2.35e-147) tmp = 2.0 / (log((exp((k_m ^ 4.0)) ^ (l ^ -2.0))) * t); else tmp = (2.0 * (((l ^ 2.0) / k_m) * (cos(k_m) / (t * (sin(k_m) ^ 2.0))))) / k_m; end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 2.35e-147], N[(2.0 / N[(N[Log[N[Power[N[Exp[N[Power[k$95$m, 4.0], $MachinePrecision]], $MachinePrecision], N[Power[l, -2.0], $MachinePrecision]], $MachinePrecision]], $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / k$95$m), $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[(t * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k_m \leq 2.35 \cdot 10^{-147}:\\
\;\;\;\;\frac{2}{\log \left({\left(e^{{k_m}^{4}}\right)}^{\left({\ell}^{-2}\right)}\right) \cdot t}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \left(\frac{{\ell}^{2}}{k_m} \cdot \frac{\cos k_m}{t \cdot {\sin k_m}^{2}}\right)}{k_m}\\
\end{array}
\end{array}
if k < 2.34999999999999994e-147Initial program 43.0%
associate-*l*43.0%
associate--l+43.0%
Simplified43.0%
unpow343.1%
times-frac51.4%
pow251.4%
Applied egg-rr51.4%
Taylor expanded in k around 0 69.7%
associate-/l*66.8%
associate-/r/69.7%
Simplified69.7%
div-inv69.7%
pow-flip69.7%
metadata-eval69.7%
add-log-exp67.6%
exp-prod72.9%
Applied egg-rr72.9%
if 2.34999999999999994e-147 < k Initial program 32.2%
associate-/r*32.1%
associate-*l*32.1%
associate-*l/32.1%
associate-/l*32.1%
+-commutative32.1%
unpow232.1%
sqr-neg32.1%
distribute-frac-neg32.1%
distribute-frac-neg32.1%
unpow232.1%
associate--l+41.1%
metadata-eval41.1%
+-rgt-identity41.1%
unpow241.1%
distribute-frac-neg41.1%
distribute-frac-neg41.1%
sqr-neg41.1%
unpow241.1%
Simplified41.1%
*-un-lft-identity41.1%
unpow241.1%
times-frac49.1%
clear-num49.1%
associate-/r/49.1%
associate-/r*49.1%
pow249.1%
Applied egg-rr49.1%
associate-*l/49.1%
associate-*r/49.1%
associate-*r/49.1%
associate-/l/50.0%
*-lft-identity50.0%
associate-*l/50.0%
associate-*r/50.0%
associate-*l/50.0%
*-lft-identity50.0%
Simplified50.0%
associate-*l/50.9%
associate-/l/50.9%
Applied egg-rr50.9%
Taylor expanded in t around 0 82.9%
*-commutative82.9%
Simplified82.9%
times-frac83.9%
*-commutative83.9%
Applied egg-rr83.9%
Final simplification77.6%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 9.2e-5)
(/ 2.0 (* (/ (pow k_m 2.0) (pow l 2.0)) (/ (* t (pow k_m 2.0)) (cos k_m))))
(/
(*
2.0
(/
(* (pow l 2.0) (cos k_m))
(* k_m (* t (- 0.5 (/ (cos (* k_m 2.0)) 2.0))))))
k_m)))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 9.2e-5) {
tmp = 2.0 / ((pow(k_m, 2.0) / pow(l, 2.0)) * ((t * pow(k_m, 2.0)) / cos(k_m)));
} else {
tmp = (2.0 * ((pow(l, 2.0) * cos(k_m)) / (k_m * (t * (0.5 - (cos((k_m * 2.0)) / 2.0)))))) / k_m;
}
return tmp;
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 9.2d-5) then
tmp = 2.0d0 / (((k_m ** 2.0d0) / (l ** 2.0d0)) * ((t * (k_m ** 2.0d0)) / cos(k_m)))
else
tmp = (2.0d0 * (((l ** 2.0d0) * cos(k_m)) / (k_m * (t * (0.5d0 - (cos((k_m * 2.0d0)) / 2.0d0)))))) / k_m
end if
code = tmp
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 9.2e-5) {
tmp = 2.0 / ((Math.pow(k_m, 2.0) / Math.pow(l, 2.0)) * ((t * Math.pow(k_m, 2.0)) / Math.cos(k_m)));
} else {
tmp = (2.0 * ((Math.pow(l, 2.0) * Math.cos(k_m)) / (k_m * (t * (0.5 - (Math.cos((k_m * 2.0)) / 2.0)))))) / k_m;
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 9.2e-5: tmp = 2.0 / ((math.pow(k_m, 2.0) / math.pow(l, 2.0)) * ((t * math.pow(k_m, 2.0)) / math.cos(k_m))) else: tmp = (2.0 * ((math.pow(l, 2.0) * math.cos(k_m)) / (k_m * (t * (0.5 - (math.cos((k_m * 2.0)) / 2.0)))))) / k_m return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 9.2e-5) tmp = Float64(2.0 / Float64(Float64((k_m ^ 2.0) / (l ^ 2.0)) * Float64(Float64(t * (k_m ^ 2.0)) / cos(k_m)))); else tmp = Float64(Float64(2.0 * Float64(Float64((l ^ 2.0) * cos(k_m)) / Float64(k_m * Float64(t * Float64(0.5 - Float64(cos(Float64(k_m * 2.0)) / 2.0)))))) / k_m); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 9.2e-5) tmp = 2.0 / (((k_m ^ 2.0) / (l ^ 2.0)) * ((t * (k_m ^ 2.0)) / cos(k_m))); else tmp = (2.0 * (((l ^ 2.0) * cos(k_m)) / (k_m * (t * (0.5 - (cos((k_m * 2.0)) / 2.0)))))) / k_m; end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 9.2e-5], N[(2.0 / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(t * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(k$95$m * N[(t * N[(0.5 - N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k_m \leq 9.2 \cdot 10^{-5}:\\
\;\;\;\;\frac{2}{\frac{{k_m}^{2}}{{\ell}^{2}} \cdot \frac{t \cdot {k_m}^{2}}{\cos k_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \frac{{\ell}^{2} \cdot \cos k_m}{k_m \cdot \left(t \cdot \left(0.5 - \frac{\cos \left(k_m \cdot 2\right)}{2}\right)\right)}}{k_m}\\
\end{array}
\end{array}
if k < 9.20000000000000001e-5Initial program 41.1%
associate-*l*41.1%
associate--l+41.1%
Simplified41.1%
unpow341.1%
times-frac49.5%
pow249.5%
Applied egg-rr49.5%
Taylor expanded in t around 0 75.4%
times-frac78.1%
*-commutative78.1%
Simplified78.1%
Taylor expanded in k around 0 75.5%
if 9.20000000000000001e-5 < k Initial program 32.3%
associate-/r*32.1%
associate-*l*32.1%
associate-*l/32.1%
associate-/l*32.1%
+-commutative32.1%
unpow232.1%
sqr-neg32.1%
distribute-frac-neg32.1%
distribute-frac-neg32.1%
unpow232.1%
associate--l+37.4%
metadata-eval37.4%
+-rgt-identity37.4%
unpow237.4%
distribute-frac-neg37.4%
distribute-frac-neg37.4%
sqr-neg37.4%
unpow237.4%
Simplified37.4%
*-un-lft-identity37.4%
unpow237.4%
times-frac47.4%
clear-num47.4%
associate-/r/47.4%
associate-/r*47.4%
pow247.4%
Applied egg-rr47.4%
associate-*l/47.4%
associate-*r/47.4%
associate-*r/47.4%
associate-/l/48.6%
*-lft-identity48.6%
associate-*l/48.6%
associate-*r/48.7%
associate-*l/48.7%
*-lft-identity48.7%
Simplified48.7%
associate-*l/48.7%
associate-/l/48.7%
Applied egg-rr48.7%
Taylor expanded in t around 0 81.3%
*-commutative81.3%
Simplified81.3%
unpow281.3%
sin-mult81.3%
Applied egg-rr81.3%
div-sub81.3%
+-inverses81.3%
cos-081.3%
metadata-eval81.3%
count-281.3%
*-commutative81.3%
Simplified81.3%
Final simplification77.2%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (/ (* 2.0 (* (/ (pow l 2.0) k_m) (/ (cos k_m) (* t (pow (sin k_m) 2.0))))) k_m))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return (2.0 * ((pow(l, 2.0) / k_m) * (cos(k_m) / (t * pow(sin(k_m), 2.0))))) / k_m;
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = (2.0d0 * (((l ** 2.0d0) / k_m) * (cos(k_m) / (t * (sin(k_m) ** 2.0d0))))) / k_m
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return (2.0 * ((Math.pow(l, 2.0) / k_m) * (Math.cos(k_m) / (t * Math.pow(Math.sin(k_m), 2.0))))) / k_m;
}
k_m = math.fabs(k) def code(t, l, k_m): return (2.0 * ((math.pow(l, 2.0) / k_m) * (math.cos(k_m) / (t * math.pow(math.sin(k_m), 2.0))))) / k_m
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(2.0 * Float64(Float64((l ^ 2.0) / k_m) * Float64(cos(k_m) / Float64(t * (sin(k_m) ^ 2.0))))) / k_m) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = (2.0 * (((l ^ 2.0) / k_m) * (cos(k_m) / (t * (sin(k_m) ^ 2.0))))) / k_m; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / k$95$m), $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[(t * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2 \cdot \left(\frac{{\ell}^{2}}{k_m} \cdot \frac{\cos k_m}{t \cdot {\sin k_m}^{2}}\right)}{k_m}
\end{array}
Initial program 38.4%
associate-/r*38.4%
associate-*l*38.4%
associate-*l/38.4%
associate-/l*38.4%
+-commutative38.4%
unpow238.4%
sqr-neg38.4%
distribute-frac-neg38.4%
distribute-frac-neg38.4%
unpow238.4%
associate--l+45.4%
metadata-eval45.4%
+-rgt-identity45.4%
unpow245.4%
distribute-frac-neg45.4%
distribute-frac-neg45.4%
sqr-neg45.4%
unpow245.4%
Simplified45.4%
*-un-lft-identity45.4%
unpow245.4%
times-frac51.2%
clear-num51.2%
associate-/r/51.2%
associate-/r*51.2%
pow251.2%
Applied egg-rr51.2%
associate-*l/51.2%
associate-*r/51.2%
associate-*r/51.2%
associate-/l/51.6%
*-lft-identity51.6%
associate-*l/51.6%
associate-*r/51.6%
associate-*l/51.6%
*-lft-identity51.6%
Simplified51.6%
associate-*l/52.0%
associate-/l/52.0%
Applied egg-rr52.0%
Taylor expanded in t around 0 81.6%
*-commutative81.6%
Simplified81.6%
times-frac82.1%
*-commutative82.1%
Applied egg-rr82.1%
Final simplification82.1%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (* (pow l 2.0) (cos k_m))))
(if (<= k_m 9.2e-5)
(/ (* 2.0 (/ t_1 (* k_m (* t (pow k_m 2.0))))) k_m)
(/ (* 2.0 (/ t_1 (* k_m (* t (- 0.5 (/ (cos (* k_m 2.0)) 2.0)))))) k_m))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = pow(l, 2.0) * cos(k_m);
double tmp;
if (k_m <= 9.2e-5) {
tmp = (2.0 * (t_1 / (k_m * (t * pow(k_m, 2.0))))) / k_m;
} else {
tmp = (2.0 * (t_1 / (k_m * (t * (0.5 - (cos((k_m * 2.0)) / 2.0)))))) / k_m;
}
return tmp;
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: t_1
real(8) :: tmp
t_1 = (l ** 2.0d0) * cos(k_m)
if (k_m <= 9.2d-5) then
tmp = (2.0d0 * (t_1 / (k_m * (t * (k_m ** 2.0d0))))) / k_m
else
tmp = (2.0d0 * (t_1 / (k_m * (t * (0.5d0 - (cos((k_m * 2.0d0)) / 2.0d0)))))) / k_m
end if
code = tmp
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double t_1 = Math.pow(l, 2.0) * Math.cos(k_m);
double tmp;
if (k_m <= 9.2e-5) {
tmp = (2.0 * (t_1 / (k_m * (t * Math.pow(k_m, 2.0))))) / k_m;
} else {
tmp = (2.0 * (t_1 / (k_m * (t * (0.5 - (Math.cos((k_m * 2.0)) / 2.0)))))) / k_m;
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): t_1 = math.pow(l, 2.0) * math.cos(k_m) tmp = 0 if k_m <= 9.2e-5: tmp = (2.0 * (t_1 / (k_m * (t * math.pow(k_m, 2.0))))) / k_m else: tmp = (2.0 * (t_1 / (k_m * (t * (0.5 - (math.cos((k_m * 2.0)) / 2.0)))))) / k_m return tmp
k_m = abs(k) function code(t, l, k_m) t_1 = Float64((l ^ 2.0) * cos(k_m)) tmp = 0.0 if (k_m <= 9.2e-5) tmp = Float64(Float64(2.0 * Float64(t_1 / Float64(k_m * Float64(t * (k_m ^ 2.0))))) / k_m); else tmp = Float64(Float64(2.0 * Float64(t_1 / Float64(k_m * Float64(t * Float64(0.5 - Float64(cos(Float64(k_m * 2.0)) / 2.0)))))) / k_m); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) t_1 = (l ^ 2.0) * cos(k_m); tmp = 0.0; if (k_m <= 9.2e-5) tmp = (2.0 * (t_1 / (k_m * (t * (k_m ^ 2.0))))) / k_m; else tmp = (2.0 * (t_1 / (k_m * (t * (0.5 - (cos((k_m * 2.0)) / 2.0)))))) / k_m; end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k$95$m, 9.2e-5], N[(N[(2.0 * N[(t$95$1 / N[(k$95$m * N[(t * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision], N[(N[(2.0 * N[(t$95$1 / N[(k$95$m * N[(t * N[(0.5 - N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := {\ell}^{2} \cdot \cos k_m\\
\mathbf{if}\;k_m \leq 9.2 \cdot 10^{-5}:\\
\;\;\;\;\frac{2 \cdot \frac{t_1}{k_m \cdot \left(t \cdot {k_m}^{2}\right)}}{k_m}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \frac{t_1}{k_m \cdot \left(t \cdot \left(0.5 - \frac{\cos \left(k_m \cdot 2\right)}{2}\right)\right)}}{k_m}\\
\end{array}
\end{array}
if k < 9.20000000000000001e-5Initial program 41.1%
associate-/r*41.1%
associate-*l*41.1%
associate-*l/41.1%
associate-/l*41.1%
+-commutative41.1%
unpow241.1%
sqr-neg41.1%
distribute-frac-neg41.1%
distribute-frac-neg41.1%
unpow241.1%
associate--l+48.8%
metadata-eval48.8%
+-rgt-identity48.8%
unpow248.8%
distribute-frac-neg48.8%
distribute-frac-neg48.8%
sqr-neg48.8%
unpow248.8%
Simplified48.8%
*-un-lft-identity48.8%
unpow248.8%
times-frac52.8%
clear-num52.8%
associate-/r/52.8%
associate-/r*52.8%
pow252.8%
Applied egg-rr52.8%
associate-*l/52.9%
associate-*r/52.9%
associate-*r/52.9%
associate-/l/52.9%
*-lft-identity52.9%
associate-*l/52.9%
associate-*r/52.9%
associate-*l/52.9%
*-lft-identity52.9%
Simplified52.9%
associate-*l/53.5%
associate-/l/53.5%
Applied egg-rr53.5%
Taylor expanded in t around 0 81.8%
*-commutative81.8%
Simplified81.8%
Taylor expanded in k around 0 75.4%
if 9.20000000000000001e-5 < k Initial program 32.3%
associate-/r*32.1%
associate-*l*32.1%
associate-*l/32.1%
associate-/l*32.1%
+-commutative32.1%
unpow232.1%
sqr-neg32.1%
distribute-frac-neg32.1%
distribute-frac-neg32.1%
unpow232.1%
associate--l+37.4%
metadata-eval37.4%
+-rgt-identity37.4%
unpow237.4%
distribute-frac-neg37.4%
distribute-frac-neg37.4%
sqr-neg37.4%
unpow237.4%
Simplified37.4%
*-un-lft-identity37.4%
unpow237.4%
times-frac47.4%
clear-num47.4%
associate-/r/47.4%
associate-/r*47.4%
pow247.4%
Applied egg-rr47.4%
associate-*l/47.4%
associate-*r/47.4%
associate-*r/47.4%
associate-/l/48.6%
*-lft-identity48.6%
associate-*l/48.6%
associate-*r/48.7%
associate-*l/48.7%
*-lft-identity48.7%
Simplified48.7%
associate-*l/48.7%
associate-/l/48.7%
Applied egg-rr48.7%
Taylor expanded in t around 0 81.3%
*-commutative81.3%
Simplified81.3%
unpow281.3%
sin-mult81.3%
Applied egg-rr81.3%
div-sub81.3%
+-inverses81.3%
cos-081.3%
metadata-eval81.3%
count-281.3%
*-commutative81.3%
Simplified81.3%
Final simplification77.2%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (/ (* 2.0 (/ (* (pow l 2.0) (cos k_m)) (* k_m (* t (pow k_m 2.0))))) k_m))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return (2.0 * ((pow(l, 2.0) * cos(k_m)) / (k_m * (t * pow(k_m, 2.0))))) / k_m;
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = (2.0d0 * (((l ** 2.0d0) * cos(k_m)) / (k_m * (t * (k_m ** 2.0d0))))) / k_m
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return (2.0 * ((Math.pow(l, 2.0) * Math.cos(k_m)) / (k_m * (t * Math.pow(k_m, 2.0))))) / k_m;
}
k_m = math.fabs(k) def code(t, l, k_m): return (2.0 * ((math.pow(l, 2.0) * math.cos(k_m)) / (k_m * (t * math.pow(k_m, 2.0))))) / k_m
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(2.0 * Float64(Float64((l ^ 2.0) * cos(k_m)) / Float64(k_m * Float64(t * (k_m ^ 2.0))))) / k_m) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = (2.0 * (((l ^ 2.0) * cos(k_m)) / (k_m * (t * (k_m ^ 2.0))))) / k_m; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(k$95$m * N[(t * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2 \cdot \frac{{\ell}^{2} \cdot \cos k_m}{k_m \cdot \left(t \cdot {k_m}^{2}\right)}}{k_m}
\end{array}
Initial program 38.4%
associate-/r*38.4%
associate-*l*38.4%
associate-*l/38.4%
associate-/l*38.4%
+-commutative38.4%
unpow238.4%
sqr-neg38.4%
distribute-frac-neg38.4%
distribute-frac-neg38.4%
unpow238.4%
associate--l+45.4%
metadata-eval45.4%
+-rgt-identity45.4%
unpow245.4%
distribute-frac-neg45.4%
distribute-frac-neg45.4%
sqr-neg45.4%
unpow245.4%
Simplified45.4%
*-un-lft-identity45.4%
unpow245.4%
times-frac51.2%
clear-num51.2%
associate-/r/51.2%
associate-/r*51.2%
pow251.2%
Applied egg-rr51.2%
associate-*l/51.2%
associate-*r/51.2%
associate-*r/51.2%
associate-/l/51.6%
*-lft-identity51.6%
associate-*l/51.6%
associate-*r/51.6%
associate-*l/51.6%
*-lft-identity51.6%
Simplified51.6%
associate-*l/52.0%
associate-/l/52.0%
Applied egg-rr52.0%
Taylor expanded in t around 0 81.6%
*-commutative81.6%
Simplified81.6%
Taylor expanded in k around 0 69.9%
Final simplification69.9%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (/ (* 2.0 (/ (pow l 2.0) (* t (pow k_m 3.0)))) k_m))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return (2.0 * (pow(l, 2.0) / (t * pow(k_m, 3.0)))) / k_m;
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = (2.0d0 * ((l ** 2.0d0) / (t * (k_m ** 3.0d0)))) / k_m
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return (2.0 * (Math.pow(l, 2.0) / (t * Math.pow(k_m, 3.0)))) / k_m;
}
k_m = math.fabs(k) def code(t, l, k_m): return (2.0 * (math.pow(l, 2.0) / (t * math.pow(k_m, 3.0)))) / k_m
k_m = abs(k) function code(t, l, k_m) return Float64(Float64(2.0 * Float64((l ^ 2.0) / Float64(t * (k_m ^ 3.0)))) / k_m) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = (2.0 * ((l ^ 2.0) / (t * (k_m ^ 3.0)))) / k_m; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] / N[(t * N[Power[k$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2 \cdot \frac{{\ell}^{2}}{t \cdot {k_m}^{3}}}{k_m}
\end{array}
Initial program 38.4%
associate-/r*38.4%
associate-*l*38.4%
associate-*l/38.4%
associate-/l*38.4%
+-commutative38.4%
unpow238.4%
sqr-neg38.4%
distribute-frac-neg38.4%
distribute-frac-neg38.4%
unpow238.4%
associate--l+45.4%
metadata-eval45.4%
+-rgt-identity45.4%
unpow245.4%
distribute-frac-neg45.4%
distribute-frac-neg45.4%
sqr-neg45.4%
unpow245.4%
Simplified45.4%
*-un-lft-identity45.4%
unpow245.4%
times-frac51.2%
clear-num51.2%
associate-/r/51.2%
associate-/r*51.2%
pow251.2%
Applied egg-rr51.2%
associate-*l/51.2%
associate-*r/51.2%
associate-*r/51.2%
associate-/l/51.6%
*-lft-identity51.6%
associate-*l/51.6%
associate-*r/51.6%
associate-*l/51.6%
*-lft-identity51.6%
Simplified51.6%
associate-*l/52.0%
associate-/l/52.0%
Applied egg-rr52.0%
Taylor expanded in k around 0 67.1%
Final simplification67.1%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (/ 2.0 (* t (* (pow k_m 4.0) (pow l -2.0)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / (t * (pow(k_m, 4.0) * pow(l, -2.0)));
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = 2.0d0 / (t * ((k_m ** 4.0d0) * (l ** (-2.0d0))))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return 2.0 / (t * (Math.pow(k_m, 4.0) * Math.pow(l, -2.0)));
}
k_m = math.fabs(k) def code(t, l, k_m): return 2.0 / (t * (math.pow(k_m, 4.0) * math.pow(l, -2.0)))
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64(t * Float64((k_m ^ 4.0) * (l ^ -2.0)))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = 2.0 / (t * ((k_m ^ 4.0) * (l ^ -2.0))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(t * N[(N[Power[k$95$m, 4.0], $MachinePrecision] * N[Power[l, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{t \cdot \left({k_m}^{4} \cdot {\ell}^{-2}\right)}
\end{array}
Initial program 38.4%
associate-*l*38.4%
associate--l+38.4%
Simplified38.4%
unpow338.4%
times-frac45.2%
pow245.2%
Applied egg-rr45.2%
Taylor expanded in k around 0 64.9%
associate-/l*63.2%
associate-/r/64.5%
Simplified64.5%
expm1-log1p-u64.5%
expm1-udef62.2%
div-inv62.2%
pow-flip62.2%
metadata-eval62.2%
Applied egg-rr62.2%
expm1-def64.5%
expm1-log1p64.5%
Simplified64.5%
Final simplification64.5%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (/ 2.0 (* (pow l -2.0) (* (pow k_m 4.0) t))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / (pow(l, -2.0) * (pow(k_m, 4.0) * t));
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = 2.0d0 / ((l ** (-2.0d0)) * ((k_m ** 4.0d0) * t))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return 2.0 / (Math.pow(l, -2.0) * (Math.pow(k_m, 4.0) * t));
}
k_m = math.fabs(k) def code(t, l, k_m): return 2.0 / (math.pow(l, -2.0) * (math.pow(k_m, 4.0) * t))
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64((l ^ -2.0) * Float64((k_m ^ 4.0) * t))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = 2.0 / ((l ^ -2.0) * ((k_m ^ 4.0) * t)); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(N[Power[l, -2.0], $MachinePrecision] * N[(N[Power[k$95$m, 4.0], $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{{\ell}^{-2} \cdot \left({k_m}^{4} \cdot t\right)}
\end{array}
Initial program 38.4%
associate-*l*38.4%
associate--l+38.4%
Simplified38.4%
Taylor expanded in k around 0 64.9%
expm1-log1p-u44.2%
expm1-udef23.0%
div-inv23.0%
*-commutative23.0%
pow-flip23.0%
metadata-eval23.0%
Applied egg-rr23.0%
expm1-def44.2%
expm1-log1p64.9%
*-commutative64.9%
*-commutative64.9%
Simplified64.9%
Final simplification64.9%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (/ 2.0 (/ (* (pow k_m 4.0) t) (pow l 2.0))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / ((pow(k_m, 4.0) * t) / pow(l, 2.0));
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = 2.0d0 / (((k_m ** 4.0d0) * t) / (l ** 2.0d0))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return 2.0 / ((Math.pow(k_m, 4.0) * t) / Math.pow(l, 2.0));
}
k_m = math.fabs(k) def code(t, l, k_m): return 2.0 / ((math.pow(k_m, 4.0) * t) / math.pow(l, 2.0))
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64(Float64((k_m ^ 4.0) * t) / (l ^ 2.0))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = 2.0 / (((k_m ^ 4.0) * t) / (l ^ 2.0)); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(N[(N[Power[k$95$m, 4.0], $MachinePrecision] * t), $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{\frac{{k_m}^{4} \cdot t}{{\ell}^{2}}}
\end{array}
Initial program 38.4%
associate-*l*38.4%
associate--l+38.4%
Simplified38.4%
Taylor expanded in k around 0 64.9%
Final simplification64.9%
herbie shell --seed 2024021
(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))))