
(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 12 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.35e-228)
(/ (* (* t_m (sqrt 2.0)) (sqrt (/ x 2.0))) l_m)
(if (<= t_m 3.2e-159)
(sqrt (/ (+ x -1.0) (+ x 1.0)))
(if (<= t_m 8e-31)
(*
t_m
(sqrt
(/
2.0
(+
(+ (* 2.0 (+ (* t_m t_m) (/ (* t_m t_m) x))) (/ (* l_m l_m) x))
(/ (+ (* l_m l_m) (* 2.0 (* t_m t_m))) x)))))
(+ 1.0 (/ (- -1.0 (/ (- (/ 0.5 x) 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) {
double tmp;
if (t_m <= 2.35e-228) {
tmp = ((t_m * sqrt(2.0)) * sqrt((x / 2.0))) / l_m;
} else if (t_m <= 3.2e-159) {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
} else if (t_m <= 8e-31) {
tmp = t_m * sqrt((2.0 / (((2.0 * ((t_m * t_m) + ((t_m * t_m) / x))) + ((l_m * l_m) / x)) + (((l_m * l_m) + (2.0 * (t_m * t_m))) / x))));
} else {
tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x);
}
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-228) then
tmp = ((t_m * sqrt(2.0d0)) * sqrt((x / 2.0d0))) / l_m
else if (t_m <= 3.2d-159) then
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
else if (t_m <= 8d-31) then
tmp = t_m * sqrt((2.0d0 / (((2.0d0 * ((t_m * t_m) + ((t_m * t_m) / x))) + ((l_m * l_m) / x)) + (((l_m * l_m) + (2.0d0 * (t_m * t_m))) / x))))
else
tmp = 1.0d0 + (((-1.0d0) - (((0.5d0 / x) - 0.5d0) / x)) / x)
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-228) {
tmp = ((t_m * Math.sqrt(2.0)) * Math.sqrt((x / 2.0))) / l_m;
} else if (t_m <= 3.2e-159) {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
} else if (t_m <= 8e-31) {
tmp = t_m * Math.sqrt((2.0 / (((2.0 * ((t_m * t_m) + ((t_m * t_m) / x))) + ((l_m * l_m) / x)) + (((l_m * l_m) + (2.0 * (t_m * t_m))) / x))));
} else {
tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x);
}
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-228: tmp = ((t_m * math.sqrt(2.0)) * math.sqrt((x / 2.0))) / l_m elif t_m <= 3.2e-159: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) elif t_m <= 8e-31: tmp = t_m * math.sqrt((2.0 / (((2.0 * ((t_m * t_m) + ((t_m * t_m) / x))) + ((l_m * l_m) / x)) + (((l_m * l_m) + (2.0 * (t_m * t_m))) / x)))) else: tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x) 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-228) tmp = Float64(Float64(Float64(t_m * sqrt(2.0)) * sqrt(Float64(x / 2.0))) / l_m); elseif (t_m <= 3.2e-159) tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); elseif (t_m <= 8e-31) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(Float64(Float64(2.0 * Float64(Float64(t_m * t_m) + Float64(Float64(t_m * t_m) / x))) + Float64(Float64(l_m * l_m) / x)) + Float64(Float64(Float64(l_m * l_m) + Float64(2.0 * Float64(t_m * t_m))) / x))))); else tmp = Float64(1.0 + Float64(Float64(-1.0 - Float64(Float64(Float64(0.5 / x) - 0.5) / x)) / x)); 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-228) tmp = ((t_m * sqrt(2.0)) * sqrt((x / 2.0))) / l_m; elseif (t_m <= 3.2e-159) tmp = sqrt(((x + -1.0) / (x + 1.0))); elseif (t_m <= 8e-31) tmp = t_m * sqrt((2.0 / (((2.0 * ((t_m * t_m) + ((t_m * t_m) / x))) + ((l_m * l_m) / x)) + (((l_m * l_m) + (2.0 * (t_m * t_m))) / x)))); else tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x); 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-228], N[(N[(N[(t$95$m * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(x / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / l$95$m), $MachinePrecision], If[LessEqual[t$95$m, 3.2e-159], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[t$95$m, 8e-31], N[(t$95$m * N[Sqrt[N[(2.0 / N[(N[(N[(2.0 * N[(N[(t$95$m * t$95$m), $MachinePrecision] + N[(N[(t$95$m * t$95$m), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(l$95$m * l$95$m), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(l$95$m * l$95$m), $MachinePrecision] + N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(N[(-1.0 - N[(N[(N[(0.5 / x), $MachinePrecision] - 0.5), $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 \begin{array}{l}
\mathbf{if}\;t\_m \leq 2.35 \cdot 10^{-228}:\\
\;\;\;\;\frac{\left(t\_m \cdot \sqrt{2}\right) \cdot \sqrt{\frac{x}{2}}}{l\_m}\\
\mathbf{elif}\;t\_m \leq 3.2 \cdot 10^{-159}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\mathbf{elif}\;t\_m \leq 8 \cdot 10^{-31}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{\left(2 \cdot \left(t\_m \cdot t\_m + \frac{t\_m \cdot t\_m}{x}\right) + \frac{l\_m \cdot l\_m}{x}\right) + \frac{l\_m \cdot l\_m + 2 \cdot \left(t\_m \cdot t\_m\right)}{x}}}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1 - \frac{\frac{0.5}{x} - 0.5}{x}}{x}\\
\end{array}
\end{array}
if t < 2.3500000000000001e-228Initial program 32.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-rr32.3%
Taylor expanded in l around inf
associate-/r*N/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
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6412.8%
Simplified12.8%
Taylor expanded in t around 0
associate-*l/N/A
/-lowering-/.f64N/A
Simplified12.7%
Taylor expanded in x around inf
*-commutativeN/A
metadata-evalN/A
associate-/l*N/A
*-commutativeN/A
mul-1-negN/A
distribute-frac-negN/A
distribute-neg-frac2N/A
metadata-evalN/A
/-lowering-/.f6417.3%
Simplified17.3%
if 2.3500000000000001e-228 < t < 3.2e-159Initial program 2.5%
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-+.f6475.9%
Simplified75.9%
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-+.f6475.9%
Simplified75.9%
if 3.2e-159 < t < 8.000000000000001e-31Initial program 50.6%
*-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-rr50.9%
Taylor expanded in x around inf
sub-negN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-lowering-+.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
Simplified90.2%
if 8.000000000000001e-31 < t Initial program 34.6%
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-+.f6491.8%
Simplified91.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-+.f6491.8%
Simplified91.8%
Taylor expanded in x around inf
associate--r+N/A
sub-negN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.8%
Simplified91.8%
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
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6491.8%
Simplified91.8%
Final simplification48.7%
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.7e-226)
(/ 1.0 (/ l_m (* t_m (sqrt x))))
(if (<= t_m 2.6e-160)
(sqrt (/ (+ x -1.0) (+ x 1.0)))
(if (<= t_m 8e-31)
(*
t_m
(sqrt
(/
2.0
(+
(+ (* 2.0 (+ (* t_m t_m) (/ (* t_m t_m) x))) (/ (* l_m l_m) x))
(/ (+ (* l_m l_m) (* 2.0 (* t_m t_m))) x)))))
(+ 1.0 (/ (- -1.0 (/ (- (/ 0.5 x) 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) {
double tmp;
if (t_m <= 1.7e-226) {
tmp = 1.0 / (l_m / (t_m * sqrt(x)));
} else if (t_m <= 2.6e-160) {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
} else if (t_m <= 8e-31) {
tmp = t_m * sqrt((2.0 / (((2.0 * ((t_m * t_m) + ((t_m * t_m) / x))) + ((l_m * l_m) / x)) + (((l_m * l_m) + (2.0 * (t_m * t_m))) / x))));
} else {
tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x);
}
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.7d-226) then
tmp = 1.0d0 / (l_m / (t_m * sqrt(x)))
else if (t_m <= 2.6d-160) then
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
else if (t_m <= 8d-31) then
tmp = t_m * sqrt((2.0d0 / (((2.0d0 * ((t_m * t_m) + ((t_m * t_m) / x))) + ((l_m * l_m) / x)) + (((l_m * l_m) + (2.0d0 * (t_m * t_m))) / x))))
else
tmp = 1.0d0 + (((-1.0d0) - (((0.5d0 / x) - 0.5d0) / x)) / x)
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.7e-226) {
tmp = 1.0 / (l_m / (t_m * Math.sqrt(x)));
} else if (t_m <= 2.6e-160) {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
} else if (t_m <= 8e-31) {
tmp = t_m * Math.sqrt((2.0 / (((2.0 * ((t_m * t_m) + ((t_m * t_m) / x))) + ((l_m * l_m) / x)) + (((l_m * l_m) + (2.0 * (t_m * t_m))) / x))));
} else {
tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x);
}
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.7e-226: tmp = 1.0 / (l_m / (t_m * math.sqrt(x))) elif t_m <= 2.6e-160: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) elif t_m <= 8e-31: tmp = t_m * math.sqrt((2.0 / (((2.0 * ((t_m * t_m) + ((t_m * t_m) / x))) + ((l_m * l_m) / x)) + (((l_m * l_m) + (2.0 * (t_m * t_m))) / x)))) else: tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x) 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.7e-226) tmp = Float64(1.0 / Float64(l_m / Float64(t_m * sqrt(x)))); elseif (t_m <= 2.6e-160) tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); elseif (t_m <= 8e-31) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(Float64(Float64(2.0 * Float64(Float64(t_m * t_m) + Float64(Float64(t_m * t_m) / x))) + Float64(Float64(l_m * l_m) / x)) + Float64(Float64(Float64(l_m * l_m) + Float64(2.0 * Float64(t_m * t_m))) / x))))); else tmp = Float64(1.0 + Float64(Float64(-1.0 - Float64(Float64(Float64(0.5 / x) - 0.5) / x)) / x)); 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.7e-226) tmp = 1.0 / (l_m / (t_m * sqrt(x))); elseif (t_m <= 2.6e-160) tmp = sqrt(((x + -1.0) / (x + 1.0))); elseif (t_m <= 8e-31) tmp = t_m * sqrt((2.0 / (((2.0 * ((t_m * t_m) + ((t_m * t_m) / x))) + ((l_m * l_m) / x)) + (((l_m * l_m) + (2.0 * (t_m * t_m))) / x)))); else tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x); 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.7e-226], N[(1.0 / N[(l$95$m / N[(t$95$m * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 2.6e-160], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[t$95$m, 8e-31], N[(t$95$m * N[Sqrt[N[(2.0 / N[(N[(N[(2.0 * N[(N[(t$95$m * t$95$m), $MachinePrecision] + N[(N[(t$95$m * t$95$m), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(l$95$m * l$95$m), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(l$95$m * l$95$m), $MachinePrecision] + N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(N[(-1.0 - N[(N[(N[(0.5 / x), $MachinePrecision] - 0.5), $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 \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.7 \cdot 10^{-226}:\\
\;\;\;\;\frac{1}{\frac{l\_m}{t\_m \cdot \sqrt{x}}}\\
\mathbf{elif}\;t\_m \leq 2.6 \cdot 10^{-160}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\mathbf{elif}\;t\_m \leq 8 \cdot 10^{-31}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{\left(2 \cdot \left(t\_m \cdot t\_m + \frac{t\_m \cdot t\_m}{x}\right) + \frac{l\_m \cdot l\_m}{x}\right) + \frac{l\_m \cdot l\_m + 2 \cdot \left(t\_m \cdot t\_m\right)}{x}}}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1 - \frac{\frac{0.5}{x} - 0.5}{x}}{x}\\
\end{array}
\end{array}
if t < 1.70000000000000004e-226Initial program 32.2%
Taylor expanded in l around inf
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-commutativeN/A
metadata-evalN/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-+.f642.9%
Simplified2.9%
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
pow1/2N/A
pow1/2N/A
pow-prod-downN/A
pow-lowering-pow.f64N/A
Applied egg-rr3.1%
Taylor expanded in x around inf
sqrt-lowering-sqrt.f6417.2%
Simplified17.2%
if 1.70000000000000004e-226 < t < 2.60000000000000003e-160Initial program 2.5%
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-+.f6475.9%
Simplified75.9%
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-+.f6475.9%
Simplified75.9%
if 2.60000000000000003e-160 < t < 8.000000000000001e-31Initial program 50.6%
*-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-rr50.9%
Taylor expanded in x around inf
sub-negN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-lowering-+.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
Simplified90.2%
if 8.000000000000001e-31 < t Initial program 34.6%
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-+.f6491.8%
Simplified91.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-+.f6491.8%
Simplified91.8%
Taylor expanded in x around inf
associate--r+N/A
sub-negN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.8%
Simplified91.8%
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
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6491.8%
Simplified91.8%
Final simplification48.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
(if (<= l_m 4e+168)
(sqrt (/ (+ x -1.0) (+ x 1.0)))
(/
1.0
(/
l_m
(*
t_m
(pow
(*
2.0
(/
1.0
(/
(+ (+ (/ 2.0 x) (/ 2.0 (* x (* x x)))) (+ 2.0 (/ 2.0 (* x x))))
x)))
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 (l_m <= 4e+168) {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
} else {
tmp = 1.0 / (l_m / (t_m * pow((2.0 * (1.0 / ((((2.0 / x) + (2.0 / (x * (x * x)))) + (2.0 + (2.0 / (x * x)))) / x))), 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 (l_m <= 4d+168) then
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
else
tmp = 1.0d0 / (l_m / (t_m * ((2.0d0 * (1.0d0 / ((((2.0d0 / x) + (2.0d0 / (x * (x * x)))) + (2.0d0 + (2.0d0 / (x * x)))) / x))) ** 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 (l_m <= 4e+168) {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
} else {
tmp = 1.0 / (l_m / (t_m * Math.pow((2.0 * (1.0 / ((((2.0 / x) + (2.0 / (x * (x * x)))) + (2.0 + (2.0 / (x * x)))) / x))), 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 l_m <= 4e+168: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) else: tmp = 1.0 / (l_m / (t_m * math.pow((2.0 * (1.0 / ((((2.0 / x) + (2.0 / (x * (x * x)))) + (2.0 + (2.0 / (x * x)))) / x))), 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 (l_m <= 4e+168) tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); else tmp = Float64(1.0 / Float64(l_m / Float64(t_m * (Float64(2.0 * Float64(1.0 / Float64(Float64(Float64(Float64(2.0 / x) + Float64(2.0 / Float64(x * Float64(x * x)))) + Float64(2.0 + Float64(2.0 / Float64(x * x)))) / x))) ^ 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 (l_m <= 4e+168) tmp = sqrt(((x + -1.0) / (x + 1.0))); else tmp = 1.0 / (l_m / (t_m * ((2.0 * (1.0 / ((((2.0 / x) + (2.0 / (x * (x * x)))) + (2.0 + (2.0 / (x * x)))) / x))) ^ 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[l$95$m, 4e+168], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[(1.0 / N[(l$95$m / N[(t$95$m * N[Power[N[(2.0 * N[(1.0 / N[(N[(N[(N[(2.0 / x), $MachinePrecision] + N[(2.0 / N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(2.0 + N[(2.0 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision]), $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}\;l\_m \leq 4 \cdot 10^{+168}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{l\_m}{t\_m \cdot {\left(2 \cdot \frac{1}{\frac{\left(\frac{2}{x} + \frac{2}{x \cdot \left(x \cdot x\right)}\right) + \left(2 + \frac{2}{x \cdot x}\right)}{x}}\right)}^{0.5}}}\\
\end{array}
\end{array}
if l < 3.9999999999999997e168Initial program 36.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-+.f6441.0%
Simplified41.0%
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-+.f6441.0%
Simplified41.0%
if 3.9999999999999997e168 < l Initial program 0.0%
Taylor expanded in l around inf
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-commutativeN/A
metadata-evalN/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-+.f646.2%
Simplified6.2%
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
pow1/2N/A
pow1/2N/A
pow-prod-downN/A
pow-lowering-pow.f64N/A
Applied egg-rr6.7%
Taylor expanded in x around inf
/-lowering-/.f64N/A
Simplified59.0%
Final simplification42.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
(if (<= l_m 2.85e+168)
(sqrt (/ (+ x -1.0) (+ x 1.0)))
(/
1.0
(/ l_m (* t_m (pow (* 2.0 (/ 1.0 (/ (+ 2.0 (/ 2.0 x)) x))) 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 (l_m <= 2.85e+168) {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
} else {
tmp = 1.0 / (l_m / (t_m * pow((2.0 * (1.0 / ((2.0 + (2.0 / x)) / x))), 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 (l_m <= 2.85d+168) then
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
else
tmp = 1.0d0 / (l_m / (t_m * ((2.0d0 * (1.0d0 / ((2.0d0 + (2.0d0 / x)) / x))) ** 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 (l_m <= 2.85e+168) {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
} else {
tmp = 1.0 / (l_m / (t_m * Math.pow((2.0 * (1.0 / ((2.0 + (2.0 / x)) / x))), 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 l_m <= 2.85e+168: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) else: tmp = 1.0 / (l_m / (t_m * math.pow((2.0 * (1.0 / ((2.0 + (2.0 / x)) / x))), 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 (l_m <= 2.85e+168) tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); else tmp = Float64(1.0 / Float64(l_m / Float64(t_m * (Float64(2.0 * Float64(1.0 / Float64(Float64(2.0 + Float64(2.0 / x)) / x))) ^ 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 (l_m <= 2.85e+168) tmp = sqrt(((x + -1.0) / (x + 1.0))); else tmp = 1.0 / (l_m / (t_m * ((2.0 * (1.0 / ((2.0 + (2.0 / x)) / x))) ^ 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[l$95$m, 2.85e+168], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[(1.0 / N[(l$95$m / N[(t$95$m * N[Power[N[(2.0 * N[(1.0 / N[(N[(2.0 + N[(2.0 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision]), $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}\;l\_m \leq 2.85 \cdot 10^{+168}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{l\_m}{t\_m \cdot {\left(2 \cdot \frac{1}{\frac{2 + \frac{2}{x}}{x}}\right)}^{0.5}}}\\
\end{array}
\end{array}
if l < 2.85000000000000019e168Initial program 36.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-+.f6441.0%
Simplified41.0%
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-+.f6441.0%
Simplified41.0%
if 2.85000000000000019e168 < l Initial program 0.0%
Taylor expanded in l around inf
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-commutativeN/A
metadata-evalN/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-+.f646.2%
Simplified6.2%
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
pow1/2N/A
pow1/2N/A
pow-prod-downN/A
pow-lowering-pow.f64N/A
Applied egg-rr6.7%
Taylor expanded in x around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6457.5%
Simplified57.5%
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.8e-227)
(/ 1.0 (/ l_m (* t_m (sqrt x))))
(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.8e-227) {
tmp = 1.0 / (l_m / (t_m * sqrt(x)));
} 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.8d-227) then
tmp = 1.0d0 / (l_m / (t_m * sqrt(x)))
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.8e-227) {
tmp = 1.0 / (l_m / (t_m * Math.sqrt(x)));
} 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.8e-227: tmp = 1.0 / (l_m / (t_m * math.sqrt(x))) 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.8e-227) tmp = Float64(1.0 / Float64(l_m / Float64(t_m * sqrt(x)))); 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.8e-227) tmp = 1.0 / (l_m / (t_m * sqrt(x))); 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.8e-227], N[(1.0 / N[(l$95$m / N[(t$95$m * N[Sqrt[x], $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.8 \cdot 10^{-227}:\\
\;\;\;\;\frac{1}{\frac{l\_m}{t\_m \cdot \sqrt{x}}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\end{array}
\end{array}
if t < 2.7999999999999998e-227Initial program 32.2%
Taylor expanded in l around inf
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-commutativeN/A
metadata-evalN/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-+.f642.9%
Simplified2.9%
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
pow1/2N/A
pow1/2N/A
pow-prod-downN/A
pow-lowering-pow.f64N/A
Applied egg-rr3.1%
Taylor expanded in x around inf
sqrt-lowering-sqrt.f6417.2%
Simplified17.2%
if 2.7999999999999998e-227 < t Initial program 34.6%
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-+.f6484.2%
Simplified84.2%
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-+.f6484.2%
Simplified84.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
(if (<= t_m 1.3e-241)
(* 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 <= 1.3e-241) {
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 <= 1.3d-241) 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 <= 1.3e-241) {
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 <= 1.3e-241: 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 <= 1.3e-241) 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 <= 1.3e-241) 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, 1.3e-241], 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 1.3 \cdot 10^{-241}:\\
\;\;\;\;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 < 1.3e-241Initial program 32.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-rr32.3%
Taylor expanded in l around inf
associate-/r*N/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
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6412.8%
Simplified12.8%
Taylor expanded in x around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6417.8%
Simplified17.8%
if 1.3e-241 < t Initial program 34.6%
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-+.f6484.2%
Simplified84.2%
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-+.f6484.2%
Simplified84.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
(if (<= t_m 9.5e-269)
(/ (+ (/ x 2.0) -0.5) (* x (* x x)))
(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 <= 9.5e-269) {
tmp = ((x / 2.0) + -0.5) / (x * (x * x));
} 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 <= 9.5d-269) then
tmp = ((x / 2.0d0) + (-0.5d0)) / (x * (x * x))
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 <= 9.5e-269) {
tmp = ((x / 2.0) + -0.5) / (x * (x * x));
} 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 <= 9.5e-269: tmp = ((x / 2.0) + -0.5) / (x * (x * x)) 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 <= 9.5e-269) tmp = Float64(Float64(Float64(x / 2.0) + -0.5) / Float64(x * Float64(x * x))); 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 <= 9.5e-269) tmp = ((x / 2.0) + -0.5) / (x * (x * x)); 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, 9.5e-269], N[(N[(N[(x / 2.0), $MachinePrecision] + -0.5), $MachinePrecision] / N[(x * N[(x * x), $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 9.5 \cdot 10^{-269}:\\
\;\;\;\;\frac{\frac{x}{2} + -0.5}{x \cdot \left(x \cdot x\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\end{array}
\end{array}
if t < 9.5000000000000006e-269Initial program 32.1%
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-+.f642.9%
Simplified2.9%
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-+.f642.9%
Simplified2.9%
Taylor expanded in x around inf
associate--r+N/A
sub-negN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f642.9%
Simplified2.9%
Taylor expanded in x around 0
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
metadata-evalN/A
associate-/l*N/A
*-commutativeN/A
mul-1-negN/A
distribute-frac-negN/A
distribute-neg-frac2N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f649.0%
Simplified9.0%
if 9.5000000000000006e-269 < t Initial program 34.6%
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.9%
Simplified82.9%
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.9%
Simplified82.9%
Final simplification41.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
(if (<= t_m 2.9e-268)
(/ (+ (/ x 2.0) -0.5) (* x (* x x)))
(+ 1.0 (/ (- -1.0 (/ (- (/ 0.5 x) 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) {
double tmp;
if (t_m <= 2.9e-268) {
tmp = ((x / 2.0) + -0.5) / (x * (x * x));
} else {
tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x);
}
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.9d-268) then
tmp = ((x / 2.0d0) + (-0.5d0)) / (x * (x * x))
else
tmp = 1.0d0 + (((-1.0d0) - (((0.5d0 / x) - 0.5d0) / x)) / x)
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.9e-268) {
tmp = ((x / 2.0) + -0.5) / (x * (x * x));
} else {
tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x);
}
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.9e-268: tmp = ((x / 2.0) + -0.5) / (x * (x * x)) else: tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x) 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.9e-268) tmp = Float64(Float64(Float64(x / 2.0) + -0.5) / Float64(x * Float64(x * x))); else tmp = Float64(1.0 + Float64(Float64(-1.0 - Float64(Float64(Float64(0.5 / x) - 0.5) / x)) / x)); 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.9e-268) tmp = ((x / 2.0) + -0.5) / (x * (x * x)); else tmp = 1.0 + ((-1.0 - (((0.5 / x) - 0.5) / x)) / x); 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.9e-268], N[(N[(N[(x / 2.0), $MachinePrecision] + -0.5), $MachinePrecision] / N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(N[(-1.0 - N[(N[(N[(0.5 / x), $MachinePrecision] - 0.5), $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 \begin{array}{l}
\mathbf{if}\;t\_m \leq 2.9 \cdot 10^{-268}:\\
\;\;\;\;\frac{\frac{x}{2} + -0.5}{x \cdot \left(x \cdot x\right)}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1 - \frac{\frac{0.5}{x} - 0.5}{x}}{x}\\
\end{array}
\end{array}
if t < 2.9000000000000002e-268Initial program 32.1%
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-+.f642.9%
Simplified2.9%
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-+.f642.9%
Simplified2.9%
Taylor expanded in x around inf
associate--r+N/A
sub-negN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f642.9%
Simplified2.9%
Taylor expanded in x around 0
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
metadata-evalN/A
associate-/l*N/A
*-commutativeN/A
mul-1-negN/A
distribute-frac-negN/A
distribute-neg-frac2N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f649.0%
Simplified9.0%
if 2.9000000000000002e-268 < t Initial program 34.6%
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.9%
Simplified82.9%
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.9%
Simplified82.9%
Taylor expanded in x around inf
associate--r+N/A
sub-negN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.9%
Simplified82.9%
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
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6482.9%
Simplified82.9%
Final simplification41.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
(if (<= t_m 8.5e-269)
(/ (+ (/ x 2.0) -0.5) (* x (* x x)))
(+ 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) {
double tmp;
if (t_m <= 8.5e-269) {
tmp = ((x / 2.0) + -0.5) / (x * (x * x));
} else {
tmp = 1.0 + ((-1.0 + (0.5 / x)) / x);
}
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 <= 8.5d-269) then
tmp = ((x / 2.0d0) + (-0.5d0)) / (x * (x * x))
else
tmp = 1.0d0 + (((-1.0d0) + (0.5d0 / x)) / x)
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 <= 8.5e-269) {
tmp = ((x / 2.0) + -0.5) / (x * (x * x));
} else {
tmp = 1.0 + ((-1.0 + (0.5 / x)) / x);
}
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 <= 8.5e-269: tmp = ((x / 2.0) + -0.5) / (x * (x * x)) else: tmp = 1.0 + ((-1.0 + (0.5 / x)) / x) 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 <= 8.5e-269) tmp = Float64(Float64(Float64(x / 2.0) + -0.5) / Float64(x * Float64(x * x))); else tmp = Float64(1.0 + Float64(Float64(-1.0 + Float64(0.5 / x)) / x)); 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 <= 8.5e-269) tmp = ((x / 2.0) + -0.5) / (x * (x * x)); else tmp = 1.0 + ((-1.0 + (0.5 / x)) / x); 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, 8.5e-269], N[(N[(N[(x / 2.0), $MachinePrecision] + -0.5), $MachinePrecision] / N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 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 \begin{array}{l}
\mathbf{if}\;t\_m \leq 8.5 \cdot 10^{-269}:\\
\;\;\;\;\frac{\frac{x}{2} + -0.5}{x \cdot \left(x \cdot x\right)}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1 + \frac{0.5}{x}}{x}\\
\end{array}
\end{array}
if t < 8.5e-269Initial program 32.1%
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-+.f642.9%
Simplified2.9%
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-+.f642.9%
Simplified2.9%
Taylor expanded in x around inf
associate--r+N/A
sub-negN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f642.9%
Simplified2.9%
Taylor expanded in x around 0
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
metadata-evalN/A
associate-/l*N/A
*-commutativeN/A
mul-1-negN/A
distribute-frac-negN/A
distribute-neg-frac2N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f649.0%
Simplified9.0%
if 8.5e-269 < t Initial program 34.6%
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.9%
Simplified82.9%
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.9%
Simplified82.9%
Taylor expanded in x around inf
associate--r+N/A
sub-negN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.9%
Simplified82.9%
Taylor expanded in x around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6482.7%
Simplified82.7%
Final simplification41.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 (+ 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 33.2%
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-+.f6438.5%
Simplified38.5%
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-+.f6438.5%
Simplified38.5%
Taylor expanded in x around inf
associate--r+N/A
sub-negN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6438.5%
Simplified38.5%
Taylor expanded in x around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6438.4%
Simplified38.4%
Final simplification38.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 (+ 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 33.2%
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-+.f6438.5%
Simplified38.5%
Taylor expanded in x around inf
--lowering--.f64N/A
/-lowering-/.f6438.3%
Simplified38.3%
Final simplification38.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))
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 33.2%
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-+.f6438.5%
Simplified38.5%
Taylor expanded in x around inf
Simplified38.0%
herbie shell --seed 2024144
(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)))))