
(FPCore (x l t) :precision binary64 (/ (* (sqrt 2.0) t) (sqrt (- (* (/ (+ x 1.0) (- x 1.0)) (+ (* l l) (* 2.0 (* t t)))) (* l l)))))
double code(double x, double l, double t) {
return (sqrt(2.0) * t) / sqrt(((((x + 1.0) / (x - 1.0)) * ((l * l) + (2.0 * (t * t)))) - (l * l)));
}
real(8) function code(x, l, t)
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t
code = (sqrt(2.0d0) * t) / sqrt(((((x + 1.0d0) / (x - 1.0d0)) * ((l * l) + (2.0d0 * (t * t)))) - (l * l)))
end function
public static double code(double x, double l, double t) {
return (Math.sqrt(2.0) * t) / Math.sqrt(((((x + 1.0) / (x - 1.0)) * ((l * l) + (2.0 * (t * t)))) - (l * l)));
}
def code(x, l, t): return (math.sqrt(2.0) * t) / math.sqrt(((((x + 1.0) / (x - 1.0)) * ((l * l) + (2.0 * (t * t)))) - (l * l)))
function code(x, l, t) return Float64(Float64(sqrt(2.0) * t) / sqrt(Float64(Float64(Float64(Float64(x + 1.0) / Float64(x - 1.0)) * Float64(Float64(l * l) + Float64(2.0 * Float64(t * t)))) - Float64(l * l)))) end
function tmp = code(x, l, t) tmp = (sqrt(2.0) * t) / sqrt(((((x + 1.0) / (x - 1.0)) * ((l * l) + (2.0 * (t * t)))) - (l * l))); end
code[x_, l_, t_] := N[(N[(N[Sqrt[2.0], $MachinePrecision] * t), $MachinePrecision] / N[Sqrt[N[(N[(N[(N[(x + 1.0), $MachinePrecision] / N[(x - 1.0), $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] + N[(2.0 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(l * l), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2} \cdot t}{\sqrt{\frac{x + 1}{x - 1} \cdot \left(\ell \cdot \ell + 2 \cdot \left(t \cdot t\right)\right) - \ell \cdot \ell}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 13 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x l t) :precision binary64 (/ (* (sqrt 2.0) t) (sqrt (- (* (/ (+ x 1.0) (- x 1.0)) (+ (* l l) (* 2.0 (* t t)))) (* l l)))))
double code(double x, double l, double t) {
return (sqrt(2.0) * t) / sqrt(((((x + 1.0) / (x - 1.0)) * ((l * l) + (2.0 * (t * t)))) - (l * l)));
}
real(8) function code(x, l, t)
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t
code = (sqrt(2.0d0) * t) / sqrt(((((x + 1.0d0) / (x - 1.0d0)) * ((l * l) + (2.0d0 * (t * t)))) - (l * l)))
end function
public static double code(double x, double l, double t) {
return (Math.sqrt(2.0) * t) / Math.sqrt(((((x + 1.0) / (x - 1.0)) * ((l * l) + (2.0 * (t * t)))) - (l * l)));
}
def code(x, l, t): return (math.sqrt(2.0) * t) / math.sqrt(((((x + 1.0) / (x - 1.0)) * ((l * l) + (2.0 * (t * t)))) - (l * l)))
function code(x, l, t) return Float64(Float64(sqrt(2.0) * t) / sqrt(Float64(Float64(Float64(Float64(x + 1.0) / Float64(x - 1.0)) * Float64(Float64(l * l) + Float64(2.0 * Float64(t * t)))) - Float64(l * l)))) end
function tmp = code(x, l, t) tmp = (sqrt(2.0) * t) / sqrt(((((x + 1.0) / (x - 1.0)) * ((l * l) + (2.0 * (t * t)))) - (l * l))); end
code[x_, l_, t_] := N[(N[(N[Sqrt[2.0], $MachinePrecision] * t), $MachinePrecision] / N[Sqrt[N[(N[(N[(N[(x + 1.0), $MachinePrecision] / N[(x - 1.0), $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] + N[(2.0 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(l * l), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2} \cdot t}{\sqrt{\frac{x + 1}{x - 1} \cdot \left(\ell \cdot \ell + 2 \cdot \left(t \cdot t\right)\right) - \ell \cdot \ell}}
\end{array}
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= t_m 2.4e-241)
(/ (* t_m (* (* (sqrt 0.5) (sqrt 2.0)) (sqrt x))) l_m)
(if (<= t_m 1.8e-165)
(- 1.0 (/ (- 1.0 (/ (+ 0.5 (/ -0.5 x)) x)) x))
(if (<= t_m 3.25e+18)
(*
t_m
(sqrt
(/
2.0
(+
(* l_m (* l_m (/ 2.0 x)))
(/ (* 2.0 (* (* t_m t_m) (+ x 1.0))) (+ x -1.0))))))
(pow (/ (+ x 1.0) (+ x -1.0)) -0.5))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 2.4e-241) {
tmp = (t_m * ((sqrt(0.5) * sqrt(2.0)) * sqrt(x))) / l_m;
} else if (t_m <= 1.8e-165) {
tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x);
} else if (t_m <= 3.25e+18) {
tmp = t_m * sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0)))));
} else {
tmp = pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
real(8) :: tmp
if (t_m <= 2.4d-241) then
tmp = (t_m * ((sqrt(0.5d0) * sqrt(2.0d0)) * sqrt(x))) / l_m
else if (t_m <= 1.8d-165) then
tmp = 1.0d0 - ((1.0d0 - ((0.5d0 + ((-0.5d0) / x)) / x)) / x)
else if (t_m <= 3.25d+18) then
tmp = t_m * sqrt((2.0d0 / ((l_m * (l_m * (2.0d0 / x))) + ((2.0d0 * ((t_m * t_m) * (x + 1.0d0))) / (x + (-1.0d0))))))
else
tmp = ((x + 1.0d0) / (x + (-1.0d0))) ** (-0.5d0)
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 2.4e-241) {
tmp = (t_m * ((Math.sqrt(0.5) * Math.sqrt(2.0)) * Math.sqrt(x))) / l_m;
} else if (t_m <= 1.8e-165) {
tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x);
} else if (t_m <= 3.25e+18) {
tmp = t_m * Math.sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0)))));
} else {
tmp = Math.pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): tmp = 0 if t_m <= 2.4e-241: tmp = (t_m * ((math.sqrt(0.5) * math.sqrt(2.0)) * math.sqrt(x))) / l_m elif t_m <= 1.8e-165: tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x) elif t_m <= 3.25e+18: tmp = t_m * math.sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0))))) else: tmp = math.pow(((x + 1.0) / (x + -1.0)), -0.5) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) tmp = 0.0 if (t_m <= 2.4e-241) tmp = Float64(Float64(t_m * Float64(Float64(sqrt(0.5) * sqrt(2.0)) * sqrt(x))) / l_m); elseif (t_m <= 1.8e-165) tmp = Float64(1.0 - Float64(Float64(1.0 - Float64(Float64(0.5 + Float64(-0.5 / x)) / x)) / x)); elseif (t_m <= 3.25e+18) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(Float64(l_m * Float64(l_m * Float64(2.0 / x))) + Float64(Float64(2.0 * Float64(Float64(t_m * t_m) * Float64(x + 1.0))) / Float64(x + -1.0)))))); else tmp = Float64(Float64(x + 1.0) / Float64(x + -1.0)) ^ -0.5; end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l_m, t_m) tmp = 0.0; if (t_m <= 2.4e-241) tmp = (t_m * ((sqrt(0.5) * sqrt(2.0)) * sqrt(x))) / l_m; elseif (t_m <= 1.8e-165) tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x); elseif (t_m <= 3.25e+18) tmp = t_m * sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0))))); else tmp = ((x + 1.0) / (x + -1.0)) ^ -0.5; end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 2.4e-241], N[(N[(t$95$m * N[(N[(N[Sqrt[0.5], $MachinePrecision] * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision] * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / l$95$m), $MachinePrecision], If[LessEqual[t$95$m, 1.8e-165], N[(1.0 - N[(N[(1.0 - N[(N[(0.5 + N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 3.25e+18], N[(t$95$m * N[Sqrt[N[(2.0 / N[(N[(l$95$m * N[(l$95$m * N[(2.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(2.0 * N[(N[(t$95$m * t$95$m), $MachinePrecision] * N[(x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(x + 1.0), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision], -0.5], $MachinePrecision]]]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
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^{-241}:\\
\;\;\;\;\frac{t\_m \cdot \left(\left(\sqrt{0.5} \cdot \sqrt{2}\right) \cdot \sqrt{x}\right)}{l\_m}\\
\mathbf{elif}\;t\_m \leq 1.8 \cdot 10^{-165}:\\
\;\;\;\;1 - \frac{1 - \frac{0.5 + \frac{-0.5}{x}}{x}}{x}\\
\mathbf{elif}\;t\_m \leq 3.25 \cdot 10^{+18}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{l\_m \cdot \left(l\_m \cdot \frac{2}{x}\right) + \frac{2 \cdot \left(\left(t\_m \cdot t\_m\right) \cdot \left(x + 1\right)\right)}{x + -1}}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{x + 1}{x + -1}\right)}^{-0.5}\\
\end{array}
\end{array}
if t < 2.4e-241Initial program 25.1%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr25.2%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified31.8%
Taylor expanded in x around inf
/-lowering-/.f6440.5%
Simplified40.5%
Taylor expanded in t around 0
associate-*l/N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f6419.5%
Simplified19.5%
if 2.4e-241 < t < 1.79999999999999992e-165Initial program 8.7%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6468.5%
Simplified68.5%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6468.5%
Applied egg-rr68.5%
Taylor expanded in x around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6468.5%
Simplified68.5%
if 1.79999999999999992e-165 < t < 3.25e18Initial program 58.2%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr58.4%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified75.2%
Taylor expanded in x around inf
/-lowering-/.f6489.1%
Simplified89.1%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6495.8%
Applied egg-rr95.8%
if 3.25e18 < t Initial program 34.7%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6498.7%
Simplified98.7%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6498.7%
Applied egg-rr98.7%
Final simplification56.3%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= t_m 1.1e-239)
(* t_m (/ (* (* (sqrt 0.5) (sqrt 2.0)) (sqrt x)) l_m))
(if (<= t_m 8e-165)
(- 1.0 (/ (- 1.0 (/ (+ 0.5 (/ -0.5 x)) x)) x))
(if (<= t_m 3.2e+18)
(*
t_m
(sqrt
(/
2.0
(+
(* l_m (* l_m (/ 2.0 x)))
(/ (* 2.0 (* (* t_m t_m) (+ x 1.0))) (+ x -1.0))))))
(pow (/ (+ x 1.0) (+ x -1.0)) -0.5))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 1.1e-239) {
tmp = t_m * (((sqrt(0.5) * sqrt(2.0)) * sqrt(x)) / l_m);
} else if (t_m <= 8e-165) {
tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x);
} else if (t_m <= 3.2e+18) {
tmp = t_m * sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0)))));
} else {
tmp = pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
real(8) :: tmp
if (t_m <= 1.1d-239) then
tmp = t_m * (((sqrt(0.5d0) * sqrt(2.0d0)) * sqrt(x)) / l_m)
else if (t_m <= 8d-165) then
tmp = 1.0d0 - ((1.0d0 - ((0.5d0 + ((-0.5d0) / x)) / x)) / x)
else if (t_m <= 3.2d+18) then
tmp = t_m * sqrt((2.0d0 / ((l_m * (l_m * (2.0d0 / x))) + ((2.0d0 * ((t_m * t_m) * (x + 1.0d0))) / (x + (-1.0d0))))))
else
tmp = ((x + 1.0d0) / (x + (-1.0d0))) ** (-0.5d0)
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 1.1e-239) {
tmp = t_m * (((Math.sqrt(0.5) * Math.sqrt(2.0)) * Math.sqrt(x)) / l_m);
} else if (t_m <= 8e-165) {
tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x);
} else if (t_m <= 3.2e+18) {
tmp = t_m * Math.sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0)))));
} else {
tmp = Math.pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): tmp = 0 if t_m <= 1.1e-239: tmp = t_m * (((math.sqrt(0.5) * math.sqrt(2.0)) * math.sqrt(x)) / l_m) elif t_m <= 8e-165: tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x) elif t_m <= 3.2e+18: tmp = t_m * math.sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0))))) else: tmp = math.pow(((x + 1.0) / (x + -1.0)), -0.5) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) tmp = 0.0 if (t_m <= 1.1e-239) tmp = Float64(t_m * Float64(Float64(Float64(sqrt(0.5) * sqrt(2.0)) * sqrt(x)) / l_m)); elseif (t_m <= 8e-165) tmp = Float64(1.0 - Float64(Float64(1.0 - Float64(Float64(0.5 + Float64(-0.5 / x)) / x)) / x)); elseif (t_m <= 3.2e+18) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(Float64(l_m * Float64(l_m * Float64(2.0 / x))) + Float64(Float64(2.0 * Float64(Float64(t_m * t_m) * Float64(x + 1.0))) / Float64(x + -1.0)))))); else tmp = Float64(Float64(x + 1.0) / Float64(x + -1.0)) ^ -0.5; end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l_m, t_m) tmp = 0.0; if (t_m <= 1.1e-239) tmp = t_m * (((sqrt(0.5) * sqrt(2.0)) * sqrt(x)) / l_m); elseif (t_m <= 8e-165) tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x); elseif (t_m <= 3.2e+18) tmp = t_m * sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0))))); else tmp = ((x + 1.0) / (x + -1.0)) ^ -0.5; end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 1.1e-239], N[(t$95$m * N[(N[(N[(N[Sqrt[0.5], $MachinePrecision] * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision] * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] / l$95$m), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 8e-165], N[(1.0 - N[(N[(1.0 - N[(N[(0.5 + N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 3.2e+18], N[(t$95$m * N[Sqrt[N[(2.0 / N[(N[(l$95$m * N[(l$95$m * N[(2.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(2.0 * N[(N[(t$95$m * t$95$m), $MachinePrecision] * N[(x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(x + 1.0), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision], -0.5], $MachinePrecision]]]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.1 \cdot 10^{-239}:\\
\;\;\;\;t\_m \cdot \frac{\left(\sqrt{0.5} \cdot \sqrt{2}\right) \cdot \sqrt{x}}{l\_m}\\
\mathbf{elif}\;t\_m \leq 8 \cdot 10^{-165}:\\
\;\;\;\;1 - \frac{1 - \frac{0.5 + \frac{-0.5}{x}}{x}}{x}\\
\mathbf{elif}\;t\_m \leq 3.2 \cdot 10^{+18}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{l\_m \cdot \left(l\_m \cdot \frac{2}{x}\right) + \frac{2 \cdot \left(\left(t\_m \cdot t\_m\right) \cdot \left(x + 1\right)\right)}{x + -1}}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{x + 1}{x + -1}\right)}^{-0.5}\\
\end{array}
\end{array}
if t < 1.09999999999999991e-239Initial program 25.1%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr25.2%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified31.8%
Taylor expanded in x around inf
/-lowering-/.f6440.5%
Simplified40.5%
Taylor expanded in l around inf
associate-*l/N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f6419.5%
Simplified19.5%
if 1.09999999999999991e-239 < t < 8.0000000000000001e-165Initial program 8.7%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6468.5%
Simplified68.5%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6468.5%
Applied egg-rr68.5%
Taylor expanded in x around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6468.5%
Simplified68.5%
if 8.0000000000000001e-165 < t < 3.2e18Initial program 58.2%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr58.4%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified75.2%
Taylor expanded in x around inf
/-lowering-/.f6489.1%
Simplified89.1%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6495.8%
Applied egg-rr95.8%
if 3.2e18 < t Initial program 34.7%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6498.7%
Simplified98.7%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6498.7%
Applied egg-rr98.7%
Final simplification56.3%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l_m t_m)
:precision binary64
(let* ((t_2 (* 2.0 (* t_m t_m))))
(*
t_s
(if (<= t_m 8e-288)
(* t_m (sqrt (/ x (* l_m l_m))))
(if (<= t_m 2.8e-165)
(- 1.0 (/ (- 1.0 (/ (+ 0.5 (/ -0.5 x)) x)) x))
(if (<= t_m 4.6e+57)
(* t_m (sqrt (/ 2.0 (- t_2 (/ (* -2.0 (+ (* l_m l_m) t_2)) x)))))
(pow (/ (+ x 1.0) (+ x -1.0)) -0.5)))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
double t_2 = 2.0 * (t_m * t_m);
double tmp;
if (t_m <= 8e-288) {
tmp = t_m * sqrt((x / (l_m * l_m)));
} else if (t_m <= 2.8e-165) {
tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x);
} else if (t_m <= 4.6e+57) {
tmp = t_m * sqrt((2.0 / (t_2 - ((-2.0 * ((l_m * l_m) + t_2)) / x))));
} else {
tmp = pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
real(8) :: t_2
real(8) :: tmp
t_2 = 2.0d0 * (t_m * t_m)
if (t_m <= 8d-288) then
tmp = t_m * sqrt((x / (l_m * l_m)))
else if (t_m <= 2.8d-165) then
tmp = 1.0d0 - ((1.0d0 - ((0.5d0 + ((-0.5d0) / x)) / x)) / x)
else if (t_m <= 4.6d+57) then
tmp = t_m * sqrt((2.0d0 / (t_2 - (((-2.0d0) * ((l_m * l_m) + t_2)) / x))))
else
tmp = ((x + 1.0d0) / (x + (-1.0d0))) ** (-0.5d0)
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
double t_2 = 2.0 * (t_m * t_m);
double tmp;
if (t_m <= 8e-288) {
tmp = t_m * Math.sqrt((x / (l_m * l_m)));
} else if (t_m <= 2.8e-165) {
tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x);
} else if (t_m <= 4.6e+57) {
tmp = t_m * Math.sqrt((2.0 / (t_2 - ((-2.0 * ((l_m * l_m) + t_2)) / x))));
} else {
tmp = Math.pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): t_2 = 2.0 * (t_m * t_m) tmp = 0 if t_m <= 8e-288: tmp = t_m * math.sqrt((x / (l_m * l_m))) elif t_m <= 2.8e-165: tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x) elif t_m <= 4.6e+57: tmp = t_m * math.sqrt((2.0 / (t_2 - ((-2.0 * ((l_m * l_m) + t_2)) / x)))) else: tmp = math.pow(((x + 1.0) / (x + -1.0)), -0.5) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) t_2 = Float64(2.0 * Float64(t_m * t_m)) tmp = 0.0 if (t_m <= 8e-288) tmp = Float64(t_m * sqrt(Float64(x / Float64(l_m * l_m)))); elseif (t_m <= 2.8e-165) tmp = Float64(1.0 - Float64(Float64(1.0 - Float64(Float64(0.5 + Float64(-0.5 / x)) / x)) / x)); elseif (t_m <= 4.6e+57) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(t_2 - Float64(Float64(-2.0 * Float64(Float64(l_m * l_m) + t_2)) / x))))); else tmp = Float64(Float64(x + 1.0) / Float64(x + -1.0)) ^ -0.5; end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l_m, t_m) t_2 = 2.0 * (t_m * t_m); tmp = 0.0; if (t_m <= 8e-288) tmp = t_m * sqrt((x / (l_m * l_m))); elseif (t_m <= 2.8e-165) tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x); elseif (t_m <= 4.6e+57) tmp = t_m * sqrt((2.0 / (t_2 - ((-2.0 * ((l_m * l_m) + t_2)) / x)))); else tmp = ((x + 1.0) / (x + -1.0)) ^ -0.5; end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := Block[{t$95$2 = N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 8e-288], N[(t$95$m * N[Sqrt[N[(x / N[(l$95$m * l$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 2.8e-165], N[(1.0 - N[(N[(1.0 - N[(N[(0.5 + N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 4.6e+57], N[(t$95$m * N[Sqrt[N[(2.0 / N[(t$95$2 - N[(N[(-2.0 * N[(N[(l$95$m * l$95$m), $MachinePrecision] + t$95$2), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(x + 1.0), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision], -0.5], $MachinePrecision]]]]), $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := 2 \cdot \left(t\_m \cdot t\_m\right)\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 8 \cdot 10^{-288}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{x}{l\_m \cdot l\_m}}\\
\mathbf{elif}\;t\_m \leq 2.8 \cdot 10^{-165}:\\
\;\;\;\;1 - \frac{1 - \frac{0.5 + \frac{-0.5}{x}}{x}}{x}\\
\mathbf{elif}\;t\_m \leq 4.6 \cdot 10^{+57}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{t\_2 - \frac{-2 \cdot \left(l\_m \cdot l\_m + t\_2\right)}{x}}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{x + 1}{x + -1}\right)}^{-0.5}\\
\end{array}
\end{array}
\end{array}
if t < 8.00000000000000046e-288Initial program 26.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr26.5%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified33.0%
Taylor expanded in x around inf
/-lowering-/.f6441.0%
Simplified41.0%
Taylor expanded in l around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6418.0%
Simplified18.0%
if 8.00000000000000046e-288 < t < 2.7999999999999999e-165Initial program 6.3%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6455.0%
Simplified55.0%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6455.0%
Applied egg-rr55.0%
Taylor expanded in x around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6455.0%
Simplified55.0%
if 2.7999999999999999e-165 < t < 4.5999999999999998e57Initial program 62.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr62.7%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified74.3%
Taylor expanded in x around -inf
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified89.2%
if 4.5999999999999998e57 < t Initial program 27.4%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64100.0%
Simplified100.0%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64100.0%
Applied egg-rr100.0%
Final simplification54.8%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= t_m 2.35e-287)
(* t_m (sqrt (/ x (* l_m l_m))))
(if (<= t_m 8e-165)
(- 1.0 (/ (- 1.0 (/ (+ 0.5 (/ -0.5 x)) x)) x))
(if (<= t_m 2.5e+57)
(*
t_m
(sqrt (/ 2.0 (+ (* 2.0 (* t_m t_m)) (* (/ 2.0 x) (* l_m l_m))))))
(pow (/ (+ x 1.0) (+ x -1.0)) -0.5))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 2.35e-287) {
tmp = t_m * sqrt((x / (l_m * l_m)));
} else if (t_m <= 8e-165) {
tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x);
} else if (t_m <= 2.5e+57) {
tmp = t_m * sqrt((2.0 / ((2.0 * (t_m * t_m)) + ((2.0 / x) * (l_m * l_m)))));
} else {
tmp = pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
real(8) :: tmp
if (t_m <= 2.35d-287) then
tmp = t_m * sqrt((x / (l_m * l_m)))
else if (t_m <= 8d-165) then
tmp = 1.0d0 - ((1.0d0 - ((0.5d0 + ((-0.5d0) / x)) / x)) / x)
else if (t_m <= 2.5d+57) then
tmp = t_m * sqrt((2.0d0 / ((2.0d0 * (t_m * t_m)) + ((2.0d0 / x) * (l_m * l_m)))))
else
tmp = ((x + 1.0d0) / (x + (-1.0d0))) ** (-0.5d0)
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 2.35e-287) {
tmp = t_m * Math.sqrt((x / (l_m * l_m)));
} else if (t_m <= 8e-165) {
tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x);
} else if (t_m <= 2.5e+57) {
tmp = t_m * Math.sqrt((2.0 / ((2.0 * (t_m * t_m)) + ((2.0 / x) * (l_m * l_m)))));
} else {
tmp = Math.pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): tmp = 0 if t_m <= 2.35e-287: tmp = t_m * math.sqrt((x / (l_m * l_m))) elif t_m <= 8e-165: tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x) elif t_m <= 2.5e+57: tmp = t_m * math.sqrt((2.0 / ((2.0 * (t_m * t_m)) + ((2.0 / x) * (l_m * l_m))))) else: tmp = math.pow(((x + 1.0) / (x + -1.0)), -0.5) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) tmp = 0.0 if (t_m <= 2.35e-287) tmp = Float64(t_m * sqrt(Float64(x / Float64(l_m * l_m)))); elseif (t_m <= 8e-165) tmp = Float64(1.0 - Float64(Float64(1.0 - Float64(Float64(0.5 + Float64(-0.5 / x)) / x)) / x)); elseif (t_m <= 2.5e+57) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(Float64(2.0 * Float64(t_m * t_m)) + Float64(Float64(2.0 / x) * Float64(l_m * l_m)))))); else tmp = Float64(Float64(x + 1.0) / Float64(x + -1.0)) ^ -0.5; end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l_m, t_m) tmp = 0.0; if (t_m <= 2.35e-287) tmp = t_m * sqrt((x / (l_m * l_m))); elseif (t_m <= 8e-165) tmp = 1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x); elseif (t_m <= 2.5e+57) tmp = t_m * sqrt((2.0 / ((2.0 * (t_m * t_m)) + ((2.0 / x) * (l_m * l_m))))); else tmp = ((x + 1.0) / (x + -1.0)) ^ -0.5; end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 2.35e-287], N[(t$95$m * N[Sqrt[N[(x / N[(l$95$m * l$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 8e-165], N[(1.0 - N[(N[(1.0 - N[(N[(0.5 + N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 2.5e+57], N[(t$95$m * N[Sqrt[N[(2.0 / N[(N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] + N[(N[(2.0 / x), $MachinePrecision] * N[(l$95$m * l$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(x + 1.0), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision], -0.5], $MachinePrecision]]]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
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.35 \cdot 10^{-287}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{x}{l\_m \cdot l\_m}}\\
\mathbf{elif}\;t\_m \leq 8 \cdot 10^{-165}:\\
\;\;\;\;1 - \frac{1 - \frac{0.5 + \frac{-0.5}{x}}{x}}{x}\\
\mathbf{elif}\;t\_m \leq 2.5 \cdot 10^{+57}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{2 \cdot \left(t\_m \cdot t\_m\right) + \frac{2}{x} \cdot \left(l\_m \cdot l\_m\right)}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{x + 1}{x + -1}\right)}^{-0.5}\\
\end{array}
\end{array}
if t < 2.3499999999999999e-287Initial program 26.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr26.5%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified33.0%
Taylor expanded in x around inf
/-lowering-/.f6441.0%
Simplified41.0%
Taylor expanded in l around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6418.0%
Simplified18.0%
if 2.3499999999999999e-287 < t < 8.0000000000000001e-165Initial program 6.3%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6455.0%
Simplified55.0%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6455.0%
Applied egg-rr55.0%
Taylor expanded in x around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6455.0%
Simplified55.0%
if 8.0000000000000001e-165 < t < 2.49999999999999986e57Initial program 62.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr62.7%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified74.3%
Taylor expanded in x around inf
/-lowering-/.f6485.4%
Simplified85.4%
Taylor expanded in x around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6489.2%
Simplified89.2%
if 2.49999999999999986e57 < t Initial program 27.4%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64100.0%
Simplified100.0%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64100.0%
Applied egg-rr100.0%
Final simplification54.8%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= t_m 6.7e+18)
(*
t_m
(sqrt
(/
2.0
(+
(* l_m (* l_m (/ 2.0 x)))
(/ (* 2.0 (* (* t_m t_m) (+ x 1.0))) (+ x -1.0))))))
(pow (/ (+ x 1.0) (+ x -1.0)) -0.5))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 6.7e+18) {
tmp = t_m * sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0)))));
} else {
tmp = pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
real(8) :: tmp
if (t_m <= 6.7d+18) then
tmp = t_m * sqrt((2.0d0 / ((l_m * (l_m * (2.0d0 / x))) + ((2.0d0 * ((t_m * t_m) * (x + 1.0d0))) / (x + (-1.0d0))))))
else
tmp = ((x + 1.0d0) / (x + (-1.0d0))) ** (-0.5d0)
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 6.7e+18) {
tmp = t_m * Math.sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0)))));
} else {
tmp = Math.pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): tmp = 0 if t_m <= 6.7e+18: tmp = t_m * math.sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0))))) else: tmp = math.pow(((x + 1.0) / (x + -1.0)), -0.5) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) tmp = 0.0 if (t_m <= 6.7e+18) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(Float64(l_m * Float64(l_m * Float64(2.0 / x))) + Float64(Float64(2.0 * Float64(Float64(t_m * t_m) * Float64(x + 1.0))) / Float64(x + -1.0)))))); else tmp = Float64(Float64(x + 1.0) / Float64(x + -1.0)) ^ -0.5; end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l_m, t_m) tmp = 0.0; if (t_m <= 6.7e+18) tmp = t_m * sqrt((2.0 / ((l_m * (l_m * (2.0 / x))) + ((2.0 * ((t_m * t_m) * (x + 1.0))) / (x + -1.0))))); else tmp = ((x + 1.0) / (x + -1.0)) ^ -0.5; end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 6.7e+18], N[(t$95$m * N[Sqrt[N[(2.0 / N[(N[(l$95$m * N[(l$95$m * N[(2.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(2.0 * N[(N[(t$95$m * t$95$m), $MachinePrecision] * N[(x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(x + 1.0), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision], -0.5], $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 6.7 \cdot 10^{+18}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{l\_m \cdot \left(l\_m \cdot \frac{2}{x}\right) + \frac{2 \cdot \left(\left(t\_m \cdot t\_m\right) \cdot \left(x + 1\right)\right)}{x + -1}}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{x + 1}{x + -1}\right)}^{-0.5}\\
\end{array}
\end{array}
if t < 6.7e18Initial program 31.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr31.5%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified40.2%
Taylor expanded in x around inf
/-lowering-/.f6450.2%
Simplified50.2%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6454.6%
Applied egg-rr54.6%
if 6.7e18 < t Initial program 34.7%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6498.7%
Simplified98.7%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6498.7%
Applied egg-rr98.7%
Final simplification67.0%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= t_m 1.9e-286)
(* t_m (sqrt (/ x (* l_m l_m))))
(pow (/ (+ x 1.0) (+ x -1.0)) -0.5))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 1.9e-286) {
tmp = t_m * sqrt((x / (l_m * l_m)));
} else {
tmp = pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
real(8) :: tmp
if (t_m <= 1.9d-286) then
tmp = t_m * sqrt((x / (l_m * l_m)))
else
tmp = ((x + 1.0d0) / (x + (-1.0d0))) ** (-0.5d0)
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 1.9e-286) {
tmp = t_m * Math.sqrt((x / (l_m * l_m)));
} else {
tmp = Math.pow(((x + 1.0) / (x + -1.0)), -0.5);
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): tmp = 0 if t_m <= 1.9e-286: tmp = t_m * math.sqrt((x / (l_m * l_m))) else: tmp = math.pow(((x + 1.0) / (x + -1.0)), -0.5) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) tmp = 0.0 if (t_m <= 1.9e-286) tmp = Float64(t_m * sqrt(Float64(x / Float64(l_m * l_m)))); else tmp = Float64(Float64(x + 1.0) / Float64(x + -1.0)) ^ -0.5; end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l_m, t_m) tmp = 0.0; if (t_m <= 1.9e-286) tmp = t_m * sqrt((x / (l_m * l_m))); else tmp = ((x + 1.0) / (x + -1.0)) ^ -0.5; end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 1.9e-286], N[(t$95$m * N[Sqrt[N[(x / N[(l$95$m * l$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(x + 1.0), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision], -0.5], $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.9 \cdot 10^{-286}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{x}{l\_m \cdot l\_m}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{x + 1}{x + -1}\right)}^{-0.5}\\
\end{array}
\end{array}
if t < 1.9000000000000001e-286Initial program 26.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr26.5%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified33.0%
Taylor expanded in x around inf
/-lowering-/.f6441.0%
Simplified41.0%
Taylor expanded in l around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6418.0%
Simplified18.0%
if 1.9000000000000001e-286 < t Initial program 37.9%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6482.8%
Simplified82.8%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6482.8%
Applied egg-rr82.8%
Final simplification51.4%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= t_m 2.7e-286)
(* t_m (sqrt (/ x (* l_m l_m))))
(sqrt (/ (+ x -1.0) (+ x 1.0))))))l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 2.7e-286) {
tmp = t_m * sqrt((x / (l_m * l_m)));
} else {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
real(8) :: tmp
if (t_m <= 2.7d-286) then
tmp = t_m * sqrt((x / (l_m * l_m)))
else
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
end if
code = t_s * tmp
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
double tmp;
if (t_m <= 2.7e-286) {
tmp = t_m * Math.sqrt((x / (l_m * l_m)));
} else {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): tmp = 0 if t_m <= 2.7e-286: tmp = t_m * math.sqrt((x / (l_m * l_m))) else: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) return t_s * tmp
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) tmp = 0.0 if (t_m <= 2.7e-286) tmp = Float64(t_m * sqrt(Float64(x / Float64(l_m * l_m)))); else tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); end return Float64(t_s * tmp) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l_m, t_m) tmp = 0.0; if (t_m <= 2.7e-286) tmp = t_m * sqrt((x / (l_m * l_m))); else tmp = sqrt(((x + -1.0) / (x + 1.0))); end tmp_2 = t_s * tmp; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 2.7e-286], N[(t$95$m * N[Sqrt[N[(x / N[(l$95$m * l$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
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.7 \cdot 10^{-286}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{x}{l\_m \cdot l\_m}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\end{array}
\end{array}
if t < 2.7000000000000002e-286Initial program 26.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
Applied egg-rr26.5%
Taylor expanded in l around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified33.0%
Taylor expanded in x around inf
/-lowering-/.f6441.0%
Simplified41.0%
Taylor expanded in l around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6418.0%
Simplified18.0%
if 2.7000000000000002e-286 < t Initial program 37.9%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6482.8%
Simplified82.8%
Taylor expanded in t around 0
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6482.8%
Simplified82.8%
l_m = (fabs.f64 l) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s x l_m t_m) :precision binary64 (* t_s (sqrt (/ (+ x -1.0) (+ x 1.0)))))
l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
return t_s * sqrt(((x + -1.0) / (x + 1.0)));
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
code = t_s * sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
return t_s * Math.sqrt(((x + -1.0) / (x + 1.0)));
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): return t_s * math.sqrt(((x + -1.0) / (x + 1.0)))
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) return Float64(t_s * sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0)))) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, x, l_m, t_m) tmp = t_s * sqrt(((x + -1.0) / (x + 1.0))); end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \sqrt{\frac{x + -1}{x + 1}}
\end{array}
Initial program 32.3%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6443.6%
Simplified43.6%
Taylor expanded in t around 0
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6443.6%
Simplified43.6%
l_m = (fabs.f64 l) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s x l_m t_m) :precision binary64 (* t_s (- 1.0 (/ (- 1.0 (/ (+ 0.5 (/ -0.5 x)) x)) x))))
l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
return t_s * (1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x));
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
code = t_s * (1.0d0 - ((1.0d0 - ((0.5d0 + ((-0.5d0) / x)) / x)) / x))
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
return t_s * (1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x));
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): return t_s * (1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x))
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) return Float64(t_s * Float64(1.0 - Float64(Float64(1.0 - Float64(Float64(0.5 + Float64(-0.5 / x)) / x)) / x))) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, x, l_m, t_m) tmp = t_s * (1.0 - ((1.0 - ((0.5 + (-0.5 / x)) / x)) / x)); end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * N[(1.0 - N[(N[(1.0 - N[(N[(0.5 + N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(1 - \frac{1 - \frac{0.5 + \frac{-0.5}{x}}{x}}{x}\right)
\end{array}
Initial program 32.3%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6443.6%
Simplified43.6%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6443.6%
Applied egg-rr43.6%
Taylor expanded in x around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6443.3%
Simplified43.3%
l_m = (fabs.f64 l) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s x l_m t_m) :precision binary64 (* t_s (+ (+ 1.0 (/ 0.5 (* x x))) (/ -1.0 x))))
l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
return t_s * ((1.0 + (0.5 / (x * x))) + (-1.0 / x));
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
code = t_s * ((1.0d0 + (0.5d0 / (x * x))) + ((-1.0d0) / x))
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
return t_s * ((1.0 + (0.5 / (x * x))) + (-1.0 / x));
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): return t_s * ((1.0 + (0.5 / (x * x))) + (-1.0 / x))
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) return Float64(t_s * Float64(Float64(1.0 + Float64(0.5 / Float64(x * x))) + Float64(-1.0 / x))) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, x, l_m, t_m) tmp = t_s * ((1.0 + (0.5 / (x * x))) + (-1.0 / x)); end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * N[(N[(1.0 + N[(0.5 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(1 + \frac{0.5}{x \cdot x}\right) + \frac{-1}{x}\right)
\end{array}
Initial program 32.3%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6443.6%
Simplified43.6%
Taylor expanded in x around inf
--lowering--.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f6443.2%
Simplified43.2%
Final simplification43.2%
l_m = (fabs.f64 l) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s x l_m t_m) :precision binary64 (* t_s (+ 1.0 (/ (- -1.0 (/ -0.5 x)) x))))
l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
return t_s * (1.0 + ((-1.0 - (-0.5 / x)) / x));
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
code = t_s * (1.0d0 + (((-1.0d0) - ((-0.5d0) / x)) / x))
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
return t_s * (1.0 + ((-1.0 - (-0.5 / x)) / x));
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): return t_s * (1.0 + ((-1.0 - (-0.5 / x)) / x))
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) return Float64(t_s * Float64(1.0 + Float64(Float64(-1.0 - Float64(-0.5 / x)) / x))) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, x, l_m, t_m) tmp = t_s * (1.0 + ((-1.0 - (-0.5 / x)) / x)); end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * N[(1.0 + N[(N[(-1.0 - N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(1 + \frac{-1 - \frac{-0.5}{x}}{x}\right)
\end{array}
Initial program 32.3%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6443.6%
Simplified43.6%
*-commutativeN/A
associate-/r*N/A
*-inversesN/A
pow1/2N/A
pow-flipN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6443.6%
Applied egg-rr43.6%
Taylor expanded in x around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6443.2%
Simplified43.2%
Final simplification43.2%
l_m = (fabs.f64 l) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s x l_m t_m) :precision binary64 (* t_s (+ 1.0 (/ -1.0 x))))
l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
return t_s * (1.0 + (-1.0 / x));
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
code = t_s * (1.0d0 + ((-1.0d0) / x))
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
return t_s * (1.0 + (-1.0 / x));
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): return t_s * (1.0 + (-1.0 / x))
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) return Float64(t_s * Float64(1.0 + Float64(-1.0 / x))) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, x, l_m, t_m) tmp = t_s * (1.0 + (-1.0 / x)); end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * N[(1.0 + N[(-1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(1 + \frac{-1}{x}\right)
\end{array}
Initial program 32.3%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6443.6%
Simplified43.6%
Taylor expanded in x around inf
--lowering--.f64N/A
/-lowering-/.f6442.9%
Simplified42.9%
Final simplification42.9%
l_m = (fabs.f64 l) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s x l_m t_m) :precision binary64 (* t_s 1.0))
l_m = fabs(l);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l_m, double t_m) {
return t_s * 1.0;
}
l_m = abs(l)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l_m, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l_m
real(8), intent (in) :: t_m
code = t_s * 1.0d0
end function
l_m = Math.abs(l);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l_m, double t_m) {
return t_s * 1.0;
}
l_m = math.fabs(l) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l_m, t_m): return t_s * 1.0
l_m = abs(l) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l_m, t_m) return Float64(t_s * 1.0) end
l_m = abs(l); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, x, l_m, t_m) tmp = t_s * 1.0; end
l_m = N[Abs[l], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l$95$m_, t$95$m_] := N[(t$95$s * 1.0), $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot 1
\end{array}
Initial program 32.3%
Taylor expanded in l around 0
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6443.6%
Simplified43.6%
Taylor expanded in x around inf
Simplified42.6%
herbie shell --seed 2024288
(FPCore (x l t)
:name "Toniolo and Linder, Equation (7)"
:precision binary64
(/ (* (sqrt 2.0) t) (sqrt (- (* (/ (+ x 1.0) (- x 1.0)) (+ (* l l) (* 2.0 (* t t)))) (* l l)))))