
(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}
(FPCore (t l k) :precision binary64 (/ 2.0 (/ (* (* k (/ (* k t) l)) (* (sin k) (* (sin k) (/ 1.0 l)))) (cos k))))
double code(double t, double l, double k) {
return 2.0 / (((k * ((k * t) / l)) * (sin(k) * (sin(k) * (1.0 / l)))) / cos(k));
}
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 / (((k * ((k * t) / l)) * (sin(k) * (sin(k) * (1.0d0 / l)))) / cos(k))
end function
public static double code(double t, double l, double k) {
return 2.0 / (((k * ((k * t) / l)) * (Math.sin(k) * (Math.sin(k) * (1.0 / l)))) / Math.cos(k));
}
def code(t, l, k): return 2.0 / (((k * ((k * t) / l)) * (math.sin(k) * (math.sin(k) * (1.0 / l)))) / math.cos(k))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(k * Float64(Float64(k * t) / l)) * Float64(sin(k) * Float64(sin(k) * Float64(1.0 / l)))) / cos(k))) end
function tmp = code(t, l, k) tmp = 2.0 / (((k * ((k * t) / l)) * (sin(k) * (sin(k) * (1.0 / l)))) / cos(k)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(k * N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(1.0 / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\frac{\left(k \cdot \frac{k \cdot t}{\ell}\right) \cdot \left(\sin k \cdot \left(\sin k \cdot \frac{1}{\ell}\right)\right)}{\cos k}}
\end{array}
Initial program 35.4%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6485.4%
Simplified85.4%
div-invN/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6492.7%
Applied egg-rr92.7%
div-invN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f6494.8%
Applied egg-rr94.8%
(FPCore (t l k) :precision binary64 (/ 2.0 (/ (* (* k (/ (* k t) l)) (/ (sin k) (/ l (sin k)))) (cos k))))
double code(double t, double l, double k) {
return 2.0 / (((k * ((k * t) / l)) * (sin(k) / (l / sin(k)))) / cos(k));
}
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 / (((k * ((k * t) / l)) * (sin(k) / (l / sin(k)))) / cos(k))
end function
public static double code(double t, double l, double k) {
return 2.0 / (((k * ((k * t) / l)) * (Math.sin(k) / (l / Math.sin(k)))) / Math.cos(k));
}
def code(t, l, k): return 2.0 / (((k * ((k * t) / l)) * (math.sin(k) / (l / math.sin(k)))) / math.cos(k))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(k * Float64(Float64(k * t) / l)) * Float64(sin(k) / Float64(l / sin(k)))) / cos(k))) end
function tmp = code(t, l, k) tmp = 2.0 / (((k * ((k * t) / l)) * (sin(k) / (l / sin(k)))) / cos(k)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(k * N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] / N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\frac{\left(k \cdot \frac{k \cdot t}{\ell}\right) \cdot \frac{\sin k}{\frac{\ell}{\sin k}}}{\cos k}}
\end{array}
Initial program 35.4%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6485.4%
Simplified85.4%
div-invN/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6492.7%
Applied egg-rr92.7%
div-invN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f6494.8%
Applied egg-rr94.8%
un-div-invN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6494.8%
Applied egg-rr94.8%
(FPCore (t l k) :precision binary64 (/ 2.0 (/ (* (* k (/ (* k t) l)) (* (sin k) (/ (sin k) l))) (cos k))))
double code(double t, double l, double k) {
return 2.0 / (((k * ((k * t) / l)) * (sin(k) * (sin(k) / l))) / cos(k));
}
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 / (((k * ((k * t) / l)) * (sin(k) * (sin(k) / l))) / cos(k))
end function
public static double code(double t, double l, double k) {
return 2.0 / (((k * ((k * t) / l)) * (Math.sin(k) * (Math.sin(k) / l))) / Math.cos(k));
}
def code(t, l, k): return 2.0 / (((k * ((k * t) / l)) * (math.sin(k) * (math.sin(k) / l))) / math.cos(k))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(k * Float64(Float64(k * t) / l)) * Float64(sin(k) * Float64(sin(k) / l))) / cos(k))) end
function tmp = code(t, l, k) tmp = 2.0 / (((k * ((k * t) / l)) * (sin(k) * (sin(k) / l))) / cos(k)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(k * N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\frac{\left(k \cdot \frac{k \cdot t}{\ell}\right) \cdot \left(\sin k \cdot \frac{\sin k}{\ell}\right)}{\cos k}}
\end{array}
Initial program 35.4%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6485.4%
Simplified85.4%
div-invN/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6492.7%
Applied egg-rr92.7%
unpow2N/A
associate-/l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6494.8%
Applied egg-rr94.8%
(FPCore (t l k)
:precision binary64
(if (<= k 4e-116)
(/
(/ 2.0 (/ (* t (* k k)) l))
(* k (* k (+ (/ 1.0 l) (/ (* (* k k) 0.16666666666666666) l)))))
(/ 2.0 (* (/ k (/ l t)) (* k (* (/ (sin k) l) (tan k)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 4e-116) {
tmp = (2.0 / ((t * (k * k)) / l)) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l))));
} else {
tmp = 2.0 / ((k / (l / t)) * (k * ((sin(k) / l) * tan(k))));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 4d-116) then
tmp = (2.0d0 / ((t * (k * k)) / l)) / (k * (k * ((1.0d0 / l) + (((k * k) * 0.16666666666666666d0) / l))))
else
tmp = 2.0d0 / ((k / (l / t)) * (k * ((sin(k) / l) * tan(k))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 4e-116) {
tmp = (2.0 / ((t * (k * k)) / l)) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l))));
} else {
tmp = 2.0 / ((k / (l / t)) * (k * ((Math.sin(k) / l) * Math.tan(k))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 4e-116: tmp = (2.0 / ((t * (k * k)) / l)) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l)))) else: tmp = 2.0 / ((k / (l / t)) * (k * ((math.sin(k) / l) * math.tan(k)))) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 4e-116) tmp = Float64(Float64(2.0 / Float64(Float64(t * Float64(k * k)) / l)) / Float64(k * Float64(k * Float64(Float64(1.0 / l) + Float64(Float64(Float64(k * k) * 0.16666666666666666) / l))))); else tmp = Float64(2.0 / Float64(Float64(k / Float64(l / t)) * Float64(k * Float64(Float64(sin(k) / l) * tan(k))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 4e-116) tmp = (2.0 / ((t * (k * k)) / l)) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l)))); else tmp = 2.0 / ((k / (l / t)) * (k * ((sin(k) / l) * tan(k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 4e-116], N[(N[(2.0 / N[(N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / N[(k * N[(k * N[(N[(1.0 / l), $MachinePrecision] + N[(N[(N[(k * k), $MachinePrecision] * 0.16666666666666666), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k / N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(k * N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 4 \cdot 10^{-116}:\\
\;\;\;\;\frac{\frac{2}{\frac{t \cdot \left(k \cdot k\right)}{\ell}}}{k \cdot \left(k \cdot \left(\frac{1}{\ell} + \frac{\left(k \cdot k\right) \cdot 0.16666666666666666}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{k}{\frac{\ell}{t}} \cdot \left(k \cdot \left(\frac{\sin k}{\ell} \cdot \tan k\right)\right)}\\
\end{array}
\end{array}
if k < 4e-116Initial program 39.7%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6487.0%
Simplified87.0%
associate-/l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
unpow2N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6476.4%
Applied egg-rr76.4%
Taylor expanded in k around 0
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
Simplified78.3%
if 4e-116 < k Initial program 27.7%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6482.5%
Simplified82.5%
div-invN/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6496.1%
Applied egg-rr96.1%
div-invN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f6496.1%
Applied egg-rr96.1%
associate-*r/N/A
associate-*l*N/A
associate-/l*N/A
associate-*l*N/A
associate-*r/N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
Applied egg-rr94.7%
Final simplification84.1%
(FPCore (t l k)
:precision binary64
(if (<= k 2.2e-116)
(/
(/ 2.0 (/ (* t (* k k)) l))
(* k (* k (+ (/ 1.0 l) (/ (* (* k k) 0.16666666666666666) l)))))
(/ 2.0 (* k (* (/ k (/ l t)) (* (/ (sin k) l) (tan k)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 2.2e-116) {
tmp = (2.0 / ((t * (k * k)) / l)) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l))));
} else {
tmp = 2.0 / (k * ((k / (l / t)) * ((sin(k) / l) * tan(k))));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 2.2d-116) then
tmp = (2.0d0 / ((t * (k * k)) / l)) / (k * (k * ((1.0d0 / l) + (((k * k) * 0.16666666666666666d0) / l))))
else
tmp = 2.0d0 / (k * ((k / (l / t)) * ((sin(k) / l) * tan(k))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 2.2e-116) {
tmp = (2.0 / ((t * (k * k)) / l)) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l))));
} else {
tmp = 2.0 / (k * ((k / (l / t)) * ((Math.sin(k) / l) * Math.tan(k))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 2.2e-116: tmp = (2.0 / ((t * (k * k)) / l)) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l)))) else: tmp = 2.0 / (k * ((k / (l / t)) * ((math.sin(k) / l) * math.tan(k)))) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 2.2e-116) tmp = Float64(Float64(2.0 / Float64(Float64(t * Float64(k * k)) / l)) / Float64(k * Float64(k * Float64(Float64(1.0 / l) + Float64(Float64(Float64(k * k) * 0.16666666666666666) / l))))); else tmp = Float64(2.0 / Float64(k * Float64(Float64(k / Float64(l / t)) * Float64(Float64(sin(k) / l) * tan(k))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 2.2e-116) tmp = (2.0 / ((t * (k * k)) / l)) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l)))); else tmp = 2.0 / (k * ((k / (l / t)) * ((sin(k) / l) * tan(k)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 2.2e-116], N[(N[(2.0 / N[(N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / N[(k * N[(k * N[(N[(1.0 / l), $MachinePrecision] + N[(N[(N[(k * k), $MachinePrecision] * 0.16666666666666666), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(k * N[(N[(k / N[(l / t), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 2.2 \cdot 10^{-116}:\\
\;\;\;\;\frac{\frac{2}{\frac{t \cdot \left(k \cdot k\right)}{\ell}}}{k \cdot \left(k \cdot \left(\frac{1}{\ell} + \frac{\left(k \cdot k\right) \cdot 0.16666666666666666}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{k \cdot \left(\frac{k}{\frac{\ell}{t}} \cdot \left(\frac{\sin k}{\ell} \cdot \tan k\right)\right)}\\
\end{array}
\end{array}
if k < 2.2000000000000001e-116Initial program 39.7%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6487.0%
Simplified87.0%
associate-/l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
unpow2N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6476.4%
Applied egg-rr76.4%
Taylor expanded in k around 0
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
Simplified78.3%
if 2.2000000000000001e-116 < k Initial program 27.7%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6482.5%
Simplified82.5%
div-invN/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6496.1%
Applied egg-rr96.1%
div-invN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f6496.1%
Applied egg-rr96.1%
associate-*r/N/A
associate-*l*N/A
associate-/l*N/A
associate-*l*N/A
associate-*r/N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/l*N/A
Applied egg-rr93.6%
Final simplification83.8%
(FPCore (t l k) :precision binary64 (* 2.0 (/ (/ l (* t (* k k))) (* (/ (sin k) l) (tan k)))))
double code(double t, double l, double k) {
return 2.0 * ((l / (t * (k * k))) / ((sin(k) / l) * tan(k)));
}
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 * ((l / (t * (k * k))) / ((sin(k) / l) * tan(k)))
end function
public static double code(double t, double l, double k) {
return 2.0 * ((l / (t * (k * k))) / ((Math.sin(k) / l) * Math.tan(k)));
}
def code(t, l, k): return 2.0 * ((l / (t * (k * k))) / ((math.sin(k) / l) * math.tan(k)))
function code(t, l, k) return Float64(2.0 * Float64(Float64(l / Float64(t * Float64(k * k))) / Float64(Float64(sin(k) / l) * tan(k)))) end
function tmp = code(t, l, k) tmp = 2.0 * ((l / (t * (k * k))) / ((sin(k) / l) * tan(k))); end
code[t_, l_, k_] := N[(2.0 * N[(N[(l / N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \frac{\frac{\ell}{t \cdot \left(k \cdot k\right)}}{\frac{\sin k}{\ell} \cdot \tan k}
\end{array}
Initial program 35.4%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6485.4%
Simplified85.4%
div-invN/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6492.7%
Applied egg-rr92.7%
div-invN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f6494.8%
Applied egg-rr94.8%
clear-numN/A
associate-/r/N/A
clear-numN/A
*-lowering-*.f64N/A
Applied egg-rr85.7%
Final simplification85.7%
(FPCore (t l k) :precision binary64 (/ 2.0 (/ (* (* k (/ (* k t) l)) (/ (* k k) l)) (cos k))))
double code(double t, double l, double k) {
return 2.0 / (((k * ((k * t) / l)) * ((k * k) / l)) / cos(k));
}
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 / (((k * ((k * t) / l)) * ((k * k) / l)) / cos(k))
end function
public static double code(double t, double l, double k) {
return 2.0 / (((k * ((k * t) / l)) * ((k * k) / l)) / Math.cos(k));
}
def code(t, l, k): return 2.0 / (((k * ((k * t) / l)) * ((k * k) / l)) / math.cos(k))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(k * Float64(Float64(k * t) / l)) * Float64(Float64(k * k) / l)) / cos(k))) end
function tmp = code(t, l, k) tmp = 2.0 / (((k * ((k * t) / l)) * ((k * k) / l)) / cos(k)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(k * N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\frac{\left(k \cdot \frac{k \cdot t}{\ell}\right) \cdot \frac{k \cdot k}{\ell}}{\cos k}}
\end{array}
Initial program 35.4%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6485.4%
Simplified85.4%
div-invN/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6492.7%
Applied egg-rr92.7%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6474.5%
Simplified74.5%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ (* t (* k k)) l)))
(if (<= (* l l) 1e-22)
(/
(/ 2.0 t_1)
(*
(* k k)
(+
(/ 1.0 l)
(*
(* k k)
(+
(/ 0.16666666666666666 l)
(*
(* k k)
(+
(* (/ (* k k) l) 0.03432539682539683)
(/ 0.08611111111111111 l))))))))
(/ 2.0 (* (* k k) (/ t_1 l))))))
double code(double t, double l, double k) {
double t_1 = (t * (k * k)) / l;
double tmp;
if ((l * l) <= 1e-22) {
tmp = (2.0 / t_1) / ((k * k) * ((1.0 / l) + ((k * k) * ((0.16666666666666666 / l) + ((k * k) * ((((k * k) / l) * 0.03432539682539683) + (0.08611111111111111 / l)))))));
} else {
tmp = 2.0 / ((k * k) * (t_1 / l));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = (t * (k * k)) / l
if ((l * l) <= 1d-22) then
tmp = (2.0d0 / t_1) / ((k * k) * ((1.0d0 / l) + ((k * k) * ((0.16666666666666666d0 / l) + ((k * k) * ((((k * k) / l) * 0.03432539682539683d0) + (0.08611111111111111d0 / l)))))))
else
tmp = 2.0d0 / ((k * k) * (t_1 / l))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = (t * (k * k)) / l;
double tmp;
if ((l * l) <= 1e-22) {
tmp = (2.0 / t_1) / ((k * k) * ((1.0 / l) + ((k * k) * ((0.16666666666666666 / l) + ((k * k) * ((((k * k) / l) * 0.03432539682539683) + (0.08611111111111111 / l)))))));
} else {
tmp = 2.0 / ((k * k) * (t_1 / l));
}
return tmp;
}
def code(t, l, k): t_1 = (t * (k * k)) / l tmp = 0 if (l * l) <= 1e-22: tmp = (2.0 / t_1) / ((k * k) * ((1.0 / l) + ((k * k) * ((0.16666666666666666 / l) + ((k * k) * ((((k * k) / l) * 0.03432539682539683) + (0.08611111111111111 / l))))))) else: tmp = 2.0 / ((k * k) * (t_1 / l)) return tmp
function code(t, l, k) t_1 = Float64(Float64(t * Float64(k * k)) / l) tmp = 0.0 if (Float64(l * l) <= 1e-22) tmp = Float64(Float64(2.0 / t_1) / Float64(Float64(k * k) * Float64(Float64(1.0 / l) + Float64(Float64(k * k) * Float64(Float64(0.16666666666666666 / l) + Float64(Float64(k * k) * Float64(Float64(Float64(Float64(k * k) / l) * 0.03432539682539683) + Float64(0.08611111111111111 / l)))))))); else tmp = Float64(2.0 / Float64(Float64(k * k) * Float64(t_1 / l))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = (t * (k * k)) / l; tmp = 0.0; if ((l * l) <= 1e-22) tmp = (2.0 / t_1) / ((k * k) * ((1.0 / l) + ((k * k) * ((0.16666666666666666 / l) + ((k * k) * ((((k * k) / l) * 0.03432539682539683) + (0.08611111111111111 / l))))))); else tmp = 2.0 / ((k * k) * (t_1 / l)); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]}, If[LessEqual[N[(l * l), $MachinePrecision], 1e-22], N[(N[(2.0 / t$95$1), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] * N[(N[(1.0 / l), $MachinePrecision] + N[(N[(k * k), $MachinePrecision] * N[(N[(0.16666666666666666 / l), $MachinePrecision] + N[(N[(k * k), $MachinePrecision] * N[(N[(N[(N[(k * k), $MachinePrecision] / l), $MachinePrecision] * 0.03432539682539683), $MachinePrecision] + N[(0.08611111111111111 / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k * k), $MachinePrecision] * N[(t$95$1 / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t \cdot \left(k \cdot k\right)}{\ell}\\
\mathbf{if}\;\ell \cdot \ell \leq 10^{-22}:\\
\;\;\;\;\frac{\frac{2}{t\_1}}{\left(k \cdot k\right) \cdot \left(\frac{1}{\ell} + \left(k \cdot k\right) \cdot \left(\frac{0.16666666666666666}{\ell} + \left(k \cdot k\right) \cdot \left(\frac{k \cdot k}{\ell} \cdot 0.03432539682539683 + \frac{0.08611111111111111}{\ell}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(k \cdot k\right) \cdot \frac{t\_1}{\ell}}\\
\end{array}
\end{array}
if (*.f64 l l) < 1e-22Initial program 30.7%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6490.6%
Simplified90.6%
associate-/l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
unpow2N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6472.6%
Applied egg-rr72.6%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
Simplified83.6%
if 1e-22 < (*.f64 l l) Initial program 40.2%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6456.7%
Simplified56.7%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-/r*N/A
associate-*r/N/A
associate-/l*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6465.7%
Applied egg-rr65.7%
Final simplification74.7%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ (* t (* k k)) l)))
(if (<= (* l l) 1e-22)
(/
(/ 2.0 t_1)
(* k (* k (+ (/ 1.0 l) (/ (* (* k k) 0.16666666666666666) l)))))
(/ 2.0 (* (* k k) (/ t_1 l))))))
double code(double t, double l, double k) {
double t_1 = (t * (k * k)) / l;
double tmp;
if ((l * l) <= 1e-22) {
tmp = (2.0 / t_1) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l))));
} else {
tmp = 2.0 / ((k * k) * (t_1 / l));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = (t * (k * k)) / l
if ((l * l) <= 1d-22) then
tmp = (2.0d0 / t_1) / (k * (k * ((1.0d0 / l) + (((k * k) * 0.16666666666666666d0) / l))))
else
tmp = 2.0d0 / ((k * k) * (t_1 / l))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = (t * (k * k)) / l;
double tmp;
if ((l * l) <= 1e-22) {
tmp = (2.0 / t_1) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l))));
} else {
tmp = 2.0 / ((k * k) * (t_1 / l));
}
return tmp;
}
def code(t, l, k): t_1 = (t * (k * k)) / l tmp = 0 if (l * l) <= 1e-22: tmp = (2.0 / t_1) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l)))) else: tmp = 2.0 / ((k * k) * (t_1 / l)) return tmp
function code(t, l, k) t_1 = Float64(Float64(t * Float64(k * k)) / l) tmp = 0.0 if (Float64(l * l) <= 1e-22) tmp = Float64(Float64(2.0 / t_1) / Float64(k * Float64(k * Float64(Float64(1.0 / l) + Float64(Float64(Float64(k * k) * 0.16666666666666666) / l))))); else tmp = Float64(2.0 / Float64(Float64(k * k) * Float64(t_1 / l))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = (t * (k * k)) / l; tmp = 0.0; if ((l * l) <= 1e-22) tmp = (2.0 / t_1) / (k * (k * ((1.0 / l) + (((k * k) * 0.16666666666666666) / l)))); else tmp = 2.0 / ((k * k) * (t_1 / l)); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]}, If[LessEqual[N[(l * l), $MachinePrecision], 1e-22], N[(N[(2.0 / t$95$1), $MachinePrecision] / N[(k * N[(k * N[(N[(1.0 / l), $MachinePrecision] + N[(N[(N[(k * k), $MachinePrecision] * 0.16666666666666666), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k * k), $MachinePrecision] * N[(t$95$1 / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t \cdot \left(k \cdot k\right)}{\ell}\\
\mathbf{if}\;\ell \cdot \ell \leq 10^{-22}:\\
\;\;\;\;\frac{\frac{2}{t\_1}}{k \cdot \left(k \cdot \left(\frac{1}{\ell} + \frac{\left(k \cdot k\right) \cdot 0.16666666666666666}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(k \cdot k\right) \cdot \frac{t\_1}{\ell}}\\
\end{array}
\end{array}
if (*.f64 l l) < 1e-22Initial program 30.7%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6490.6%
Simplified90.6%
associate-/l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
unpow2N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6472.6%
Applied egg-rr72.6%
Taylor expanded in k around 0
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
Simplified83.6%
if 1e-22 < (*.f64 l l) Initial program 40.2%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6456.7%
Simplified56.7%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-/r*N/A
associate-*r/N/A
associate-/l*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6465.7%
Applied egg-rr65.7%
Final simplification74.7%
(FPCore (t l k) :precision binary64 (if (<= l 2e-161) (* l (/ (/ 2.0 (* k (* k (* k k)))) (/ t l))) (* (/ (/ 2.0 t) (* k k)) (/ (* l l) (* k k)))))
double code(double t, double l, double k) {
double tmp;
if (l <= 2e-161) {
tmp = l * ((2.0 / (k * (k * (k * k)))) / (t / l));
} else {
tmp = ((2.0 / t) / (k * k)) * ((l * l) / (k * k));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (l <= 2d-161) then
tmp = l * ((2.0d0 / (k * (k * (k * k)))) / (t / l))
else
tmp = ((2.0d0 / t) / (k * k)) * ((l * l) / (k * k))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (l <= 2e-161) {
tmp = l * ((2.0 / (k * (k * (k * k)))) / (t / l));
} else {
tmp = ((2.0 / t) / (k * k)) * ((l * l) / (k * k));
}
return tmp;
}
def code(t, l, k): tmp = 0 if l <= 2e-161: tmp = l * ((2.0 / (k * (k * (k * k)))) / (t / l)) else: tmp = ((2.0 / t) / (k * k)) * ((l * l) / (k * k)) return tmp
function code(t, l, k) tmp = 0.0 if (l <= 2e-161) tmp = Float64(l * Float64(Float64(2.0 / Float64(k * Float64(k * Float64(k * k)))) / Float64(t / l))); else tmp = Float64(Float64(Float64(2.0 / t) / Float64(k * k)) * Float64(Float64(l * l) / Float64(k * k))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (l <= 2e-161) tmp = l * ((2.0 / (k * (k * (k * k)))) / (t / l)); else tmp = ((2.0 / t) / (k * k)) * ((l * l) / (k * k)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[l, 2e-161], N[(l * N[(N[(2.0 / N[(k * N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / t), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \leq 2 \cdot 10^{-161}:\\
\;\;\;\;\ell \cdot \frac{\frac{2}{k \cdot \left(k \cdot \left(k \cdot k\right)\right)}}{\frac{t}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{t}}{k \cdot k} \cdot \frac{\ell \cdot \ell}{k \cdot k}\\
\end{array}
\end{array}
if l < 2.00000000000000006e-161Initial program 34.1%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6458.0%
Simplified58.0%
associate-/r*N/A
associate-/r*N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
cube-unmultN/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6465.9%
Applied egg-rr65.9%
if 2.00000000000000006e-161 < l Initial program 37.9%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6456.5%
Simplified56.5%
*-commutativeN/A
associate-/r*N/A
un-div-invN/A
clear-numN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r/N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
cube-unmultN/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f6456.4%
Applied egg-rr56.4%
associate-/r/N/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6467.4%
Applied egg-rr67.4%
Final simplification66.4%
(FPCore (t l k) :precision binary64 (/ (/ 2.0 (/ (* t (* k k)) l)) (/ (* k k) l)))
double code(double t, double l, double k) {
return (2.0 / ((t * (k * k)) / l)) / ((k * k) / l);
}
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 * (k * k)) / l)) / ((k * k) / l)
end function
public static double code(double t, double l, double k) {
return (2.0 / ((t * (k * k)) / l)) / ((k * k) / l);
}
def code(t, l, k): return (2.0 / ((t * (k * k)) / l)) / ((k * k) / l)
function code(t, l, k) return Float64(Float64(2.0 / Float64(Float64(t * Float64(k * k)) / l)) / Float64(Float64(k * k) / l)) end
function tmp = code(t, l, k) tmp = (2.0 / ((t * (k * k)) / l)) / ((k * k) / l); end
code[t_, l_, k_] := N[(N[(2.0 / N[(N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{2}{\frac{t \cdot \left(k \cdot k\right)}{\ell}}}{\frac{k \cdot k}{\ell}}
\end{array}
Initial program 35.4%
Taylor expanded in t around 0
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
unpow2N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6485.4%
Simplified85.4%
associate-/l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
unpow2N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6473.4%
Applied egg-rr73.4%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6473.0%
Simplified73.0%
Final simplification73.0%
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* k k) (/ (/ (* t (* k k)) l) l))))
double code(double t, double l, double k) {
return 2.0 / ((k * k) * (((t * (k * k)) / l) / l));
}
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 / ((k * k) * (((t * (k * k)) / l) / l))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((k * k) * (((t * (k * k)) / l) / l));
}
def code(t, l, k): return 2.0 / ((k * k) * (((t * (k * k)) / l) / l))
function code(t, l, k) return Float64(2.0 / Float64(Float64(k * k) * Float64(Float64(Float64(t * Float64(k * k)) / l) / l))) end
function tmp = code(t, l, k) tmp = 2.0 / ((k * k) * (((t * (k * k)) / l) / l)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(k * k), $MachinePrecision] * N[(N[(N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(k \cdot k\right) \cdot \frac{\frac{t \cdot \left(k \cdot k\right)}{\ell}}{\ell}}
\end{array}
Initial program 35.4%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.5%
Simplified57.5%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-/r*N/A
associate-*r/N/A
associate-/l*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.0%
Applied egg-rr73.0%
Final simplification73.0%
(FPCore (t l k) :precision binary64 (* (/ (/ 2.0 t) (* k k)) (/ (* l l) (* k k))))
double code(double t, double l, double k) {
return ((2.0 / t) / (k * k)) * ((l * l) / (k * k));
}
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) / (k * k)) * ((l * l) / (k * k))
end function
public static double code(double t, double l, double k) {
return ((2.0 / t) / (k * k)) * ((l * l) / (k * k));
}
def code(t, l, k): return ((2.0 / t) / (k * k)) * ((l * l) / (k * k))
function code(t, l, k) return Float64(Float64(Float64(2.0 / t) / Float64(k * k)) * Float64(Float64(l * l) / Float64(k * k))) end
function tmp = code(t, l, k) tmp = ((2.0 / t) / (k * k)) * ((l * l) / (k * k)); end
code[t_, l_, k_] := N[(N[(N[(2.0 / t), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{2}{t}}{k \cdot k} \cdot \frac{\ell \cdot \ell}{k \cdot k}
\end{array}
Initial program 35.4%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.5%
Simplified57.5%
*-commutativeN/A
associate-/r*N/A
un-div-invN/A
clear-numN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r/N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
cube-unmultN/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f6457.4%
Applied egg-rr57.4%
associate-/r/N/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.3%
Applied egg-rr63.3%
(FPCore (t l k) :precision binary64 (* (/ (/ 2.0 t) k) (/ (* l l) (* k (* k k)))))
double code(double t, double l, double k) {
return ((2.0 / t) / k) * ((l * l) / (k * (k * k)));
}
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) / k) * ((l * l) / (k * (k * k)))
end function
public static double code(double t, double l, double k) {
return ((2.0 / t) / k) * ((l * l) / (k * (k * k)));
}
def code(t, l, k): return ((2.0 / t) / k) * ((l * l) / (k * (k * k)))
function code(t, l, k) return Float64(Float64(Float64(2.0 / t) / k) * Float64(Float64(l * l) / Float64(k * Float64(k * k)))) end
function tmp = code(t, l, k) tmp = ((2.0 / t) / k) * ((l * l) / (k * (k * k))); end
code[t_, l_, k_] := N[(N[(N[(2.0 / t), $MachinePrecision] / k), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] / N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{2}{t}}{k} \cdot \frac{\ell \cdot \ell}{k \cdot \left(k \cdot k\right)}
\end{array}
Initial program 35.4%
Taylor expanded in k around 0
associate-/l*N/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.5%
Simplified57.5%
*-commutativeN/A
associate-/r*N/A
un-div-invN/A
clear-numN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r/N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
cube-unmultN/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f6457.4%
Applied egg-rr57.4%
associate-/r/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6460.6%
Applied egg-rr60.6%
(FPCore (t l k) :precision binary64 (/ (- 0.0 0.11666666666666667) (/ t (* l l))))
double code(double t, double l, double k) {
return (0.0 - 0.11666666666666667) / (t / (l * l));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = (0.0d0 - 0.11666666666666667d0) / (t / (l * l))
end function
public static double code(double t, double l, double k) {
return (0.0 - 0.11666666666666667) / (t / (l * l));
}
def code(t, l, k): return (0.0 - 0.11666666666666667) / (t / (l * l))
function code(t, l, k) return Float64(Float64(0.0 - 0.11666666666666667) / Float64(t / Float64(l * l))) end
function tmp = code(t, l, k) tmp = (0.0 - 0.11666666666666667) / (t / (l * l)); end
code[t_, l_, k_] := N[(N[(0.0 - 0.11666666666666667), $MachinePrecision] / N[(t / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0 - 0.11666666666666667}{\frac{t}{\ell \cdot \ell}}
\end{array}
Initial program 35.4%
Taylor expanded in k around 0
Simplified34.5%
Taylor expanded in k around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6420.5%
Simplified20.5%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6420.5%
Applied egg-rr20.5%
clear-numN/A
un-div-invN/A
frac-2negN/A
/-lowering-/.f64N/A
metadata-evalN/A
neg-sub0N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6428.8%
Applied egg-rr28.8%
Final simplification20.5%
(FPCore (t l k) :precision binary64 (* -0.11666666666666667 (/ (* l l) t)))
double code(double t, double l, double k) {
return -0.11666666666666667 * ((l * l) / t);
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = (-0.11666666666666667d0) * ((l * l) / t)
end function
public static double code(double t, double l, double k) {
return -0.11666666666666667 * ((l * l) / t);
}
def code(t, l, k): return -0.11666666666666667 * ((l * l) / t)
function code(t, l, k) return Float64(-0.11666666666666667 * Float64(Float64(l * l) / t)) end
function tmp = code(t, l, k) tmp = -0.11666666666666667 * ((l * l) / t); end
code[t_, l_, k_] := N[(-0.11666666666666667 * N[(N[(l * l), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.11666666666666667 \cdot \frac{\ell \cdot \ell}{t}
\end{array}
Initial program 35.4%
Taylor expanded in k around 0
Simplified34.5%
Taylor expanded in k around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6420.5%
Simplified20.5%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6420.5%
Applied egg-rr20.5%
(FPCore (t l k) :precision binary64 (* -0.11666666666666667 (/ l (/ t l))))
double code(double t, double l, double k) {
return -0.11666666666666667 * (l / (t / l));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = (-0.11666666666666667d0) * (l / (t / l))
end function
public static double code(double t, double l, double k) {
return -0.11666666666666667 * (l / (t / l));
}
def code(t, l, k): return -0.11666666666666667 * (l / (t / l))
function code(t, l, k) return Float64(-0.11666666666666667 * Float64(l / Float64(t / l))) end
function tmp = code(t, l, k) tmp = -0.11666666666666667 * (l / (t / l)); end
code[t_, l_, k_] := N[(-0.11666666666666667 * N[(l / N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.11666666666666667 \cdot \frac{\ell}{\frac{t}{\ell}}
\end{array}
Initial program 35.4%
Taylor expanded in k around 0
Simplified34.5%
Taylor expanded in k around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6420.5%
Simplified20.5%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6420.5%
Applied egg-rr20.5%
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6417.2%
Applied egg-rr17.2%
herbie shell --seed 2024139
(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))))