
(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 (/ (* (sin k) (* t (tan k))) (/ l k))) (/ k l)))
double code(double t, double l, double k) {
return (2.0 / ((sin(k) * (t * tan(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 / ((sin(k) * (t * tan(k))) / (l / k))) / (k / l)
end function
public static double code(double t, double l, double k) {
return (2.0 / ((Math.sin(k) * (t * Math.tan(k))) / (l / k))) / (k / l);
}
def code(t, l, k): return (2.0 / ((math.sin(k) * (t * math.tan(k))) / (l / k))) / (k / l)
function code(t, l, k) return Float64(Float64(2.0 / Float64(Float64(sin(k) * Float64(t * tan(k))) / Float64(l / k))) / Float64(k / l)) end
function tmp = code(t, l, k) tmp = (2.0 / ((sin(k) * (t * tan(k))) / (l / k))) / (k / l); end
code[t_, l_, k_] := N[(N[(2.0 / N[(N[(N[Sin[k], $MachinePrecision] * N[(t * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k / l), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{2}{\frac{\sin k \cdot \left(t \cdot \tan k\right)}{\frac{\ell}{k}}}}{\frac{k}{\ell}}
\end{array}
Initial program 34.5%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6482.6%
Simplified82.6%
associate-/l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6490.4%
Applied egg-rr90.4%
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6491.8%
Applied egg-rr91.8%
associate-*l*N/A
associate-/r*N/A
associate-*l*N/A
associate-*r*N/A
associate-/r*N/A
sqr-sin-aN/A
associate-/r*N/A
div-invN/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr97.5%
(FPCore (t l k) :precision binary64 (* 2.0 (/ (/ l k) (/ (* (sin k) (* t (tan k))) (/ l k)))))
double code(double t, double l, double k) {
return 2.0 * ((l / k) / ((sin(k) * (t * tan(k))) / (l / 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 / k) / ((sin(k) * (t * tan(k))) / (l / k)))
end function
public static double code(double t, double l, double k) {
return 2.0 * ((l / k) / ((Math.sin(k) * (t * Math.tan(k))) / (l / k)));
}
def code(t, l, k): return 2.0 * ((l / k) / ((math.sin(k) * (t * math.tan(k))) / (l / k)))
function code(t, l, k) return Float64(2.0 * Float64(Float64(l / k) / Float64(Float64(sin(k) * Float64(t * tan(k))) / Float64(l / k)))) end
function tmp = code(t, l, k) tmp = 2.0 * ((l / k) / ((sin(k) * (t * tan(k))) / (l / k))); end
code[t_, l_, k_] := N[(2.0 * N[(N[(l / k), $MachinePrecision] / N[(N[(N[Sin[k], $MachinePrecision] * N[(t * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \frac{\frac{\ell}{k}}{\frac{\sin k \cdot \left(t \cdot \tan k\right)}{\frac{\ell}{k}}}
\end{array}
Initial program 34.5%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6482.6%
Simplified82.6%
associate-/l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6490.4%
Applied egg-rr90.4%
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6491.8%
Applied egg-rr91.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr97.4%
Final simplification97.4%
(FPCore (t l k) :precision binary64 (* (/ 2.0 k) (/ l (/ (* (sin k) (* t (tan k))) (/ l k)))))
double code(double t, double l, double k) {
return (2.0 / k) * (l / ((sin(k) * (t * tan(k))) / (l / 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) * (l / ((sin(k) * (t * tan(k))) / (l / k)))
end function
public static double code(double t, double l, double k) {
return (2.0 / k) * (l / ((Math.sin(k) * (t * Math.tan(k))) / (l / k)));
}
def code(t, l, k): return (2.0 / k) * (l / ((math.sin(k) * (t * math.tan(k))) / (l / k)))
function code(t, l, k) return Float64(Float64(2.0 / k) * Float64(l / Float64(Float64(sin(k) * Float64(t * tan(k))) / Float64(l / k)))) end
function tmp = code(t, l, k) tmp = (2.0 / k) * (l / ((sin(k) * (t * tan(k))) / (l / k))); end
code[t_, l_, k_] := N[(N[(2.0 / k), $MachinePrecision] * N[(l / N[(N[(N[Sin[k], $MachinePrecision] * N[(t * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{k} \cdot \frac{\ell}{\frac{\sin k \cdot \left(t \cdot \tan k\right)}{\frac{\ell}{k}}}
\end{array}
Initial program 34.5%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6482.6%
Simplified82.6%
associate-/l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6490.4%
Applied egg-rr90.4%
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6491.8%
Applied egg-rr91.8%
associate-*l*N/A
associate-/r*N/A
associate-*l*N/A
associate-*r*N/A
associate-/r*N/A
sqr-sin-aN/A
associate-/r*N/A
associate-/r/N/A
associate-/l*N/A
un-div-invN/A
*-lowering-*.f64N/A
Applied egg-rr95.5%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ 2.0 (/ k l))))
(if (<= k 1.6e-105)
(/
(/ (/ t_1 (/ k l)) k)
(* k (+ t (* t (* (* k k) 0.16666666666666666)))))
(if (<= k 9.5e-7)
(/
t_1
(*
(* k (* k k))
(+ (/ t l) (* k (* k (/ (* t 0.16666666666666666) l))))))
(/ t_1 (* (/ (* k t) l) (- 0.5 (* 0.5 (cos (* 2.0 k))))))))))
double code(double t, double l, double k) {
double t_1 = 2.0 / (k / l);
double tmp;
if (k <= 1.6e-105) {
tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666))));
} else if (k <= 9.5e-7) {
tmp = t_1 / ((k * (k * k)) * ((t / l) + (k * (k * ((t * 0.16666666666666666) / l)))));
} else {
tmp = t_1 / (((k * t) / l) * (0.5 - (0.5 * cos((2.0 * 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) :: t_1
real(8) :: tmp
t_1 = 2.0d0 / (k / l)
if (k <= 1.6d-105) then
tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666d0))))
else if (k <= 9.5d-7) then
tmp = t_1 / ((k * (k * k)) * ((t / l) + (k * (k * ((t * 0.16666666666666666d0) / l)))))
else
tmp = t_1 / (((k * t) / l) * (0.5d0 - (0.5d0 * cos((2.0d0 * k)))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = 2.0 / (k / l);
double tmp;
if (k <= 1.6e-105) {
tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666))));
} else if (k <= 9.5e-7) {
tmp = t_1 / ((k * (k * k)) * ((t / l) + (k * (k * ((t * 0.16666666666666666) / l)))));
} else {
tmp = t_1 / (((k * t) / l) * (0.5 - (0.5 * Math.cos((2.0 * k)))));
}
return tmp;
}
def code(t, l, k): t_1 = 2.0 / (k / l) tmp = 0 if k <= 1.6e-105: tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666)))) elif k <= 9.5e-7: tmp = t_1 / ((k * (k * k)) * ((t / l) + (k * (k * ((t * 0.16666666666666666) / l))))) else: tmp = t_1 / (((k * t) / l) * (0.5 - (0.5 * math.cos((2.0 * k))))) return tmp
function code(t, l, k) t_1 = Float64(2.0 / Float64(k / l)) tmp = 0.0 if (k <= 1.6e-105) tmp = Float64(Float64(Float64(t_1 / Float64(k / l)) / k) / Float64(k * Float64(t + Float64(t * Float64(Float64(k * k) * 0.16666666666666666))))); elseif (k <= 9.5e-7) tmp = Float64(t_1 / Float64(Float64(k * Float64(k * k)) * Float64(Float64(t / l) + Float64(k * Float64(k * Float64(Float64(t * 0.16666666666666666) / l)))))); else tmp = Float64(t_1 / Float64(Float64(Float64(k * t) / l) * Float64(0.5 - Float64(0.5 * cos(Float64(2.0 * k)))))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = 2.0 / (k / l); tmp = 0.0; if (k <= 1.6e-105) tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666)))); elseif (k <= 9.5e-7) tmp = t_1 / ((k * (k * k)) * ((t / l) + (k * (k * ((t * 0.16666666666666666) / l))))); else tmp = t_1 / (((k * t) / l) * (0.5 - (0.5 * cos((2.0 * k))))); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(2.0 / N[(k / l), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, 1.6e-105], N[(N[(N[(t$95$1 / N[(k / l), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision] / N[(k * N[(t + N[(t * N[(N[(k * k), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 9.5e-7], N[(t$95$1 / N[(N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] + N[(k * N[(k * N[(N[(t * 0.16666666666666666), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 / N[(N[(N[(k * t), $MachinePrecision] / l), $MachinePrecision] * N[(0.5 - N[(0.5 * N[Cos[N[(2.0 * k), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{\frac{k}{\ell}}\\
\mathbf{if}\;k \leq 1.6 \cdot 10^{-105}:\\
\;\;\;\;\frac{\frac{\frac{t\_1}{\frac{k}{\ell}}}{k}}{k \cdot \left(t + t \cdot \left(\left(k \cdot k\right) \cdot 0.16666666666666666\right)\right)}\\
\mathbf{elif}\;k \leq 9.5 \cdot 10^{-7}:\\
\;\;\;\;\frac{t\_1}{\left(k \cdot \left(k \cdot k\right)\right) \cdot \left(\frac{t}{\ell} + k \cdot \left(k \cdot \frac{t \cdot 0.16666666666666666}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_1}{\frac{k \cdot t}{\ell} \cdot \left(0.5 - 0.5 \cdot \cos \left(2 \cdot k\right)\right)}\\
\end{array}
\end{array}
if k < 1.59999999999999991e-105Initial program 40.1%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6485.6%
Simplified85.6%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6477.1%
Simplified77.1%
associate-/r*N/A
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-/l/N/A
frac-timesN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6480.4%
Applied egg-rr80.4%
if 1.59999999999999991e-105 < k < 9.5000000000000001e-7Initial program 24.4%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6486.4%
Simplified86.4%
associate-/l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6486.4%
Applied egg-rr86.4%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f6459.8%
Applied egg-rr59.8%
Taylor expanded in k around 0
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*l/N/A
metadata-evalN/A
distribute-rgt-out--N/A
/-lowering-/.f64N/A
distribute-rgt-out--N/A
metadata-evalN/A
*-lowering-*.f6499.8%
Simplified99.8%
if 9.5000000000000001e-7 < k Initial program 23.1%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6471.6%
Simplified71.6%
associate-/l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6493.5%
Applied egg-rr93.5%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f6499.5%
Applied egg-rr99.5%
Taylor expanded in k around 0
/-lowering-/.f64N/A
*-lowering-*.f6462.5%
Simplified62.5%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ 2.0 (/ k l))))
(if (<= t 1.7e+43)
(/
t_1
(*
k
(*
(* k k)
(+
(/ t l)
(*
(* k k)
(+
(/ (* t 0.16666666666666666) l)
(*
(* k k)
(+
(-
(* (/ t l) 0.044444444444444446)
(* (/ t l) -0.08333333333333333))
(* (/ t l) -0.041666666666666664)))))))))
(/
(/ (/ t_1 (/ k l)) k)
(* k (+ t (* t (* (* k k) 0.16666666666666666))))))))
double code(double t, double l, double k) {
double t_1 = 2.0 / (k / l);
double tmp;
if (t <= 1.7e+43) {
tmp = t_1 / (k * ((k * k) * ((t / l) + ((k * k) * (((t * 0.16666666666666666) / l) + ((k * k) * ((((t / l) * 0.044444444444444446) - ((t / l) * -0.08333333333333333)) + ((t / l) * -0.041666666666666664))))))));
} else {
tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666))));
}
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 = 2.0d0 / (k / l)
if (t <= 1.7d+43) then
tmp = t_1 / (k * ((k * k) * ((t / l) + ((k * k) * (((t * 0.16666666666666666d0) / l) + ((k * k) * ((((t / l) * 0.044444444444444446d0) - ((t / l) * (-0.08333333333333333d0))) + ((t / l) * (-0.041666666666666664d0)))))))))
else
tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666d0))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = 2.0 / (k / l);
double tmp;
if (t <= 1.7e+43) {
tmp = t_1 / (k * ((k * k) * ((t / l) + ((k * k) * (((t * 0.16666666666666666) / l) + ((k * k) * ((((t / l) * 0.044444444444444446) - ((t / l) * -0.08333333333333333)) + ((t / l) * -0.041666666666666664))))))));
} else {
tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666))));
}
return tmp;
}
def code(t, l, k): t_1 = 2.0 / (k / l) tmp = 0 if t <= 1.7e+43: tmp = t_1 / (k * ((k * k) * ((t / l) + ((k * k) * (((t * 0.16666666666666666) / l) + ((k * k) * ((((t / l) * 0.044444444444444446) - ((t / l) * -0.08333333333333333)) + ((t / l) * -0.041666666666666664)))))))) else: tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666)))) return tmp
function code(t, l, k) t_1 = Float64(2.0 / Float64(k / l)) tmp = 0.0 if (t <= 1.7e+43) tmp = Float64(t_1 / Float64(k * Float64(Float64(k * k) * Float64(Float64(t / l) + Float64(Float64(k * k) * Float64(Float64(Float64(t * 0.16666666666666666) / l) + Float64(Float64(k * k) * Float64(Float64(Float64(Float64(t / l) * 0.044444444444444446) - Float64(Float64(t / l) * -0.08333333333333333)) + Float64(Float64(t / l) * -0.041666666666666664))))))))); else tmp = Float64(Float64(Float64(t_1 / Float64(k / l)) / k) / Float64(k * Float64(t + Float64(t * Float64(Float64(k * k) * 0.16666666666666666))))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = 2.0 / (k / l); tmp = 0.0; if (t <= 1.7e+43) tmp = t_1 / (k * ((k * k) * ((t / l) + ((k * k) * (((t * 0.16666666666666666) / l) + ((k * k) * ((((t / l) * 0.044444444444444446) - ((t / l) * -0.08333333333333333)) + ((t / l) * -0.041666666666666664)))))))); else tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666)))); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(2.0 / N[(k / l), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, 1.7e+43], N[(t$95$1 / N[(k * N[(N[(k * k), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] + N[(N[(k * k), $MachinePrecision] * N[(N[(N[(t * 0.16666666666666666), $MachinePrecision] / l), $MachinePrecision] + N[(N[(k * k), $MachinePrecision] * N[(N[(N[(N[(t / l), $MachinePrecision] * 0.044444444444444446), $MachinePrecision] - N[(N[(t / l), $MachinePrecision] * -0.08333333333333333), $MachinePrecision]), $MachinePrecision] + N[(N[(t / l), $MachinePrecision] * -0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t$95$1 / N[(k / l), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision] / N[(k * N[(t + N[(t * N[(N[(k * k), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{\frac{k}{\ell}}\\
\mathbf{if}\;t \leq 1.7 \cdot 10^{+43}:\\
\;\;\;\;\frac{t\_1}{k \cdot \left(\left(k \cdot k\right) \cdot \left(\frac{t}{\ell} + \left(k \cdot k\right) \cdot \left(\frac{t \cdot 0.16666666666666666}{\ell} + \left(k \cdot k\right) \cdot \left(\left(\frac{t}{\ell} \cdot 0.044444444444444446 - \frac{t}{\ell} \cdot -0.08333333333333333\right) + \frac{t}{\ell} \cdot -0.041666666666666664\right)\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{t\_1}{\frac{k}{\ell}}}{k}}{k \cdot \left(t + t \cdot \left(\left(k \cdot k\right) \cdot 0.16666666666666666\right)\right)}\\
\end{array}
\end{array}
if t < 1.70000000000000006e43Initial program 38.0%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6480.6%
Simplified80.6%
associate-/l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6490.0%
Applied egg-rr90.0%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f6487.6%
Applied egg-rr87.6%
Taylor expanded in k around 0
Simplified67.5%
if 1.70000000000000006e43 < t Initial program 21.2%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6490.1%
Simplified90.1%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6475.8%
Simplified75.8%
associate-/r*N/A
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-/l/N/A
frac-timesN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6479.6%
Applied egg-rr79.6%
Final simplification70.0%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ 2.0 (/ k l))))
(if (<= k 1.6e-105)
(/
(/ (/ t_1 (/ k l)) k)
(* k (+ t (* t (* (* k k) 0.16666666666666666)))))
(/
t_1
(*
(* k (* k k))
(+ (/ t l) (* k (* k (/ (* t 0.16666666666666666) l)))))))))
double code(double t, double l, double k) {
double t_1 = 2.0 / (k / l);
double tmp;
if (k <= 1.6e-105) {
tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666))));
} else {
tmp = t_1 / ((k * (k * k)) * ((t / l) + (k * (k * ((t * 0.16666666666666666) / 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 = 2.0d0 / (k / l)
if (k <= 1.6d-105) then
tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666d0))))
else
tmp = t_1 / ((k * (k * k)) * ((t / l) + (k * (k * ((t * 0.16666666666666666d0) / l)))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = 2.0 / (k / l);
double tmp;
if (k <= 1.6e-105) {
tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666))));
} else {
tmp = t_1 / ((k * (k * k)) * ((t / l) + (k * (k * ((t * 0.16666666666666666) / l)))));
}
return tmp;
}
def code(t, l, k): t_1 = 2.0 / (k / l) tmp = 0 if k <= 1.6e-105: tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666)))) else: tmp = t_1 / ((k * (k * k)) * ((t / l) + (k * (k * ((t * 0.16666666666666666) / l))))) return tmp
function code(t, l, k) t_1 = Float64(2.0 / Float64(k / l)) tmp = 0.0 if (k <= 1.6e-105) tmp = Float64(Float64(Float64(t_1 / Float64(k / l)) / k) / Float64(k * Float64(t + Float64(t * Float64(Float64(k * k) * 0.16666666666666666))))); else tmp = Float64(t_1 / Float64(Float64(k * Float64(k * k)) * Float64(Float64(t / l) + Float64(k * Float64(k * Float64(Float64(t * 0.16666666666666666) / l)))))); end return tmp end
function tmp_2 = code(t, l, k) t_1 = 2.0 / (k / l); tmp = 0.0; if (k <= 1.6e-105) tmp = ((t_1 / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666)))); else tmp = t_1 / ((k * (k * k)) * ((t / l) + (k * (k * ((t * 0.16666666666666666) / l))))); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(2.0 / N[(k / l), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, 1.6e-105], N[(N[(N[(t$95$1 / N[(k / l), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision] / N[(k * N[(t + N[(t * N[(N[(k * k), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 / N[(N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] + N[(k * N[(k * N[(N[(t * 0.16666666666666666), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{\frac{k}{\ell}}\\
\mathbf{if}\;k \leq 1.6 \cdot 10^{-105}:\\
\;\;\;\;\frac{\frac{\frac{t\_1}{\frac{k}{\ell}}}{k}}{k \cdot \left(t + t \cdot \left(\left(k \cdot k\right) \cdot 0.16666666666666666\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_1}{\left(k \cdot \left(k \cdot k\right)\right) \cdot \left(\frac{t}{\ell} + k \cdot \left(k \cdot \frac{t \cdot 0.16666666666666666}{\ell}\right)\right)}\\
\end{array}
\end{array}
if k < 1.59999999999999991e-105Initial program 40.1%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6485.6%
Simplified85.6%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6477.1%
Simplified77.1%
associate-/r*N/A
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-/l/N/A
frac-timesN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6480.4%
Applied egg-rr80.4%
if 1.59999999999999991e-105 < k Initial program 23.6%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6476.6%
Simplified76.6%
associate-/l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6491.1%
Applied egg-rr91.1%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f6486.1%
Applied egg-rr86.1%
Taylor expanded in k around 0
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*l/N/A
metadata-evalN/A
distribute-rgt-out--N/A
/-lowering-/.f64N/A
distribute-rgt-out--N/A
metadata-evalN/A
*-lowering-*.f6468.6%
Simplified68.6%
(FPCore (t l k)
:precision binary64
(if (<= (* l l) 1e+248)
(/
(/ (/ (/ 2.0 (/ k l)) (/ k l)) k)
(* k (+ t (* t (* (* k k) 0.16666666666666666)))))
(/ (/ 2.0 (* k k)) (/ (* k k) (/ (* l l) t)))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e+248) {
tmp = (((2.0 / (k / l)) / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666))));
} else {
tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t));
}
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 * l) <= 1d+248) then
tmp = (((2.0d0 / (k / l)) / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666d0))))
else
tmp = (2.0d0 / (k * k)) / ((k * k) / ((l * l) / t))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e+248) {
tmp = (((2.0 / (k / l)) / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666))));
} else {
tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (l * l) <= 1e+248: tmp = (((2.0 / (k / l)) / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666)))) else: tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t)) return tmp
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e+248) tmp = Float64(Float64(Float64(Float64(2.0 / Float64(k / l)) / Float64(k / l)) / k) / Float64(k * Float64(t + Float64(t * Float64(Float64(k * k) * 0.16666666666666666))))); else tmp = Float64(Float64(2.0 / Float64(k * k)) / Float64(Float64(k * k) / Float64(Float64(l * l) / t))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((l * l) <= 1e+248) tmp = (((2.0 / (k / l)) / (k / l)) / k) / (k * (t + (t * ((k * k) * 0.16666666666666666)))); else tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e+248], N[(N[(N[(N[(2.0 / N[(k / l), $MachinePrecision]), $MachinePrecision] / N[(k / l), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision] / N[(k * N[(t + N[(t * N[(N[(k * k), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] / N[(N[(l * l), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{+248}:\\
\;\;\;\;\frac{\frac{\frac{\frac{2}{\frac{k}{\ell}}}{\frac{k}{\ell}}}{k}}{k \cdot \left(t + t \cdot \left(\left(k \cdot k\right) \cdot 0.16666666666666666\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{k \cdot k}}{\frac{k \cdot k}{\frac{\ell \cdot \ell}{t}}}\\
\end{array}
\end{array}
if (*.f64 l l) < 1.00000000000000005e248Initial program 30.7%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6485.0%
Simplified85.0%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6476.6%
Simplified76.6%
associate-/r*N/A
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-/l/N/A
frac-timesN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6480.4%
Applied egg-rr80.4%
if 1.00000000000000005e248 < (*.f64 l l) Initial program 43.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-*.f6464.7%
Simplified64.7%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6469.8%
Applied egg-rr69.8%
(FPCore (t l k)
:precision binary64
(if (<= (* l l) 1e-123)
(/
(/ 2.0 (/ (* k k) (/ l (/ k (/ l k)))))
(+ t (* t (* (* k k) 0.16666666666666666))))
(/ (/ 2.0 (* k k)) (/ (* k k) (/ (* l l) t)))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-123) {
tmp = (2.0 / ((k * k) / (l / (k / (l / k))))) / (t + (t * ((k * k) * 0.16666666666666666)));
} else {
tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t));
}
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 * l) <= 1d-123) then
tmp = (2.0d0 / ((k * k) / (l / (k / (l / k))))) / (t + (t * ((k * k) * 0.16666666666666666d0)))
else
tmp = (2.0d0 / (k * k)) / ((k * k) / ((l * l) / t))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-123) {
tmp = (2.0 / ((k * k) / (l / (k / (l / k))))) / (t + (t * ((k * k) * 0.16666666666666666)));
} else {
tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (l * l) <= 1e-123: tmp = (2.0 / ((k * k) / (l / (k / (l / k))))) / (t + (t * ((k * k) * 0.16666666666666666))) else: tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t)) return tmp
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-123) tmp = Float64(Float64(2.0 / Float64(Float64(k * k) / Float64(l / Float64(k / Float64(l / k))))) / Float64(t + Float64(t * Float64(Float64(k * k) * 0.16666666666666666)))); else tmp = Float64(Float64(2.0 / Float64(k * k)) / Float64(Float64(k * k) / Float64(Float64(l * l) / t))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((l * l) <= 1e-123) tmp = (2.0 / ((k * k) / (l / (k / (l / k))))) / (t + (t * ((k * k) * 0.16666666666666666))); else tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e-123], N[(N[(2.0 / N[(N[(k * k), $MachinePrecision] / N[(l / N[(k / N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t + N[(t * N[(N[(k * k), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] / N[(N[(l * l), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-123}:\\
\;\;\;\;\frac{\frac{2}{\frac{k \cdot k}{\frac{\ell}{\frac{k}{\frac{\ell}{k}}}}}}{t + t \cdot \left(\left(k \cdot k\right) \cdot 0.16666666666666666\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{k \cdot k}}{\frac{k \cdot k}{\frac{\ell \cdot \ell}{t}}}\\
\end{array}
\end{array}
if (*.f64 l l) < 1.0000000000000001e-123Initial program 24.9%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6483.0%
Simplified83.0%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6481.0%
Simplified81.0%
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Applied egg-rr86.0%
if 1.0000000000000001e-123 < (*.f64 l l) Initial program 40.3%
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-*.f6467.5%
Simplified67.5%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6470.6%
Applied egg-rr70.6%
(FPCore (t l k)
:precision binary64
(if (<= t 6.4e-86)
(/ 2.0 (* (/ t l) (/ (* k (* k (* k k))) l)))
(/
2.0
(*
(* k (+ t (* t (* (* k k) 0.16666666666666666))))
(* k (/ k (/ l (/ k l))))))))
double code(double t, double l, double k) {
double tmp;
if (t <= 6.4e-86) {
tmp = 2.0 / ((t / l) * ((k * (k * (k * k))) / l));
} else {
tmp = 2.0 / ((k * (t + (t * ((k * k) * 0.16666666666666666)))) * (k * (k / (l / (k / 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) :: tmp
if (t <= 6.4d-86) then
tmp = 2.0d0 / ((t / l) * ((k * (k * (k * k))) / l))
else
tmp = 2.0d0 / ((k * (t + (t * ((k * k) * 0.16666666666666666d0)))) * (k * (k / (l / (k / l)))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (t <= 6.4e-86) {
tmp = 2.0 / ((t / l) * ((k * (k * (k * k))) / l));
} else {
tmp = 2.0 / ((k * (t + (t * ((k * k) * 0.16666666666666666)))) * (k * (k / (l / (k / l)))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if t <= 6.4e-86: tmp = 2.0 / ((t / l) * ((k * (k * (k * k))) / l)) else: tmp = 2.0 / ((k * (t + (t * ((k * k) * 0.16666666666666666)))) * (k * (k / (l / (k / l))))) return tmp
function code(t, l, k) tmp = 0.0 if (t <= 6.4e-86) tmp = Float64(2.0 / Float64(Float64(t / l) * Float64(Float64(k * Float64(k * Float64(k * k))) / l))); else tmp = Float64(2.0 / Float64(Float64(k * Float64(t + Float64(t * Float64(Float64(k * k) * 0.16666666666666666)))) * Float64(k * Float64(k / Float64(l / Float64(k / l)))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (t <= 6.4e-86) tmp = 2.0 / ((t / l) * ((k * (k * (k * k))) / l)); else tmp = 2.0 / ((k * (t + (t * ((k * k) * 0.16666666666666666)))) * (k * (k / (l / (k / l))))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[t, 6.4e-86], N[(2.0 / N[(N[(t / l), $MachinePrecision] * N[(N[(k * N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k * N[(t + N[(t * N[(N[(k * k), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(k * N[(k / N[(l / N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 6.4 \cdot 10^{-86}:\\
\;\;\;\;\frac{2}{\frac{t}{\ell} \cdot \frac{k \cdot \left(k \cdot \left(k \cdot k\right)\right)}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(k \cdot \left(t + t \cdot \left(\left(k \cdot k\right) \cdot 0.16666666666666666\right)\right)\right) \cdot \left(k \cdot \frac{k}{\frac{\ell}{\frac{k}{\ell}}}\right)}\\
\end{array}
\end{array}
if t < 6.40000000000000011e-86Initial program 34.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-*.f6465.5%
Simplified65.5%
associate-*r/N/A
times-fracN/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-/.f6472.5%
Applied egg-rr72.5%
if 6.40000000000000011e-86 < t Initial program 33.8%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6482.2%
Simplified82.2%
Taylor expanded in k around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.9%
Simplified67.9%
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-/l/N/A
frac-timesN/A
clear-numN/A
frac-timesN/A
*-rgt-identityN/A
/-lowering-/.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6469.8%
Applied egg-rr69.8%
Final simplification71.7%
(FPCore (t l k) :precision binary64 (if (<= (* l l) 1e-139) (/ (/ 2.0 (/ k l)) (/ (* k (* t (* k k))) l)) (/ (/ 2.0 (* k k)) (/ (* k k) (/ (* l l) t)))))
double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-139) {
tmp = (2.0 / (k / l)) / ((k * (t * (k * k))) / l);
} else {
tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t));
}
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 * l) <= 1d-139) then
tmp = (2.0d0 / (k / l)) / ((k * (t * (k * k))) / l)
else
tmp = (2.0d0 / (k * k)) / ((k * k) / ((l * l) / t))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if ((l * l) <= 1e-139) {
tmp = (2.0 / (k / l)) / ((k * (t * (k * k))) / l);
} else {
tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t));
}
return tmp;
}
def code(t, l, k): tmp = 0 if (l * l) <= 1e-139: tmp = (2.0 / (k / l)) / ((k * (t * (k * k))) / l) else: tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t)) return tmp
function code(t, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-139) tmp = Float64(Float64(2.0 / Float64(k / l)) / Float64(Float64(k * Float64(t * Float64(k * k))) / l)); else tmp = Float64(Float64(2.0 / Float64(k * k)) / Float64(Float64(k * k) / Float64(Float64(l * l) / t))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if ((l * l) <= 1e-139) tmp = (2.0 / (k / l)) / ((k * (t * (k * k))) / l); else tmp = (2.0 / (k * k)) / ((k * k) / ((l * l) / t)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[N[(l * l), $MachinePrecision], 1e-139], N[(N[(2.0 / N[(k / l), $MachinePrecision]), $MachinePrecision] / N[(N[(k * N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] / N[(N[(l * l), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-139}:\\
\;\;\;\;\frac{\frac{2}{\frac{k}{\ell}}}{\frac{k \cdot \left(t \cdot \left(k \cdot k\right)\right)}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{k \cdot k}}{\frac{k \cdot k}{\frac{\ell \cdot \ell}{t}}}\\
\end{array}
\end{array}
if (*.f64 l l) < 1.00000000000000003e-139Initial program 24.2%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6481.9%
Simplified81.9%
associate-/l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6481.8%
Applied egg-rr81.8%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f6473.0%
Applied egg-rr73.0%
Taylor expanded in k around 0
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6481.6%
Simplified81.6%
if 1.00000000000000003e-139 < (*.f64 l l) Initial program 40.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-*.f6468.1%
Simplified68.1%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6471.6%
Applied egg-rr71.6%
(FPCore (t l k) :precision binary64 (if (<= t 11600.0) (* l (/ (* (/ 2.0 k) (/ l t)) (* k (* k k)))) (/ 2.0 (* (* t (* k k)) (/ (/ (* k k) l) l)))))
double code(double t, double l, double k) {
double tmp;
if (t <= 11600.0) {
tmp = l * (((2.0 / k) * (l / t)) / (k * (k * k)));
} else {
tmp = 2.0 / ((t * (k * k)) * (((k * k) / l) / 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) :: tmp
if (t <= 11600.0d0) then
tmp = l * (((2.0d0 / k) * (l / t)) / (k * (k * k)))
else
tmp = 2.0d0 / ((t * (k * k)) * (((k * k) / l) / l))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (t <= 11600.0) {
tmp = l * (((2.0 / k) * (l / t)) / (k * (k * k)));
} else {
tmp = 2.0 / ((t * (k * k)) * (((k * k) / l) / l));
}
return tmp;
}
def code(t, l, k): tmp = 0 if t <= 11600.0: tmp = l * (((2.0 / k) * (l / t)) / (k * (k * k))) else: tmp = 2.0 / ((t * (k * k)) * (((k * k) / l) / l)) return tmp
function code(t, l, k) tmp = 0.0 if (t <= 11600.0) tmp = Float64(l * Float64(Float64(Float64(2.0 / k) * Float64(l / t)) / Float64(k * Float64(k * k)))); else tmp = Float64(2.0 / Float64(Float64(t * Float64(k * k)) * Float64(Float64(Float64(k * k) / l) / l))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (t <= 11600.0) tmp = l * (((2.0 / k) * (l / t)) / (k * (k * k))); else tmp = 2.0 / ((t * (k * k)) * (((k * k) / l) / l)); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[t, 11600.0], N[(l * N[(N[(N[(2.0 / k), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision] / N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t * N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(k * k), $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 11600:\\
\;\;\;\;\ell \cdot \frac{\frac{2}{k} \cdot \frac{\ell}{t}}{k \cdot \left(k \cdot k\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(t \cdot \left(k \cdot k\right)\right) \cdot \frac{\frac{k \cdot k}{\ell}}{\ell}}\\
\end{array}
\end{array}
if t < 11600Initial program 36.8%
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-*.f6462.7%
Simplified62.7%
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-/.f6469.0%
Applied egg-rr69.0%
div-invN/A
associate-/r*N/A
clear-numN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6470.5%
Applied egg-rr70.5%
if 11600 < t Initial program 27.4%
Taylor expanded in t around 0
times-fracN/A
*-lowering-*.f64N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6488.6%
Simplified88.6%
Taylor expanded in k around 0
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6476.8%
Simplified76.8%
Final simplification72.1%
(FPCore (t l k) :precision binary64 (if (<= t 1e-6) (/ 1.0 (/ t (* l (* l -0.11666666666666667)))) (* t (* -0.11666666666666667 (* l (/ l (* t t)))))))
double code(double t, double l, double k) {
double tmp;
if (t <= 1e-6) {
tmp = 1.0 / (t / (l * (l * -0.11666666666666667)));
} else {
tmp = t * (-0.11666666666666667 * (l * (l / (t * t))));
}
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 (t <= 1d-6) then
tmp = 1.0d0 / (t / (l * (l * (-0.11666666666666667d0))))
else
tmp = t * ((-0.11666666666666667d0) * (l * (l / (t * t))))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (t <= 1e-6) {
tmp = 1.0 / (t / (l * (l * -0.11666666666666667)));
} else {
tmp = t * (-0.11666666666666667 * (l * (l / (t * t))));
}
return tmp;
}
def code(t, l, k): tmp = 0 if t <= 1e-6: tmp = 1.0 / (t / (l * (l * -0.11666666666666667))) else: tmp = t * (-0.11666666666666667 * (l * (l / (t * t)))) return tmp
function code(t, l, k) tmp = 0.0 if (t <= 1e-6) tmp = Float64(1.0 / Float64(t / Float64(l * Float64(l * -0.11666666666666667)))); else tmp = Float64(t * Float64(-0.11666666666666667 * Float64(l * Float64(l / Float64(t * t))))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (t <= 1e-6) tmp = 1.0 / (t / (l * (l * -0.11666666666666667))); else tmp = t * (-0.11666666666666667 * (l * (l / (t * t)))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[t, 1e-6], N[(1.0 / N[(t / N[(l * N[(l * -0.11666666666666667), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t * N[(-0.11666666666666667 * N[(l * N[(l / N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 10^{-6}:\\
\;\;\;\;\frac{1}{\frac{t}{\ell \cdot \left(\ell \cdot -0.11666666666666667\right)}}\\
\mathbf{else}:\\
\;\;\;\;t \cdot \left(-0.11666666666666667 \cdot \left(\ell \cdot \frac{\ell}{t \cdot t}\right)\right)\\
\end{array}
\end{array}
if t < 9.99999999999999955e-7Initial program 36.9%
Applied egg-rr62.7%
Taylor expanded in k around 0
Simplified8.1%
Taylor expanded in k around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6413.0%
Simplified13.0%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6413.0%
Applied egg-rr13.0%
if 9.99999999999999955e-7 < t Initial program 27.8%
Applied egg-rr59.6%
Taylor expanded in k around 0
Simplified23.9%
Taylor expanded in k around inf
*-lowering-*.f64N/A
unpow2N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6420.6%
Simplified20.6%
Final simplification15.0%
(FPCore (t l k) :precision binary64 (* l (/ (* l (/ 2.0 k)) (* t (* k (* k k))))))
double code(double t, double l, double k) {
return l * ((l * (2.0 / k)) / (t * (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 = l * ((l * (2.0d0 / k)) / (t * (k * (k * k))))
end function
public static double code(double t, double l, double k) {
return l * ((l * (2.0 / k)) / (t * (k * (k * k))));
}
def code(t, l, k): return l * ((l * (2.0 / k)) / (t * (k * (k * k))))
function code(t, l, k) return Float64(l * Float64(Float64(l * Float64(2.0 / k)) / Float64(t * Float64(k * Float64(k * k))))) end
function tmp = code(t, l, k) tmp = l * ((l * (2.0 / k)) / (t * (k * (k * k)))); end
code[t_, l_, k_] := N[(l * N[(N[(l * N[(2.0 / k), $MachinePrecision]), $MachinePrecision] / N[(t * N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\ell \cdot \frac{\ell \cdot \frac{2}{k}}{t \cdot \left(k \cdot \left(k \cdot k\right)\right)}
\end{array}
Initial program 34.5%
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-*.f6463.0%
Simplified63.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-/.f6468.9%
Applied egg-rr68.9%
div-invN/A
associate-/r*N/A
clear-numN/A
frac-timesN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6471.5%
Applied egg-rr71.5%
Final simplification71.5%
(FPCore (t l k) :precision binary64 (* (/ l t) (/ (* 2.0 l) (* k (* k (* k k))))))
double code(double t, double l, double k) {
return (l / t) * ((2.0 * l) / (k * (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 = (l / t) * ((2.0d0 * l) / (k * (k * (k * k))))
end function
public static double code(double t, double l, double k) {
return (l / t) * ((2.0 * l) / (k * (k * (k * k))));
}
def code(t, l, k): return (l / t) * ((2.0 * l) / (k * (k * (k * k))))
function code(t, l, k) return Float64(Float64(l / t) * Float64(Float64(2.0 * l) / Float64(k * Float64(k * Float64(k * k))))) end
function tmp = code(t, l, k) tmp = (l / t) * ((2.0 * l) / (k * (k * (k * k)))); end
code[t_, l_, k_] := N[(N[(l / t), $MachinePrecision] * N[(N[(2.0 * l), $MachinePrecision] / N[(k * N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell}{t} \cdot \frac{2 \cdot \ell}{k \cdot \left(k \cdot \left(k \cdot k\right)\right)}
\end{array}
Initial program 34.5%
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-*.f6463.0%
Simplified63.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-/.f6468.9%
Applied egg-rr68.9%
associate-*l/N/A
div-invN/A
clear-numN/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6470.0%
Applied egg-rr70.0%
Final simplification70.0%
(FPCore (t l k) :precision binary64 (/ 1.0 (/ t (* l (* l -0.11666666666666667)))))
double code(double t, double l, double k) {
return 1.0 / (t / (l * (l * -0.11666666666666667)));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 1.0d0 / (t / (l * (l * (-0.11666666666666667d0))))
end function
public static double code(double t, double l, double k) {
return 1.0 / (t / (l * (l * -0.11666666666666667)));
}
def code(t, l, k): return 1.0 / (t / (l * (l * -0.11666666666666667)))
function code(t, l, k) return Float64(1.0 / Float64(t / Float64(l * Float64(l * -0.11666666666666667)))) end
function tmp = code(t, l, k) tmp = 1.0 / (t / (l * (l * -0.11666666666666667))); end
code[t_, l_, k_] := N[(1.0 / N[(t / N[(l * N[(l * -0.11666666666666667), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{t}{\ell \cdot \left(\ell \cdot -0.11666666666666667\right)}}
\end{array}
Initial program 34.5%
Applied egg-rr61.9%
Taylor expanded in k around 0
Simplified12.1%
Taylor expanded in k around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6414.4%
Simplified14.4%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6414.4%
Applied egg-rr14.4%
(FPCore (t l k) :precision binary64 (/ (* l (* l -0.11666666666666667)) t))
double code(double t, double l, double k) {
return (l * (l * -0.11666666666666667)) / t;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = (l * (l * (-0.11666666666666667d0))) / t
end function
public static double code(double t, double l, double k) {
return (l * (l * -0.11666666666666667)) / t;
}
def code(t, l, k): return (l * (l * -0.11666666666666667)) / t
function code(t, l, k) return Float64(Float64(l * Float64(l * -0.11666666666666667)) / t) end
function tmp = code(t, l, k) tmp = (l * (l * -0.11666666666666667)) / t; end
code[t_, l_, k_] := N[(N[(l * N[(l * -0.11666666666666667), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell \cdot \left(\ell \cdot -0.11666666666666667\right)}{t}
\end{array}
Initial program 34.5%
Applied egg-rr61.9%
Taylor expanded in k around 0
Simplified12.1%
Taylor expanded in k around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6414.4%
Simplified14.4%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6414.4%
Applied egg-rr14.4%
Final simplification14.4%
(FPCore (t l k) :precision binary64 (* (/ (* l l) t) -0.11666666666666667))
double code(double t, double l, double k) {
return ((l * l) / t) * -0.11666666666666667;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = ((l * l) / t) * (-0.11666666666666667d0)
end function
public static double code(double t, double l, double k) {
return ((l * l) / t) * -0.11666666666666667;
}
def code(t, l, k): return ((l * l) / t) * -0.11666666666666667
function code(t, l, k) return Float64(Float64(Float64(l * l) / t) * -0.11666666666666667) end
function tmp = code(t, l, k) tmp = ((l * l) / t) * -0.11666666666666667; end
code[t_, l_, k_] := N[(N[(N[(l * l), $MachinePrecision] / t), $MachinePrecision] * -0.11666666666666667), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell \cdot \ell}{t} \cdot -0.11666666666666667
\end{array}
Initial program 34.5%
Applied egg-rr61.9%
Taylor expanded in k around 0
Simplified12.1%
Taylor expanded in k around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6414.4%
Simplified14.4%
div-invN/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
*-lowering-*.f6414.4%
Applied egg-rr14.4%
Final simplification14.4%
herbie shell --seed 2024161
(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))))