
(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 9 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 1 t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= l_m 2.4e+121)
(sqrt (/ (+ x -1.0) (+ x 1.0)))
(if (<= l_m 4.5e+147)
(* (/ t_m l_m) (sqrt (* 2.0 (fma x 0.5 -0.5))))
(if (<= l_m 8.8e+189)
1.0
(/ (sqrt 2.0) (/ (* l_m (/ (sqrt 2.0) (sqrt x))) t_m)))))))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.4e+121) {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
} else if (l_m <= 4.5e+147) {
tmp = (t_m / l_m) * sqrt((2.0 * fma(x, 0.5, -0.5)));
} else if (l_m <= 8.8e+189) {
tmp = 1.0;
} else {
tmp = sqrt(2.0) / ((l_m * (sqrt(2.0) / sqrt(x))) / t_m);
}
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.4e+121) tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); elseif (l_m <= 4.5e+147) tmp = Float64(Float64(t_m / l_m) * sqrt(Float64(2.0 * fma(x, 0.5, -0.5)))); elseif (l_m <= 8.8e+189) tmp = 1.0; else tmp = Float64(sqrt(2.0) / Float64(Float64(l_m * Float64(sqrt(2.0) / sqrt(x))) / t_m)); end return Float64(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.4e+121], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[l$95$m, 4.5e+147], N[(N[(t$95$m / l$95$m), $MachinePrecision] * N[Sqrt[N[(2.0 * N[(x * 0.5 + -0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[l$95$m, 8.8e+189], 1.0, N[(N[Sqrt[2.0], $MachinePrecision] / N[(N[(l$95$m * N[(N[Sqrt[2.0], $MachinePrecision] / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$m), $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.4 \cdot 10^{+121}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\mathbf{elif}\;l\_m \leq 4.5 \cdot 10^{+147}:\\
\;\;\;\;\frac{t\_m}{l\_m} \cdot \sqrt{2 \cdot \mathsf{fma}\left(x, 0.5, -0.5\right)}\\
\mathbf{elif}\;l\_m \leq 8.8 \cdot 10^{+189}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2}}{\frac{l\_m \cdot \frac{\sqrt{2}}{\sqrt{x}}}{t\_m}}\\
\end{array}
\end{array}
if l < 2.4e121Initial program 41.7%
Simplified41.6%
Taylor expanded in l around 0 47.4%
associate-*l*47.5%
+-commutative47.5%
sub-neg47.5%
metadata-eval47.5%
+-commutative47.5%
Simplified47.5%
Taylor expanded in t around 0 47.5%
if 2.4e121 < l < 4.50000000000000008e147Initial program 13.8%
Simplified13.8%
Taylor expanded in l around inf 14.0%
*-commutative14.0%
associate--l+29.0%
sub-neg29.0%
metadata-eval29.0%
+-commutative29.0%
sub-neg29.0%
metadata-eval29.0%
+-commutative29.0%
Simplified29.0%
Taylor expanded in x around 0 98.5%
associate-*r*99.0%
clear-num99.0%
un-div-inv99.0%
sqrt-unprod99.2%
*-commutative99.2%
fma-neg99.2%
metadata-eval99.2%
Applied egg-rr99.2%
associate-/r/98.7%
associate-*l/99.5%
associate-*r/99.5%
*-commutative99.5%
Simplified99.5%
if 4.50000000000000008e147 < l < 8.8000000000000002e189Initial program 0.0%
Simplified0.0%
Taylor expanded in l around 0 41.3%
associate-*l*41.3%
+-commutative41.3%
sub-neg41.3%
metadata-eval41.3%
+-commutative41.3%
Simplified41.3%
Taylor expanded in x around inf 41.3%
if 8.8000000000000002e189 < l Initial program 0.0%
Simplified0.0%
Taylor expanded in l around inf 0.0%
associate--l+34.5%
sub-neg34.5%
metadata-eval34.5%
+-commutative34.5%
sub-neg34.5%
metadata-eval34.5%
+-commutative34.5%
Simplified34.5%
Taylor expanded in x around inf 82.2%
clear-num82.0%
un-div-inv82.1%
associate-*l*82.3%
sqrt-div82.3%
metadata-eval82.3%
un-div-inv82.4%
Applied egg-rr82.4%
Final simplification51.1%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= l_m 2.4e+121)
(sqrt (/ (+ x -1.0) (+ x 1.0)))
(if (<= l_m 4e+147)
(* (/ t_m l_m) (sqrt (* 2.0 (fma x 0.5 -0.5))))
(if (<= l_m 1e+190)
1.0
(* (/ (sqrt 2.0) l_m) (/ t_m (sqrt (/ 2.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) {
double tmp;
if (l_m <= 2.4e+121) {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
} else if (l_m <= 4e+147) {
tmp = (t_m / l_m) * sqrt((2.0 * fma(x, 0.5, -0.5)));
} else if (l_m <= 1e+190) {
tmp = 1.0;
} else {
tmp = (sqrt(2.0) / l_m) * (t_m / sqrt((2.0 / 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 (l_m <= 2.4e+121) tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); elseif (l_m <= 4e+147) tmp = Float64(Float64(t_m / l_m) * sqrt(Float64(2.0 * fma(x, 0.5, -0.5)))); elseif (l_m <= 1e+190) tmp = 1.0; else tmp = Float64(Float64(sqrt(2.0) / l_m) * Float64(t_m / sqrt(Float64(2.0 / x)))); end return Float64(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.4e+121], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[l$95$m, 4e+147], N[(N[(t$95$m / l$95$m), $MachinePrecision] * N[Sqrt[N[(2.0 * N[(x * 0.5 + -0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[l$95$m, 1e+190], 1.0, N[(N[(N[Sqrt[2.0], $MachinePrecision] / l$95$m), $MachinePrecision] * N[(t$95$m / N[Sqrt[N[(2.0 / x), $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.4 \cdot 10^{+121}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\mathbf{elif}\;l\_m \leq 4 \cdot 10^{+147}:\\
\;\;\;\;\frac{t\_m}{l\_m} \cdot \sqrt{2 \cdot \mathsf{fma}\left(x, 0.5, -0.5\right)}\\
\mathbf{elif}\;l\_m \leq 10^{+190}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2}}{l\_m} \cdot \frac{t\_m}{\sqrt{\frac{2}{x}}}\\
\end{array}
\end{array}
if l < 2.4e121Initial program 41.7%
Simplified41.6%
Taylor expanded in l around 0 47.4%
associate-*l*47.5%
+-commutative47.5%
sub-neg47.5%
metadata-eval47.5%
+-commutative47.5%
Simplified47.5%
Taylor expanded in t around 0 47.5%
if 2.4e121 < l < 3.9999999999999999e147Initial program 13.8%
Simplified13.8%
Taylor expanded in l around inf 14.0%
*-commutative14.0%
associate--l+29.0%
sub-neg29.0%
metadata-eval29.0%
+-commutative29.0%
sub-neg29.0%
metadata-eval29.0%
+-commutative29.0%
Simplified29.0%
Taylor expanded in x around 0 98.5%
associate-*r*99.0%
clear-num99.0%
un-div-inv99.0%
sqrt-unprod99.2%
*-commutative99.2%
fma-neg99.2%
metadata-eval99.2%
Applied egg-rr99.2%
associate-/r/98.7%
associate-*l/99.5%
associate-*r/99.5%
*-commutative99.5%
Simplified99.5%
if 3.9999999999999999e147 < l < 1.0000000000000001e190Initial program 0.0%
Simplified0.0%
Taylor expanded in l around 0 41.3%
associate-*l*41.3%
+-commutative41.3%
sub-neg41.3%
metadata-eval41.3%
+-commutative41.3%
Simplified41.3%
Taylor expanded in x around inf 41.3%
if 1.0000000000000001e190 < l Initial program 0.0%
Simplified0.0%
Taylor expanded in l around inf 0.0%
associate--l+34.5%
sub-neg34.5%
metadata-eval34.5%
+-commutative34.5%
sub-neg34.5%
metadata-eval34.5%
+-commutative34.5%
Simplified34.5%
Taylor expanded in x around inf 82.2%
clear-num82.0%
un-div-inv82.1%
associate-*l*82.3%
sqrt-div82.3%
metadata-eval82.3%
un-div-inv82.4%
Applied egg-rr82.4%
*-un-lft-identity82.4%
associate-/r/82.4%
sqrt-undiv82.4%
Applied egg-rr82.4%
*-lft-identity82.4%
associate-*l/82.3%
times-frac82.4%
Simplified82.4%
Final simplification51.1%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= l_m 2.3e+121)
(sqrt (/ (+ x -1.0) (+ x 1.0)))
(if (<= l_m 4.2e+147)
(* (/ t_m l_m) (sqrt (* 2.0 (fma x 0.5 -0.5))))
(if (<= l_m 1.22e+190) 1.0 (/ (* t_m (sqrt x)) l_m))))))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.3e+121) {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
} else if (l_m <= 4.2e+147) {
tmp = (t_m / l_m) * sqrt((2.0 * fma(x, 0.5, -0.5)));
} else if (l_m <= 1.22e+190) {
tmp = 1.0;
} else {
tmp = (t_m * sqrt(x)) / l_m;
}
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.3e+121) tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); elseif (l_m <= 4.2e+147) tmp = Float64(Float64(t_m / l_m) * sqrt(Float64(2.0 * fma(x, 0.5, -0.5)))); elseif (l_m <= 1.22e+190) tmp = 1.0; else tmp = Float64(Float64(t_m * sqrt(x)) / l_m); end return Float64(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.3e+121], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[l$95$m, 4.2e+147], N[(N[(t$95$m / l$95$m), $MachinePrecision] * N[Sqrt[N[(2.0 * N[(x * 0.5 + -0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[l$95$m, 1.22e+190], 1.0, N[(N[(t$95$m * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] / l$95$m), $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.3 \cdot 10^{+121}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\mathbf{elif}\;l\_m \leq 4.2 \cdot 10^{+147}:\\
\;\;\;\;\frac{t\_m}{l\_m} \cdot \sqrt{2 \cdot \mathsf{fma}\left(x, 0.5, -0.5\right)}\\
\mathbf{elif}\;l\_m \leq 1.22 \cdot 10^{+190}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_m \cdot \sqrt{x}}{l\_m}\\
\end{array}
\end{array}
if l < 2.2999999999999999e121Initial program 41.7%
Simplified41.6%
Taylor expanded in l around 0 47.4%
associate-*l*47.5%
+-commutative47.5%
sub-neg47.5%
metadata-eval47.5%
+-commutative47.5%
Simplified47.5%
Taylor expanded in t around 0 47.5%
if 2.2999999999999999e121 < l < 4.20000000000000012e147Initial program 13.8%
Simplified13.8%
Taylor expanded in l around inf 14.0%
*-commutative14.0%
associate--l+29.0%
sub-neg29.0%
metadata-eval29.0%
+-commutative29.0%
sub-neg29.0%
metadata-eval29.0%
+-commutative29.0%
Simplified29.0%
Taylor expanded in x around 0 98.5%
associate-*r*99.0%
clear-num99.0%
un-div-inv99.0%
sqrt-unprod99.2%
*-commutative99.2%
fma-neg99.2%
metadata-eval99.2%
Applied egg-rr99.2%
associate-/r/98.7%
associate-*l/99.5%
associate-*r/99.5%
*-commutative99.5%
Simplified99.5%
if 4.20000000000000012e147 < l < 1.21999999999999995e190Initial program 0.0%
Simplified0.0%
Taylor expanded in l around 0 41.3%
associate-*l*41.3%
+-commutative41.3%
sub-neg41.3%
metadata-eval41.3%
+-commutative41.3%
Simplified41.3%
Taylor expanded in x around inf 41.3%
if 1.21999999999999995e190 < l Initial program 0.0%
Simplified0.0%
Taylor expanded in l around inf 0.0%
associate--l+34.5%
sub-neg34.5%
metadata-eval34.5%
+-commutative34.5%
sub-neg34.5%
metadata-eval34.5%
+-commutative34.5%
Simplified34.5%
Taylor expanded in x around inf 82.2%
clear-num82.0%
un-div-inv82.1%
associate-*l*82.3%
sqrt-div82.3%
metadata-eval82.3%
un-div-inv82.4%
Applied egg-rr82.4%
Taylor expanded in l around 0 71.3%
associate-*l/82.4%
Simplified82.4%
Final simplification51.1%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= l_m 1.9e+121)
(+ 1.0 (/ (- -1.0 (/ -0.5 x)) x))
(if (or (<= l_m 4.2e+147) (not (<= l_m 1.1e+190)))
(* (/ t_m l_m) (sqrt 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 (l_m <= 1.9e+121) {
tmp = 1.0 + ((-1.0 - (-0.5 / x)) / x);
} else if ((l_m <= 4.2e+147) || !(l_m <= 1.1e+190)) {
tmp = (t_m / l_m) * sqrt(x);
} else {
tmp = 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 (l_m <= 1.9d+121) then
tmp = 1.0d0 + (((-1.0d0) - ((-0.5d0) / x)) / x)
else if ((l_m <= 4.2d+147) .or. (.not. (l_m <= 1.1d+190))) then
tmp = (t_m / l_m) * sqrt(x)
else
tmp = 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 (l_m <= 1.9e+121) {
tmp = 1.0 + ((-1.0 - (-0.5 / x)) / x);
} else if ((l_m <= 4.2e+147) || !(l_m <= 1.1e+190)) {
tmp = (t_m / l_m) * Math.sqrt(x);
} else {
tmp = 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 l_m <= 1.9e+121: tmp = 1.0 + ((-1.0 - (-0.5 / x)) / x) elif (l_m <= 4.2e+147) or not (l_m <= 1.1e+190): tmp = (t_m / l_m) * math.sqrt(x) else: tmp = 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 (l_m <= 1.9e+121) tmp = Float64(1.0 + Float64(Float64(-1.0 - Float64(-0.5 / x)) / x)); elseif ((l_m <= 4.2e+147) || !(l_m <= 1.1e+190)) tmp = Float64(Float64(t_m / l_m) * sqrt(x)); else tmp = 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 (l_m <= 1.9e+121) tmp = 1.0 + ((-1.0 - (-0.5 / x)) / x); elseif ((l_m <= 4.2e+147) || ~((l_m <= 1.1e+190))) tmp = (t_m / l_m) * sqrt(x); else tmp = 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[l$95$m, 1.9e+121], N[(1.0 + N[(N[(-1.0 - N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[l$95$m, 4.2e+147], N[Not[LessEqual[l$95$m, 1.1e+190]], $MachinePrecision]], N[(N[(t$95$m / l$95$m), $MachinePrecision] * N[Sqrt[x], $MachinePrecision]), $MachinePrecision], 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 \begin{array}{l}
\mathbf{if}\;l\_m \leq 1.9 \cdot 10^{+121}:\\
\;\;\;\;1 + \frac{-1 - \frac{-0.5}{x}}{x}\\
\mathbf{elif}\;l\_m \leq 4.2 \cdot 10^{+147} \lor \neg \left(l\_m \leq 1.1 \cdot 10^{+190}\right):\\
\;\;\;\;\frac{t\_m}{l\_m} \cdot \sqrt{x}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if l < 1.9e121Initial program 41.7%
Simplified41.6%
Taylor expanded in l around 0 47.4%
associate-*l*47.5%
+-commutative47.5%
sub-neg47.5%
metadata-eval47.5%
+-commutative47.5%
Simplified47.5%
Taylor expanded in x around -inf 0.0%
mul-1-neg0.0%
unsub-neg0.0%
Simplified46.8%
if 1.9e121 < l < 4.20000000000000012e147 or 1.1e190 < l Initial program 3.3%
Simplified3.3%
Taylor expanded in l around inf 3.4%
associate--l+33.2%
sub-neg33.2%
metadata-eval33.2%
+-commutative33.2%
sub-neg33.2%
metadata-eval33.2%
+-commutative33.2%
Simplified33.2%
Taylor expanded in x around inf 84.5%
Taylor expanded in t around 0 76.2%
if 4.20000000000000012e147 < l < 1.1e190Initial program 0.0%
Simplified0.0%
Taylor expanded in l around 0 41.3%
associate-*l*41.3%
+-commutative41.3%
sub-neg41.3%
metadata-eval41.3%
+-commutative41.3%
Simplified41.3%
Taylor expanded in x around inf 41.3%
Final simplification49.4%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= l_m 2.4e+121)
(+ 1.0 (/ (- -1.0 (/ -0.5 x)) x))
(if (<= l_m 4.2e+147)
(* (/ t_m l_m) (sqrt x))
(if (<= l_m 1.35e+190) 1.0 (/ (* t_m (sqrt x)) l_m))))))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.4e+121) {
tmp = 1.0 + ((-1.0 - (-0.5 / x)) / x);
} else if (l_m <= 4.2e+147) {
tmp = (t_m / l_m) * sqrt(x);
} else if (l_m <= 1.35e+190) {
tmp = 1.0;
} else {
tmp = (t_m * sqrt(x)) / l_m;
}
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.4d+121) then
tmp = 1.0d0 + (((-1.0d0) - ((-0.5d0) / x)) / x)
else if (l_m <= 4.2d+147) then
tmp = (t_m / l_m) * sqrt(x)
else if (l_m <= 1.35d+190) then
tmp = 1.0d0
else
tmp = (t_m * sqrt(x)) / l_m
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.4e+121) {
tmp = 1.0 + ((-1.0 - (-0.5 / x)) / x);
} else if (l_m <= 4.2e+147) {
tmp = (t_m / l_m) * Math.sqrt(x);
} else if (l_m <= 1.35e+190) {
tmp = 1.0;
} else {
tmp = (t_m * Math.sqrt(x)) / l_m;
}
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.4e+121: tmp = 1.0 + ((-1.0 - (-0.5 / x)) / x) elif l_m <= 4.2e+147: tmp = (t_m / l_m) * math.sqrt(x) elif l_m <= 1.35e+190: tmp = 1.0 else: tmp = (t_m * math.sqrt(x)) / l_m 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.4e+121) tmp = Float64(1.0 + Float64(Float64(-1.0 - Float64(-0.5 / x)) / x)); elseif (l_m <= 4.2e+147) tmp = Float64(Float64(t_m / l_m) * sqrt(x)); elseif (l_m <= 1.35e+190) tmp = 1.0; else tmp = Float64(Float64(t_m * sqrt(x)) / l_m); 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.4e+121) tmp = 1.0 + ((-1.0 - (-0.5 / x)) / x); elseif (l_m <= 4.2e+147) tmp = (t_m / l_m) * sqrt(x); elseif (l_m <= 1.35e+190) tmp = 1.0; else tmp = (t_m * sqrt(x)) / l_m; 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.4e+121], N[(1.0 + N[(N[(-1.0 - N[(-0.5 / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[l$95$m, 4.2e+147], N[(N[(t$95$m / l$95$m), $MachinePrecision] * N[Sqrt[x], $MachinePrecision]), $MachinePrecision], If[LessEqual[l$95$m, 1.35e+190], 1.0, N[(N[(t$95$m * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] / l$95$m), $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.4 \cdot 10^{+121}:\\
\;\;\;\;1 + \frac{-1 - \frac{-0.5}{x}}{x}\\
\mathbf{elif}\;l\_m \leq 4.2 \cdot 10^{+147}:\\
\;\;\;\;\frac{t\_m}{l\_m} \cdot \sqrt{x}\\
\mathbf{elif}\;l\_m \leq 1.35 \cdot 10^{+190}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_m \cdot \sqrt{x}}{l\_m}\\
\end{array}
\end{array}
if l < 2.4e121Initial program 41.7%
Simplified41.6%
Taylor expanded in l around 0 47.4%
associate-*l*47.5%
+-commutative47.5%
sub-neg47.5%
metadata-eval47.5%
+-commutative47.5%
Simplified47.5%
Taylor expanded in x around -inf 0.0%
mul-1-neg0.0%
unsub-neg0.0%
Simplified46.8%
if 2.4e121 < l < 4.20000000000000012e147Initial program 13.8%
Simplified13.8%
Taylor expanded in l around inf 14.0%
associate--l+29.0%
sub-neg29.0%
metadata-eval29.0%
+-commutative29.0%
sub-neg29.0%
metadata-eval29.0%
+-commutative29.0%
Simplified29.0%
Taylor expanded in x around inf 91.7%
Taylor expanded in t around 0 91.7%
if 4.20000000000000012e147 < l < 1.35000000000000002e190Initial program 0.0%
Simplified0.0%
Taylor expanded in l around 0 41.3%
associate-*l*41.3%
+-commutative41.3%
sub-neg41.3%
metadata-eval41.3%
+-commutative41.3%
Simplified41.3%
Taylor expanded in x around inf 41.3%
if 1.35000000000000002e190 < l Initial program 0.0%
Simplified0.0%
Taylor expanded in l around inf 0.0%
associate--l+34.5%
sub-neg34.5%
metadata-eval34.5%
+-commutative34.5%
sub-neg34.5%
metadata-eval34.5%
+-commutative34.5%
Simplified34.5%
Taylor expanded in x around inf 82.2%
clear-num82.0%
un-div-inv82.1%
associate-*l*82.3%
sqrt-div82.3%
metadata-eval82.3%
un-div-inv82.4%
Applied egg-rr82.4%
Taylor expanded in l around 0 71.3%
associate-*l/82.4%
Simplified82.4%
Final simplification50.3%
l_m = (fabs.f64 l)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 1 t)
(FPCore (t_s x l_m t_m)
:precision binary64
(*
t_s
(if (<= l_m 2.3e+121)
(sqrt (/ (+ x -1.0) (+ x 1.0)))
(if (<= l_m 4.2e+147)
(* (/ t_m l_m) (sqrt x))
(if (<= l_m 8.5e+189) 1.0 (/ (* t_m (sqrt x)) l_m))))))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.3e+121) {
tmp = sqrt(((x + -1.0) / (x + 1.0)));
} else if (l_m <= 4.2e+147) {
tmp = (t_m / l_m) * sqrt(x);
} else if (l_m <= 8.5e+189) {
tmp = 1.0;
} else {
tmp = (t_m * sqrt(x)) / l_m;
}
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.3d+121) then
tmp = sqrt(((x + (-1.0d0)) / (x + 1.0d0)))
else if (l_m <= 4.2d+147) then
tmp = (t_m / l_m) * sqrt(x)
else if (l_m <= 8.5d+189) then
tmp = 1.0d0
else
tmp = (t_m * sqrt(x)) / l_m
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.3e+121) {
tmp = Math.sqrt(((x + -1.0) / (x + 1.0)));
} else if (l_m <= 4.2e+147) {
tmp = (t_m / l_m) * Math.sqrt(x);
} else if (l_m <= 8.5e+189) {
tmp = 1.0;
} else {
tmp = (t_m * Math.sqrt(x)) / l_m;
}
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.3e+121: tmp = math.sqrt(((x + -1.0) / (x + 1.0))) elif l_m <= 4.2e+147: tmp = (t_m / l_m) * math.sqrt(x) elif l_m <= 8.5e+189: tmp = 1.0 else: tmp = (t_m * math.sqrt(x)) / l_m 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.3e+121) tmp = sqrt(Float64(Float64(x + -1.0) / Float64(x + 1.0))); elseif (l_m <= 4.2e+147) tmp = Float64(Float64(t_m / l_m) * sqrt(x)); elseif (l_m <= 8.5e+189) tmp = 1.0; else tmp = Float64(Float64(t_m * sqrt(x)) / l_m); 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.3e+121) tmp = sqrt(((x + -1.0) / (x + 1.0))); elseif (l_m <= 4.2e+147) tmp = (t_m / l_m) * sqrt(x); elseif (l_m <= 8.5e+189) tmp = 1.0; else tmp = (t_m * sqrt(x)) / l_m; 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.3e+121], N[Sqrt[N[(N[(x + -1.0), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[l$95$m, 4.2e+147], N[(N[(t$95$m / l$95$m), $MachinePrecision] * N[Sqrt[x], $MachinePrecision]), $MachinePrecision], If[LessEqual[l$95$m, 8.5e+189], 1.0, N[(N[(t$95$m * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] / l$95$m), $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.3 \cdot 10^{+121}:\\
\;\;\;\;\sqrt{\frac{x + -1}{x + 1}}\\
\mathbf{elif}\;l\_m \leq 4.2 \cdot 10^{+147}:\\
\;\;\;\;\frac{t\_m}{l\_m} \cdot \sqrt{x}\\
\mathbf{elif}\;l\_m \leq 8.5 \cdot 10^{+189}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_m \cdot \sqrt{x}}{l\_m}\\
\end{array}
\end{array}
if l < 2.2999999999999999e121Initial program 41.7%
Simplified41.6%
Taylor expanded in l around 0 47.4%
associate-*l*47.5%
+-commutative47.5%
sub-neg47.5%
metadata-eval47.5%
+-commutative47.5%
Simplified47.5%
Taylor expanded in t around 0 47.5%
if 2.2999999999999999e121 < l < 4.20000000000000012e147Initial program 13.8%
Simplified13.8%
Taylor expanded in l around inf 14.0%
associate--l+29.0%
sub-neg29.0%
metadata-eval29.0%
+-commutative29.0%
sub-neg29.0%
metadata-eval29.0%
+-commutative29.0%
Simplified29.0%
Taylor expanded in x around inf 91.7%
Taylor expanded in t around 0 91.7%
if 4.20000000000000012e147 < l < 8.4999999999999998e189Initial program 0.0%
Simplified0.0%
Taylor expanded in l around 0 41.3%
associate-*l*41.3%
+-commutative41.3%
sub-neg41.3%
metadata-eval41.3%
+-commutative41.3%
Simplified41.3%
Taylor expanded in x around inf 41.3%
if 8.4999999999999998e189 < l Initial program 0.0%
Simplified0.0%
Taylor expanded in l around inf 0.0%
associate--l+34.5%
sub-neg34.5%
metadata-eval34.5%
+-commutative34.5%
sub-neg34.5%
metadata-eval34.5%
+-commutative34.5%
Simplified34.5%
Taylor expanded in x around inf 82.2%
clear-num82.0%
un-div-inv82.1%
associate-*l*82.3%
sqrt-div82.3%
metadata-eval82.3%
un-div-inv82.4%
Applied egg-rr82.4%
Taylor expanded in l around 0 71.3%
associate-*l/82.4%
Simplified82.4%
Final simplification50.9%
l_m = (fabs.f64 l) t\_m = (fabs.f64 t) t\_s = (copysign.f64 1 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 36.3%
Simplified36.2%
Taylor expanded in l around 0 44.1%
associate-*l*44.1%
+-commutative44.1%
sub-neg44.1%
metadata-eval44.1%
+-commutative44.1%
Simplified44.1%
Taylor expanded in x around -inf 0.0%
mul-1-neg0.0%
unsub-neg0.0%
Simplified43.5%
Final simplification43.5%
l_m = (fabs.f64 l) t\_m = (fabs.f64 t) t\_s = (copysign.f64 1 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 36.3%
Simplified36.2%
Taylor expanded in l around 0 44.1%
associate-*l*44.1%
+-commutative44.1%
sub-neg44.1%
metadata-eval44.1%
+-commutative44.1%
Simplified44.1%
Taylor expanded in x around inf 43.4%
Final simplification43.4%
l_m = (fabs.f64 l) t\_m = (fabs.f64 t) t\_s = (copysign.f64 1 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 36.3%
Simplified36.2%
Taylor expanded in l around 0 44.1%
associate-*l*44.1%
+-commutative44.1%
sub-neg44.1%
metadata-eval44.1%
+-commutative44.1%
Simplified44.1%
Taylor expanded in x around inf 42.9%
Final simplification42.9%
herbie shell --seed 2024051
(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)))))