
(FPCore (x y z t a b c i) :precision binary64 (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i);
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
code = ((((((((x * y) + z) * y) + 27464.7644705d0) * y) + 230661.510616d0) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i);
}
def code(x, y, z, t, a, b, c, i): return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(y + a) * y) + b) * y) + c) * y) + i)) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(N[(N[(N[(N[(N[(N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision] * y), $MachinePrecision] + 27464.7644705), $MachinePrecision] * y), $MachinePrecision] + 230661.510616), $MachinePrecision] * y), $MachinePrecision] + t), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(y + a), $MachinePrecision] * y), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] + c), $MachinePrecision] * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 22 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c i) :precision binary64 (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i);
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
code = ((((((((x * y) + z) * y) + 27464.7644705d0) * y) + 230661.510616d0) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i);
}
def code(x, y, z, t, a, b, c, i): return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(y + a) * y) + b) * y) + c) * y) + i)) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(N[(N[(N[(N[(N[(N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision] * y), $MachinePrecision] + 27464.7644705), $MachinePrecision] * y), $MachinePrecision] + 230661.510616), $MachinePrecision] * y), $MachinePrecision] + t), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(y + a), $MachinePrecision] * y), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] + c), $MachinePrecision] * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}
\end{array}
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (fma y (fma y (fma y (+ y a) b) c) i))
(t_2 (+ (/ z y) (/ 27464.7644705 (* y y)))))
(if (<= y -1.12e+77)
(+ t_2 (- x (/ (* x b) (* y y))))
(if (<= y 1.2e+43)
(fma
(fma y (fma y z 27464.7644705) 230661.510616)
(/ y t_1)
(fma x (/ (* y (* y (* y y))) t_1) (/ t t_1)))
(+ t_2 (- x (* a (+ (/ x y) (/ (- z (* x a)) (* y y))))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = fma(y, fma(y, fma(y, (y + a), b), c), i);
double t_2 = (z / y) + (27464.7644705 / (y * y));
double tmp;
if (y <= -1.12e+77) {
tmp = t_2 + (x - ((x * b) / (y * y)));
} else if (y <= 1.2e+43) {
tmp = fma(fma(y, fma(y, z, 27464.7644705), 230661.510616), (y / t_1), fma(x, ((y * (y * (y * y))) / t_1), (t / t_1)));
} else {
tmp = t_2 + (x - (a * ((x / y) + ((z - (x * a)) / (y * y)))));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(y, fma(y, fma(y, Float64(y + a), b), c), i) t_2 = Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) tmp = 0.0 if (y <= -1.12e+77) tmp = Float64(t_2 + Float64(x - Float64(Float64(x * b) / Float64(y * y)))); elseif (y <= 1.2e+43) tmp = fma(fma(y, fma(y, z, 27464.7644705), 230661.510616), Float64(y / t_1), fma(x, Float64(Float64(y * Float64(y * Float64(y * y))) / t_1), Float64(t / t_1))); else tmp = Float64(t_2 + Float64(x - Float64(a * Float64(Float64(x / y) + Float64(Float64(z - Float64(x * a)) / Float64(y * y)))))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.12e+77], N[(t$95$2 + N[(x - N[(N[(x * b), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.2e+43], N[(N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] * N[(y / t$95$1), $MachinePrecision] + N[(x * N[(N[(y * N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision] + N[(t / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$2 + N[(x - N[(a * N[(N[(x / y), $MachinePrecision] + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)\\
t_2 := \frac{z}{y} + \frac{27464.7644705}{y \cdot y}\\
\mathbf{if}\;y \leq -1.12 \cdot 10^{+77}:\\
\;\;\;\;t\_2 + \left(x - \frac{x \cdot b}{y \cdot y}\right)\\
\mathbf{elif}\;y \leq 1.2 \cdot 10^{+43}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right), \frac{y}{t\_1}, \mathsf{fma}\left(x, \frac{y \cdot \left(y \cdot \left(y \cdot y\right)\right)}{t\_1}, \frac{t}{t\_1}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2 + \left(x - a \cdot \left(\frac{x}{y} + \frac{z - x \cdot a}{y \cdot y}\right)\right)\\
\end{array}
\end{array}
if y < -1.1199999999999999e77Initial program 0.5%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites75.6%
Taylor expanded in a around 0
Applied rewrites79.6%
if -1.1199999999999999e77 < y < 1.20000000000000012e43Initial program 91.0%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites95.0%
Applied rewrites95.7%
if 1.20000000000000012e43 < y Initial program 5.5%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites63.0%
Taylor expanded in b around 0
Applied rewrites71.3%
Final simplification86.6%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (fma y (fma y (fma y (+ y a) b) c) i))
(t_2 (+ (/ z y) (/ 27464.7644705 (* y y)))))
(if (<= y -1.6e+91)
(+ t_2 (- x (/ (* x a) y)))
(if (<= y 1.2e+43)
(fma
y
(/ (fma y (fma y z 27464.7644705) 230661.510616) t_1)
(fma x (* (* y y) (/ (* y y) t_1)) (/ t t_1)))
(+ t_2 (- x (* a (+ (/ x y) (/ (- z (* x a)) (* y y))))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = fma(y, fma(y, fma(y, (y + a), b), c), i);
double t_2 = (z / y) + (27464.7644705 / (y * y));
double tmp;
if (y <= -1.6e+91) {
tmp = t_2 + (x - ((x * a) / y));
} else if (y <= 1.2e+43) {
tmp = fma(y, (fma(y, fma(y, z, 27464.7644705), 230661.510616) / t_1), fma(x, ((y * y) * ((y * y) / t_1)), (t / t_1)));
} else {
tmp = t_2 + (x - (a * ((x / y) + ((z - (x * a)) / (y * y)))));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(y, fma(y, fma(y, Float64(y + a), b), c), i) t_2 = Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) tmp = 0.0 if (y <= -1.6e+91) tmp = Float64(t_2 + Float64(x - Float64(Float64(x * a) / y))); elseif (y <= 1.2e+43) tmp = fma(y, Float64(fma(y, fma(y, z, 27464.7644705), 230661.510616) / t_1), fma(x, Float64(Float64(y * y) * Float64(Float64(y * y) / t_1)), Float64(t / t_1))); else tmp = Float64(t_2 + Float64(x - Float64(a * Float64(Float64(x / y) + Float64(Float64(z - Float64(x * a)) / Float64(y * y)))))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.6e+91], N[(t$95$2 + N[(x - N[(N[(x * a), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.2e+43], N[(y * N[(N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] / t$95$1), $MachinePrecision] + N[(x * N[(N[(y * y), $MachinePrecision] * N[(N[(y * y), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision] + N[(t / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$2 + N[(x - N[(a * N[(N[(x / y), $MachinePrecision] + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)\\
t_2 := \frac{z}{y} + \frac{27464.7644705}{y \cdot y}\\
\mathbf{if}\;y \leq -1.6 \cdot 10^{+91}:\\
\;\;\;\;t\_2 + \left(x - \frac{x \cdot a}{y}\right)\\
\mathbf{elif}\;y \leq 1.2 \cdot 10^{+43}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right)}{t\_1}, \mathsf{fma}\left(x, \left(y \cdot y\right) \cdot \frac{y \cdot y}{t\_1}, \frac{t}{t\_1}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2 + \left(x - a \cdot \left(\frac{x}{y} + \frac{z - x \cdot a}{y \cdot y}\right)\right)\\
\end{array}
\end{array}
if y < -1.59999999999999995e91Initial program 0.3%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites79.9%
Taylor expanded in y around inf
Applied rewrites84.2%
if -1.59999999999999995e91 < y < 1.20000000000000012e43Initial program 89.2%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites93.0%
Applied rewrites93.1%
if 1.20000000000000012e43 < y Initial program 5.5%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites63.0%
Taylor expanded in b around 0
Applied rewrites71.3%
Final simplification86.3%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (fma y (fma y (fma y (+ y a) b) c) i))
(t_2 (+ (/ z y) (/ 27464.7644705 (* y y)))))
(if (<= y -4.8e+62)
(+ t_2 (- x (/ (* x b) (* y y))))
(if (<= y 1.12e+43)
(fma
y
(/ (fma y (fma y (fma y x z) 27464.7644705) 230661.510616) t_1)
(/ t t_1))
(+ t_2 (- x (* a (+ (/ x y) (/ (- z (* x a)) (* y y))))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = fma(y, fma(y, fma(y, (y + a), b), c), i);
double t_2 = (z / y) + (27464.7644705 / (y * y));
double tmp;
if (y <= -4.8e+62) {
tmp = t_2 + (x - ((x * b) / (y * y)));
} else if (y <= 1.12e+43) {
tmp = fma(y, (fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616) / t_1), (t / t_1));
} else {
tmp = t_2 + (x - (a * ((x / y) + ((z - (x * a)) / (y * y)))));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(y, fma(y, fma(y, Float64(y + a), b), c), i) t_2 = Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) tmp = 0.0 if (y <= -4.8e+62) tmp = Float64(t_2 + Float64(x - Float64(Float64(x * b) / Float64(y * y)))); elseif (y <= 1.12e+43) tmp = fma(y, Float64(fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616) / t_1), Float64(t / t_1)); else tmp = Float64(t_2 + Float64(x - Float64(a * Float64(Float64(x / y) + Float64(Float64(z - Float64(x * a)) / Float64(y * y)))))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4.8e+62], N[(t$95$2 + N[(x - N[(N[(x * b), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.12e+43], N[(y * N[(N[(y * N[(y * N[(y * x + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] / t$95$1), $MachinePrecision] + N[(t / t$95$1), $MachinePrecision]), $MachinePrecision], N[(t$95$2 + N[(x - N[(a * N[(N[(x / y), $MachinePrecision] + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)\\
t_2 := \frac{z}{y} + \frac{27464.7644705}{y \cdot y}\\
\mathbf{if}\;y \leq -4.8 \cdot 10^{+62}:\\
\;\;\;\;t\_2 + \left(x - \frac{x \cdot b}{y \cdot y}\right)\\
\mathbf{elif}\;y \leq 1.12 \cdot 10^{+43}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, x, z\right), 27464.7644705\right), 230661.510616\right)}{t\_1}, \frac{t}{t\_1}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2 + \left(x - a \cdot \left(\frac{x}{y} + \frac{z - x \cdot a}{y \cdot y}\right)\right)\\
\end{array}
\end{array}
if y < -4.8e62Initial program 0.5%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites75.6%
Taylor expanded in a around 0
Applied rewrites79.6%
if -4.8e62 < y < 1.12e43Initial program 91.0%
Taylor expanded in t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites91.7%
if 1.12e43 < y Initial program 5.5%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites63.0%
Taylor expanded in b around 0
Applied rewrites71.3%
Final simplification84.4%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (/ z y) (/ 27464.7644705 (* y y)))))
(if (<= y -4e+60)
(+ t_1 (- x (/ (* x b) (* y y))))
(if (<= y 1.25e+20)
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (fma (* y (fma y (+ y a) b)) y (* y c))))
(+ t_1 (- x (* a (+ (/ x y) (/ (- z (* x a)) (* y y))))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (z / y) + (27464.7644705 / (y * y));
double tmp;
if (y <= -4e+60) {
tmp = t_1 + (x - ((x * b) / (y * y)));
} else if (y <= 1.25e+20) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + fma((y * fma(y, (y + a), b)), y, (y * c)));
} else {
tmp = t_1 + (x - (a * ((x / y) + ((z - (x * a)) / (y * y)))));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) tmp = 0.0 if (y <= -4e+60) tmp = Float64(t_1 + Float64(x - Float64(Float64(x * b) / Float64(y * y)))); elseif (y <= 1.25e+20) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + fma(Float64(y * fma(y, Float64(y + a), b)), y, Float64(y * c)))); else tmp = Float64(t_1 + Float64(x - Float64(a * Float64(Float64(x / y) + Float64(Float64(z - Float64(x * a)) / Float64(y * y)))))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4e+60], N[(t$95$1 + N[(x - N[(N[(x * b), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.25e+20], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] * y + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 + N[(x - N[(a * N[(N[(x / y), $MachinePrecision] + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z}{y} + \frac{27464.7644705}{y \cdot y}\\
\mathbf{if}\;y \leq -4 \cdot 10^{+60}:\\
\;\;\;\;t\_1 + \left(x - \frac{x \cdot b}{y \cdot y}\right)\\
\mathbf{elif}\;y \leq 1.25 \cdot 10^{+20}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + \mathsf{fma}\left(y \cdot \mathsf{fma}\left(y, y + a, b\right), y, y \cdot c\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1 + \left(x - a \cdot \left(\frac{x}{y} + \frac{z - x \cdot a}{y \cdot y}\right)\right)\\
\end{array}
\end{array}
if y < -3.9999999999999998e60Initial program 0.5%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites75.6%
Taylor expanded in a around 0
Applied rewrites79.6%
if -3.9999999999999998e60 < y < 1.25e20Initial program 92.8%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f6492.8
Applied rewrites92.8%
if 1.25e20 < y Initial program 6.9%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites62.3%
Taylor expanded in b around 0
Applied rewrites70.1%
Final simplification84.4%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
INFINITY)
(/ t i)
(/ z y)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= ((double) INFINITY)) {
tmp = t / i;
} else {
tmp = z / y;
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= Double.POSITIVE_INFINITY) {
tmp = t / i;
} else {
tmp = z / y;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if ((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= math.inf: tmp = t / i else: tmp = z / y return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))) <= Inf) tmp = Float64(t / i); else tmp = Float64(z / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (((t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))) <= Inf) tmp = t / i; else tmp = z / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(t / i), $MachinePrecision], N[(z / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)} \leq \infty:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;\frac{z}{y}\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < +inf.0Initial program 86.2%
Taylor expanded in y around 0
lower-/.f6444.0
Applied rewrites44.0%
if +inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 0.0%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f641.2
Applied rewrites1.2%
Taylor expanded in y around inf
Applied rewrites28.3%
Final simplification37.8%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (/ z y) (/ 27464.7644705 (* y y)))))
(if (<= y -4e+60)
(+ t_1 (- x (/ (* x b) (* y y))))
(if (<= y 1.25e+20)
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (fma (* y (fma y (+ y a) b)) y (* y c))))
(+ t_1 (- x (/ (* x a) y)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (z / y) + (27464.7644705 / (y * y));
double tmp;
if (y <= -4e+60) {
tmp = t_1 + (x - ((x * b) / (y * y)));
} else if (y <= 1.25e+20) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + fma((y * fma(y, (y + a), b)), y, (y * c)));
} else {
tmp = t_1 + (x - ((x * a) / y));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) tmp = 0.0 if (y <= -4e+60) tmp = Float64(t_1 + Float64(x - Float64(Float64(x * b) / Float64(y * y)))); elseif (y <= 1.25e+20) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + fma(Float64(y * fma(y, Float64(y + a), b)), y, Float64(y * c)))); else tmp = Float64(t_1 + Float64(x - Float64(Float64(x * a) / y))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4e+60], N[(t$95$1 + N[(x - N[(N[(x * b), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.25e+20], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] * y + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 + N[(x - N[(N[(x * a), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z}{y} + \frac{27464.7644705}{y \cdot y}\\
\mathbf{if}\;y \leq -4 \cdot 10^{+60}:\\
\;\;\;\;t\_1 + \left(x - \frac{x \cdot b}{y \cdot y}\right)\\
\mathbf{elif}\;y \leq 1.25 \cdot 10^{+20}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + \mathsf{fma}\left(y \cdot \mathsf{fma}\left(y, y + a, b\right), y, y \cdot c\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1 + \left(x - \frac{x \cdot a}{y}\right)\\
\end{array}
\end{array}
if y < -3.9999999999999998e60Initial program 0.5%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites75.6%
Taylor expanded in a around 0
Applied rewrites79.6%
if -3.9999999999999998e60 < y < 1.25e20Initial program 92.8%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f6492.8
Applied rewrites92.8%
if 1.25e20 < y Initial program 6.9%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites62.3%
Taylor expanded in y around inf
Applied rewrites70.0%
Final simplification84.4%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (/ z y) (/ 27464.7644705 (* y y)))))
(if (<= y -4e+60)
(+ t_1 (- x (/ (* x b) (* y y))))
(if (<= y 1.25e+20)
(/
(+
t
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ z (* y x))))))))
(+ i (* y (+ c (* y (+ b (* y (+ y a))))))))
(+ t_1 (- x (/ (* x a) y)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (z / y) + (27464.7644705 / (y * y));
double tmp;
if (y <= -4e+60) {
tmp = t_1 + (x - ((x * b) / (y * y)));
} else if (y <= 1.25e+20) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))));
} else {
tmp = t_1 + (x - ((x * a) / y));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = (z / y) + (27464.7644705d0 / (y * y))
if (y <= (-4d+60)) then
tmp = t_1 + (x - ((x * b) / (y * y)))
else if (y <= 1.25d+20) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))))
else
tmp = t_1 + (x - ((x * a) / y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (z / y) + (27464.7644705 / (y * y));
double tmp;
if (y <= -4e+60) {
tmp = t_1 + (x - ((x * b) / (y * y)));
} else if (y <= 1.25e+20) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a)))))));
} else {
tmp = t_1 + (x - ((x * a) / y));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (z / y) + (27464.7644705 / (y * y)) tmp = 0 if y <= -4e+60: tmp = t_1 + (x - ((x * b) / (y * y))) elif y <= 1.25e+20: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a))))))) else: tmp = t_1 + (x - ((x * a) / y)) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) tmp = 0.0 if (y <= -4e+60) tmp = Float64(t_1 + Float64(x - Float64(Float64(x * b) / Float64(y * y)))); elseif (y <= 1.25e+20) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(z + Float64(y * x)))))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * Float64(y + a)))))))); else tmp = Float64(t_1 + Float64(x - Float64(Float64(x * a) / y))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = (z / y) + (27464.7644705 / (y * y)); tmp = 0.0; if (y <= -4e+60) tmp = t_1 + (x - ((x * b) / (y * y))); elseif (y <= 1.25e+20) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * (z + (y * x)))))))) / (i + (y * (c + (y * (b + (y * (y + a))))))); else tmp = t_1 + (x - ((x * a) / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4e+60], N[(t$95$1 + N[(x - N[(N[(x * b), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.25e+20], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(z + N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 + N[(x - N[(N[(x * a), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z}{y} + \frac{27464.7644705}{y \cdot y}\\
\mathbf{if}\;y \leq -4 \cdot 10^{+60}:\\
\;\;\;\;t\_1 + \left(x - \frac{x \cdot b}{y \cdot y}\right)\\
\mathbf{elif}\;y \leq 1.25 \cdot 10^{+20}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(z + y \cdot x\right)\right)\right)}{i + y \cdot \left(c + y \cdot \left(b + y \cdot \left(y + a\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1 + \left(x - \frac{x \cdot a}{y}\right)\\
\end{array}
\end{array}
if y < -3.9999999999999998e60Initial program 0.5%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites75.6%
Taylor expanded in a around 0
Applied rewrites79.6%
if -3.9999999999999998e60 < y < 1.25e20Initial program 92.8%
if 1.25e20 < y Initial program 6.9%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites62.3%
Taylor expanded in y around inf
Applied rewrites70.0%
Final simplification84.4%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (/ z y) (/ 27464.7644705 (* y y)))))
(if (<= y -4e+60)
(+ t_1 (- x (/ (* x b) (* y y))))
(if (<= y 1.25e+20)
(*
(fma y (fma y (fma y (fma x y z) 27464.7644705) 230661.510616) t)
(/ 1.0 (fma y (fma y (fma y (+ y a) b) c) i)))
(+ t_1 (- x (/ (* x a) y)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (z / y) + (27464.7644705 / (y * y));
double tmp;
if (y <= -4e+60) {
tmp = t_1 + (x - ((x * b) / (y * y)));
} else if (y <= 1.25e+20) {
tmp = fma(y, fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616), t) * (1.0 / fma(y, fma(y, fma(y, (y + a), b), c), i));
} else {
tmp = t_1 + (x - ((x * a) / y));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) tmp = 0.0 if (y <= -4e+60) tmp = Float64(t_1 + Float64(x - Float64(Float64(x * b) / Float64(y * y)))); elseif (y <= 1.25e+20) tmp = Float64(fma(y, fma(y, fma(y, fma(x, y, z), 27464.7644705), 230661.510616), t) * Float64(1.0 / fma(y, fma(y, fma(y, Float64(y + a), b), c), i))); else tmp = Float64(t_1 + Float64(x - Float64(Float64(x * a) / y))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4e+60], N[(t$95$1 + N[(x - N[(N[(x * b), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.25e+20], N[(N[(y * N[(y * N[(y * N[(x * y + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] * N[(1.0 / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 + N[(x - N[(N[(x * a), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z}{y} + \frac{27464.7644705}{y \cdot y}\\
\mathbf{if}\;y \leq -4 \cdot 10^{+60}:\\
\;\;\;\;t\_1 + \left(x - \frac{x \cdot b}{y \cdot y}\right)\\
\mathbf{elif}\;y \leq 1.25 \cdot 10^{+20}:\\
\;\;\;\;\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(x, y, z\right), 27464.7644705\right), 230661.510616\right), t\right) \cdot \frac{1}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1 + \left(x - \frac{x \cdot a}{y}\right)\\
\end{array}
\end{array}
if y < -3.9999999999999998e60Initial program 0.5%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites75.6%
Taylor expanded in a around 0
Applied rewrites79.6%
if -3.9999999999999998e60 < y < 1.25e20Initial program 92.8%
lift-/.f64N/A
div-invN/A
lower-*.f64N/A
Applied rewrites92.7%
if 1.25e20 < y Initial program 6.9%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites62.3%
Taylor expanded in y around inf
Applied rewrites70.0%
Final simplification84.3%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (/ z y) (/ 27464.7644705 (* y y)))))
(if (<= y -4e+60)
(+ t_1 (- x (/ (* x b) (* y y))))
(if (<= y -5e-44)
(/
(fma y (fma y (fma y (fma y x z) 27464.7644705) 230661.510616) t)
(fma y (fma (* y y) (+ y a) c) i))
(if (<= y 4.8e+14)
(/
(fma y (fma y (fma y z 27464.7644705) 230661.510616) t)
(fma y (fma y (fma y (+ y a) b) c) i))
(+ t_1 (- x (/ (* x a) y))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (z / y) + (27464.7644705 / (y * y));
double tmp;
if (y <= -4e+60) {
tmp = t_1 + (x - ((x * b) / (y * y)));
} else if (y <= -5e-44) {
tmp = fma(y, fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616), t) / fma(y, fma((y * y), (y + a), c), i);
} else if (y <= 4.8e+14) {
tmp = fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, (y + a), b), c), i);
} else {
tmp = t_1 + (x - ((x * a) / y));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) tmp = 0.0 if (y <= -4e+60) tmp = Float64(t_1 + Float64(x - Float64(Float64(x * b) / Float64(y * y)))); elseif (y <= -5e-44) tmp = Float64(fma(y, fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616), t) / fma(y, fma(Float64(y * y), Float64(y + a), c), i)); elseif (y <= 4.8e+14) tmp = Float64(fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)); else tmp = Float64(t_1 + Float64(x - Float64(Float64(x * a) / y))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4e+60], N[(t$95$1 + N[(x - N[(N[(x * b), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -5e-44], N[(N[(y * N[(y * N[(y * N[(y * x + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(N[(y * y), $MachinePrecision] * N[(y + a), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.8e+14], N[(N[(y * N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], N[(t$95$1 + N[(x - N[(N[(x * a), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z}{y} + \frac{27464.7644705}{y \cdot y}\\
\mathbf{if}\;y \leq -4 \cdot 10^{+60}:\\
\;\;\;\;t\_1 + \left(x - \frac{x \cdot b}{y \cdot y}\right)\\
\mathbf{elif}\;y \leq -5 \cdot 10^{-44}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, x, z\right), 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y \cdot y, y + a, c\right), i\right)}\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{+14}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1 + \left(x - \frac{x \cdot a}{y}\right)\\
\end{array}
\end{array}
if y < -3.9999999999999998e60Initial program 0.5%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites75.6%
Taylor expanded in a around 0
Applied rewrites79.6%
if -3.9999999999999998e60 < y < -5.00000000000000039e-44Initial program 68.8%
Taylor expanded in b around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f6456.4
Applied rewrites56.4%
if -5.00000000000000039e-44 < y < 4.8e14Initial program 99.7%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6498.8
Applied rewrites98.8%
if 4.8e14 < y Initial program 8.3%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites61.3%
Taylor expanded in y around inf
Applied rewrites69.1%
Final simplification82.2%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (+ (/ z y) (/ 27464.7644705 (* y y))) (- x (/ (* x a) y)))))
(if (<= y -4e+60)
t_1
(if (<= y -5e-44)
(/
(fma y (fma y (fma y (fma y x z) 27464.7644705) 230661.510616) t)
(fma y (fma (* y y) (+ y a) c) i))
(if (<= y 4.8e+14)
(/
(fma y (fma y (fma y z 27464.7644705) 230661.510616) t)
(fma y (fma y (fma y (+ y a) b) c) i))
t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = ((z / y) + (27464.7644705 / (y * y))) + (x - ((x * a) / y));
double tmp;
if (y <= -4e+60) {
tmp = t_1;
} else if (y <= -5e-44) {
tmp = fma(y, fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616), t) / fma(y, fma((y * y), (y + a), c), i);
} else if (y <= 4.8e+14) {
tmp = fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, (y + a), b), c), i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) + Float64(x - Float64(Float64(x * a) / y))) tmp = 0.0 if (y <= -4e+60) tmp = t_1; elseif (y <= -5e-44) tmp = Float64(fma(y, fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616), t) / fma(y, fma(Float64(y * y), Float64(y + a), c), i)); elseif (y <= 4.8e+14) tmp = Float64(fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x - N[(N[(x * a), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4e+60], t$95$1, If[LessEqual[y, -5e-44], N[(N[(y * N[(y * N[(y * N[(y * x + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(N[(y * y), $MachinePrecision] * N[(y + a), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.8e+14], N[(N[(y * N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\frac{z}{y} + \frac{27464.7644705}{y \cdot y}\right) + \left(x - \frac{x \cdot a}{y}\right)\\
\mathbf{if}\;y \leq -4 \cdot 10^{+60}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -5 \cdot 10^{-44}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, x, z\right), 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y \cdot y, y + a, c\right), i\right)}\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{+14}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -3.9999999999999998e60 or 4.8e14 < y Initial program 4.8%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites67.7%
Taylor expanded in y around inf
Applied rewrites73.8%
if -3.9999999999999998e60 < y < -5.00000000000000039e-44Initial program 68.8%
Taylor expanded in b around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f6456.4
Applied rewrites56.4%
if -5.00000000000000039e-44 < y < 4.8e14Initial program 99.7%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6498.8
Applied rewrites98.8%
Final simplification82.1%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (+ (/ z y) (/ 27464.7644705 (* y y))) (- x (/ (* x a) y)))))
(if (<= y -3e+60)
t_1
(if (<= y -2.25e-42)
(/
(fma y (fma y (fma y (fma y x z) 27464.7644705) 230661.510616) t)
(fma y (fma y (fma y y b) c) i))
(if (<= y 4.8e+14)
(/
(fma y (fma y (fma y z 27464.7644705) 230661.510616) t)
(fma y (fma y (fma y (+ y a) b) c) i))
t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = ((z / y) + (27464.7644705 / (y * y))) + (x - ((x * a) / y));
double tmp;
if (y <= -3e+60) {
tmp = t_1;
} else if (y <= -2.25e-42) {
tmp = fma(y, fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, y, b), c), i);
} else if (y <= 4.8e+14) {
tmp = fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, (y + a), b), c), i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) + Float64(x - Float64(Float64(x * a) / y))) tmp = 0.0 if (y <= -3e+60) tmp = t_1; elseif (y <= -2.25e-42) tmp = Float64(fma(y, fma(y, fma(y, fma(y, x, z), 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, y, b), c), i)); elseif (y <= 4.8e+14) tmp = Float64(fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x - N[(N[(x * a), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3e+60], t$95$1, If[LessEqual[y, -2.25e-42], N[(N[(y * N[(y * N[(y * N[(y * x + z), $MachinePrecision] + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(y * N[(y * y + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.8e+14], N[(N[(y * N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\frac{z}{y} + \frac{27464.7644705}{y \cdot y}\right) + \left(x - \frac{x \cdot a}{y}\right)\\
\mathbf{if}\;y \leq -3 \cdot 10^{+60}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -2.25 \cdot 10^{-42}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, x, z\right), 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y, b\right), c\right), i\right)}\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{+14}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.9999999999999998e60 or 4.8e14 < y Initial program 4.8%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites67.7%
Taylor expanded in y around inf
Applied rewrites73.8%
if -2.9999999999999998e60 < y < -2.25e-42Initial program 68.8%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
unpow2N/A
lower-fma.f6456.0
Applied rewrites56.0%
if -2.25e-42 < y < 4.8e14Initial program 99.7%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6498.8
Applied rewrites98.8%
Final simplification82.1%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ (- z (* x a)) y))))
(if (<= y -1.6e+91)
t_1
(if (<= y 6.4e-21)
(/
(fma y (fma y (fma y z 27464.7644705) 230661.510616) t)
(fma y (+ c (* y b)) i))
(if (<= y 4.8e+14)
(/ (fma y (* z (* y y)) t) (fma y (fma y (fma y (+ y a) b) c) i))
t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -1.6e+91) {
tmp = t_1;
} else if (y <= 6.4e-21) {
tmp = fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, (c + (y * b)), i);
} else if (y <= 4.8e+14) {
tmp = fma(y, (z * (y * y)), t) / fma(y, fma(y, fma(y, (y + a), b), c), i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -1.6e+91) tmp = t_1; elseif (y <= 6.4e-21) tmp = Float64(fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, Float64(c + Float64(y * b)), i)); elseif (y <= 4.8e+14) tmp = Float64(fma(y, Float64(z * Float64(y * y)), t) / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.6e+91], t$95$1, If[LessEqual[y, 6.4e-21], N[(N[(y * N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.8e+14], N[(N[(y * N[(z * N[(y * y), $MachinePrecision]), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -1.6 \cdot 10^{+91}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 6.4 \cdot 10^{-21}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, c + y \cdot b, i\right)}\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{+14}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, z \cdot \left(y \cdot y\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.59999999999999995e91 or 4.8e14 < y Initial program 4.8%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
associate-*r*N/A
mul-1-negN/A
cancel-sign-subN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6473.6
Applied rewrites73.6%
if -1.59999999999999995e91 < y < 6.4000000000000003e-21Initial program 90.2%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6483.5
Applied rewrites83.5%
Taylor expanded in y around 0
Applied rewrites79.5%
if 6.4000000000000003e-21 < y < 4.8e14Initial program 99.7%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6488.7
Applied rewrites88.7%
Taylor expanded in y around inf
Applied rewrites88.5%
Final simplification77.1%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (+ (/ z y) (/ 27464.7644705 (* y y))) (- x (/ (* x a) y)))))
(if (<= y -1.6e+91)
t_1
(if (<= y 4.8e+14)
(/
(fma y (fma y (fma y z 27464.7644705) 230661.510616) t)
(fma y (fma y (fma y (+ y a) b) c) i))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = ((z / y) + (27464.7644705 / (y * y))) + (x - ((x * a) / y));
double tmp;
if (y <= -1.6e+91) {
tmp = t_1;
} else if (y <= 4.8e+14) {
tmp = fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, (y + a), b), c), i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(Float64(z / y) + Float64(27464.7644705 / Float64(y * y))) + Float64(x - Float64(Float64(x * a) / y))) tmp = 0.0 if (y <= -1.6e+91) tmp = t_1; elseif (y <= 4.8e+14) tmp = Float64(fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(N[(z / y), $MachinePrecision] + N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x - N[(N[(x * a), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.6e+91], t$95$1, If[LessEqual[y, 4.8e+14], N[(N[(y * N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\frac{z}{y} + \frac{27464.7644705}{y \cdot y}\right) + \left(x - \frac{x \cdot a}{y}\right)\\
\mathbf{if}\;y \leq -1.6 \cdot 10^{+91}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{+14}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.59999999999999995e91 or 4.8e14 < y Initial program 4.8%
Taylor expanded in y around inf
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower--.f64N/A
associate-+r+N/A
Applied rewrites69.4%
Taylor expanded in y around inf
Applied rewrites75.7%
if -1.59999999999999995e91 < y < 4.8e14Initial program 90.8%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6483.8
Applied rewrites83.8%
Final simplification80.2%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ (- z (* x a)) y))))
(if (<= y -1.6e+91)
t_1
(if (<= y 4.8e+14)
(/
(fma y (fma y (fma y z 27464.7644705) 230661.510616) t)
(fma y (fma y (fma y (+ y a) b) c) i))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -1.6e+91) {
tmp = t_1;
} else if (y <= 4.8e+14) {
tmp = fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, (y + a), b), c), i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -1.6e+91) tmp = t_1; elseif (y <= 4.8e+14) tmp = Float64(fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.6e+91], t$95$1, If[LessEqual[y, 4.8e+14], N[(N[(y * N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -1.6 \cdot 10^{+91}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{+14}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.59999999999999995e91 or 4.8e14 < y Initial program 4.8%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
associate-*r*N/A
mul-1-negN/A
cancel-sign-subN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6473.6
Applied rewrites73.6%
if -1.59999999999999995e91 < y < 4.8e14Initial program 90.8%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6483.8
Applied rewrites83.8%
Final simplification79.2%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ (- z (* x a)) y))))
(if (<= y -1.6e+91)
t_1
(if (<= y 4.8e+14)
(/
(fma y (fma y (fma y z 27464.7644705) 230661.510616) t)
(fma y (+ c (* y b)) i))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -1.6e+91) {
tmp = t_1;
} else if (y <= 4.8e+14) {
tmp = fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, (c + (y * b)), i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -1.6e+91) tmp = t_1; elseif (y <= 4.8e+14) tmp = Float64(fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / fma(y, Float64(c + Float64(y * b)), i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.6e+91], t$95$1, If[LessEqual[y, 4.8e+14], N[(N[(y * N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -1.6 \cdot 10^{+91}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{+14}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right), t\right)}{\mathsf{fma}\left(y, c + y \cdot b, i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.59999999999999995e91 or 4.8e14 < y Initial program 4.8%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
associate-*r*N/A
mul-1-negN/A
cancel-sign-subN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6473.6
Applied rewrites73.6%
if -1.59999999999999995e91 < y < 4.8e14Initial program 90.8%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6483.8
Applied rewrites83.8%
Taylor expanded in y around 0
Applied rewrites76.6%
Final simplification75.3%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ (- z (* x a)) y))))
(if (<= y -8.8e+43)
t_1
(if (<= y 3.6e+14)
(/ (fma y 230661.510616 t) (fma y (fma y (fma y (+ y a) b) c) i))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -8.8e+43) {
tmp = t_1;
} else if (y <= 3.6e+14) {
tmp = fma(y, 230661.510616, t) / fma(y, fma(y, fma(y, (y + a), b), c), i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -8.8e+43) tmp = t_1; elseif (y <= 3.6e+14) tmp = Float64(fma(y, 230661.510616, t) / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -8.8e+43], t$95$1, If[LessEqual[y, 3.6e+14], N[(N[(y * 230661.510616 + t), $MachinePrecision] / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -8.8 \cdot 10^{+43}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.6 \cdot 10^{+14}:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, 230661.510616, t\right)}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -8.80000000000000002e43 or 3.6e14 < y Initial program 6.3%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
associate-*r*N/A
mul-1-negN/A
cancel-sign-subN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6470.3
Applied rewrites70.3%
if -8.80000000000000002e43 < y < 3.6e14Initial program 94.0%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6487.3
Applied rewrites87.3%
Taylor expanded in y around 0
Applied rewrites76.8%
Final simplification73.7%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ (- z (* x a)) y))))
(if (<= y -8.2e+43)
t_1
(if (<= y 0.00095)
(/
(fma y (fma y (fma y z 27464.7644705) 230661.510616) t)
(+ i (* y c)))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -8.2e+43) {
tmp = t_1;
} else if (y <= 0.00095) {
tmp = fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / (i + (y * c));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -8.2e+43) tmp = t_1; elseif (y <= 0.00095) tmp = Float64(fma(y, fma(y, fma(y, z, 27464.7644705), 230661.510616), t) / Float64(i + Float64(y * c))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -8.2e+43], t$95$1, If[LessEqual[y, 0.00095], N[(N[(y * N[(y * N[(y * z + 27464.7644705), $MachinePrecision] + 230661.510616), $MachinePrecision] + t), $MachinePrecision] / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -8.2 \cdot 10^{+43}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 0.00095:\\
\;\;\;\;\frac{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, z, 27464.7644705\right), 230661.510616\right), t\right)}{i + y \cdot c}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -8.2000000000000001e43 or 9.49999999999999998e-4 < y Initial program 11.4%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
associate-*r*N/A
mul-1-negN/A
cancel-sign-subN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6466.7
Applied rewrites66.7%
if -8.2000000000000001e43 < y < 9.49999999999999998e-4Initial program 93.7%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6487.4
Applied rewrites87.4%
Taylor expanded in y around 0
Applied rewrites74.4%
Final simplification70.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ (- z (* x a)) y))))
(if (<= y -8.8e+43)
t_1
(if (<= y 3.6e+14) (/ t (fma y (fma y (fma y (+ y a) b) c) i)) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -8.8e+43) {
tmp = t_1;
} else if (y <= 3.6e+14) {
tmp = t / fma(y, fma(y, fma(y, (y + a), b), c), i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -8.8e+43) tmp = t_1; elseif (y <= 3.6e+14) tmp = Float64(t / fma(y, fma(y, fma(y, Float64(y + a), b), c), i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -8.8e+43], t$95$1, If[LessEqual[y, 3.6e+14], N[(t / N[(y * N[(y * N[(y * N[(y + a), $MachinePrecision] + b), $MachinePrecision] + c), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -8.8 \cdot 10^{+43}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.6 \cdot 10^{+14}:\\
\;\;\;\;\frac{t}{\mathsf{fma}\left(y, \mathsf{fma}\left(y, \mathsf{fma}\left(y, y + a, b\right), c\right), i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -8.80000000000000002e43 or 3.6e14 < y Initial program 6.3%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
associate-*r*N/A
mul-1-negN/A
cancel-sign-subN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6470.3
Applied rewrites70.3%
if -8.80000000000000002e43 < y < 3.6e14Initial program 94.0%
Taylor expanded in t around inf
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6467.2
Applied rewrites67.2%
Final simplification68.7%
(FPCore (x y z t a b c i) :precision binary64 (let* ((t_1 (+ x (/ (- z (* x a)) y)))) (if (<= y -5.8e+43) t_1 (if (<= y 1.15e-17) (/ t i) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -5.8e+43) {
tmp = t_1;
} else if (y <= 1.15e-17) {
tmp = t / i;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z - (x * a)) / y)
if (y <= (-5.8d+43)) then
tmp = t_1
else if (y <= 1.15d-17) then
tmp = t / i
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -5.8e+43) {
tmp = t_1;
} else if (y <= 1.15e-17) {
tmp = t / i;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z - (x * a)) / y) tmp = 0 if y <= -5.8e+43: tmp = t_1 elif y <= 1.15e-17: tmp = t / i else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -5.8e+43) tmp = t_1; elseif (y <= 1.15e-17) tmp = Float64(t / i); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z - (x * a)) / y); tmp = 0.0; if (y <= -5.8e+43) tmp = t_1; elseif (y <= 1.15e-17) tmp = t / i; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -5.8e+43], t$95$1, If[LessEqual[y, 1.15e-17], N[(t / i), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -5.8 \cdot 10^{+43}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.15 \cdot 10^{-17}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.8000000000000004e43 or 1.15000000000000004e-17 < y Initial program 12.7%
Taylor expanded in y around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
associate-*r*N/A
mul-1-negN/A
cancel-sign-subN/A
lower-+.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6465.7
Applied rewrites65.7%
if -5.8000000000000004e43 < y < 1.15000000000000004e-17Initial program 93.6%
Taylor expanded in y around 0
lower-/.f6453.8
Applied rewrites53.8%
Final simplification59.9%
(FPCore (x y z t a b c i) :precision binary64 (let* ((t_1 (/ 1.0 (/ 1.0 x)))) (if (<= y -3.55e+44) t_1 (if (<= y 1.75e-6) (/ t i) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = 1.0 / (1.0 / x);
double tmp;
if (y <= -3.55e+44) {
tmp = t_1;
} else if (y <= 1.75e-6) {
tmp = t / i;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = 1.0d0 / (1.0d0 / x)
if (y <= (-3.55d+44)) then
tmp = t_1
else if (y <= 1.75d-6) then
tmp = t / i
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = 1.0 / (1.0 / x);
double tmp;
if (y <= -3.55e+44) {
tmp = t_1;
} else if (y <= 1.75e-6) {
tmp = t / i;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = 1.0 / (1.0 / x) tmp = 0 if y <= -3.55e+44: tmp = t_1 elif y <= 1.75e-6: tmp = t / i else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(1.0 / Float64(1.0 / x)) tmp = 0.0 if (y <= -3.55e+44) tmp = t_1; elseif (y <= 1.75e-6) tmp = Float64(t / i); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = 1.0 / (1.0 / x); tmp = 0.0; if (y <= -3.55e+44) tmp = t_1; elseif (y <= 1.75e-6) tmp = t / i; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(1.0 / N[(1.0 / x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.55e+44], t$95$1, If[LessEqual[y, 1.75e-6], N[(t / i), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{1}{\frac{1}{x}}\\
\mathbf{if}\;y \leq -3.55 \cdot 10^{+44}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.75 \cdot 10^{-6}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -3.55e44 or 1.74999999999999997e-6 < y Initial program 11.4%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6411.4
Applied rewrites11.4%
Taylor expanded in y around inf
lower-/.f6448.0
Applied rewrites48.0%
if -3.55e44 < y < 1.74999999999999997e-6Initial program 93.7%
Taylor expanded in y around 0
lower-/.f6453.0
Applied rewrites53.0%
(FPCore (x y z t a b c i) :precision binary64 (/ z y))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return z / y;
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
code = z / y
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return z / y;
}
def code(x, y, z, t, a, b, c, i): return z / y
function code(x, y, z, t, a, b, c, i) return Float64(z / y) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = z / y; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(z / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{z}{y}
\end{array}
Initial program 52.2%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6447.5
Applied rewrites47.5%
Taylor expanded in y around inf
Applied rewrites14.9%
(FPCore (x y z t a b c i) :precision binary64 (/ 27464.7644705 b))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return 27464.7644705 / b;
}
real(8) function code(x, y, z, t, a, b, c, i)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
code = 27464.7644705d0 / b
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return 27464.7644705 / b;
}
def code(x, y, z, t, a, b, c, i): return 27464.7644705 / b
function code(x, y, z, t, a, b, c, i) return Float64(27464.7644705 / b) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = 27464.7644705 / b; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(27464.7644705 / b), $MachinePrecision]
\begin{array}{l}
\\
\frac{27464.7644705}{b}
\end{array}
Initial program 52.2%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites54.6%
Taylor expanded in b around inf
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f649.1
Applied rewrites9.1%
Taylor expanded in y around 0
Applied rewrites6.2%
Taylor expanded in y around inf
Applied rewrites3.3%
herbie shell --seed 2024221
(FPCore (x y z t a b c i)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2"
:precision binary64
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))