
(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}
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l t_m)
:precision binary64
(let* ((t_2 (* 2.0 (* t_m t_m)))
(t_3 (+ (* l l) t_2))
(t_4 (* t_m (sqrt 2.0))))
(*
t_s
(if (<= t_m 7e-156)
(/ t_4 (+ t_4 (/ (* 0.5 (* 2.0 t_3)) (* t_m (* (sqrt 2.0) x)))))
(if (<= t_m 3.3e+51)
(*
t_m
(sqrt
(/
2.0
(+
t_2
(/
(+
(+ t_2 (+ (* l l) t_3))
(/ (+ (+ t_3 t_3) (+ (+ (/ t_2 x) (/ (* l l) x)) (/ t_3 x))) x))
x)))))
(sqrt (/ (+ x -1.0) (+ x 1.0))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
double t_2 = 2.0 * (t_m * t_m);
double t_3 = (l * l) + t_2;
double t_4 = t_m * sqrt(2.0);
double tmp;
if (t_m <= 7e-156) {
tmp = t_4 / (t_4 + ((0.5 * (2.0 * t_3)) / (t_m * (sqrt(2.0) * x))));
} else if (t_m <= 3.3e+51) {
tmp = t_m * sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x))));
} else {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
real(8) :: t_2
real(8) :: t_3
real(8) :: t_4
real(8) :: tmp
t_2 = 2.0d0 * (t_m * t_m)
t_3 = (l * l) + t_2
t_4 = t_m * sqrt(2.0d0)
if (t_m <= 7d-156) then
tmp = t_4 / (t_4 + ((0.5d0 * (2.0d0 * t_3)) / (t_m * (sqrt(2.0d0) * x))))
else if (t_m <= 3.3d+51) then
tmp = t_m * sqrt((2.0d0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x))))
else
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
double t_2 = 2.0 * (t_m * t_m);
double t_3 = (l * l) + t_2;
double t_4 = t_m * Math.sqrt(2.0);
double tmp;
if (t_m <= 7e-156) {
tmp = t_4 / (t_4 + ((0.5 * (2.0 * t_3)) / (t_m * (Math.sqrt(2.0) * x))));
} else if (t_m <= 3.3e+51) {
tmp = t_m * Math.sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x))));
} else {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): t_2 = 2.0 * (t_m * t_m) t_3 = (l * l) + t_2 t_4 = t_m * math.sqrt(2.0) tmp = 0 if t_m <= 7e-156: tmp = t_4 / (t_4 + ((0.5 * (2.0 * t_3)) / (t_m * (math.sqrt(2.0) * x)))) elif t_m <= 3.3e+51: tmp = t_m * math.sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x)))) else: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) t_2 = Float64(2.0 * Float64(t_m * t_m)) t_3 = Float64(Float64(l * l) + t_2) t_4 = Float64(t_m * sqrt(2.0)) tmp = 0.0 if (t_m <= 7e-156) tmp = Float64(t_4 / Float64(t_4 + Float64(Float64(0.5 * Float64(2.0 * t_3)) / Float64(t_m * Float64(sqrt(2.0) * x))))); elseif (t_m <= 3.3e+51) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(t_2 + Float64(Float64(Float64(t_2 + Float64(Float64(l * l) + t_3)) + Float64(Float64(Float64(t_3 + t_3) + Float64(Float64(Float64(t_2 / x) + Float64(Float64(l * l) / x)) + Float64(t_3 / x))) / x)) / x))))); else tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l, t_m) t_2 = 2.0 * (t_m * t_m); t_3 = (l * l) + t_2; t_4 = t_m * sqrt(2.0); tmp = 0.0; if (t_m <= 7e-156) tmp = t_4 / (t_4 + ((0.5 * (2.0 * t_3)) / (t_m * (sqrt(2.0) * x)))); elseif (t_m <= 3.3e+51) tmp = t_m * sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x)))); else tmp = sqrt(((x + -1.0) / (x + 1.0))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := Block[{t$95$2 = N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(l * l), $MachinePrecision] + t$95$2), $MachinePrecision]}, Block[{t$95$4 = N[(t$95$m * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 7e-156], N[(t$95$4 / N[(t$95$4 + N[(N[(0.5 * N[(2.0 * t$95$3), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[(N[Sqrt[2.0], $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 3.3e+51], N[(t$95$m * N[Sqrt[N[(2.0 / N[(t$95$2 + N[(N[(N[(t$95$2 + N[(N[(l * l), $MachinePrecision] + t$95$3), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(t$95$3 + t$95$3), $MachinePrecision] + N[(N[(N[(t$95$2 / x), $MachinePrecision] + N[(N[(l * l), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] + N[(t$95$3 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]), $MachinePrecision]]]]
\begin{array}{l}
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_3 := \ell \cdot \ell + t\_2\\
t_4 := t\_m \cdot \sqrt{2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 7 \cdot 10^{-156}:\\
\;\;\;\;\frac{t\_4}{t\_4 + \frac{0.5 \cdot \left(2 \cdot t\_3\right)}{t\_m \cdot \left(\sqrt{2} \cdot x\right)}}\\
\mathbf{elif}\;t\_m \leq 3.3 \cdot 10^{+51}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{t\_2 + \frac{\left(t\_2 + \left(\ell \cdot \ell + t\_3\right)\right) + \frac{\left(t\_3 + t\_3\right) + \left(\left(\frac{t\_2}{x} + \frac{\ell \cdot \ell}{x}\right) + \frac{t\_3}{x}\right)}{x}}{x}}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\end{array}
\end{array}
\end{array}
if t < 6.9999999999999999e-156Initial program 28.5%
Taylor expanded in x around inf
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
Simplified16.4%
if 6.9999999999999999e-156 < t < 3.2999999999999997e51Initial program 46.8%
*-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-rr46.9%
Taylor expanded in x around -inf
Simplified85.9%
if 3.2999999999999997e51 < t Initial program 30.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-+.f6494.5%
Simplified94.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-+.f6494.5%
Simplified94.5%
Final simplification48.5%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l t_m)
:precision binary64
(let* ((t_2 (* 2.0 (* t_m t_m))) (t_3 (+ (* l l) t_2)))
(*
t_s
(if (<= t_m 9.8e-171)
(* t_m (* (/ (sqrt 2.0) l) (sqrt (/ 1.0 (+ (/ 2.0 x) (/ 2.0 (* x x)))))))
(if (<= t_m 8.5e+49)
(*
t_m
(sqrt
(/
2.0
(+
t_2
(/
(+
(+ t_2 (+ (* l l) t_3))
(/ (+ (+ t_3 t_3) (+ (+ (/ t_2 x) (/ (* l l) x)) (/ t_3 x))) x))
x)))))
(sqrt (/ (+ x -1.0) (+ x 1.0))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
double t_2 = 2.0 * (t_m * t_m);
double t_3 = (l * l) + t_2;
double tmp;
if (t_m <= 9.8e-171) {
tmp = t_m * ((sqrt(2.0) / l) * sqrt((1.0 / ((2.0 / x) + (2.0 / (x * x))))));
} else if (t_m <= 8.5e+49) {
tmp = t_m * sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x))));
} else {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_2 = 2.0d0 * (t_m * t_m)
t_3 = (l * l) + t_2
if (t_m <= 9.8d-171) then
tmp = t_m * ((sqrt(2.0d0) / l) * sqrt((1.0d0 / ((2.0d0 / x) + (2.0d0 / (x * x))))))
else if (t_m <= 8.5d+49) then
tmp = t_m * sqrt((2.0d0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x))))
else
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
double t_2 = 2.0 * (t_m * t_m);
double t_3 = (l * l) + t_2;
double tmp;
if (t_m <= 9.8e-171) {
tmp = t_m * ((Math.sqrt(2.0) / l) * Math.sqrt((1.0 / ((2.0 / x) + (2.0 / (x * x))))));
} else if (t_m <= 8.5e+49) {
tmp = t_m * Math.sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x))));
} else {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): t_2 = 2.0 * (t_m * t_m) t_3 = (l * l) + t_2 tmp = 0 if t_m <= 9.8e-171: tmp = t_m * ((math.sqrt(2.0) / l) * math.sqrt((1.0 / ((2.0 / x) + (2.0 / (x * x)))))) elif t_m <= 8.5e+49: tmp = t_m * math.sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x)))) else: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) t_2 = Float64(2.0 * Float64(t_m * t_m)) t_3 = Float64(Float64(l * l) + t_2) tmp = 0.0 if (t_m <= 9.8e-171) tmp = Float64(t_m * Float64(Float64(sqrt(2.0) / l) * sqrt(Float64(1.0 / Float64(Float64(2.0 / x) + Float64(2.0 / Float64(x * x))))))); elseif (t_m <= 8.5e+49) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(t_2 + Float64(Float64(Float64(t_2 + Float64(Float64(l * l) + t_3)) + Float64(Float64(Float64(t_3 + t_3) + Float64(Float64(Float64(t_2 / x) + Float64(Float64(l * l) / x)) + Float64(t_3 / x))) / x)) / x))))); else tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l, t_m) t_2 = 2.0 * (t_m * t_m); t_3 = (l * l) + t_2; tmp = 0.0; if (t_m <= 9.8e-171) tmp = t_m * ((sqrt(2.0) / l) * sqrt((1.0 / ((2.0 / x) + (2.0 / (x * x)))))); elseif (t_m <= 8.5e+49) tmp = t_m * sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x)))); else tmp = sqrt(((x + -1.0) / (x + 1.0))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := Block[{t$95$2 = N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(l * l), $MachinePrecision] + t$95$2), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 9.8e-171], N[(t$95$m * N[(N[(N[Sqrt[2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[N[(1.0 / N[(N[(2.0 / x), $MachinePrecision] + N[(2.0 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 8.5e+49], N[(t$95$m * N[Sqrt[N[(2.0 / N[(t$95$2 + N[(N[(N[(t$95$2 + N[(N[(l * l), $MachinePrecision] + t$95$3), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(t$95$3 + t$95$3), $MachinePrecision] + N[(N[(N[(t$95$2 / x), $MachinePrecision] + N[(N[(l * l), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] + N[(t$95$3 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]), $MachinePrecision]]]
\begin{array}{l}
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_3 := \ell \cdot \ell + t\_2\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 9.8 \cdot 10^{-171}:\\
\;\;\;\;t\_m \cdot \left(\frac{\sqrt{2}}{\ell} \cdot \sqrt{\frac{1}{\frac{2}{x} + \frac{2}{x \cdot x}}}\right)\\
\mathbf{elif}\;t\_m \leq 8.5 \cdot 10^{+49}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{t\_2 + \frac{\left(t\_2 + \left(\ell \cdot \ell + t\_3\right)\right) + \frac{\left(t\_3 + t\_3\right) + \left(\left(\frac{t\_2}{x} + \frac{\ell \cdot \ell}{x}\right) + \frac{t\_3}{x}\right)}{x}}{x}}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\end{array}
\end{array}
\end{array}
if t < 9.79999999999999965e-171Initial program 29.0%
*-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-rr29.1%
Taylor expanded in x around inf
--lowering--.f64N/A
Simplified46.4%
Taylor expanded in l around inf
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
distribute-rgt1-inN/A
metadata-evalN/A
distribute-lft-outN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Simplified18.8%
if 9.79999999999999965e-171 < t < 8.4999999999999996e49Initial program 44.5%
*-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-rr44.6%
Taylor expanded in x around -inf
Simplified81.4%
if 8.4999999999999996e49 < t Initial program 30.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-+.f6494.5%
Simplified94.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-+.f6494.5%
Simplified94.5%
Final simplification49.6%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l t_m)
:precision binary64
(let* ((t_2 (* 2.0 (* t_m t_m))) (t_3 (+ (* l l) t_2)))
(*
t_s
(if (<= t_m 1.7e-171)
(/ (* t_m (sqrt x)) l)
(if (<= t_m 4e+49)
(*
t_m
(sqrt
(/
2.0
(+
t_2
(/
(+
(+ t_2 (+ (* l l) t_3))
(/ (+ (+ t_3 t_3) (+ (+ (/ t_2 x) (/ (* l l) x)) (/ t_3 x))) x))
x)))))
(sqrt (/ (+ x -1.0) (+ x 1.0))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
double t_2 = 2.0 * (t_m * t_m);
double t_3 = (l * l) + t_2;
double tmp;
if (t_m <= 1.7e-171) {
tmp = (t_m * sqrt(x)) / l;
} else if (t_m <= 4e+49) {
tmp = t_m * sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x))));
} else {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_2 = 2.0d0 * (t_m * t_m)
t_3 = (l * l) + t_2
if (t_m <= 1.7d-171) then
tmp = (t_m * sqrt(x)) / l
else if (t_m <= 4d+49) then
tmp = t_m * sqrt((2.0d0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x))))
else
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
double t_2 = 2.0 * (t_m * t_m);
double t_3 = (l * l) + t_2;
double tmp;
if (t_m <= 1.7e-171) {
tmp = (t_m * Math.sqrt(x)) / l;
} else if (t_m <= 4e+49) {
tmp = t_m * Math.sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x))));
} else {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): t_2 = 2.0 * (t_m * t_m) t_3 = (l * l) + t_2 tmp = 0 if t_m <= 1.7e-171: tmp = (t_m * math.sqrt(x)) / l elif t_m <= 4e+49: tmp = t_m * math.sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x)))) else: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) t_2 = Float64(2.0 * Float64(t_m * t_m)) t_3 = Float64(Float64(l * l) + t_2) tmp = 0.0 if (t_m <= 1.7e-171) tmp = Float64(Float64(t_m * sqrt(x)) / l); elseif (t_m <= 4e+49) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(t_2 + Float64(Float64(Float64(t_2 + Float64(Float64(l * l) + t_3)) + Float64(Float64(Float64(t_3 + t_3) + Float64(Float64(Float64(t_2 / x) + Float64(Float64(l * l) / x)) + Float64(t_3 / x))) / x)) / x))))); else tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l, t_m) t_2 = 2.0 * (t_m * t_m); t_3 = (l * l) + t_2; tmp = 0.0; if (t_m <= 1.7e-171) tmp = (t_m * sqrt(x)) / l; elseif (t_m <= 4e+49) tmp = t_m * sqrt((2.0 / (t_2 + (((t_2 + ((l * l) + t_3)) + (((t_3 + t_3) + (((t_2 / x) + ((l * l) / x)) + (t_3 / x))) / x)) / x)))); else tmp = sqrt(((x + -1.0) / (x + 1.0))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := Block[{t$95$2 = N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(l * l), $MachinePrecision] + t$95$2), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 1.7e-171], N[(N[(t$95$m * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision], If[LessEqual[t$95$m, 4e+49], N[(t$95$m * N[Sqrt[N[(2.0 / N[(t$95$2 + N[(N[(N[(t$95$2 + N[(N[(l * l), $MachinePrecision] + t$95$3), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(t$95$3 + t$95$3), $MachinePrecision] + N[(N[(N[(t$95$2 / x), $MachinePrecision] + N[(N[(l * l), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] + N[(t$95$3 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]), $MachinePrecision]]]
\begin{array}{l}
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_3 := \ell \cdot \ell + t\_2\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.7 \cdot 10^{-171}:\\
\;\;\;\;\frac{t\_m \cdot \sqrt{x}}{\ell}\\
\mathbf{elif}\;t\_m \leq 4 \cdot 10^{+49}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{t\_2 + \frac{\left(t\_2 + \left(\ell \cdot \ell + t\_3\right)\right) + \frac{\left(t\_3 + t\_3\right) + \left(\left(\frac{t\_2}{x} + \frac{\ell \cdot \ell}{x}\right) + \frac{t\_3}{x}\right)}{x}}{x}}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\end{array}
\end{array}
\end{array}
if t < 1.69999999999999993e-171Initial program 29.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
--lowering--.f64N/A
Simplified19.8%
Taylor expanded in l around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f643.5%
Simplified3.5%
Taylor expanded in x around inf
/-lowering-/.f64N/A
cancel-sign-sub-invN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6417.1%
Simplified17.1%
*-lft-identityN/A
distribute-rgt-outN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f6417.1%
Applied egg-rr17.1%
Taylor expanded in t around 0
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6418.5%
Simplified18.5%
if 1.69999999999999993e-171 < t < 3.99999999999999979e49Initial program 44.5%
*-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-rr44.6%
Taylor expanded in x around -inf
Simplified81.4%
if 3.99999999999999979e49 < t Initial program 30.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-+.f6494.5%
Simplified94.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-+.f6494.5%
Simplified94.5%
Final simplification49.5%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l t_m)
:precision binary64
(let* ((t_2 (* 2.0 (* t_m t_m))))
(*
t_s
(if (<= t_m 1.65e-171)
(/ (* t_m (sqrt x)) l)
(if (<= t_m 1.05e+49)
(*
t_m
(/
1.0
(sqrt
(/
(+
(+
t_2
(+
(* (/ l x) (+ l (/ l x)))
(* (/ (* t_m t_m) x) (+ 2.0 (/ 2.0 x)))))
(* (+ (* l l) t_2) (+ (/ 1.0 x) (/ 1.0 (* x x)))))
2.0))))
(sqrt (/ (+ x -1.0) (+ x 1.0))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
double t_2 = 2.0 * (t_m * t_m);
double tmp;
if (t_m <= 1.65e-171) {
tmp = (t_m * sqrt(x)) / l;
} else if (t_m <= 1.05e+49) {
tmp = t_m * (1.0 / sqrt((((t_2 + (((l / x) * (l + (l / x))) + (((t_m * t_m) / x) * (2.0 + (2.0 / x))))) + (((l * l) + t_2) * ((1.0 / x) + (1.0 / (x * x))))) / 2.0)));
} else {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
real(8) :: t_2
real(8) :: tmp
t_2 = 2.0d0 * (t_m * t_m)
if (t_m <= 1.65d-171) then
tmp = (t_m * sqrt(x)) / l
else if (t_m <= 1.05d+49) then
tmp = t_m * (1.0d0 / sqrt((((t_2 + (((l / x) * (l + (l / x))) + (((t_m * t_m) / x) * (2.0d0 + (2.0d0 / x))))) + (((l * l) + t_2) * ((1.0d0 / x) + (1.0d0 / (x * x))))) / 2.0d0)))
else
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
double t_2 = 2.0 * (t_m * t_m);
double tmp;
if (t_m <= 1.65e-171) {
tmp = (t_m * Math.sqrt(x)) / l;
} else if (t_m <= 1.05e+49) {
tmp = t_m * (1.0 / Math.sqrt((((t_2 + (((l / x) * (l + (l / x))) + (((t_m * t_m) / x) * (2.0 + (2.0 / x))))) + (((l * l) + t_2) * ((1.0 / x) + (1.0 / (x * x))))) / 2.0)));
} else {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): t_2 = 2.0 * (t_m * t_m) tmp = 0 if t_m <= 1.65e-171: tmp = (t_m * math.sqrt(x)) / l elif t_m <= 1.05e+49: tmp = t_m * (1.0 / math.sqrt((((t_2 + (((l / x) * (l + (l / x))) + (((t_m * t_m) / x) * (2.0 + (2.0 / x))))) + (((l * l) + t_2) * ((1.0 / x) + (1.0 / (x * x))))) / 2.0))) else: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) t_2 = Float64(2.0 * Float64(t_m * t_m)) tmp = 0.0 if (t_m <= 1.65e-171) tmp = Float64(Float64(t_m * sqrt(x)) / l); elseif (t_m <= 1.05e+49) tmp = Float64(t_m * Float64(1.0 / sqrt(Float64(Float64(Float64(t_2 + Float64(Float64(Float64(l / x) * Float64(l + Float64(l / x))) + Float64(Float64(Float64(t_m * t_m) / x) * Float64(2.0 + Float64(2.0 / x))))) + Float64(Float64(Float64(l * l) + t_2) * Float64(Float64(1.0 / x) + Float64(1.0 / Float64(x * x))))) / 2.0)))); else tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l, t_m) t_2 = 2.0 * (t_m * t_m); tmp = 0.0; if (t_m <= 1.65e-171) tmp = (t_m * sqrt(x)) / l; elseif (t_m <= 1.05e+49) tmp = t_m * (1.0 / sqrt((((t_2 + (((l / x) * (l + (l / x))) + (((t_m * t_m) / x) * (2.0 + (2.0 / x))))) + (((l * l) + t_2) * ((1.0 / x) + (1.0 / (x * x))))) / 2.0))); else tmp = sqrt(((x + -1.0) / (x + 1.0))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, 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, 1.65e-171], N[(N[(t$95$m * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision], If[LessEqual[t$95$m, 1.05e+49], N[(t$95$m * N[(1.0 / N[Sqrt[N[(N[(N[(t$95$2 + N[(N[(N[(l / x), $MachinePrecision] * N[(l + N[(l / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(t$95$m * t$95$m), $MachinePrecision] / x), $MachinePrecision] * N[(2.0 + N[(2.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(l * l), $MachinePrecision] + t$95$2), $MachinePrecision] * N[(N[(1.0 / x), $MachinePrecision] + N[(1.0 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
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 1.65 \cdot 10^{-171}:\\
\;\;\;\;\frac{t\_m \cdot \sqrt{x}}{\ell}\\
\mathbf{elif}\;t\_m \leq 1.05 \cdot 10^{+49}:\\
\;\;\;\;t\_m \cdot \frac{1}{\sqrt{\frac{\left(t\_2 + \left(\frac{\ell}{x} \cdot \left(\ell + \frac{\ell}{x}\right) + \frac{t\_m \cdot t\_m}{x} \cdot \left(2 + \frac{2}{x}\right)\right)\right) + \left(\ell \cdot \ell + t\_2\right) \cdot \left(\frac{1}{x} + \frac{1}{x \cdot x}\right)}{2}}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\end{array}
\end{array}
\end{array}
if t < 1.6500000000000001e-171Initial program 29.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
--lowering--.f64N/A
Simplified19.8%
Taylor expanded in l around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f643.5%
Simplified3.5%
Taylor expanded in x around inf
/-lowering-/.f64N/A
cancel-sign-sub-invN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6417.1%
Simplified17.1%
*-lft-identityN/A
distribute-rgt-outN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f6417.1%
Applied egg-rr17.1%
Taylor expanded in t around 0
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6418.5%
Simplified18.5%
if 1.6500000000000001e-171 < t < 1.05000000000000005e49Initial program 44.5%
*-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-rr44.6%
Taylor expanded in x around inf
--lowering--.f64N/A
Simplified81.0%
Applied egg-rr81.3%
if 1.05000000000000005e49 < t Initial program 30.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-+.f6494.5%
Simplified94.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-+.f6494.5%
Simplified94.5%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l t_m)
:precision binary64
(*
t_s
(if (<= t_m 1.12e-170)
(/ (* t_m (sqrt x)) l)
(if (<= t_m 3.5e+47)
(*
t_m
(sqrt
(/
2.0
(+
(/ (+ (* l l) (* 2.0 (* t_m t_m))) x)
(+ (/ (* l l) x) (* 2.0 (+ (* t_m t_m) (/ (* t_m t_m) x))))))))
(sqrt (/ (+ x -1.0) (+ x 1.0)))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
double tmp;
if (t_m <= 1.12e-170) {
tmp = (t_m * sqrt(x)) / l;
} else if (t_m <= 3.5e+47) {
tmp = t_m * sqrt((2.0 / ((((l * l) + (2.0 * (t_m * t_m))) / x) + (((l * l) / x) + (2.0 * ((t_m * t_m) + ((t_m * t_m) / x)))))));
} else {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
real(8) :: tmp
if (t_m <= 1.12d-170) then
tmp = (t_m * sqrt(x)) / l
else if (t_m <= 3.5d+47) then
tmp = t_m * sqrt((2.0d0 / ((((l * l) + (2.0d0 * (t_m * t_m))) / x) + (((l * l) / x) + (2.0d0 * ((t_m * t_m) + ((t_m * t_m) / x)))))))
else
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
double tmp;
if (t_m <= 1.12e-170) {
tmp = (t_m * Math.sqrt(x)) / l;
} else if (t_m <= 3.5e+47) {
tmp = t_m * Math.sqrt((2.0 / ((((l * l) + (2.0 * (t_m * t_m))) / x) + (((l * l) / x) + (2.0 * ((t_m * t_m) + ((t_m * t_m) / x)))))));
} else {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): tmp = 0 if t_m <= 1.12e-170: tmp = (t_m * math.sqrt(x)) / l elif t_m <= 3.5e+47: tmp = t_m * math.sqrt((2.0 / ((((l * l) + (2.0 * (t_m * t_m))) / x) + (((l * l) / x) + (2.0 * ((t_m * t_m) + ((t_m * t_m) / x))))))) else: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) tmp = 0.0 if (t_m <= 1.12e-170) tmp = Float64(Float64(t_m * sqrt(x)) / l); elseif (t_m <= 3.5e+47) tmp = Float64(t_m * sqrt(Float64(2.0 / Float64(Float64(Float64(Float64(l * l) + Float64(2.0 * Float64(t_m * t_m))) / x) + Float64(Float64(Float64(l * l) / x) + Float64(2.0 * Float64(Float64(t_m * t_m) + Float64(Float64(t_m * t_m) / x)))))))); else tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l, t_m) tmp = 0.0; if (t_m <= 1.12e-170) tmp = (t_m * sqrt(x)) / l; elseif (t_m <= 3.5e+47) tmp = t_m * sqrt((2.0 / ((((l * l) + (2.0 * (t_m * t_m))) / x) + (((l * l) / x) + (2.0 * ((t_m * t_m) + ((t_m * t_m) / x))))))); else tmp = sqrt(((x + -1.0) / (x + 1.0))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 1.12e-170], N[(N[(t$95$m * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision], If[LessEqual[t$95$m, 3.5e+47], N[(t$95$m * N[Sqrt[N[(2.0 / N[(N[(N[(N[(l * l), $MachinePrecision] + N[(2.0 * N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] + N[(N[(N[(l * l), $MachinePrecision] / x), $MachinePrecision] + 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]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.12 \cdot 10^{-170}:\\
\;\;\;\;\frac{t\_m \cdot \sqrt{x}}{\ell}\\
\mathbf{elif}\;t\_m \leq 3.5 \cdot 10^{+47}:\\
\;\;\;\;t\_m \cdot \sqrt{\frac{2}{\frac{\ell \cdot \ell + 2 \cdot \left(t\_m \cdot t\_m\right)}{x} + \left(\frac{\ell \cdot \ell}{x} + 2 \cdot \left(t\_m \cdot t\_m + \frac{t\_m \cdot t\_m}{x}\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\end{array}
\end{array}
if t < 1.12000000000000009e-170Initial program 29.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
--lowering--.f64N/A
Simplified19.8%
Taylor expanded in l around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f643.5%
Simplified3.5%
Taylor expanded in x around inf
/-lowering-/.f64N/A
cancel-sign-sub-invN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6417.1%
Simplified17.1%
*-lft-identityN/A
distribute-rgt-outN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f6417.1%
Applied egg-rr17.1%
Taylor expanded in t around 0
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6418.5%
Simplified18.5%
if 1.12000000000000009e-170 < t < 3.50000000000000015e47Initial program 44.5%
*-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-rr44.6%
Taylor expanded in x around inf
cancel-sign-sub-invN/A
+-lowering-+.f64N/A
Simplified80.2%
if 3.50000000000000015e47 < t Initial program 30.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-+.f6494.5%
Simplified94.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-+.f6494.5%
Simplified94.5%
Final simplification49.2%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l t_m)
:precision binary64
(*
t_s
(if (<= t_m 1.32e-194)
(/ (* t_m (sqrt x)) l)
(sqrt (/ (+ x -1.0) (+ x 1.0))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
double tmp;
if (t_m <= 1.32e-194) {
tmp = (t_m * sqrt(x)) / l;
} else {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
real(8) :: tmp
if (t_m <= 1.32d-194) then
tmp = (t_m * sqrt(x)) / l
else
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
double tmp;
if (t_m <= 1.32e-194) {
tmp = (t_m * Math.sqrt(x)) / l;
} else {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): tmp = 0 if t_m <= 1.32e-194: tmp = (t_m * math.sqrt(x)) / l else: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) tmp = 0.0 if (t_m <= 1.32e-194) tmp = Float64(Float64(t_m * sqrt(x)) / l); else tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l, t_m) tmp = 0.0; if (t_m <= 1.32e-194) tmp = (t_m * sqrt(x)) / l; else tmp = sqrt(((x + -1.0) / (x + 1.0))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 1.32e-194], N[(N[(t$95$m * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.32 \cdot 10^{-194}:\\
\;\;\;\;\frac{t\_m \cdot \sqrt{x}}{\ell}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\end{array}
\end{array}
if t < 1.32e-194Initial program 29.4%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
--lowering--.f64N/A
Simplified20.0%
Taylor expanded in l around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f643.5%
Simplified3.5%
Taylor expanded in x around inf
/-lowering-/.f64N/A
cancel-sign-sub-invN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6417.3%
Simplified17.3%
*-lft-identityN/A
distribute-rgt-outN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f6417.3%
Applied egg-rr17.3%
Taylor expanded in t around 0
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6418.8%
Simplified18.8%
if 1.32e-194 < t Initial program 36.8%
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-+.f6478.8%
Simplified78.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-+.f6478.8%
Simplified78.8%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l t_m)
:precision binary64
(*
t_s
(if (<= t_m 8e-191)
(/ (* t_m (sqrt x)) l)
(+ (/ (+ -1.0 (/ (+ 0.5 (/ -0.5 x)) x)) x) 1.0))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
double tmp;
if (t_m <= 8e-191) {
tmp = (t_m * sqrt(x)) / l;
} else {
tmp = ((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0;
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
real(8) :: tmp
if (t_m <= 8d-191) then
tmp = (t_m * sqrt(x)) / l
else
tmp = (((-1.0d0) + ((0.5d0 + ((-0.5d0) / x)) / x)) / x) + 1.0d0
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
double tmp;
if (t_m <= 8e-191) {
tmp = (t_m * Math.sqrt(x)) / l;
} else {
tmp = ((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0;
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): tmp = 0 if t_m <= 8e-191: tmp = (t_m * math.sqrt(x)) / l else: tmp = ((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0 return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) tmp = 0.0 if (t_m <= 8e-191) tmp = Float64(Float64(t_m * sqrt(x)) / l); else tmp = Float64(Float64(Float64(-1.0 + Float64(Float64(0.5 + Float64(-0.5 / x)) / x)) / x) + 1.0); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l, t_m) tmp = 0.0; if (t_m <= 8e-191) tmp = (t_m * sqrt(x)) / l; else tmp = ((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0; end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 8e-191], N[(N[(t$95$m * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision], N[(N[(N[(-1.0 + N[(N[(0.5 + N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] + 1.0), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 8 \cdot 10^{-191}:\\
\;\;\;\;\frac{t\_m \cdot \sqrt{x}}{\ell}\\
\mathbf{else}:\\
\;\;\;\;\frac{-1 + \frac{0.5 + \frac{-0.5}{x}}{x}}{x} + 1\\
\end{array}
\end{array}
if t < 8.0000000000000002e-191Initial program 29.4%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
--lowering--.f64N/A
Simplified20.0%
Taylor expanded in l around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f643.5%
Simplified3.5%
Taylor expanded in x around inf
/-lowering-/.f64N/A
cancel-sign-sub-invN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6417.3%
Simplified17.3%
*-lft-identityN/A
distribute-rgt-outN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f6417.3%
Applied egg-rr17.3%
Taylor expanded in t around 0
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6418.8%
Simplified18.8%
if 8.0000000000000002e-191 < t Initial program 36.8%
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-+.f6478.8%
Simplified78.8%
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
pow1/2N/A
pow1/2N/A
pow-divN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6478.6%
Applied egg-rr78.6%
Taylor expanded in x around inf
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/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-*.f6478.4%
Simplified78.4%
Taylor expanded in x around -inf
Simplified78.6%
Final simplification45.9%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s x l t_m)
:precision binary64
(*
t_s
(if (<= l 1.25e+276)
(+ (/ (+ -1.0 (/ (+ 0.5 (/ -0.5 x)) x)) x) 1.0)
(* (sqrt x) (/ t_m l)))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
double tmp;
if (l <= 1.25e+276) {
tmp = ((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0;
} else {
tmp = sqrt(x) * (t_m / l);
}
return t_s * tmp;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
real(8) :: tmp
if (l <= 1.25d+276) then
tmp = (((-1.0d0) + ((0.5d0 + ((-0.5d0) / x)) / x)) / x) + 1.0d0
else
tmp = sqrt(x) * (t_m / l)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
double tmp;
if (l <= 1.25e+276) {
tmp = ((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0;
} else {
tmp = Math.sqrt(x) * (t_m / l);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): tmp = 0 if l <= 1.25e+276: tmp = ((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0 else: tmp = math.sqrt(x) * (t_m / l) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) tmp = 0.0 if (l <= 1.25e+276) tmp = Float64(Float64(Float64(-1.0 + Float64(Float64(0.5 + Float64(-0.5 / x)) / x)) / x) + 1.0); else tmp = Float64(sqrt(x) * Float64(t_m / l)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, x, l, t_m) tmp = 0.0; if (l <= 1.25e+276) tmp = ((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0; else tmp = sqrt(x) * (t_m / l); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := N[(t$95$s * If[LessEqual[l, 1.25e+276], N[(N[(N[(-1.0 + N[(N[(0.5 + N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] + 1.0), $MachinePrecision], N[(N[Sqrt[x], $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \leq 1.25 \cdot 10^{+276}:\\
\;\;\;\;\frac{-1 + \frac{0.5 + \frac{-0.5}{x}}{x}}{x} + 1\\
\mathbf{else}:\\
\;\;\;\;\sqrt{x} \cdot \frac{t\_m}{\ell}\\
\end{array}
\end{array}
if l < 1.25e276Initial program 33.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-+.f6439.0%
Simplified39.0%
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
pow1/2N/A
pow1/2N/A
pow-divN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6438.9%
Applied egg-rr38.9%
Taylor expanded in x around inf
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/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.8%
Simplified38.8%
Taylor expanded in x around -inf
Simplified38.9%
if 1.25e276 < l Initial program 0.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
--lowering--.f64N/A
Simplified0.0%
Taylor expanded in l around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f640.0%
Simplified0.0%
Taylor expanded in x around inf
/-lowering-/.f64N/A
cancel-sign-sub-invN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6461.9%
Simplified61.9%
Taylor expanded in t around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f6480.6%
Simplified80.6%
Final simplification39.7%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s x l t_m) :precision binary64 (* t_s (+ (/ (+ -1.0 (/ (+ 0.5 (/ -0.5 x)) x)) x) 1.0)))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
return t_s * (((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0);
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
code = t_s * ((((-1.0d0) + ((0.5d0 + ((-0.5d0) / x)) / x)) / x) + 1.0d0)
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
return t_s * (((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0);
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): return t_s * (((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0)
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) return Float64(t_s * Float64(Float64(Float64(-1.0 + Float64(Float64(0.5 + Float64(-0.5 / x)) / x)) / x) + 1.0)) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, x, l, t_m) tmp = t_s * (((-1.0 + ((0.5 + (-0.5 / x)) / x)) / x) + 1.0); end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := N[(t$95$s * N[(N[(N[(-1.0 + N[(N[(0.5 + N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\frac{-1 + \frac{0.5 + \frac{-0.5}{x}}{x}}{x} + 1\right)
\end{array}
Initial program 32.8%
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.6%
Simplified38.6%
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
pow1/2N/A
pow1/2N/A
pow-divN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6438.5%
Applied egg-rr38.5%
Taylor expanded in x around inf
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/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
Simplified38.6%
Final simplification38.6%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s x l t_m) :precision binary64 (* t_s (+ (/ (- -1.0 (/ -0.5 x)) x) 1.0)))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
return t_s * (((-1.0 - (-0.5 / x)) / x) + 1.0);
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
code = t_s * ((((-1.0d0) - ((-0.5d0) / x)) / x) + 1.0d0)
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
return t_s * (((-1.0 - (-0.5 / x)) / x) + 1.0);
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): return t_s * (((-1.0 - (-0.5 / x)) / x) + 1.0)
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) return Float64(t_s * Float64(Float64(Float64(-1.0 - Float64(-0.5 / x)) / x) + 1.0)) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, x, l, t_m) tmp = t_s * (((-1.0 - (-0.5 / x)) / x) + 1.0); end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := N[(t$95$s * N[(N[(N[(-1.0 - N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\frac{-1 - \frac{-0.5}{x}}{x} + 1\right)
\end{array}
Initial program 32.8%
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.6%
Simplified38.6%
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
pow1/2N/A
pow1/2N/A
pow-divN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f6438.5%
Applied egg-rr38.5%
Taylor expanded in x around -inf
Simplified38.5%
Final simplification38.5%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s x l t_m) :precision binary64 (* t_s (+ (/ -1.0 x) 1.0)))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
return t_s * ((-1.0 / x) + 1.0);
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
code = t_s * (((-1.0d0) / x) + 1.0d0)
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
return t_s * ((-1.0 / x) + 1.0);
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): return t_s * ((-1.0 / x) + 1.0)
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) return Float64(t_s * Float64(Float64(-1.0 / x) + 1.0)) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, x, l, t_m) tmp = t_s * ((-1.0 / x) + 1.0); end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := N[(t$95$s * N[(N[(-1.0 / x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\frac{-1}{x} + 1\right)
\end{array}
Initial program 32.8%
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.6%
Simplified38.6%
Taylor expanded in x around inf
--lowering--.f64N/A
/-lowering-/.f6438.5%
Simplified38.5%
Final simplification38.5%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s x l t_m) :precision binary64 (* t_s 1.0))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double x, double l, double t_m) {
return t_s * 1.0;
}
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, x, l, t_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: x
real(8), intent (in) :: l
real(8), intent (in) :: t_m
code = t_s * 1.0d0
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double x, double l, double t_m) {
return t_s * 1.0;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, x, l, t_m): return t_s * 1.0
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, x, l, t_m) return Float64(t_s * 1.0) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, x, l, t_m) tmp = t_s * 1.0; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, x_, l_, t$95$m_] := N[(t$95$s * 1.0), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot 1
\end{array}
Initial program 32.8%
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.6%
Simplified38.6%
Taylor expanded in x around inf
Simplified38.4%
herbie shell --seed 2024162
(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)))))