
(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 20 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 (fma (fma (+ a y) y b) y c) y i)))
(if (<= y -4.2e+74)
(- (+ (/ z y) (+ (/ 27464.7644705 (* y y)) x)) (* (/ b y) (/ x y)))
(if (<= y 4.7e+73)
(fma
y
(/ (fma (fma z y 27464.7644705) y 230661.510616) t_1)
(fma x (/ (pow y 4.0) t_1) (/ t t_1)))
(/ 1.0 (/ (+ (/ a y) 1.0) x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = fma(fma(fma((a + y), y, b), y, c), y, i);
double tmp;
if (y <= -4.2e+74) {
tmp = ((z / y) + ((27464.7644705 / (y * y)) + x)) - ((b / y) * (x / y));
} else if (y <= 4.7e+73) {
tmp = fma(y, (fma(fma(z, y, 27464.7644705), y, 230661.510616) / t_1), fma(x, (pow(y, 4.0) / t_1), (t / t_1)));
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(fma(fma(Float64(a + y), y, b), y, c), y, i) tmp = 0.0 if (y <= -4.2e+74) tmp = Float64(Float64(Float64(z / y) + Float64(Float64(27464.7644705 / Float64(y * y)) + x)) - Float64(Float64(b / y) * Float64(x / y))); elseif (y <= 4.7e+73) tmp = fma(y, Float64(fma(fma(z, y, 27464.7644705), y, 230661.510616) / t_1), fma(x, Float64((y ^ 4.0) / t_1), Float64(t / t_1))); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(N[(N[(a + y), $MachinePrecision] * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]}, If[LessEqual[y, -4.2e+74], N[(N[(N[(z / y), $MachinePrecision] + N[(N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision] - N[(N[(b / y), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.7e+73], N[(y * N[(N[(N[(z * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision] / t$95$1), $MachinePrecision] + N[(x * N[(N[Power[y, 4.0], $MachinePrecision] / t$95$1), $MachinePrecision] + N[(t / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(a + y, y, b\right), y, c\right), y, i\right)\\
\mathbf{if}\;y \leq -4.2 \cdot 10^{+74}:\\
\;\;\;\;\left(\frac{z}{y} + \left(\frac{27464.7644705}{y \cdot y} + x\right)\right) - \frac{b}{y} \cdot \frac{x}{y}\\
\mathbf{elif}\;y \leq 4.7 \cdot 10^{+73}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{\mathsf{fma}\left(\mathsf{fma}\left(z, y, 27464.7644705\right), y, 230661.510616\right)}{t\_1}, \mathsf{fma}\left(x, \frac{{y}^{4}}{t\_1}, \frac{t}{t\_1}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -4.1999999999999998e74Initial program 0.3%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.3%
Taylor expanded in y around inf
Applied rewrites79.7%
Taylor expanded in z around inf
Applied rewrites30.8%
Taylor expanded in y around inf
Applied rewrites81.4%
if -4.1999999999999998e74 < y < 4.7000000000000002e73Initial program 91.6%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites93.9%
if 4.7000000000000002e73 < y Initial program 2.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f642.5
Applied rewrites2.5%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6465.6
Applied rewrites65.6%
Taylor expanded in x around inf
Applied rewrites80.4%
Final simplification89.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (fma (fma (fma (+ a y) y b) y c) y i)))
(if (<= y -4.2e+74)
(- (+ (/ z y) (+ (/ 27464.7644705 (* y y)) x)) (* (/ b y) (/ x y)))
(if (<= y 4.7e+73)
(fma
y
(/ (fma (fma z y 27464.7644705) y 230661.510616) t_1)
(fma x (* (/ (* y y) t_1) (* y y)) (/ t t_1)))
(/ 1.0 (/ (+ (/ a y) 1.0) x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = fma(fma(fma((a + y), y, b), y, c), y, i);
double tmp;
if (y <= -4.2e+74) {
tmp = ((z / y) + ((27464.7644705 / (y * y)) + x)) - ((b / y) * (x / y));
} else if (y <= 4.7e+73) {
tmp = fma(y, (fma(fma(z, y, 27464.7644705), y, 230661.510616) / t_1), fma(x, (((y * y) / t_1) * (y * y)), (t / t_1)));
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(fma(fma(Float64(a + y), y, b), y, c), y, i) tmp = 0.0 if (y <= -4.2e+74) tmp = Float64(Float64(Float64(z / y) + Float64(Float64(27464.7644705 / Float64(y * y)) + x)) - Float64(Float64(b / y) * Float64(x / y))); elseif (y <= 4.7e+73) tmp = fma(y, Float64(fma(fma(z, y, 27464.7644705), y, 230661.510616) / t_1), fma(x, Float64(Float64(Float64(y * y) / t_1) * Float64(y * y)), Float64(t / t_1))); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(N[(N[(a + y), $MachinePrecision] * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]}, If[LessEqual[y, -4.2e+74], N[(N[(N[(z / y), $MachinePrecision] + N[(N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision] - N[(N[(b / y), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.7e+73], N[(y * N[(N[(N[(z * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision] / t$95$1), $MachinePrecision] + N[(x * N[(N[(N[(y * y), $MachinePrecision] / t$95$1), $MachinePrecision] * N[(y * y), $MachinePrecision]), $MachinePrecision] + N[(t / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(a + y, y, b\right), y, c\right), y, i\right)\\
\mathbf{if}\;y \leq -4.2 \cdot 10^{+74}:\\
\;\;\;\;\left(\frac{z}{y} + \left(\frac{27464.7644705}{y \cdot y} + x\right)\right) - \frac{b}{y} \cdot \frac{x}{y}\\
\mathbf{elif}\;y \leq 4.7 \cdot 10^{+73}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{\mathsf{fma}\left(\mathsf{fma}\left(z, y, 27464.7644705\right), y, 230661.510616\right)}{t\_1}, \mathsf{fma}\left(x, \frac{y \cdot y}{t\_1} \cdot \left(y \cdot y\right), \frac{t}{t\_1}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -4.1999999999999998e74Initial program 0.3%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.3%
Taylor expanded in y around inf
Applied rewrites79.7%
Taylor expanded in z around inf
Applied rewrites30.8%
Taylor expanded in y around inf
Applied rewrites81.4%
if -4.1999999999999998e74 < y < 4.7000000000000002e73Initial program 91.6%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
Applied rewrites93.9%
Applied rewrites93.9%
if 4.7000000000000002e73 < y Initial program 2.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f642.5
Applied rewrites2.5%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6465.6
Applied rewrites65.6%
Taylor expanded in x around inf
Applied rewrites80.4%
Final simplification89.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (fma (fma (fma (+ a y) y b) y c) y i)))
(if (<= y -5.6e+64)
(- (+ (/ z y) (+ (/ 27464.7644705 (* y y)) x)) (* (/ b y) (/ x y)))
(if (<= y 4.7e+73)
(fma
y
(/ (fma (fma (fma y x z) y 27464.7644705) y 230661.510616) t_1)
(/ t t_1))
(/ 1.0 (/ (+ (/ a y) 1.0) x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = fma(fma(fma((a + y), y, b), y, c), y, i);
double tmp;
if (y <= -5.6e+64) {
tmp = ((z / y) + ((27464.7644705 / (y * y)) + x)) - ((b / y) * (x / y));
} else if (y <= 4.7e+73) {
tmp = fma(y, (fma(fma(fma(y, x, z), y, 27464.7644705), y, 230661.510616) / t_1), (t / t_1));
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = fma(fma(fma(Float64(a + y), y, b), y, c), y, i) tmp = 0.0 if (y <= -5.6e+64) tmp = Float64(Float64(Float64(z / y) + Float64(Float64(27464.7644705 / Float64(y * y)) + x)) - Float64(Float64(b / y) * Float64(x / y))); elseif (y <= 4.7e+73) tmp = fma(y, Float64(fma(fma(fma(y, x, z), y, 27464.7644705), y, 230661.510616) / t_1), Float64(t / t_1)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(N[(N[(a + y), $MachinePrecision] * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]}, If[LessEqual[y, -5.6e+64], N[(N[(N[(z / y), $MachinePrecision] + N[(N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision] - N[(N[(b / y), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.7e+73], N[(y * N[(N[(N[(N[(y * x + z), $MachinePrecision] * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision] / t$95$1), $MachinePrecision] + N[(t / t$95$1), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(a + y, y, b\right), y, c\right), y, i\right)\\
\mathbf{if}\;y \leq -5.6 \cdot 10^{+64}:\\
\;\;\;\;\left(\frac{z}{y} + \left(\frac{27464.7644705}{y \cdot y} + x\right)\right) - \frac{b}{y} \cdot \frac{x}{y}\\
\mathbf{elif}\;y \leq 4.7 \cdot 10^{+73}:\\
\;\;\;\;\mathsf{fma}\left(y, \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, x, z\right), y, 27464.7644705\right), y, 230661.510616\right)}{t\_1}, \frac{t}{t\_1}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -5.60000000000000047e64Initial program 0.4%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.6%
Taylor expanded in y around inf
Applied rewrites75.6%
Taylor expanded in z around inf
Applied rewrites30.6%
Taylor expanded in y around inf
Applied rewrites77.2%
if -5.60000000000000047e64 < y < 4.7000000000000002e73Initial program 93.8%
Taylor expanded in t around 0
+-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites93.8%
if 4.7000000000000002e73 < y Initial program 2.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f642.5
Applied rewrites2.5%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6465.6
Applied rewrites65.6%
Taylor expanded in x around inf
Applied rewrites80.4%
Final simplification88.4%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -2.45e+64)
(- (+ (/ z y) (+ (/ 27464.7644705 (* y y)) x)) (* (/ b y) (/ x y)))
(if (<= y 1.32e+73)
(/
(+
(/ y (/ 1.0 (fma (fma (fma y x z) y 27464.7644705) y 230661.510616)))
t)
(+ (* (+ (* (+ (* (+ a y) y) b) y) c) y) i))
(/ 1.0 (/ (+ (/ a y) 1.0) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -2.45e+64) {
tmp = ((z / y) + ((27464.7644705 / (y * y)) + x)) - ((b / y) * (x / y));
} else if (y <= 1.32e+73) {
tmp = ((y / (1.0 / fma(fma(fma(y, x, z), y, 27464.7644705), y, 230661.510616))) + t) / (((((((a + y) * y) + b) * y) + c) * y) + i);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -2.45e+64) tmp = Float64(Float64(Float64(z / y) + Float64(Float64(27464.7644705 / Float64(y * y)) + x)) - Float64(Float64(b / y) * Float64(x / y))); elseif (y <= 1.32e+73) tmp = Float64(Float64(Float64(y / Float64(1.0 / fma(fma(fma(y, x, z), y, 27464.7644705), y, 230661.510616))) + t) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(a + y) * y) + b) * y) + c) * y) + i)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -2.45e+64], N[(N[(N[(z / y), $MachinePrecision] + N[(N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision] - N[(N[(b / y), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.32e+73], N[(N[(N[(y / N[(1.0 / N[(N[(N[(y * x + z), $MachinePrecision] * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(a + y), $MachinePrecision] * y), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] + c), $MachinePrecision] * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.45 \cdot 10^{+64}:\\
\;\;\;\;\left(\frac{z}{y} + \left(\frac{27464.7644705}{y \cdot y} + x\right)\right) - \frac{b}{y} \cdot \frac{x}{y}\\
\mathbf{elif}\;y \leq 1.32 \cdot 10^{+73}:\\
\;\;\;\;\frac{\frac{y}{\frac{1}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, x, z\right), y, 27464.7644705\right), y, 230661.510616\right)}} + t}{\left(\left(\left(a + y\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -2.4500000000000001e64Initial program 0.4%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.6%
Taylor expanded in y around inf
Applied rewrites75.6%
Taylor expanded in z around inf
Applied rewrites30.6%
Taylor expanded in y around inf
Applied rewrites77.2%
if -2.4500000000000001e64 < y < 1.32e73Initial program 93.8%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
clear-numN/A
flip-+N/A
Applied rewrites93.8%
if 1.32e73 < y Initial program 2.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f642.5
Applied rewrites2.5%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6465.6
Applied rewrites65.6%
Taylor expanded in x around inf
Applied rewrites80.4%
Final simplification88.3%
(FPCore (x y z t a b c i)
:precision binary64
(if (<=
(/
(+
(* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y)
t)
(+ (* (+ (* (+ (* (+ a y) y) b) y) c) y) i))
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 ((((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((a + y) * y) + b) * y) + c) * y) + i)) <= ((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 ((((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((a + y) * y) + b) * y) + c) * y) + i)) <= 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 (((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((a + y) * y) + b) * y) + c) * y) + i)) <= 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(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(a + y) * y) + b) * y) + c) * y) + i)) <= 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 ((((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((a + y) * y) + b) * y) + c) * y) + i)) <= 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[(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[(a + y), $MachinePrecision] * y), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] + c), $MachinePrecision] * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], Infinity], N[(t / i), $MachinePrecision], N[(z / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\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(a + y\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i} \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 91.7%
Taylor expanded in y around 0
lower-/.f6441.1
Applied rewrites41.1%
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 a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.1%
Taylor expanded in y around inf
Applied rewrites74.2%
Taylor expanded in z around inf
Applied rewrites24.4%
Final simplification35.5%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -2.45e+64)
(- (+ (/ z y) (+ (/ 27464.7644705 (* y y)) x)) (* (/ b y) (/ x y)))
(if (<= y 1.32e+73)
(/
(+ (fma (* (fma (fma y x z) y 27464.7644705) y) y (* 230661.510616 y)) t)
(+ (* (+ (* (+ (* (+ a y) y) b) y) c) y) i))
(/ 1.0 (/ (+ (/ a y) 1.0) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -2.45e+64) {
tmp = ((z / y) + ((27464.7644705 / (y * y)) + x)) - ((b / y) * (x / y));
} else if (y <= 1.32e+73) {
tmp = (fma((fma(fma(y, x, z), y, 27464.7644705) * y), y, (230661.510616 * y)) + t) / (((((((a + y) * y) + b) * y) + c) * y) + i);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -2.45e+64) tmp = Float64(Float64(Float64(z / y) + Float64(Float64(27464.7644705 / Float64(y * y)) + x)) - Float64(Float64(b / y) * Float64(x / y))); elseif (y <= 1.32e+73) tmp = Float64(Float64(fma(Float64(fma(fma(y, x, z), y, 27464.7644705) * y), y, Float64(230661.510616 * y)) + t) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(a + y) * y) + b) * y) + c) * y) + i)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -2.45e+64], N[(N[(N[(z / y), $MachinePrecision] + N[(N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision] - N[(N[(b / y), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.32e+73], N[(N[(N[(N[(N[(N[(y * x + z), $MachinePrecision] * y + 27464.7644705), $MachinePrecision] * y), $MachinePrecision] * y + N[(230661.510616 * y), $MachinePrecision]), $MachinePrecision] + t), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(a + y), $MachinePrecision] * y), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] + c), $MachinePrecision] * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.45 \cdot 10^{+64}:\\
\;\;\;\;\left(\frac{z}{y} + \left(\frac{27464.7644705}{y \cdot y} + x\right)\right) - \frac{b}{y} \cdot \frac{x}{y}\\
\mathbf{elif}\;y \leq 1.32 \cdot 10^{+73}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, x, z\right), y, 27464.7644705\right) \cdot y, y, 230661.510616 \cdot y\right) + t}{\left(\left(\left(a + y\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -2.4500000000000001e64Initial program 0.4%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.6%
Taylor expanded in y around inf
Applied rewrites75.6%
Taylor expanded in z around inf
Applied rewrites30.6%
Taylor expanded in y around inf
Applied rewrites77.2%
if -2.4500000000000001e64 < y < 1.32e73Initial program 93.8%
lift-*.f64N/A
*-commutativeN/A
lift-+.f64N/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lower-fma.f64N/A
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f6493.8
Applied rewrites93.8%
if 1.32e73 < y Initial program 2.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f642.5
Applied rewrites2.5%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6465.6
Applied rewrites65.6%
Taylor expanded in x around inf
Applied rewrites80.4%
Final simplification88.3%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -2.45e+64)
(- (+ (/ z y) (+ (/ 27464.7644705 (* y y)) x)) (* (/ b y) (/ x y)))
(if (<= y 1.32e+73)
(*
(/ -1.0 (fma (fma (fma (+ a y) y b) y c) y i))
(- (fma (fma (fma (fma y x z) y 27464.7644705) y 230661.510616) y t)))
(/ 1.0 (/ (+ (/ a y) 1.0) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -2.45e+64) {
tmp = ((z / y) + ((27464.7644705 / (y * y)) + x)) - ((b / y) * (x / y));
} else if (y <= 1.32e+73) {
tmp = (-1.0 / fma(fma(fma((a + y), y, b), y, c), y, i)) * -fma(fma(fma(fma(y, x, z), y, 27464.7644705), y, 230661.510616), y, t);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -2.45e+64) tmp = Float64(Float64(Float64(z / y) + Float64(Float64(27464.7644705 / Float64(y * y)) + x)) - Float64(Float64(b / y) * Float64(x / y))); elseif (y <= 1.32e+73) tmp = Float64(Float64(-1.0 / fma(fma(fma(Float64(a + y), y, b), y, c), y, i)) * Float64(-fma(fma(fma(fma(y, x, z), y, 27464.7644705), y, 230661.510616), y, t))); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -2.45e+64], N[(N[(N[(z / y), $MachinePrecision] + N[(N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision] - N[(N[(b / y), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.32e+73], N[(N[(-1.0 / N[(N[(N[(N[(a + y), $MachinePrecision] * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision] * (-N[(N[(N[(N[(y * x + z), $MachinePrecision] * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision] * y + t), $MachinePrecision])), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.45 \cdot 10^{+64}:\\
\;\;\;\;\left(\frac{z}{y} + \left(\frac{27464.7644705}{y \cdot y} + x\right)\right) - \frac{b}{y} \cdot \frac{x}{y}\\
\mathbf{elif}\;y \leq 1.32 \cdot 10^{+73}:\\
\;\;\;\;\frac{-1}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(a + y, y, b\right), y, c\right), y, i\right)} \cdot \left(-\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, x, z\right), y, 27464.7644705\right), y, 230661.510616\right), y, t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -2.4500000000000001e64Initial program 0.4%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.6%
Taylor expanded in y around inf
Applied rewrites75.6%
Taylor expanded in z around inf
Applied rewrites30.6%
Taylor expanded in y around inf
Applied rewrites77.2%
if -2.4500000000000001e64 < y < 1.32e73Initial program 93.8%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-*.f64N/A
Applied rewrites93.6%
if 1.32e73 < y Initial program 2.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f642.5
Applied rewrites2.5%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6465.6
Applied rewrites65.6%
Taylor expanded in x around inf
Applied rewrites80.4%
Final simplification88.2%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -9.5e+58)
(- (+ (/ z y) (+ (/ 27464.7644705 (* y y)) x)) (* (/ b y) (/ x y)))
(if (<= y 1.2e-52)
(/
(fma (fma (fma (fma y x z) y 27464.7644705) y 230661.510616) y t)
(fma (fma (fma y y b) y c) y i))
(if (<= y 1.32e+73)
(/
(fma (fma (fma (fma x y z) y 27464.7644705) y 230661.510616) y t)
(fma (fma (* y y) (+ a y) c) y i))
(/ 1.0 (/ (+ (/ a y) 1.0) x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -9.5e+58) {
tmp = ((z / y) + ((27464.7644705 / (y * y)) + x)) - ((b / y) * (x / y));
} else if (y <= 1.2e-52) {
tmp = fma(fma(fma(fma(y, x, z), y, 27464.7644705), y, 230661.510616), y, t) / fma(fma(fma(y, y, b), y, c), y, i);
} else if (y <= 1.32e+73) {
tmp = fma(fma(fma(fma(x, y, z), y, 27464.7644705), y, 230661.510616), y, t) / fma(fma((y * y), (a + y), c), y, i);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -9.5e+58) tmp = Float64(Float64(Float64(z / y) + Float64(Float64(27464.7644705 / Float64(y * y)) + x)) - Float64(Float64(b / y) * Float64(x / y))); elseif (y <= 1.2e-52) tmp = Float64(fma(fma(fma(fma(y, x, z), y, 27464.7644705), y, 230661.510616), y, t) / fma(fma(fma(y, y, b), y, c), y, i)); elseif (y <= 1.32e+73) tmp = Float64(fma(fma(fma(fma(x, y, z), y, 27464.7644705), y, 230661.510616), y, t) / fma(fma(Float64(y * y), Float64(a + y), c), y, i)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -9.5e+58], N[(N[(N[(z / y), $MachinePrecision] + N[(N[(27464.7644705 / N[(y * y), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision] - N[(N[(b / y), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.2e-52], N[(N[(N[(N[(N[(y * x + z), $MachinePrecision] * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision] * y + t), $MachinePrecision] / N[(N[(N[(y * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.32e+73], N[(N[(N[(N[(N[(x * y + z), $MachinePrecision] * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision] * y + t), $MachinePrecision] / N[(N[(N[(y * y), $MachinePrecision] * N[(a + y), $MachinePrecision] + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.5 \cdot 10^{+58}:\\
\;\;\;\;\left(\frac{z}{y} + \left(\frac{27464.7644705}{y \cdot y} + x\right)\right) - \frac{b}{y} \cdot \frac{x}{y}\\
\mathbf{elif}\;y \leq 1.2 \cdot 10^{-52}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, x, z\right), y, 27464.7644705\right), y, 230661.510616\right), y, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, y, b\right), y, c\right), y, i\right)}\\
\mathbf{elif}\;y \leq 1.32 \cdot 10^{+73}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(x, y, z\right), y, 27464.7644705\right), y, 230661.510616\right), y, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(y \cdot y, a + y, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -9.5000000000000002e58Initial program 2.5%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.8%
Taylor expanded in y around inf
Applied rewrites74.5%
Taylor expanded in z around inf
Applied rewrites29.6%
Taylor expanded in y around inf
Applied rewrites76.2%
if -9.5000000000000002e58 < y < 1.2000000000000001e-52Initial program 96.8%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites94.6%
if 1.2000000000000001e-52 < y < 1.32e73Initial program 82.8%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites60.3%
Taylor expanded in y around inf
Applied rewrites3.9%
Taylor expanded in z around inf
Applied rewrites4.0%
Taylor expanded in b around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites79.6%
if 1.32e73 < y Initial program 2.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f642.5
Applied rewrites2.5%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6465.6
Applied rewrites65.6%
Taylor expanded in x around inf
Applied rewrites80.4%
Final simplification86.8%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -1e+59)
(+ (/ z y) x)
(if (<= y 1.2e-52)
(/
(fma (fma (fma (fma y x z) y 27464.7644705) y 230661.510616) y t)
(fma (fma (fma y y b) y c) y i))
(if (<= y 1.32e+73)
(/
(fma (fma (fma (fma x y z) y 27464.7644705) y 230661.510616) y t)
(fma (fma (* y y) (+ a y) c) y i))
(/ 1.0 (/ (+ (/ a y) 1.0) x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -1e+59) {
tmp = (z / y) + x;
} else if (y <= 1.2e-52) {
tmp = fma(fma(fma(fma(y, x, z), y, 27464.7644705), y, 230661.510616), y, t) / fma(fma(fma(y, y, b), y, c), y, i);
} else if (y <= 1.32e+73) {
tmp = fma(fma(fma(fma(x, y, z), y, 27464.7644705), y, 230661.510616), y, t) / fma(fma((y * y), (a + y), c), y, i);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -1e+59) tmp = Float64(Float64(z / y) + x); elseif (y <= 1.2e-52) tmp = Float64(fma(fma(fma(fma(y, x, z), y, 27464.7644705), y, 230661.510616), y, t) / fma(fma(fma(y, y, b), y, c), y, i)); elseif (y <= 1.32e+73) tmp = Float64(fma(fma(fma(fma(x, y, z), y, 27464.7644705), y, 230661.510616), y, t) / fma(fma(Float64(y * y), Float64(a + y), c), y, i)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -1e+59], N[(N[(z / y), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[y, 1.2e-52], N[(N[(N[(N[(N[(y * x + z), $MachinePrecision] * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision] * y + t), $MachinePrecision] / N[(N[(N[(y * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.32e+73], N[(N[(N[(N[(N[(x * y + z), $MachinePrecision] * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision] * y + t), $MachinePrecision] / N[(N[(N[(y * y), $MachinePrecision] * N[(a + y), $MachinePrecision] + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1 \cdot 10^{+59}:\\
\;\;\;\;\frac{z}{y} + x\\
\mathbf{elif}\;y \leq 1.2 \cdot 10^{-52}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, x, z\right), y, 27464.7644705\right), y, 230661.510616\right), y, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, y, b\right), y, c\right), y, i\right)}\\
\mathbf{elif}\;y \leq 1.32 \cdot 10^{+73}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(x, y, z\right), y, 27464.7644705\right), y, 230661.510616\right), y, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(y \cdot y, a + y, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -9.99999999999999972e58Initial program 2.5%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.8%
Taylor expanded in y around inf
Applied rewrites74.5%
if -9.99999999999999972e58 < y < 1.2000000000000001e-52Initial program 96.8%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites94.6%
if 1.2000000000000001e-52 < y < 1.32e73Initial program 82.8%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites60.3%
Taylor expanded in y around inf
Applied rewrites3.9%
Taylor expanded in z around inf
Applied rewrites4.0%
Taylor expanded in b around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites79.6%
if 1.32e73 < y Initial program 2.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f642.5
Applied rewrites2.5%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6465.6
Applied rewrites65.6%
Taylor expanded in x around inf
Applied rewrites80.4%
Final simplification86.5%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -1.65e+62)
(+ (/ z y) x)
(if (<= y 1.1e+73)
(/
(fma (fma (fma z y 27464.7644705) y 230661.510616) y t)
(fma (fma (fma (+ a y) y b) y c) y i))
(/ 1.0 (/ (+ (/ a y) 1.0) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -1.65e+62) {
tmp = (z / y) + x;
} else if (y <= 1.1e+73) {
tmp = fma(fma(fma(z, y, 27464.7644705), y, 230661.510616), y, t) / fma(fma(fma((a + y), y, b), y, c), y, i);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -1.65e+62) tmp = Float64(Float64(z / y) + x); elseif (y <= 1.1e+73) tmp = Float64(fma(fma(fma(z, y, 27464.7644705), y, 230661.510616), y, t) / fma(fma(fma(Float64(a + y), y, b), y, c), y, i)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -1.65e+62], N[(N[(z / y), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[y, 1.1e+73], N[(N[(N[(N[(z * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision] * y + t), $MachinePrecision] / N[(N[(N[(N[(a + y), $MachinePrecision] * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.65 \cdot 10^{+62}:\\
\;\;\;\;\frac{z}{y} + x\\
\mathbf{elif}\;y \leq 1.1 \cdot 10^{+73}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z, y, 27464.7644705\right), y, 230661.510616\right), y, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(a + y, y, b\right), y, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -1.65e62Initial program 0.4%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.6%
Taylor expanded in y around inf
Applied rewrites75.6%
if -1.65e62 < y < 1.1e73Initial program 93.8%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites88.3%
if 1.1e73 < y Initial program 2.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f642.5
Applied rewrites2.5%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6465.6
Applied rewrites65.6%
Taylor expanded in x around inf
Applied rewrites80.4%
Final simplification84.5%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -2.25e+58)
(+ (/ z y) x)
(if (<= y 1.22e+48)
(/
(fma (fma (fma z y 27464.7644705) y 230661.510616) y t)
(fma (fma (fma y y b) y c) y i))
(/ 1.0 (/ (+ (/ a y) 1.0) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -2.25e+58) {
tmp = (z / y) + x;
} else if (y <= 1.22e+48) {
tmp = fma(fma(fma(z, y, 27464.7644705), y, 230661.510616), y, t) / fma(fma(fma(y, y, b), y, c), y, i);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -2.25e+58) tmp = Float64(Float64(z / y) + x); elseif (y <= 1.22e+48) tmp = Float64(fma(fma(fma(z, y, 27464.7644705), y, 230661.510616), y, t) / fma(fma(fma(y, y, b), y, c), y, i)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -2.25e+58], N[(N[(z / y), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[y, 1.22e+48], N[(N[(N[(N[(z * y + 27464.7644705), $MachinePrecision] * y + 230661.510616), $MachinePrecision] * y + t), $MachinePrecision] / N[(N[(N[(y * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.25 \cdot 10^{+58}:\\
\;\;\;\;\frac{z}{y} + x\\
\mathbf{elif}\;y \leq 1.22 \cdot 10^{+48}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z, y, 27464.7644705\right), y, 230661.510616\right), y, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, y, b\right), y, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -2.2499999999999999e58Initial program 2.5%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.8%
Taylor expanded in y around inf
Applied rewrites74.5%
if -2.2499999999999999e58 < y < 1.22000000000000004e48Initial program 95.9%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites90.6%
Taylor expanded in y around 0
Applied rewrites74.8%
Taylor expanded in x around 0
Applied rewrites86.2%
if 1.22000000000000004e48 < y Initial program 4.8%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f644.8
Applied rewrites4.8%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6460.2
Applied rewrites60.2%
Taylor expanded in x around inf
Applied rewrites73.8%
Final simplification81.6%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -5.8e+38)
(+ (/ z y) x)
(if (<= y 1.65e+60)
(/
(fma (fma (* z y) y 230661.510616) y t)
(fma (fma (fma y y b) y c) y i))
(/ 1.0 (/ (+ (/ a y) 1.0) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -5.8e+38) {
tmp = (z / y) + x;
} else if (y <= 1.65e+60) {
tmp = fma(fma((z * y), y, 230661.510616), y, t) / fma(fma(fma(y, y, b), y, c), y, i);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -5.8e+38) tmp = Float64(Float64(z / y) + x); elseif (y <= 1.65e+60) tmp = Float64(fma(fma(Float64(z * y), y, 230661.510616), y, t) / fma(fma(fma(y, y, b), y, c), y, i)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -5.8e+38], N[(N[(z / y), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[y, 1.65e+60], N[(N[(N[(N[(z * y), $MachinePrecision] * y + 230661.510616), $MachinePrecision] * y + t), $MachinePrecision] / N[(N[(N[(y * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.8 \cdot 10^{+38}:\\
\;\;\;\;\frac{z}{y} + x\\
\mathbf{elif}\;y \leq 1.65 \cdot 10^{+60}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(z \cdot y, y, 230661.510616\right), y, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, y, b\right), y, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -5.80000000000000013e38Initial program 4.1%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites2.7%
Taylor expanded in y around inf
Applied rewrites68.2%
if -5.80000000000000013e38 < y < 1.6499999999999999e60Initial program 97.2%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites91.1%
Taylor expanded in z around inf
Applied rewrites85.3%
if 1.6499999999999999e60 < y Initial program 2.6%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f642.6
Applied rewrites2.6%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6464.3
Applied rewrites64.3%
Taylor expanded in x around inf
Applied rewrites78.7%
Final simplification80.4%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -7.2e+17)
(+ (/ z y) x)
(if (<= y 0.23)
(/ (fma 230661.510616 y t) (+ (* (+ (* (+ (* (+ a y) y) b) y) c) y) i))
(/ 1.0 (/ (+ (/ a y) 1.0) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -7.2e+17) {
tmp = (z / y) + x;
} else if (y <= 0.23) {
tmp = fma(230661.510616, y, t) / (((((((a + y) * y) + b) * y) + c) * y) + i);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -7.2e+17) tmp = Float64(Float64(z / y) + x); elseif (y <= 0.23) tmp = Float64(fma(230661.510616, y, t) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(a + y) * y) + b) * y) + c) * y) + i)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -7.2e+17], N[(N[(z / y), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[y, 0.23], N[(N[(230661.510616 * y + t), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(a + y), $MachinePrecision] * y), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] + c), $MachinePrecision] * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -7.2 \cdot 10^{+17}:\\
\;\;\;\;\frac{z}{y} + x\\
\mathbf{elif}\;y \leq 0.23:\\
\;\;\;\;\frac{\mathsf{fma}\left(230661.510616, y, t\right)}{\left(\left(\left(a + y\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -7.2e17Initial program 7.3%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites5.9%
Taylor expanded in y around inf
Applied rewrites66.1%
if -7.2e17 < y < 0.23000000000000001Initial program 99.6%
Taylor expanded in y around 0
+-commutativeN/A
lower-fma.f6487.1
Applied rewrites87.1%
if 0.23000000000000001 < y Initial program 18.2%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6418.3
Applied rewrites18.3%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6450.8
Applied rewrites50.8%
Taylor expanded in x around inf
Applied rewrites62.1%
Final simplification76.8%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -2.25e+58)
(+ (/ z y) x)
(if (<= y 1.8e+55)
(/
(fma (fma 27464.7644705 y 230661.510616) y t)
(fma (fma (fma y y b) y c) y i))
(/ 1.0 (/ (+ (/ a y) 1.0) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -2.25e+58) {
tmp = (z / y) + x;
} else if (y <= 1.8e+55) {
tmp = fma(fma(27464.7644705, y, 230661.510616), y, t) / fma(fma(fma(y, y, b), y, c), y, i);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -2.25e+58) tmp = Float64(Float64(z / y) + x); elseif (y <= 1.8e+55) tmp = Float64(fma(fma(27464.7644705, y, 230661.510616), y, t) / fma(fma(fma(y, y, b), y, c), y, i)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -2.25e+58], N[(N[(z / y), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[y, 1.8e+55], N[(N[(N[(27464.7644705 * y + 230661.510616), $MachinePrecision] * y + t), $MachinePrecision] / N[(N[(N[(y * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.25 \cdot 10^{+58}:\\
\;\;\;\;\frac{z}{y} + x\\
\mathbf{elif}\;y \leq 1.8 \cdot 10^{+55}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(27464.7644705, y, 230661.510616\right), y, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, y, b\right), y, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -2.2499999999999999e58Initial program 2.5%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites0.8%
Taylor expanded in y around inf
Applied rewrites74.5%
if -2.2499999999999999e58 < y < 1.79999999999999994e55Initial program 95.4%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites90.1%
Taylor expanded in y around 0
Applied rewrites76.2%
if 1.79999999999999994e55 < y Initial program 4.8%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f644.8
Applied rewrites4.8%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6461.6
Applied rewrites61.6%
Taylor expanded in x around inf
Applied rewrites75.4%
Final simplification75.7%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -1.75e+17)
(+ (/ z y) x)
(if (<= y 0.105)
(/ (fma 230661.510616 y t) (fma (fma (fma y y b) y c) y i))
(/ 1.0 (/ (+ (/ a y) 1.0) x)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -1.75e+17) {
tmp = (z / y) + x;
} else if (y <= 0.105) {
tmp = fma(230661.510616, y, t) / fma(fma(fma(y, y, b), y, c), y, i);
} else {
tmp = 1.0 / (((a / y) + 1.0) / x);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -1.75e+17) tmp = Float64(Float64(z / y) + x); elseif (y <= 0.105) tmp = Float64(fma(230661.510616, y, t) / fma(fma(fma(y, y, b), y, c), y, i)); else tmp = Float64(1.0 / Float64(Float64(Float64(a / y) + 1.0) / x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -1.75e+17], N[(N[(z / y), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[y, 0.105], N[(N[(230661.510616 * y + t), $MachinePrecision] / N[(N[(N[(y * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[(a / y), $MachinePrecision] + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.75 \cdot 10^{+17}:\\
\;\;\;\;\frac{z}{y} + x\\
\mathbf{elif}\;y \leq 0.105:\\
\;\;\;\;\frac{\mathsf{fma}\left(230661.510616, y, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, y, b\right), y, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a}{y} + 1}{x}}\\
\end{array}
\end{array}
if y < -1.75e17Initial program 7.3%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites5.9%
Taylor expanded in y around inf
Applied rewrites66.1%
if -1.75e17 < y < 0.104999999999999996Initial program 99.6%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites94.3%
Taylor expanded in y around 0
Applied rewrites83.3%
if 0.104999999999999996 < y Initial program 18.2%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6418.3
Applied rewrites18.3%
Taylor expanded in y around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f6450.8
Applied rewrites50.8%
Taylor expanded in x around inf
Applied rewrites62.1%
Final simplification74.7%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (/ z y) x)))
(if (<= y -1.75e+17)
t_1
(if (<= y 1.8e+29)
(/ (fma 230661.510616 y t) (fma (fma (fma y y b) y c) y 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) + x;
double tmp;
if (y <= -1.75e+17) {
tmp = t_1;
} else if (y <= 1.8e+29) {
tmp = fma(230661.510616, y, t) / fma(fma(fma(y, y, b), y, c), y, i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(z / y) + x) tmp = 0.0 if (y <= -1.75e+17) tmp = t_1; elseif (y <= 1.8e+29) tmp = Float64(fma(230661.510616, y, t) / fma(fma(fma(y, y, b), y, c), y, i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(z / y), $MachinePrecision] + x), $MachinePrecision]}, If[LessEqual[y, -1.75e+17], t$95$1, If[LessEqual[y, 1.8e+29], N[(N[(230661.510616 * y + t), $MachinePrecision] / N[(N[(N[(y * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z}{y} + x\\
\mathbf{if}\;y \leq -1.75 \cdot 10^{+17}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.8 \cdot 10^{+29}:\\
\;\;\;\;\frac{\mathsf{fma}\left(230661.510616, y, t\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, y, b\right), y, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.75e17 or 1.79999999999999988e29 < y Initial program 7.9%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites6.3%
Taylor expanded in y around inf
Applied rewrites65.2%
if -1.75e17 < y < 1.79999999999999988e29Initial program 99.6%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites93.8%
Taylor expanded in y around 0
Applied rewrites80.2%
Final simplification73.8%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (/ z y) x)))
(if (<= y -2.05e+61)
t_1
(if (<= y 1.7e+30) (/ t (fma (fma (fma (+ a y) y b) y c) y 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) + x;
double tmp;
if (y <= -2.05e+61) {
tmp = t_1;
} else if (y <= 1.7e+30) {
tmp = t / fma(fma(fma((a + y), y, b), y, c), y, i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(z / y) + x) tmp = 0.0 if (y <= -2.05e+61) tmp = t_1; elseif (y <= 1.7e+30) tmp = Float64(t / fma(fma(fma(Float64(a + y), y, b), y, c), y, i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(z / y), $MachinePrecision] + x), $MachinePrecision]}, If[LessEqual[y, -2.05e+61], t$95$1, If[LessEqual[y, 1.7e+30], N[(t / N[(N[(N[(N[(a + y), $MachinePrecision] * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z}{y} + x\\
\mathbf{if}\;y \leq -2.05 \cdot 10^{+61}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.7 \cdot 10^{+30}:\\
\;\;\;\;\frac{t}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(a + y, y, b\right), y, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.04999999999999986e61 or 1.7000000000000001e30 < y Initial program 4.5%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites3.6%
Taylor expanded in y around inf
Applied rewrites69.9%
if -2.04999999999999986e61 < y < 1.7000000000000001e30Initial program 96.5%
Taylor expanded in t around inf
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-+.f6465.1
Applied rewrites65.1%
Final simplification67.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (/ z y) x)))
(if (<= y -2.15e+17)
t_1
(if (<= y 2.6e+28) (/ t (fma (fma (fma y y b) y c) y 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) + x;
double tmp;
if (y <= -2.15e+17) {
tmp = t_1;
} else if (y <= 2.6e+28) {
tmp = t / fma(fma(fma(y, y, b), y, c), y, i);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(z / y) + x) tmp = 0.0 if (y <= -2.15e+17) tmp = t_1; elseif (y <= 2.6e+28) tmp = Float64(t / fma(fma(fma(y, y, b), y, c), y, i)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(z / y), $MachinePrecision] + x), $MachinePrecision]}, If[LessEqual[y, -2.15e+17], t$95$1, If[LessEqual[y, 2.6e+28], N[(t / N[(N[(N[(y * y + b), $MachinePrecision] * y + c), $MachinePrecision] * y + i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z}{y} + x\\
\mathbf{if}\;y \leq -2.15 \cdot 10^{+17}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2.6 \cdot 10^{+28}:\\
\;\;\;\;\frac{t}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(y, y, b\right), y, c\right), y, i\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.15e17 or 2.6000000000000002e28 < y Initial program 7.9%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites6.3%
Taylor expanded in y around inf
Applied rewrites65.2%
if -2.15e17 < y < 2.6000000000000002e28Initial program 99.6%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites93.8%
Taylor expanded in y around inf
Applied rewrites2.6%
Taylor expanded in z around inf
Applied rewrites2.8%
Taylor expanded in t around inf
Applied rewrites66.4%
Final simplification65.9%
(FPCore (x y z t a b c i) :precision binary64 (let* ((t_1 (+ (/ z y) x))) (if (<= y -9.4e+56) t_1 (if (<= y 6.5e+27) (/ 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 = (z / y) + x;
double tmp;
if (y <= -9.4e+56) {
tmp = t_1;
} else if (y <= 6.5e+27) {
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 = (z / y) + x
if (y <= (-9.4d+56)) then
tmp = t_1
else if (y <= 6.5d+27) 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 = (z / y) + x;
double tmp;
if (y <= -9.4e+56) {
tmp = t_1;
} else if (y <= 6.5e+27) {
tmp = t / i;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (z / y) + x tmp = 0 if y <= -9.4e+56: tmp = t_1 elif y <= 6.5e+27: tmp = t / i else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(z / y) + x) tmp = 0.0 if (y <= -9.4e+56) tmp = t_1; elseif (y <= 6.5e+27) 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 = (z / y) + x; tmp = 0.0; if (y <= -9.4e+56) tmp = t_1; elseif (y <= 6.5e+27) 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[(N[(z / y), $MachinePrecision] + x), $MachinePrecision]}, If[LessEqual[y, -9.4e+56], t$95$1, If[LessEqual[y, 6.5e+27], N[(t / i), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z}{y} + x\\
\mathbf{if}\;y \leq -9.4 \cdot 10^{+56}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 6.5 \cdot 10^{+27}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -9.400000000000001e56 or 6.5000000000000005e27 < y Initial program 5.4%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites3.6%
Taylor expanded in y around inf
Applied rewrites69.5%
if -9.400000000000001e56 < y < 6.5000000000000005e27Initial program 97.1%
Taylor expanded in y around 0
lower-/.f6444.9
Applied rewrites44.9%
Final simplification54.7%
(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 60.6%
Taylor expanded in a around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
Applied rewrites56.6%
Taylor expanded in y around inf
Applied rewrites29.3%
Taylor expanded in z around inf
Applied rewrites10.7%
herbie shell --seed 2024276
(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)))