
(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 19 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}
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (/ (sin k) l)))
(*
t_s
(if (<= t_m 1.5e-101)
(/ 2.0 (* t_2 (/ (* k (* k (* t_m (sin k)))) (* l (cos k)))))
(if (<= t_m 7.6e+192)
(/
(/ (/ 2.0 t_2) t_m)
(* t_m (* (+ 2.0 (/ (* k (/ k t_m)) t_m)) (* t_m (/ (tan k) l)))))
(* (/ (/ (/ l (* t_m k)) t_m) k) (/ l t_m)))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = sin(k) / l;
double tmp;
if (t_m <= 1.5e-101) {
tmp = 2.0 / (t_2 * ((k * (k * (t_m * sin(k)))) / (l * cos(k))));
} else if (t_m <= 7.6e+192) {
tmp = ((2.0 / t_2) / t_m) / (t_m * ((2.0 + ((k * (k / t_m)) / t_m)) * (t_m * (tan(k) / l))));
} else {
tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = sin(k) / l
if (t_m <= 1.5d-101) then
tmp = 2.0d0 / (t_2 * ((k * (k * (t_m * sin(k)))) / (l * cos(k))))
else if (t_m <= 7.6d+192) then
tmp = ((2.0d0 / t_2) / t_m) / (t_m * ((2.0d0 + ((k * (k / t_m)) / t_m)) * (t_m * (tan(k) / l))))
else
tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double t_2 = Math.sin(k) / l;
double tmp;
if (t_m <= 1.5e-101) {
tmp = 2.0 / (t_2 * ((k * (k * (t_m * Math.sin(k)))) / (l * Math.cos(k))));
} else if (t_m <= 7.6e+192) {
tmp = ((2.0 / t_2) / t_m) / (t_m * ((2.0 + ((k * (k / t_m)) / t_m)) * (t_m * (Math.tan(k) / l))));
} else {
tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): t_2 = math.sin(k) / l tmp = 0 if t_m <= 1.5e-101: tmp = 2.0 / (t_2 * ((k * (k * (t_m * math.sin(k)))) / (l * math.cos(k)))) elif t_m <= 7.6e+192: tmp = ((2.0 / t_2) / t_m) / (t_m * ((2.0 + ((k * (k / t_m)) / t_m)) * (t_m * (math.tan(k) / l)))) else: tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(sin(k) / l) tmp = 0.0 if (t_m <= 1.5e-101) tmp = Float64(2.0 / Float64(t_2 * Float64(Float64(k * Float64(k * Float64(t_m * sin(k)))) / Float64(l * cos(k))))); elseif (t_m <= 7.6e+192) tmp = Float64(Float64(Float64(2.0 / t_2) / t_m) / Float64(t_m * Float64(Float64(2.0 + Float64(Float64(k * Float64(k / t_m)) / t_m)) * Float64(t_m * Float64(tan(k) / l))))); else tmp = Float64(Float64(Float64(Float64(l / Float64(t_m * k)) / t_m) / k) * Float64(l / t_m)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) t_2 = sin(k) / l; tmp = 0.0; if (t_m <= 1.5e-101) tmp = 2.0 / (t_2 * ((k * (k * (t_m * sin(k)))) / (l * cos(k)))); elseif (t_m <= 7.6e+192) tmp = ((2.0 / t_2) / t_m) / (t_m * ((2.0 + ((k * (k / t_m)) / t_m)) * (t_m * (tan(k) / l)))); else tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := Block[{t$95$2 = N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 1.5e-101], N[(2.0 / N[(t$95$2 * N[(N[(k * N[(k * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 7.6e+192], N[(N[(N[(2.0 / t$95$2), $MachinePrecision] / t$95$m), $MachinePrecision] / N[(t$95$m * N[(N[(2.0 + N[(N[(k * N[(k / t$95$m), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(N[Tan[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / k), $MachinePrecision] * N[(l / t$95$m), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{\sin k}{\ell}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.5 \cdot 10^{-101}:\\
\;\;\;\;\frac{2}{t\_2 \cdot \frac{k \cdot \left(k \cdot \left(t\_m \cdot \sin k\right)\right)}{\ell \cdot \cos k}}\\
\mathbf{elif}\;t\_m \leq 7.6 \cdot 10^{+192}:\\
\;\;\;\;\frac{\frac{\frac{2}{t\_2}}{t\_m}}{t\_m \cdot \left(\left(2 + \frac{k \cdot \frac{k}{t\_m}}{t\_m}\right) \cdot \left(t\_m \cdot \frac{\tan k}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{\ell}{t\_m \cdot k}}{t\_m}}{k} \cdot \frac{\ell}{t\_m}\\
\end{array}
\end{array}
\end{array}
if t < 1.5000000000000002e-101Initial program 60.6%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.0%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
Applied egg-rr66.0%
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr77.9%
Taylor expanded in t around 0
/-lowering-/.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6476.9%
Simplified76.9%
if 1.5000000000000002e-101 < t < 7.5999999999999999e192Initial program 59.7%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified52.9%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
Applied egg-rr69.5%
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr84.1%
associate-/r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
Applied egg-rr87.2%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6491.8%
Applied egg-rr91.8%
if 7.5999999999999999e192 < t Initial program 47.5%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified36.7%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6437.0%
Simplified37.0%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6457.8%
Applied egg-rr57.8%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6487.4%
Applied egg-rr87.4%
associate-/r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6490.7%
Applied egg-rr90.7%
Final simplification82.1%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 2.1e-52)
(* l (* l (/ (/ (/ 2.0 t_m) k) (* (sin k) (* k (tan k))))))
(if (<= t_m 5.6e+193)
(/
(/ (/ 2.0 (/ (sin k) l)) t_m)
(* t_m (* (+ 2.0 (/ (* k (/ k t_m)) t_m)) (* t_m (/ (tan k) l)))))
(* (/ (/ (/ l (* t_m k)) t_m) k) (/ l t_m))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 2.1e-52) {
tmp = l * (l * (((2.0 / t_m) / k) / (sin(k) * (k * tan(k)))));
} else if (t_m <= 5.6e+193) {
tmp = ((2.0 / (sin(k) / l)) / t_m) / (t_m * ((2.0 + ((k * (k / t_m)) / t_m)) * (t_m * (tan(k) / l))));
} else {
tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 2.1d-52) then
tmp = l * (l * (((2.0d0 / t_m) / k) / (sin(k) * (k * tan(k)))))
else if (t_m <= 5.6d+193) then
tmp = ((2.0d0 / (sin(k) / l)) / t_m) / (t_m * ((2.0d0 + ((k * (k / t_m)) / t_m)) * (t_m * (tan(k) / l))))
else
tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 2.1e-52) {
tmp = l * (l * (((2.0 / t_m) / k) / (Math.sin(k) * (k * Math.tan(k)))));
} else if (t_m <= 5.6e+193) {
tmp = ((2.0 / (Math.sin(k) / l)) / t_m) / (t_m * ((2.0 + ((k * (k / t_m)) / t_m)) * (t_m * (Math.tan(k) / l))));
} else {
tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 2.1e-52: tmp = l * (l * (((2.0 / t_m) / k) / (math.sin(k) * (k * math.tan(k))))) elif t_m <= 5.6e+193: tmp = ((2.0 / (math.sin(k) / l)) / t_m) / (t_m * ((2.0 + ((k * (k / t_m)) / t_m)) * (t_m * (math.tan(k) / l)))) else: tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 2.1e-52) tmp = Float64(l * Float64(l * Float64(Float64(Float64(2.0 / t_m) / k) / Float64(sin(k) * Float64(k * tan(k)))))); elseif (t_m <= 5.6e+193) tmp = Float64(Float64(Float64(2.0 / Float64(sin(k) / l)) / t_m) / Float64(t_m * Float64(Float64(2.0 + Float64(Float64(k * Float64(k / t_m)) / t_m)) * Float64(t_m * Float64(tan(k) / l))))); else tmp = Float64(Float64(Float64(Float64(l / Float64(t_m * k)) / t_m) / k) * Float64(l / t_m)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 2.1e-52) tmp = l * (l * (((2.0 / t_m) / k) / (sin(k) * (k * tan(k))))); elseif (t_m <= 5.6e+193) tmp = ((2.0 / (sin(k) / l)) / t_m) / (t_m * ((2.0 + ((k * (k / t_m)) / t_m)) * (t_m * (tan(k) / l)))); else tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 2.1e-52], N[(l * N[(l * N[(N[(N[(2.0 / t$95$m), $MachinePrecision] / k), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[(k * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 5.6e+193], N[(N[(N[(2.0 / N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / N[(t$95$m * N[(N[(2.0 + N[(N[(k * N[(k / t$95$m), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision] * N[(t$95$m * N[(N[Tan[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / k), $MachinePrecision] * N[(l / t$95$m), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 2.1 \cdot 10^{-52}:\\
\;\;\;\;\ell \cdot \left(\ell \cdot \frac{\frac{\frac{2}{t\_m}}{k}}{\sin k \cdot \left(k \cdot \tan k\right)}\right)\\
\mathbf{elif}\;t\_m \leq 5.6 \cdot 10^{+193}:\\
\;\;\;\;\frac{\frac{\frac{2}{\frac{\sin k}{\ell}}}{t\_m}}{t\_m \cdot \left(\left(2 + \frac{k \cdot \frac{k}{t\_m}}{t\_m}\right) \cdot \left(t\_m \cdot \frac{\tan k}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{\ell}{t\_m \cdot k}}{t\_m}}{k} \cdot \frac{\ell}{t\_m}\\
\end{array}
\end{array}
if t < 2.0999999999999999e-52Initial program 60.8%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.4%
Taylor expanded in k around inf
associate-/l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6430.8%
Simplified30.8%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr32.6%
associate-/r*N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr74.3%
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr77.1%
if 2.0999999999999999e-52 < t < 5.59999999999999972e193Initial program 58.7%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified50.4%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
Applied egg-rr68.6%
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr86.3%
associate-/r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
Applied egg-rr90.2%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6495.8%
Applied egg-rr95.8%
if 5.59999999999999972e193 < t Initial program 47.5%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified36.7%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6437.0%
Simplified37.0%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6457.8%
Applied egg-rr57.8%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6487.4%
Applied egg-rr87.4%
associate-/r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6490.7%
Applied egg-rr90.7%
Final simplification82.3%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 2.4e-52)
(* l (* l (/ (/ (/ 2.0 t_m) k) (* (sin k) (* k (tan k))))))
(if (<= t_m 5e+152)
(/
2.0
(*
(* t_m (/ (* (tan k) (+ 2.0 (/ (/ k t_m) (/ t_m k)))) l))
(* (/ (sin k) l) (* t_m t_m))))
(* (/ l t_m) (/ (/ (/ l t_m) (* t_m k)) k))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 2.4e-52) {
tmp = l * (l * (((2.0 / t_m) / k) / (sin(k) * (k * tan(k)))));
} else if (t_m <= 5e+152) {
tmp = 2.0 / ((t_m * ((tan(k) * (2.0 + ((k / t_m) / (t_m / k)))) / l)) * ((sin(k) / l) * (t_m * t_m)));
} else {
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 2.4d-52) then
tmp = l * (l * (((2.0d0 / t_m) / k) / (sin(k) * (k * tan(k)))))
else if (t_m <= 5d+152) then
tmp = 2.0d0 / ((t_m * ((tan(k) * (2.0d0 + ((k / t_m) / (t_m / k)))) / l)) * ((sin(k) / l) * (t_m * t_m)))
else
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 2.4e-52) {
tmp = l * (l * (((2.0 / t_m) / k) / (Math.sin(k) * (k * Math.tan(k)))));
} else if (t_m <= 5e+152) {
tmp = 2.0 / ((t_m * ((Math.tan(k) * (2.0 + ((k / t_m) / (t_m / k)))) / l)) * ((Math.sin(k) / l) * (t_m * t_m)));
} else {
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 2.4e-52: tmp = l * (l * (((2.0 / t_m) / k) / (math.sin(k) * (k * math.tan(k))))) elif t_m <= 5e+152: tmp = 2.0 / ((t_m * ((math.tan(k) * (2.0 + ((k / t_m) / (t_m / k)))) / l)) * ((math.sin(k) / l) * (t_m * t_m))) else: tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 2.4e-52) tmp = Float64(l * Float64(l * Float64(Float64(Float64(2.0 / t_m) / k) / Float64(sin(k) * Float64(k * tan(k)))))); elseif (t_m <= 5e+152) tmp = Float64(2.0 / Float64(Float64(t_m * Float64(Float64(tan(k) * Float64(2.0 + Float64(Float64(k / t_m) / Float64(t_m / k)))) / l)) * Float64(Float64(sin(k) / l) * Float64(t_m * t_m)))); else tmp = Float64(Float64(l / t_m) * Float64(Float64(Float64(l / t_m) / Float64(t_m * k)) / k)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 2.4e-52) tmp = l * (l * (((2.0 / t_m) / k) / (sin(k) * (k * tan(k))))); elseif (t_m <= 5e+152) tmp = 2.0 / ((t_m * ((tan(k) * (2.0 + ((k / t_m) / (t_m / k)))) / l)) * ((sin(k) / l) * (t_m * t_m))); else tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 2.4e-52], N[(l * N[(l * N[(N[(N[(2.0 / t$95$m), $MachinePrecision] / k), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[(k * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 5e+152], N[(2.0 / N[(N[(t$95$m * N[(N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[(N[(k / t$95$m), $MachinePrecision] / N[(t$95$m / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / t$95$m), $MachinePrecision] * N[(N[(N[(l / t$95$m), $MachinePrecision] / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 2.4 \cdot 10^{-52}:\\
\;\;\;\;\ell \cdot \left(\ell \cdot \frac{\frac{\frac{2}{t\_m}}{k}}{\sin k \cdot \left(k \cdot \tan k\right)}\right)\\
\mathbf{elif}\;t\_m \leq 5 \cdot 10^{+152}:\\
\;\;\;\;\frac{2}{\left(t\_m \cdot \frac{\tan k \cdot \left(2 + \frac{\frac{k}{t\_m}}{\frac{t\_m}{k}}\right)}{\ell}\right) \cdot \left(\frac{\sin k}{\ell} \cdot \left(t\_m \cdot t\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\frac{\frac{\ell}{t\_m}}{t\_m \cdot k}}{k}\\
\end{array}
\end{array}
if t < 2.4000000000000002e-52Initial program 60.8%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.4%
Taylor expanded in k around inf
associate-/l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6430.8%
Simplified30.8%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr32.6%
associate-/r*N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr74.3%
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr77.1%
if 2.4000000000000002e-52 < t < 5e152Initial program 59.8%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified50.5%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
Applied egg-rr70.9%
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
Applied egg-rr95.1%
if 5e152 < t Initial program 48.1%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified38.9%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6439.3%
Simplified39.3%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6459.5%
Applied egg-rr59.5%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6486.9%
Applied egg-rr86.9%
Final simplification81.5%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 4.7e-53)
(/ (/ l (* t_m k)) (* (* t_m k) (/ t_m l)))
(if (<= k 4.0)
(/
(/ (/ 2.0 (/ (sin k) l)) t_m)
(*
k
(+
(/ (* 2.0 (* t_m t_m)) l)
(/ (* k (* k (+ 1.0 (* (* t_m t_m) 0.6666666666666666)))) l))))
(* l (* l (/ (/ (/ 2.0 t_m) k) (* (sin k) (* k (tan k))))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 4.7e-53) {
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l));
} else if (k <= 4.0) {
tmp = ((2.0 / (sin(k) / l)) / t_m) / (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l)));
} else {
tmp = l * (l * (((2.0 / t_m) / k) / (sin(k) * (k * tan(k)))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 4.7d-53) then
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l))
else if (k <= 4.0d0) then
tmp = ((2.0d0 / (sin(k) / l)) / t_m) / (k * (((2.0d0 * (t_m * t_m)) / l) + ((k * (k * (1.0d0 + ((t_m * t_m) * 0.6666666666666666d0)))) / l)))
else
tmp = l * (l * (((2.0d0 / t_m) / k) / (sin(k) * (k * tan(k)))))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 4.7e-53) {
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l));
} else if (k <= 4.0) {
tmp = ((2.0 / (Math.sin(k) / l)) / t_m) / (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l)));
} else {
tmp = l * (l * (((2.0 / t_m) / k) / (Math.sin(k) * (k * Math.tan(k)))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if k <= 4.7e-53: tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l)) elif k <= 4.0: tmp = ((2.0 / (math.sin(k) / l)) / t_m) / (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l))) else: tmp = l * (l * (((2.0 / t_m) / k) / (math.sin(k) * (k * math.tan(k))))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 4.7e-53) tmp = Float64(Float64(l / Float64(t_m * k)) / Float64(Float64(t_m * k) * Float64(t_m / l))); elseif (k <= 4.0) tmp = Float64(Float64(Float64(2.0 / Float64(sin(k) / l)) / t_m) / Float64(k * Float64(Float64(Float64(2.0 * Float64(t_m * t_m)) / l) + Float64(Float64(k * Float64(k * Float64(1.0 + Float64(Float64(t_m * t_m) * 0.6666666666666666)))) / l)))); else tmp = Float64(l * Float64(l * Float64(Float64(Float64(2.0 / t_m) / k) / Float64(sin(k) * Float64(k * tan(k)))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (k <= 4.7e-53) tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l)); elseif (k <= 4.0) tmp = ((2.0 / (sin(k) / l)) / t_m) / (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l))); else tmp = l * (l * (((2.0 / t_m) / k) / (sin(k) * (k * tan(k))))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[k, 4.7e-53], N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$m * k), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 4.0], N[(N[(N[(2.0 / N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / N[(k * N[(N[(N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] + N[(N[(k * N[(k * N[(1.0 + N[(N[(t$95$m * t$95$m), $MachinePrecision] * 0.6666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l * N[(l * N[(N[(N[(2.0 / t$95$m), $MachinePrecision] / k), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[(k * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 4.7 \cdot 10^{-53}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot k}}{\left(t\_m \cdot k\right) \cdot \frac{t\_m}{\ell}}\\
\mathbf{elif}\;k \leq 4:\\
\;\;\;\;\frac{\frac{\frac{2}{\frac{\sin k}{\ell}}}{t\_m}}{k \cdot \left(\frac{2 \cdot \left(t\_m \cdot t\_m\right)}{\ell} + \frac{k \cdot \left(k \cdot \left(1 + \left(t\_m \cdot t\_m\right) \cdot 0.6666666666666666\right)\right)}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \left(\ell \cdot \frac{\frac{\frac{2}{t\_m}}{k}}{\sin k \cdot \left(k \cdot \tan k\right)}\right)\\
\end{array}
\end{array}
if k < 4.7e-53Initial program 57.7%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified51.3%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6456.1%
Simplified56.1%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6466.8%
Applied egg-rr66.8%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6474.9%
Applied egg-rr74.9%
*-commutativeN/A
clear-numN/A
associate-/r*N/A
frac-timesN/A
clear-numN/A
div-invN/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6481.1%
Applied egg-rr81.1%
if 4.7e-53 < k < 4Initial program 53.7%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified53.6%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
Applied egg-rr53.6%
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr54.5%
associate-/r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
Applied egg-rr66.3%
Taylor expanded in k around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified99.7%
if 4 < k Initial program 63.1%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified63.2%
Taylor expanded in k around inf
associate-/l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6429.2%
Simplified29.2%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr29.4%
associate-/r*N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr79.9%
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr83.1%
Final simplification82.5%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 4e+133)
(/
(/ (/ 2.0 (/ (sin k) l)) t_m)
(*
k
(+
(/ (* 2.0 (* t_m t_m)) l)
(/ (* k (* k (+ 1.0 (* (* t_m t_m) 0.6666666666666666)))) l))))
(* (/ l t_m) (/ (/ (/ l t_m) (* t_m k)) k)))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 4e+133) {
tmp = ((2.0 / (sin(k) / l)) / t_m) / (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l)));
} else {
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 4d+133) then
tmp = ((2.0d0 / (sin(k) / l)) / t_m) / (k * (((2.0d0 * (t_m * t_m)) / l) + ((k * (k * (1.0d0 + ((t_m * t_m) * 0.6666666666666666d0)))) / l)))
else
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 4e+133) {
tmp = ((2.0 / (Math.sin(k) / l)) / t_m) / (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l)));
} else {
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 4e+133: tmp = ((2.0 / (math.sin(k) / l)) / t_m) / (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l))) else: tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 4e+133) tmp = Float64(Float64(Float64(2.0 / Float64(sin(k) / l)) / t_m) / Float64(k * Float64(Float64(Float64(2.0 * Float64(t_m * t_m)) / l) + Float64(Float64(k * Float64(k * Float64(1.0 + Float64(Float64(t_m * t_m) * 0.6666666666666666)))) / l)))); else tmp = Float64(Float64(l / t_m) * Float64(Float64(Float64(l / t_m) / Float64(t_m * k)) / k)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 4e+133) tmp = ((2.0 / (sin(k) / l)) / t_m) / (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l))); else tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 4e+133], N[(N[(N[(2.0 / N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / N[(k * N[(N[(N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] + N[(N[(k * N[(k * N[(1.0 + N[(N[(t$95$m * t$95$m), $MachinePrecision] * 0.6666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / t$95$m), $MachinePrecision] * N[(N[(N[(l / t$95$m), $MachinePrecision] / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 4 \cdot 10^{+133}:\\
\;\;\;\;\frac{\frac{\frac{2}{\frac{\sin k}{\ell}}}{t\_m}}{k \cdot \left(\frac{2 \cdot \left(t\_m \cdot t\_m\right)}{\ell} + \frac{k \cdot \left(k \cdot \left(1 + \left(t\_m \cdot t\_m\right) \cdot 0.6666666666666666\right)\right)}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\frac{\frac{\ell}{t\_m}}{t\_m \cdot k}}{k}\\
\end{array}
\end{array}
if t < 4.0000000000000001e133Initial program 61.5%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified57.6%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
Applied egg-rr68.4%
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr80.3%
associate-/r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
Applied egg-rr82.6%
Taylor expanded in k around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified82.1%
if 4.0000000000000001e133 < t Initial program 44.9%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified36.7%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6439.4%
Simplified39.4%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6459.7%
Applied egg-rr59.7%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6483.8%
Applied egg-rr83.8%
Final simplification82.3%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 8e+132)
(/
2.0
(*
(/ (sin k) l)
(*
t_m
(*
k
(+
(/ (* 2.0 (* t_m t_m)) l)
(/ (* k (* k (+ 1.0 (* (* t_m t_m) 0.6666666666666666)))) l))))))
(* (/ l t_m) (/ (/ (/ l t_m) (* t_m k)) k)))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 8e+132) {
tmp = 2.0 / ((sin(k) / l) * (t_m * (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l)))));
} else {
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 8d+132) then
tmp = 2.0d0 / ((sin(k) / l) * (t_m * (k * (((2.0d0 * (t_m * t_m)) / l) + ((k * (k * (1.0d0 + ((t_m * t_m) * 0.6666666666666666d0)))) / l)))))
else
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 8e+132) {
tmp = 2.0 / ((Math.sin(k) / l) * (t_m * (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l)))));
} else {
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 8e+132: tmp = 2.0 / ((math.sin(k) / l) * (t_m * (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l))))) else: tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 8e+132) tmp = Float64(2.0 / Float64(Float64(sin(k) / l) * Float64(t_m * Float64(k * Float64(Float64(Float64(2.0 * Float64(t_m * t_m)) / l) + Float64(Float64(k * Float64(k * Float64(1.0 + Float64(Float64(t_m * t_m) * 0.6666666666666666)))) / l)))))); else tmp = Float64(Float64(l / t_m) * Float64(Float64(Float64(l / t_m) / Float64(t_m * k)) / k)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 8e+132) tmp = 2.0 / ((sin(k) / l) * (t_m * (k * (((2.0 * (t_m * t_m)) / l) + ((k * (k * (1.0 + ((t_m * t_m) * 0.6666666666666666)))) / l))))); else tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 8e+132], N[(2.0 / N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[(t$95$m * N[(k * N[(N[(N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] + N[(N[(k * N[(k * N[(1.0 + N[(N[(t$95$m * t$95$m), $MachinePrecision] * 0.6666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / t$95$m), $MachinePrecision] * N[(N[(N[(l / t$95$m), $MachinePrecision] / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 8 \cdot 10^{+132}:\\
\;\;\;\;\frac{2}{\frac{\sin k}{\ell} \cdot \left(t\_m \cdot \left(k \cdot \left(\frac{2 \cdot \left(t\_m \cdot t\_m\right)}{\ell} + \frac{k \cdot \left(k \cdot \left(1 + \left(t\_m \cdot t\_m\right) \cdot 0.6666666666666666\right)\right)}{\ell}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\frac{\frac{\ell}{t\_m}}{t\_m \cdot k}}{k}\\
\end{array}
\end{array}
if t < 7.99999999999999993e132Initial program 61.5%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified57.6%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
Applied egg-rr68.4%
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr80.3%
Taylor expanded in k around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified81.4%
if 7.99999999999999993e132 < t Initial program 44.9%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified36.7%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6439.4%
Simplified39.4%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6459.7%
Applied egg-rr59.7%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6483.8%
Applied egg-rr83.8%
Final simplification81.8%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= (* l l) 1e+78)
(/
(/ (/ (* 2.0 l) k) t_m)
(* t_m (* (* t_m (/ (tan k) l)) (+ 2.0 (/ (/ (* k k) t_m) t_m)))))
(* (/ l t_m) (/ (/ (/ l t_m) (* t_m k)) k)))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if ((l * l) <= 1e+78) {
tmp = (((2.0 * l) / k) / t_m) / (t_m * ((t_m * (tan(k) / l)) * (2.0 + (((k * k) / t_m) / t_m))));
} else {
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if ((l * l) <= 1d+78) then
tmp = (((2.0d0 * l) / k) / t_m) / (t_m * ((t_m * (tan(k) / l)) * (2.0d0 + (((k * k) / t_m) / t_m))))
else
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if ((l * l) <= 1e+78) {
tmp = (((2.0 * l) / k) / t_m) / (t_m * ((t_m * (Math.tan(k) / l)) * (2.0 + (((k * k) / t_m) / t_m))));
} else {
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if (l * l) <= 1e+78: tmp = (((2.0 * l) / k) / t_m) / (t_m * ((t_m * (math.tan(k) / l)) * (2.0 + (((k * k) / t_m) / t_m)))) else: tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (Float64(l * l) <= 1e+78) tmp = Float64(Float64(Float64(Float64(2.0 * l) / k) / t_m) / Float64(t_m * Float64(Float64(t_m * Float64(tan(k) / l)) * Float64(2.0 + Float64(Float64(Float64(k * k) / t_m) / t_m))))); else tmp = Float64(Float64(l / t_m) * Float64(Float64(Float64(l / t_m) / Float64(t_m * k)) / k)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if ((l * l) <= 1e+78) tmp = (((2.0 * l) / k) / t_m) / (t_m * ((t_m * (tan(k) / l)) * (2.0 + (((k * k) / t_m) / t_m)))); else tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 1e+78], N[(N[(N[(N[(2.0 * l), $MachinePrecision] / k), $MachinePrecision] / t$95$m), $MachinePrecision] / N[(t$95$m * N[(N[(t$95$m * N[(N[Tan[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(2.0 + N[(N[(N[(k * k), $MachinePrecision] / t$95$m), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / t$95$m), $MachinePrecision] * N[(N[(N[(l / t$95$m), $MachinePrecision] / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{+78}:\\
\;\;\;\;\frac{\frac{\frac{2 \cdot \ell}{k}}{t\_m}}{t\_m \cdot \left(\left(t\_m \cdot \frac{\tan k}{\ell}\right) \cdot \left(2 + \frac{\frac{k \cdot k}{t\_m}}{t\_m}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\frac{\frac{\ell}{t\_m}}{t\_m \cdot k}}{k}\\
\end{array}
\end{array}
if (*.f64 l l) < 1.00000000000000001e78Initial program 64.5%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.0%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
Applied egg-rr74.2%
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr87.2%
associate-/r*N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/l*N/A
Applied egg-rr90.2%
Taylor expanded in k around 0
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6488.9%
Simplified88.9%
if 1.00000000000000001e78 < (*.f64 l l) Initial program 51.9%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified49.7%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.6%
Simplified59.6%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6468.7%
Applied egg-rr68.7%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6477.2%
Applied egg-rr77.2%
Final simplification83.6%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= (* l l) 1e-124)
(/
2.0
(*
(/ k l)
(*
t_m
(* t_m (* t_m (/ (* (tan k) (+ 2.0 (/ (/ k t_m) (/ t_m k)))) l))))))
(/ (/ l (* t_m k)) (* (* t_m k) (/ t_m l))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if ((l * l) <= 1e-124) {
tmp = 2.0 / ((k / l) * (t_m * (t_m * (t_m * ((tan(k) * (2.0 + ((k / t_m) / (t_m / k)))) / l)))));
} else {
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if ((l * l) <= 1d-124) then
tmp = 2.0d0 / ((k / l) * (t_m * (t_m * (t_m * ((tan(k) * (2.0d0 + ((k / t_m) / (t_m / k)))) / l)))))
else
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if ((l * l) <= 1e-124) {
tmp = 2.0 / ((k / l) * (t_m * (t_m * (t_m * ((Math.tan(k) * (2.0 + ((k / t_m) / (t_m / k)))) / l)))));
} else {
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if (l * l) <= 1e-124: tmp = 2.0 / ((k / l) * (t_m * (t_m * (t_m * ((math.tan(k) * (2.0 + ((k / t_m) / (t_m / k)))) / l))))) else: tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l)) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (Float64(l * l) <= 1e-124) tmp = Float64(2.0 / Float64(Float64(k / l) * Float64(t_m * Float64(t_m * Float64(t_m * Float64(Float64(tan(k) * Float64(2.0 + Float64(Float64(k / t_m) / Float64(t_m / k)))) / l)))))); else tmp = Float64(Float64(l / Float64(t_m * k)) / Float64(Float64(t_m * k) * Float64(t_m / l))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if ((l * l) <= 1e-124) tmp = 2.0 / ((k / l) * (t_m * (t_m * (t_m * ((tan(k) * (2.0 + ((k / t_m) / (t_m / k)))) / l))))); else tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l)); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 1e-124], N[(2.0 / N[(N[(k / l), $MachinePrecision] * N[(t$95$m * N[(t$95$m * N[(t$95$m * N[(N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[(N[(k / t$95$m), $MachinePrecision] / N[(t$95$m / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$m * k), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-124}:\\
\;\;\;\;\frac{2}{\frac{k}{\ell} \cdot \left(t\_m \cdot \left(t\_m \cdot \left(t\_m \cdot \frac{\tan k \cdot \left(2 + \frac{\frac{k}{t\_m}}{\frac{t\_m}{k}}\right)}{\ell}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot k}}{\left(t\_m \cdot k\right) \cdot \frac{t\_m}{\ell}}\\
\end{array}
\end{array}
if (*.f64 l l) < 9.99999999999999933e-125Initial program 65.7%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.8%
*-commutativeN/A
times-fracN/A
associate-*l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
Applied egg-rr78.1%
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr92.7%
Taylor expanded in k around 0
/-lowering-/.f6491.9%
Simplified91.9%
if 9.99999999999999933e-125 < (*.f64 l l) Initial program 54.2%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified51.2%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.2%
Simplified58.2%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6465.7%
Applied egg-rr65.7%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.9%
Applied egg-rr73.9%
*-commutativeN/A
clear-numN/A
associate-/r*N/A
frac-timesN/A
clear-numN/A
div-invN/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6475.8%
Applied egg-rr75.8%
Final simplification82.3%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 7.5e-125)
(* l (/ (/ 2.0 t_m) (/ (/ (* k k) (/ (/ 1.0 (tan k)) k)) l)))
(/ (/ l (* t_m k)) (* (* t_m k) (/ t_m l))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 7.5e-125) {
tmp = l * ((2.0 / t_m) / (((k * k) / ((1.0 / tan(k)) / k)) / l));
} else {
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 7.5d-125) then
tmp = l * ((2.0d0 / t_m) / (((k * k) / ((1.0d0 / tan(k)) / k)) / l))
else
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 7.5e-125) {
tmp = l * ((2.0 / t_m) / (((k * k) / ((1.0 / Math.tan(k)) / k)) / l));
} else {
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 7.5e-125: tmp = l * ((2.0 / t_m) / (((k * k) / ((1.0 / math.tan(k)) / k)) / l)) else: tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l)) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 7.5e-125) tmp = Float64(l * Float64(Float64(2.0 / t_m) / Float64(Float64(Float64(k * k) / Float64(Float64(1.0 / tan(k)) / k)) / l))); else tmp = Float64(Float64(l / Float64(t_m * k)) / Float64(Float64(t_m * k) * Float64(t_m / l))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 7.5e-125) tmp = l * ((2.0 / t_m) / (((k * k) / ((1.0 / tan(k)) / k)) / l)); else tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l)); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 7.5e-125], N[(l * N[(N[(2.0 / t$95$m), $MachinePrecision] / N[(N[(N[(k * k), $MachinePrecision] / N[(N[(1.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$m * k), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 7.5 \cdot 10^{-125}:\\
\;\;\;\;\ell \cdot \frac{\frac{2}{t\_m}}{\frac{\frac{k \cdot k}{\frac{\frac{1}{\tan k}}{k}}}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot k}}{\left(t\_m \cdot k\right) \cdot \frac{t\_m}{\ell}}\\
\end{array}
\end{array}
if t < 7.5e-125Initial program 61.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.8%
Taylor expanded in k around inf
associate-/l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6427.4%
Simplified27.4%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr31.2%
associate-/r*N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr73.7%
Taylor expanded in k around 0
Simplified67.1%
if 7.5e-125 < t Initial program 54.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified46.6%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6449.1%
Simplified49.1%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6460.2%
Applied egg-rr60.2%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6472.7%
Applied egg-rr72.7%
*-commutativeN/A
clear-numN/A
associate-/r*N/A
frac-timesN/A
clear-numN/A
div-invN/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6479.7%
Applied egg-rr79.7%
Final simplification71.7%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 7.8e-128)
(* (/ l t_m) (/ l (* t_m (* t_m (* k k)))))
(if (<= t_m 1.9e+145)
(* (/ l (* t_m k)) (/ l (* k (* t_m t_m))))
(* (/ l t_m) (/ l (* (* t_m k) (* t_m k))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 7.8e-128) {
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))));
} else if (t_m <= 1.9e+145) {
tmp = (l / (t_m * k)) * (l / (k * (t_m * t_m)));
} else {
tmp = (l / t_m) * (l / ((t_m * k) * (t_m * k)));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 7.8d-128) then
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))))
else if (t_m <= 1.9d+145) then
tmp = (l / (t_m * k)) * (l / (k * (t_m * t_m)))
else
tmp = (l / t_m) * (l / ((t_m * k) * (t_m * k)))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 7.8e-128) {
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))));
} else if (t_m <= 1.9e+145) {
tmp = (l / (t_m * k)) * (l / (k * (t_m * t_m)));
} else {
tmp = (l / t_m) * (l / ((t_m * k) * (t_m * k)));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 7.8e-128: tmp = (l / t_m) * (l / (t_m * (t_m * (k * k)))) elif t_m <= 1.9e+145: tmp = (l / (t_m * k)) * (l / (k * (t_m * t_m))) else: tmp = (l / t_m) * (l / ((t_m * k) * (t_m * k))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 7.8e-128) tmp = Float64(Float64(l / t_m) * Float64(l / Float64(t_m * Float64(t_m * Float64(k * k))))); elseif (t_m <= 1.9e+145) tmp = Float64(Float64(l / Float64(t_m * k)) * Float64(l / Float64(k * Float64(t_m * t_m)))); else tmp = Float64(Float64(l / t_m) * Float64(l / Float64(Float64(t_m * k) * Float64(t_m * k)))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 7.8e-128) tmp = (l / t_m) * (l / (t_m * (t_m * (k * k)))); elseif (t_m <= 1.9e+145) tmp = (l / (t_m * k)) * (l / (k * (t_m * t_m))); else tmp = (l / t_m) * (l / ((t_m * k) * (t_m * k))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 7.8e-128], N[(N[(l / t$95$m), $MachinePrecision] * N[(l / N[(t$95$m * N[(t$95$m * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 1.9e+145], N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] * N[(l / N[(k * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / t$95$m), $MachinePrecision] * N[(l / N[(N[(t$95$m * k), $MachinePrecision] * N[(t$95$m * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 7.8 \cdot 10^{-128}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\ell}{t\_m \cdot \left(t\_m \cdot \left(k \cdot k\right)\right)}\\
\mathbf{elif}\;t\_m \leq 1.9 \cdot 10^{+145}:\\
\;\;\;\;\frac{\ell}{t\_m \cdot k} \cdot \frac{\ell}{k \cdot \left(t\_m \cdot t\_m\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\ell}{\left(t\_m \cdot k\right) \cdot \left(t\_m \cdot k\right)}\\
\end{array}
\end{array}
if t < 7.79999999999999993e-128Initial program 61.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.8%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6464.7%
Simplified64.7%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6473.4%
Applied egg-rr73.4%
if 7.79999999999999993e-128 < t < 1.90000000000000006e145Initial program 61.3%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified54.1%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6456.3%
Simplified56.3%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6462.0%
Applied egg-rr62.0%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6467.2%
Applied egg-rr67.2%
associate-*l/N/A
associate-*r*N/A
times-fracN/A
associate-/r*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6475.8%
Applied egg-rr75.8%
if 1.90000000000000006e145 < t Initial program 44.5%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified35.9%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6438.7%
Simplified38.7%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6457.7%
Applied egg-rr57.7%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6480.5%
Applied egg-rr80.5%
Final simplification75.0%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 7.8e-125)
(* l (/ (/ 2.0 t_m) (/ (/ (* k k) (/ 1.0 (* k k))) l)))
(/ (/ l (* t_m k)) (* (* t_m k) (/ t_m l))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 7.8e-125) {
tmp = l * ((2.0 / t_m) / (((k * k) / (1.0 / (k * k))) / l));
} else {
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 7.8d-125) then
tmp = l * ((2.0d0 / t_m) / (((k * k) / (1.0d0 / (k * k))) / l))
else
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 7.8e-125) {
tmp = l * ((2.0 / t_m) / (((k * k) / (1.0 / (k * k))) / l));
} else {
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 7.8e-125: tmp = l * ((2.0 / t_m) / (((k * k) / (1.0 / (k * k))) / l)) else: tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l)) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 7.8e-125) tmp = Float64(l * Float64(Float64(2.0 / t_m) / Float64(Float64(Float64(k * k) / Float64(1.0 / Float64(k * k))) / l))); else tmp = Float64(Float64(l / Float64(t_m * k)) / Float64(Float64(t_m * k) * Float64(t_m / l))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 7.8e-125) tmp = l * ((2.0 / t_m) / (((k * k) / (1.0 / (k * k))) / l)); else tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l)); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 7.8e-125], N[(l * N[(N[(2.0 / t$95$m), $MachinePrecision] / N[(N[(N[(k * k), $MachinePrecision] / N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$m * k), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 7.8 \cdot 10^{-125}:\\
\;\;\;\;\ell \cdot \frac{\frac{2}{t\_m}}{\frac{\frac{k \cdot k}{\frac{1}{k \cdot k}}}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot k}}{\left(t\_m \cdot k\right) \cdot \frac{t\_m}{\ell}}\\
\end{array}
\end{array}
if t < 7.79999999999999965e-125Initial program 61.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.8%
Taylor expanded in k around inf
associate-/l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6427.4%
Simplified27.4%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr31.2%
associate-/r*N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr73.7%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6467.1%
Simplified67.1%
if 7.79999999999999965e-125 < t Initial program 54.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified46.6%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6449.1%
Simplified49.1%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6460.2%
Applied egg-rr60.2%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6472.7%
Applied egg-rr72.7%
*-commutativeN/A
clear-numN/A
associate-/r*N/A
frac-timesN/A
clear-numN/A
div-invN/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6479.7%
Applied egg-rr79.7%
Final simplification71.8%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 1e-66)
(/ (/ l (* t_m k)) (* (* t_m k) (/ t_m l)))
(/ 1.0 (/ t_m (/ (/ l (/ t_m l)) (* t_m (* k k))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 1e-66) {
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l));
} else {
tmp = 1.0 / (t_m / ((l / (t_m / l)) / (t_m * (k * k))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 1d-66) then
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l))
else
tmp = 1.0d0 / (t_m / ((l / (t_m / l)) / (t_m * (k * k))))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 1e-66) {
tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l));
} else {
tmp = 1.0 / (t_m / ((l / (t_m / l)) / (t_m * (k * k))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if k <= 1e-66: tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l)) else: tmp = 1.0 / (t_m / ((l / (t_m / l)) / (t_m * (k * k)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 1e-66) tmp = Float64(Float64(l / Float64(t_m * k)) / Float64(Float64(t_m * k) * Float64(t_m / l))); else tmp = Float64(1.0 / Float64(t_m / Float64(Float64(l / Float64(t_m / l)) / Float64(t_m * Float64(k * k))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (k <= 1e-66) tmp = (l / (t_m * k)) / ((t_m * k) * (t_m / l)); else tmp = 1.0 / (t_m / ((l / (t_m / l)) / (t_m * (k * k)))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[k, 1e-66], N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$m * k), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(t$95$m / N[(N[(l / N[(t$95$m / l), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 10^{-66}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot k}}{\left(t\_m \cdot k\right) \cdot \frac{t\_m}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{t\_m}{\frac{\frac{\ell}{\frac{t\_m}{\ell}}}{t\_m \cdot \left(k \cdot k\right)}}}\\
\end{array}
\end{array}
if k < 9.9999999999999998e-67Initial program 58.0%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified51.6%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6455.9%
Simplified55.9%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6466.6%
Applied egg-rr66.6%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6474.7%
Applied egg-rr74.7%
*-commutativeN/A
clear-numN/A
associate-/r*N/A
frac-timesN/A
clear-numN/A
div-invN/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6481.0%
Applied egg-rr81.0%
if 9.9999999999999998e-67 < k Initial program 60.7%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified60.8%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6466.3%
Simplified66.3%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6473.0%
Applied egg-rr73.0%
associate-*r/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-/r*N/A
associate-*l/N/A
*-commutativeN/A
/-lowering-/.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.4%
Applied egg-rr73.4%
Final simplification78.8%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 3.5e+91)
(/ (* (/ l (* t_m k)) (/ l t_m)) (* t_m k))
(/ (/ (* l l) (* t_m (* t_m (* k k)))) t_m))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 3.5e+91) {
tmp = ((l / (t_m * k)) * (l / t_m)) / (t_m * k);
} else {
tmp = ((l * l) / (t_m * (t_m * (k * k)))) / t_m;
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 3.5d+91) then
tmp = ((l / (t_m * k)) * (l / t_m)) / (t_m * k)
else
tmp = ((l * l) / (t_m * (t_m * (k * k)))) / t_m
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 3.5e+91) {
tmp = ((l / (t_m * k)) * (l / t_m)) / (t_m * k);
} else {
tmp = ((l * l) / (t_m * (t_m * (k * k)))) / t_m;
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if k <= 3.5e+91: tmp = ((l / (t_m * k)) * (l / t_m)) / (t_m * k) else: tmp = ((l * l) / (t_m * (t_m * (k * k)))) / t_m return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 3.5e+91) tmp = Float64(Float64(Float64(l / Float64(t_m * k)) * Float64(l / t_m)) / Float64(t_m * k)); else tmp = Float64(Float64(Float64(l * l) / Float64(t_m * Float64(t_m * Float64(k * k)))) / t_m); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (k <= 3.5e+91) tmp = ((l / (t_m * k)) * (l / t_m)) / (t_m * k); else tmp = ((l * l) / (t_m * (t_m * (k * k)))) / t_m; end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[k, 3.5e+91], N[(N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] * N[(l / t$95$m), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision], N[(N[(N[(l * l), $MachinePrecision] / N[(t$95$m * N[(t$95$m * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 3.5 \cdot 10^{+91}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot k} \cdot \frac{\ell}{t\_m}}{t\_m \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell \cdot \ell}{t\_m \cdot \left(t\_m \cdot \left(k \cdot k\right)\right)}}{t\_m}\\
\end{array}
\end{array}
if k < 3.50000000000000001e91Initial program 58.9%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified53.2%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.0%
Simplified59.0%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6468.7%
Applied egg-rr68.7%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.7%
Applied egg-rr75.7%
associate-/r*N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6481.5%
Applied egg-rr81.5%
if 3.50000000000000001e91 < k Initial program 58.6%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.7%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.7%
Simplified58.7%
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6467.2%
Applied egg-rr67.2%
Final simplification78.8%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 6.2e+121)
(/ (* (/ l (* t_m k)) (/ l t_m)) (* t_m k))
(* (/ l t_m) (/ l (* t_m (* t_m (* k k))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 6.2e+121) {
tmp = ((l / (t_m * k)) * (l / t_m)) / (t_m * k);
} else {
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 6.2d+121) then
tmp = ((l / (t_m * k)) * (l / t_m)) / (t_m * k)
else
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 6.2e+121) {
tmp = ((l / (t_m * k)) * (l / t_m)) / (t_m * k);
} else {
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if k <= 6.2e+121: tmp = ((l / (t_m * k)) * (l / t_m)) / (t_m * k) else: tmp = (l / t_m) * (l / (t_m * (t_m * (k * k)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 6.2e+121) tmp = Float64(Float64(Float64(l / Float64(t_m * k)) * Float64(l / t_m)) / Float64(t_m * k)); else tmp = Float64(Float64(l / t_m) * Float64(l / Float64(t_m * Float64(t_m * Float64(k * k))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (k <= 6.2e+121) tmp = ((l / (t_m * k)) * (l / t_m)) / (t_m * k); else tmp = (l / t_m) * (l / (t_m * (t_m * (k * k)))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[k, 6.2e+121], N[(N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] * N[(l / t$95$m), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision], N[(N[(l / t$95$m), $MachinePrecision] * N[(l / N[(t$95$m * N[(t$95$m * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 6.2 \cdot 10^{+121}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot k} \cdot \frac{\ell}{t\_m}}{t\_m \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\ell}{t\_m \cdot \left(t\_m \cdot \left(k \cdot k\right)\right)}\\
\end{array}
\end{array}
if k < 6.20000000000000016e121Initial program 58.9%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified53.4%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.1%
Simplified58.1%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6467.6%
Applied egg-rr67.6%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6474.4%
Applied egg-rr74.4%
associate-/r*N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6480.1%
Applied egg-rr80.1%
if 6.20000000000000016e121 < k Initial program 58.3%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.4%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6463.1%
Simplified63.1%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6473.2%
Applied egg-rr73.2%
Final simplification78.9%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1e-126)
(* (/ l t_m) (/ l (* t_m (* t_m (* k k)))))
(* (/ (/ (/ l (* t_m k)) t_m) k) (/ l t_m)))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1e-126) {
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))));
} else {
tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1d-126) then
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))))
else
tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1e-126) {
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))));
} else {
tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1e-126: tmp = (l / t_m) * (l / (t_m * (t_m * (k * k)))) else: tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1e-126) tmp = Float64(Float64(l / t_m) * Float64(l / Float64(t_m * Float64(t_m * Float64(k * k))))); else tmp = Float64(Float64(Float64(Float64(l / Float64(t_m * k)) / t_m) / k) * Float64(l / t_m)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1e-126) tmp = (l / t_m) * (l / (t_m * (t_m * (k * k)))); else tmp = (((l / (t_m * k)) / t_m) / k) * (l / t_m); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 1e-126], N[(N[(l / t$95$m), $MachinePrecision] * N[(l / N[(t$95$m * N[(t$95$m * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / k), $MachinePrecision] * N[(l / t$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 10^{-126}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\ell}{t\_m \cdot \left(t\_m \cdot \left(k \cdot k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{\ell}{t\_m \cdot k}}{t\_m}}{k} \cdot \frac{\ell}{t\_m}\\
\end{array}
\end{array}
if t < 9.9999999999999995e-127Initial program 61.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.8%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6464.7%
Simplified64.7%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6473.4%
Applied egg-rr73.4%
if 9.9999999999999995e-127 < t Initial program 54.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified46.6%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6449.1%
Simplified49.1%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6460.2%
Applied egg-rr60.2%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6472.7%
Applied egg-rr72.7%
associate-/r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6478.7%
Applied egg-rr78.7%
Final simplification75.4%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 7.8e-128)
(* (/ l t_m) (/ l (* t_m (* t_m (* k k)))))
(* (/ l t_m) (/ (/ (/ l t_m) (* t_m k)) k)))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 7.8e-128) {
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))));
} else {
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 7.8d-128) then
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))))
else
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 7.8e-128) {
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))));
} else {
tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 7.8e-128: tmp = (l / t_m) * (l / (t_m * (t_m * (k * k)))) else: tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 7.8e-128) tmp = Float64(Float64(l / t_m) * Float64(l / Float64(t_m * Float64(t_m * Float64(k * k))))); else tmp = Float64(Float64(l / t_m) * Float64(Float64(Float64(l / t_m) / Float64(t_m * k)) / k)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 7.8e-128) tmp = (l / t_m) * (l / (t_m * (t_m * (k * k)))); else tmp = (l / t_m) * (((l / t_m) / (t_m * k)) / k); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 7.8e-128], N[(N[(l / t$95$m), $MachinePrecision] * N[(l / N[(t$95$m * N[(t$95$m * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / t$95$m), $MachinePrecision] * N[(N[(N[(l / t$95$m), $MachinePrecision] / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 7.8 \cdot 10^{-128}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\ell}{t\_m \cdot \left(t\_m \cdot \left(k \cdot k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\frac{\frac{\ell}{t\_m}}{t\_m \cdot k}}{k}\\
\end{array}
\end{array}
if t < 7.79999999999999993e-128Initial program 61.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.8%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6464.7%
Simplified64.7%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6473.4%
Applied egg-rr73.4%
if 7.79999999999999993e-128 < t Initial program 54.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified46.6%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6449.1%
Simplified49.1%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6460.2%
Applied egg-rr60.2%
associate-/r*N/A
associate-*r*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6478.7%
Applied egg-rr78.7%
Final simplification75.4%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 1e-156)
(* (/ l (* t_m k)) (/ l (* k (* t_m t_m))))
(* (/ l t_m) (/ l (* t_m (* t_m (* k k))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 1e-156) {
tmp = (l / (t_m * k)) * (l / (k * (t_m * t_m)));
} else {
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 1d-156) then
tmp = (l / (t_m * k)) * (l / (k * (t_m * t_m)))
else
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 1e-156) {
tmp = (l / (t_m * k)) * (l / (k * (t_m * t_m)));
} else {
tmp = (l / t_m) * (l / (t_m * (t_m * (k * k))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if k <= 1e-156: tmp = (l / (t_m * k)) * (l / (k * (t_m * t_m))) else: tmp = (l / t_m) * (l / (t_m * (t_m * (k * k)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 1e-156) tmp = Float64(Float64(l / Float64(t_m * k)) * Float64(l / Float64(k * Float64(t_m * t_m)))); else tmp = Float64(Float64(l / t_m) * Float64(l / Float64(t_m * Float64(t_m * Float64(k * k))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (k <= 1e-156) tmp = (l / (t_m * k)) * (l / (k * (t_m * t_m))); else tmp = (l / t_m) * (l / (t_m * (t_m * (k * k)))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[k, 1e-156], N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] * N[(l / N[(k * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / t$95$m), $MachinePrecision] * N[(l / N[(t$95$m * N[(t$95$m * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 10^{-156}:\\
\;\;\;\;\frac{\ell}{t\_m \cdot k} \cdot \frac{\ell}{k \cdot \left(t\_m \cdot t\_m\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{t\_m} \cdot \frac{\ell}{t\_m \cdot \left(t\_m \cdot \left(k \cdot k\right)\right)}\\
\end{array}
\end{array}
if k < 1.00000000000000004e-156Initial program 58.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified51.8%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6455.9%
Simplified55.9%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6464.0%
Applied egg-rr64.0%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.3%
Applied egg-rr73.3%
associate-*l/N/A
associate-*r*N/A
times-fracN/A
associate-/r*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6469.6%
Applied egg-rr69.6%
if 1.00000000000000004e-156 < k Initial program 59.4%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified58.2%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6463.7%
Simplified63.7%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6475.6%
Applied egg-rr75.6%
Final simplification71.9%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k) :precision binary64 (* t_s (* (/ l (* t_m k)) (/ l (* k (* t_m t_m))))))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
return t_s * ((l / (t_m * k)) * (l / (k * (t_m * t_m))));
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
code = t_s * ((l / (t_m * k)) * (l / (k * (t_m * t_m))))
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
return t_s * ((l / (t_m * k)) * (l / (k * (t_m * t_m))));
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): return t_s * ((l / (t_m * k)) * (l / (k * (t_m * t_m))))
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) return Float64(t_s * Float64(Float64(l / Float64(t_m * k)) * Float64(l / Float64(k * Float64(t_m * t_m))))) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k) tmp = t_s * ((l / (t_m * k)) * (l / (k * (t_m * t_m)))); end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * N[(N[(l / N[(t$95$m * k), $MachinePrecision]), $MachinePrecision] * N[(l / N[(k * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\frac{\ell}{t\_m \cdot k} \cdot \frac{\ell}{k \cdot \left(t\_m \cdot t\_m\right)}\right)
\end{array}
Initial program 58.8%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified54.3%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.9%
Simplified58.9%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6468.5%
Applied egg-rr68.5%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.5%
Applied egg-rr73.5%
associate-*l/N/A
associate-*r*N/A
times-fracN/A
associate-/r*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6469.0%
Applied egg-rr69.0%
Final simplification69.0%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k) :precision binary64 (* t_s (* l (/ l (* t_m (* k (* k (* t_m t_m))))))))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
return t_s * (l * (l / (t_m * (k * (k * (t_m * t_m))))));
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
code = t_s * (l * (l / (t_m * (k * (k * (t_m * t_m))))))
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
return t_s * (l * (l / (t_m * (k * (k * (t_m * t_m))))));
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): return t_s * (l * (l / (t_m * (k * (k * (t_m * t_m))))))
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) return Float64(t_s * Float64(l * Float64(l / Float64(t_m * Float64(k * Float64(k * Float64(t_m * t_m))))))) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k) tmp = t_s * (l * (l / (t_m * (k * (k * (t_m * t_m)))))); end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * N[(l * N[(l / N[(t$95$m * N[(k * N[(k * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\ell \cdot \frac{\ell}{t\_m \cdot \left(k \cdot \left(k \cdot \left(t\_m \cdot t\_m\right)\right)\right)}\right)
\end{array}
Initial program 58.8%
/-lowering-/.f64N/A
associate-*l*N/A
associate-*l/N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
Simplified54.3%
Taylor expanded in k around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.9%
Simplified58.9%
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6468.5%
Applied egg-rr68.5%
associate-*r*N/A
unswap-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.5%
Applied egg-rr73.5%
frac-timesN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
swap-sqrN/A
associate-*l*N/A
/-lowering-/.f64N/A
associate-*l*N/A
swap-sqrN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6465.7%
Applied egg-rr65.7%
herbie shell --seed 2024163
(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))))