
(FPCore (x y z t a b) :precision binary64 (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))
double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
real(8) function code(x, y, z, t, a, b)
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
code = ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a)) + (((y + t) - 2.0d0) * b)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
def code(x, y, z, t, a, b): return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b)) end
function tmp = code(x, y, z, t, a, b) tmp = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 27 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b) :precision binary64 (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))
double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
real(8) function code(x, y, z, t, a, b)
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
code = ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a)) + (((y + t) - 2.0d0) * b)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
def code(x, y, z, t, a, b): return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b)) end
function tmp = code(x, y, z, t, a, b) tmp = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b
\end{array}
(FPCore (x y z t a b) :precision binary64 (fma (- 1.0 t) a (+ (fma (- b z) y (fma (- t 2.0) b x)) z)))
double code(double x, double y, double z, double t, double a, double b) {
return fma((1.0 - t), a, (fma((b - z), y, fma((t - 2.0), b, x)) + z));
}
function code(x, y, z, t, a, b) return fma(Float64(1.0 - t), a, Float64(fma(Float64(b - z), y, fma(Float64(t - 2.0), b, x)) + z)) end
code[x_, y_, z_, t_, a_, b_] := N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(N[(b - z), $MachinePrecision] * y + N[(N[(t - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(b - z, y, \mathsf{fma}\left(t - 2, b, x\right)\right) + z\right)
\end{array}
Initial program 96.8%
Taylor expanded in x around 0
Applied rewrites98.0%
Final simplification98.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+ (* (- (+ t y) 2.0) b) (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)))))
(if (<= t_1 -2e+303) (* b t) (if (<= t_1 5e+300) (+ (+ z x) a) (* b y)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (((t + y) - 2.0) * b) + ((x - ((y - 1.0) * z)) - ((t - 1.0) * a));
double tmp;
if (t_1 <= -2e+303) {
tmp = b * t;
} else if (t_1 <= 5e+300) {
tmp = (z + x) + a;
} else {
tmp = b * y;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = (((t + y) - 2.0d0) * b) + ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a))
if (t_1 <= (-2d+303)) then
tmp = b * t
else if (t_1 <= 5d+300) then
tmp = (z + x) + a
else
tmp = b * y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (((t + y) - 2.0) * b) + ((x - ((y - 1.0) * z)) - ((t - 1.0) * a));
double tmp;
if (t_1 <= -2e+303) {
tmp = b * t;
} else if (t_1 <= 5e+300) {
tmp = (z + x) + a;
} else {
tmp = b * y;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (((t + y) - 2.0) * b) + ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) tmp = 0 if t_1 <= -2e+303: tmp = b * t elif t_1 <= 5e+300: tmp = (z + x) + a else: tmp = b * y return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(Float64(t + y) - 2.0) * b) + Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a))) tmp = 0.0 if (t_1 <= -2e+303) tmp = Float64(b * t); elseif (t_1 <= 5e+300) tmp = Float64(Float64(z + x) + a); else tmp = Float64(b * y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (((t + y) - 2.0) * b) + ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)); tmp = 0.0; if (t_1 <= -2e+303) tmp = b * t; elseif (t_1 <= 5e+300) tmp = (z + x) + a; else tmp = b * y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision] + N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e+303], N[(b * t), $MachinePrecision], If[LessEqual[t$95$1, 5e+300], N[(N[(z + x), $MachinePrecision] + a), $MachinePrecision], N[(b * y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(t + y\right) - 2\right) \cdot b + \left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right)\\
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{+303}:\\
\;\;\;\;b \cdot t\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{+300}:\\
\;\;\;\;\left(z + x\right) + a\\
\mathbf{else}:\\
\;\;\;\;b \cdot y\\
\end{array}
\end{array}
if (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y #s(literal 1 binary64)) z)) (*.f64 (-.f64 t #s(literal 1 binary64)) a)) (*.f64 (-.f64 (+.f64 y t) #s(literal 2 binary64)) b)) < -2e303Initial program 100.0%
Taylor expanded in x around 0
Applied rewrites98.0%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6484.4
Applied rewrites84.4%
Taylor expanded in t around inf
Applied rewrites35.6%
if -2e303 < (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y #s(literal 1 binary64)) z)) (*.f64 (-.f64 t #s(literal 1 binary64)) a)) (*.f64 (-.f64 (+.f64 y t) #s(literal 2 binary64)) b)) < 5.00000000000000026e300Initial program 99.9%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites85.2%
Taylor expanded in b around 0
Applied rewrites62.6%
Taylor expanded in y around 0
Applied rewrites57.1%
if 5.00000000000000026e300 < (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y #s(literal 1 binary64)) z)) (*.f64 (-.f64 t #s(literal 1 binary64)) a)) (*.f64 (-.f64 (+.f64 y t) #s(literal 2 binary64)) b)) Initial program 83.3%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6467.8
Applied rewrites67.8%
Taylor expanded in y around inf
Applied rewrites38.6%
Final simplification49.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+ (* (- (+ t y) 2.0) b) (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)))))
(if (<= t_1 -2e+303) (* b t) (if (<= t_1 5e+300) (+ a x) (* b y)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (((t + y) - 2.0) * b) + ((x - ((y - 1.0) * z)) - ((t - 1.0) * a));
double tmp;
if (t_1 <= -2e+303) {
tmp = b * t;
} else if (t_1 <= 5e+300) {
tmp = a + x;
} else {
tmp = b * y;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = (((t + y) - 2.0d0) * b) + ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a))
if (t_1 <= (-2d+303)) then
tmp = b * t
else if (t_1 <= 5d+300) then
tmp = a + x
else
tmp = b * y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (((t + y) - 2.0) * b) + ((x - ((y - 1.0) * z)) - ((t - 1.0) * a));
double tmp;
if (t_1 <= -2e+303) {
tmp = b * t;
} else if (t_1 <= 5e+300) {
tmp = a + x;
} else {
tmp = b * y;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (((t + y) - 2.0) * b) + ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) tmp = 0 if t_1 <= -2e+303: tmp = b * t elif t_1 <= 5e+300: tmp = a + x else: tmp = b * y return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(Float64(t + y) - 2.0) * b) + Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a))) tmp = 0.0 if (t_1 <= -2e+303) tmp = Float64(b * t); elseif (t_1 <= 5e+300) tmp = Float64(a + x); else tmp = Float64(b * y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (((t + y) - 2.0) * b) + ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)); tmp = 0.0; if (t_1 <= -2e+303) tmp = b * t; elseif (t_1 <= 5e+300) tmp = a + x; else tmp = b * y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision] + N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e+303], N[(b * t), $MachinePrecision], If[LessEqual[t$95$1, 5e+300], N[(a + x), $MachinePrecision], N[(b * y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(t + y\right) - 2\right) \cdot b + \left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right)\\
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{+303}:\\
\;\;\;\;b \cdot t\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{+300}:\\
\;\;\;\;a + x\\
\mathbf{else}:\\
\;\;\;\;b \cdot y\\
\end{array}
\end{array}
if (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y #s(literal 1 binary64)) z)) (*.f64 (-.f64 t #s(literal 1 binary64)) a)) (*.f64 (-.f64 (+.f64 y t) #s(literal 2 binary64)) b)) < -2e303Initial program 100.0%
Taylor expanded in x around 0
Applied rewrites98.0%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6484.4
Applied rewrites84.4%
Taylor expanded in t around inf
Applied rewrites35.6%
if -2e303 < (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y #s(literal 1 binary64)) z)) (*.f64 (-.f64 t #s(literal 1 binary64)) a)) (*.f64 (-.f64 (+.f64 y t) #s(literal 2 binary64)) b)) < 5.00000000000000026e300Initial program 99.9%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites85.2%
Taylor expanded in b around 0
Applied rewrites62.6%
Taylor expanded in z around 0
Applied rewrites37.6%
if 5.00000000000000026e300 < (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y #s(literal 1 binary64)) z)) (*.f64 (-.f64 t #s(literal 1 binary64)) a)) (*.f64 (-.f64 (+.f64 y t) #s(literal 2 binary64)) b)) Initial program 83.3%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6467.8
Applied rewrites67.8%
Taylor expanded in y around inf
Applied rewrites38.6%
Final simplification37.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma (- b a) t (fma (- y 2.0) b (fma (- 1.0 y) z a)))))
(if (<= z -5e-69)
t_1
(if (<= z 155000.0) (+ (fma (- b a) t a) (fma (- y 2.0) b x)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma((b - a), t, fma((y - 2.0), b, fma((1.0 - y), z, a)));
double tmp;
if (z <= -5e-69) {
tmp = t_1;
} else if (z <= 155000.0) {
tmp = fma((b - a), t, a) + fma((y - 2.0), b, x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(Float64(b - a), t, fma(Float64(y - 2.0), b, fma(Float64(1.0 - y), z, a))) tmp = 0.0 if (z <= -5e-69) tmp = t_1; elseif (z <= 155000.0) tmp = Float64(fma(Float64(b - a), t, a) + fma(Float64(y - 2.0), b, x)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - a), $MachinePrecision] * t + N[(N[(y - 2.0), $MachinePrecision] * b + N[(N[(1.0 - y), $MachinePrecision] * z + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -5e-69], t$95$1, If[LessEqual[z, 155000.0], N[(N[(N[(b - a), $MachinePrecision] * t + a), $MachinePrecision] + N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(b - a, t, \mathsf{fma}\left(y - 2, b, \mathsf{fma}\left(1 - y, z, a\right)\right)\right)\\
\mathbf{if}\;z \leq -5 \cdot 10^{-69}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 155000:\\
\;\;\;\;\mathsf{fma}\left(b - a, t, a\right) + \mathsf{fma}\left(y - 2, b, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -5.00000000000000033e-69 or 155000 < z Initial program 95.6%
Taylor expanded in x around 0
Applied rewrites94.1%
if -5.00000000000000033e-69 < z < 155000Initial program 98.3%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6497.5
Applied rewrites97.5%
Taylor expanded in t around 0
Applied rewrites98.3%
Final simplification96.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -1.46e+144)
(* (- (+ t y) 2.0) b)
(if (<= b -2.55e+41)
(fma (- 1.0 y) z (fma (- y 2.0) b x))
(if (<= b 1.02e+93)
(fma (- 1.0 t) a (+ (fma (- z) y x) z))
(fma (- (+ t y)) (- b) (* (- b) 2.0))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.46e+144) {
tmp = ((t + y) - 2.0) * b;
} else if (b <= -2.55e+41) {
tmp = fma((1.0 - y), z, fma((y - 2.0), b, x));
} else if (b <= 1.02e+93) {
tmp = fma((1.0 - t), a, (fma(-z, y, x) + z));
} else {
tmp = fma(-(t + y), -b, (-b * 2.0));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -1.46e+144) tmp = Float64(Float64(Float64(t + y) - 2.0) * b); elseif (b <= -2.55e+41) tmp = fma(Float64(1.0 - y), z, fma(Float64(y - 2.0), b, x)); elseif (b <= 1.02e+93) tmp = fma(Float64(1.0 - t), a, Float64(fma(Float64(-z), y, x) + z)); else tmp = fma(Float64(-Float64(t + y)), Float64(-b), Float64(Float64(-b) * 2.0)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -1.46e+144], N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision], If[LessEqual[b, -2.55e+41], N[(N[(1.0 - y), $MachinePrecision] * z + N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.02e+93], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[((-z) * y + x), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision], N[((-N[(t + y), $MachinePrecision]) * (-b) + N[((-b) * 2.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.46 \cdot 10^{+144}:\\
\;\;\;\;\left(\left(t + y\right) - 2\right) \cdot b\\
\mathbf{elif}\;b \leq -2.55 \cdot 10^{+41}:\\
\;\;\;\;\mathsf{fma}\left(1 - y, z, \mathsf{fma}\left(y - 2, b, x\right)\right)\\
\mathbf{elif}\;b \leq 1.02 \cdot 10^{+93}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(-z, y, x\right) + z\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-\left(t + y\right), -b, \left(-b\right) \cdot 2\right)\\
\end{array}
\end{array}
if b < -1.46e144Initial program 91.9%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites71.9%
Taylor expanded in x around 0
Applied rewrites71.9%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
lower-+.f6493.9
Applied rewrites93.9%
if -1.46e144 < b < -2.54999999999999989e41Initial program 99.8%
Taylor expanded in x around 0
Applied rewrites99.9%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6499.8
Applied rewrites99.8%
Taylor expanded in t around 0
Applied rewrites89.7%
if -2.54999999999999989e41 < b < 1.0200000000000001e93Initial program 99.3%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in b around 0
Applied rewrites90.4%
if 1.0200000000000001e93 < b Initial program 91.5%
Taylor expanded in x around 0
Applied rewrites89.3%
Taylor expanded in b around -inf
Applied rewrites82.7%
Applied rewrites82.7%
Final simplification89.4%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -1.46e+144)
(* (- (+ t y) 2.0) b)
(if (<= b -2.55e+41)
(fma (- 1.0 y) z (fma (- y 2.0) b x))
(if (<= b 1.02e+93)
(fma (- 1.0 t) a (fma (- 1.0 y) z x))
(fma (- (+ t y)) (- b) (* (- b) 2.0))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.46e+144) {
tmp = ((t + y) - 2.0) * b;
} else if (b <= -2.55e+41) {
tmp = fma((1.0 - y), z, fma((y - 2.0), b, x));
} else if (b <= 1.02e+93) {
tmp = fma((1.0 - t), a, fma((1.0 - y), z, x));
} else {
tmp = fma(-(t + y), -b, (-b * 2.0));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -1.46e+144) tmp = Float64(Float64(Float64(t + y) - 2.0) * b); elseif (b <= -2.55e+41) tmp = fma(Float64(1.0 - y), z, fma(Float64(y - 2.0), b, x)); elseif (b <= 1.02e+93) tmp = fma(Float64(1.0 - t), a, fma(Float64(1.0 - y), z, x)); else tmp = fma(Float64(-Float64(t + y)), Float64(-b), Float64(Float64(-b) * 2.0)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -1.46e+144], N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision], If[LessEqual[b, -2.55e+41], N[(N[(1.0 - y), $MachinePrecision] * z + N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.02e+93], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(1.0 - y), $MachinePrecision] * z + x), $MachinePrecision]), $MachinePrecision], N[((-N[(t + y), $MachinePrecision]) * (-b) + N[((-b) * 2.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.46 \cdot 10^{+144}:\\
\;\;\;\;\left(\left(t + y\right) - 2\right) \cdot b\\
\mathbf{elif}\;b \leq -2.55 \cdot 10^{+41}:\\
\;\;\;\;\mathsf{fma}\left(1 - y, z, \mathsf{fma}\left(y - 2, b, x\right)\right)\\
\mathbf{elif}\;b \leq 1.02 \cdot 10^{+93}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(1 - y, z, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-\left(t + y\right), -b, \left(-b\right) \cdot 2\right)\\
\end{array}
\end{array}
if b < -1.46e144Initial program 91.9%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites71.9%
Taylor expanded in x around 0
Applied rewrites71.9%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
lower-+.f6493.9
Applied rewrites93.9%
if -1.46e144 < b < -2.54999999999999989e41Initial program 99.8%
Taylor expanded in x around 0
Applied rewrites99.9%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6499.8
Applied rewrites99.8%
Taylor expanded in t around 0
Applied rewrites89.7%
if -2.54999999999999989e41 < b < 1.0200000000000001e93Initial program 99.3%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in b around 0
Applied rewrites90.4%
if 1.0200000000000001e93 < b Initial program 91.5%
Taylor expanded in x around 0
Applied rewrites89.3%
Taylor expanded in b around -inf
Applied rewrites82.7%
Applied rewrites82.7%
Final simplification89.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- (+ t y) 2.0) b)))
(if (<= b -1.46e+144)
t_1
(if (<= b -2.55e+41)
(fma (- 1.0 y) z (fma (- y 2.0) b x))
(if (<= b 1.02e+93) (fma (- 1.0 t) a (fma (- 1.0 y) z x)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((t + y) - 2.0) * b;
double tmp;
if (b <= -1.46e+144) {
tmp = t_1;
} else if (b <= -2.55e+41) {
tmp = fma((1.0 - y), z, fma((y - 2.0), b, x));
} else if (b <= 1.02e+93) {
tmp = fma((1.0 - t), a, fma((1.0 - y), z, x));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(t + y) - 2.0) * b) tmp = 0.0 if (b <= -1.46e+144) tmp = t_1; elseif (b <= -2.55e+41) tmp = fma(Float64(1.0 - y), z, fma(Float64(y - 2.0), b, x)); elseif (b <= 1.02e+93) tmp = fma(Float64(1.0 - t), a, fma(Float64(1.0 - y), z, x)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[b, -1.46e+144], t$95$1, If[LessEqual[b, -2.55e+41], N[(N[(1.0 - y), $MachinePrecision] * z + N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.02e+93], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(1.0 - y), $MachinePrecision] * z + x), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(t + y\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -1.46 \cdot 10^{+144}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -2.55 \cdot 10^{+41}:\\
\;\;\;\;\mathsf{fma}\left(1 - y, z, \mathsf{fma}\left(y - 2, b, x\right)\right)\\
\mathbf{elif}\;b \leq 1.02 \cdot 10^{+93}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(1 - y, z, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.46e144 or 1.0200000000000001e93 < b Initial program 91.6%
Taylor expanded in x around 0
Applied rewrites94.0%
Taylor expanded in t around 0
Applied rewrites68.0%
Taylor expanded in x around 0
Applied rewrites66.3%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
lower-+.f6487.6
Applied rewrites87.6%
if -1.46e144 < b < -2.54999999999999989e41Initial program 99.8%
Taylor expanded in x around 0
Applied rewrites99.9%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6499.8
Applied rewrites99.8%
Taylor expanded in t around 0
Applied rewrites89.7%
if -2.54999999999999989e41 < b < 1.0200000000000001e93Initial program 99.3%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in b around 0
Applied rewrites90.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- (+ t y) 2.0) b)))
(if (<= b -1.25e+117)
t_1
(if (<= b -2.55e+41)
(+ (fma -2.0 b (fma (- b z) y z)) a)
(if (<= b 1.02e+93) (fma (- 1.0 t) a (fma (- 1.0 y) z x)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((t + y) - 2.0) * b;
double tmp;
if (b <= -1.25e+117) {
tmp = t_1;
} else if (b <= -2.55e+41) {
tmp = fma(-2.0, b, fma((b - z), y, z)) + a;
} else if (b <= 1.02e+93) {
tmp = fma((1.0 - t), a, fma((1.0 - y), z, x));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(t + y) - 2.0) * b) tmp = 0.0 if (b <= -1.25e+117) tmp = t_1; elseif (b <= -2.55e+41) tmp = Float64(fma(-2.0, b, fma(Float64(b - z), y, z)) + a); elseif (b <= 1.02e+93) tmp = fma(Float64(1.0 - t), a, fma(Float64(1.0 - y), z, x)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[b, -1.25e+117], t$95$1, If[LessEqual[b, -2.55e+41], N[(N[(-2.0 * b + N[(N[(b - z), $MachinePrecision] * y + z), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision], If[LessEqual[b, 1.02e+93], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(1.0 - y), $MachinePrecision] * z + x), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(t + y\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -1.25 \cdot 10^{+117}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -2.55 \cdot 10^{+41}:\\
\;\;\;\;\mathsf{fma}\left(-2, b, \mathsf{fma}\left(b - z, y, z\right)\right) + a\\
\mathbf{elif}\;b \leq 1.02 \cdot 10^{+93}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(1 - y, z, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.24999999999999996e117 or 1.0200000000000001e93 < b Initial program 92.1%
Taylor expanded in x around 0
Applied rewrites94.3%
Taylor expanded in t around 0
Applied rewrites68.6%
Taylor expanded in x around 0
Applied rewrites66.0%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
lower-+.f6486.8
Applied rewrites86.8%
if -1.24999999999999996e117 < b < -2.54999999999999989e41Initial program 99.9%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites93.2%
Taylor expanded in x around 0
Applied rewrites86.4%
if -2.54999999999999989e41 < b < 1.0200000000000001e93Initial program 99.3%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in b around 0
Applied rewrites90.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma (- 1.0 y) z a)) (t_2 (* (- (+ t y) 2.0) b)))
(if (<= b -6e+35)
t_2
(if (<= b -1.85e-230)
t_1
(if (<= b 5.2e-43) (fma (- 1.0 t) a x) (if (<= b 7e+92) t_1 t_2))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma((1.0 - y), z, a);
double t_2 = ((t + y) - 2.0) * b;
double tmp;
if (b <= -6e+35) {
tmp = t_2;
} else if (b <= -1.85e-230) {
tmp = t_1;
} else if (b <= 5.2e-43) {
tmp = fma((1.0 - t), a, x);
} else if (b <= 7e+92) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(Float64(1.0 - y), z, a) t_2 = Float64(Float64(Float64(t + y) - 2.0) * b) tmp = 0.0 if (b <= -6e+35) tmp = t_2; elseif (b <= -1.85e-230) tmp = t_1; elseif (b <= 5.2e-43) tmp = fma(Float64(1.0 - t), a, x); elseif (b <= 7e+92) tmp = t_1; else tmp = t_2; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(1.0 - y), $MachinePrecision] * z + a), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[b, -6e+35], t$95$2, If[LessEqual[b, -1.85e-230], t$95$1, If[LessEqual[b, 5.2e-43], N[(N[(1.0 - t), $MachinePrecision] * a + x), $MachinePrecision], If[LessEqual[b, 7e+92], t$95$1, t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(1 - y, z, a\right)\\
t_2 := \left(\left(t + y\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -6 \cdot 10^{+35}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;b \leq -1.85 \cdot 10^{-230}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 5.2 \cdot 10^{-43}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, x\right)\\
\mathbf{elif}\;b \leq 7 \cdot 10^{+92}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if b < -5.99999999999999981e35 or 6.99999999999999972e92 < b Initial program 93.2%
Taylor expanded in x around 0
Applied rewrites95.2%
Taylor expanded in t around 0
Applied rewrites71.3%
Taylor expanded in x around 0
Applied rewrites68.2%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
lower-+.f6481.4
Applied rewrites81.4%
if -5.99999999999999981e35 < b < -1.84999999999999991e-230 or 5.2e-43 < b < 6.99999999999999972e92Initial program 98.8%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites80.5%
Taylor expanded in x around 0
Applied rewrites71.3%
Taylor expanded in b around 0
Applied rewrites65.2%
if -1.84999999999999991e-230 < b < 5.2e-43Initial program 100.0%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6470.0
Applied rewrites70.0%
Taylor expanded in b around 0
Applied rewrites64.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- b z) y)))
(if (<= y -240000000.0)
t_1
(if (<= y -1e-193)
(+ (+ z x) a)
(if (<= y -5e-265)
(* (- b a) t)
(if (<= y 17.5) (+ (fma -2.0 b z) a) t_1))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b - z) * y;
double tmp;
if (y <= -240000000.0) {
tmp = t_1;
} else if (y <= -1e-193) {
tmp = (z + x) + a;
} else if (y <= -5e-265) {
tmp = (b - a) * t;
} else if (y <= 17.5) {
tmp = fma(-2.0, b, z) + a;
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(b - z) * y) tmp = 0.0 if (y <= -240000000.0) tmp = t_1; elseif (y <= -1e-193) tmp = Float64(Float64(z + x) + a); elseif (y <= -5e-265) tmp = Float64(Float64(b - a) * t); elseif (y <= 17.5) tmp = Float64(fma(-2.0, b, z) + a); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - z), $MachinePrecision] * y), $MachinePrecision]}, If[LessEqual[y, -240000000.0], t$95$1, If[LessEqual[y, -1e-193], N[(N[(z + x), $MachinePrecision] + a), $MachinePrecision], If[LessEqual[y, -5e-265], N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision], If[LessEqual[y, 17.5], N[(N[(-2.0 * b + z), $MachinePrecision] + a), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(b - z\right) \cdot y\\
\mathbf{if}\;y \leq -240000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1 \cdot 10^{-193}:\\
\;\;\;\;\left(z + x\right) + a\\
\mathbf{elif}\;y \leq -5 \cdot 10^{-265}:\\
\;\;\;\;\left(b - a\right) \cdot t\\
\mathbf{elif}\;y \leq 17.5:\\
\;\;\;\;\mathsf{fma}\left(-2, b, z\right) + a\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.4e8 or 17.5 < y Initial program 95.8%
Taylor expanded in y around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6469.4
Applied rewrites69.4%
if -2.4e8 < y < -1e-193Initial program 98.1%
Taylor expanded in x around 0
Applied rewrites98.1%
Taylor expanded in t around 0
Applied rewrites73.7%
Taylor expanded in b around 0
Applied rewrites63.7%
Taylor expanded in y around 0
Applied rewrites62.3%
if -1e-193 < y < -5.0000000000000001e-265Initial program 93.7%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6470.1
Applied rewrites70.1%
if -5.0000000000000001e-265 < y < 17.5Initial program 98.4%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites70.6%
Taylor expanded in x around 0
Applied rewrites58.3%
Taylor expanded in y around 0
Applied rewrites57.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma (- y 2.0) b a)) (t_2 (* (- b a) t)))
(if (<= t -1.08e+40)
t_2
(if (<= t -1.65e-115)
t_1
(if (<= t 5.2e-190) (fma (- 1.0 y) z a) (if (<= t 14600.0) t_1 t_2))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma((y - 2.0), b, a);
double t_2 = (b - a) * t;
double tmp;
if (t <= -1.08e+40) {
tmp = t_2;
} else if (t <= -1.65e-115) {
tmp = t_1;
} else if (t <= 5.2e-190) {
tmp = fma((1.0 - y), z, a);
} else if (t <= 14600.0) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(Float64(y - 2.0), b, a) t_2 = Float64(Float64(b - a) * t) tmp = 0.0 if (t <= -1.08e+40) tmp = t_2; elseif (t <= -1.65e-115) tmp = t_1; elseif (t <= 5.2e-190) tmp = fma(Float64(1.0 - y), z, a); elseif (t <= 14600.0) tmp = t_1; else tmp = t_2; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(y - 2.0), $MachinePrecision] * b + a), $MachinePrecision]}, Block[{t$95$2 = N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[t, -1.08e+40], t$95$2, If[LessEqual[t, -1.65e-115], t$95$1, If[LessEqual[t, 5.2e-190], N[(N[(1.0 - y), $MachinePrecision] * z + a), $MachinePrecision], If[LessEqual[t, 14600.0], t$95$1, t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y - 2, b, a\right)\\
t_2 := \left(b - a\right) \cdot t\\
\mathbf{if}\;t \leq -1.08 \cdot 10^{+40}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t \leq -1.65 \cdot 10^{-115}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 5.2 \cdot 10^{-190}:\\
\;\;\;\;\mathsf{fma}\left(1 - y, z, a\right)\\
\mathbf{elif}\;t \leq 14600:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if t < -1.08000000000000001e40 or 14600 < t Initial program 94.4%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6464.0
Applied rewrites64.0%
if -1.08000000000000001e40 < t < -1.64999999999999995e-115 or 5.1999999999999996e-190 < t < 14600Initial program 97.3%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites94.5%
Taylor expanded in x around 0
Applied rewrites84.2%
Taylor expanded in z around 0
Applied rewrites66.5%
if -1.64999999999999995e-115 < t < 5.1999999999999996e-190Initial program 100.0%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites100.0%
Taylor expanded in x around 0
Applied rewrites85.0%
Taylor expanded in b around 0
Applied rewrites63.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma (- (+ t y) 2.0) b x)))
(if (<= z -1.45e+72)
(fma (- 1.0 y) z t_1)
(if (<= z 2.6e+20)
(fma (- 1.0 t) a t_1)
(fma (- 1.0 t) a (+ (fma (- z) y x) z))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma(((t + y) - 2.0), b, x);
double tmp;
if (z <= -1.45e+72) {
tmp = fma((1.0 - y), z, t_1);
} else if (z <= 2.6e+20) {
tmp = fma((1.0 - t), a, t_1);
} else {
tmp = fma((1.0 - t), a, (fma(-z, y, x) + z));
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(Float64(Float64(t + y) - 2.0), b, x) tmp = 0.0 if (z <= -1.45e+72) tmp = fma(Float64(1.0 - y), z, t_1); elseif (z <= 2.6e+20) tmp = fma(Float64(1.0 - t), a, t_1); else tmp = fma(Float64(1.0 - t), a, Float64(fma(Float64(-z), y, x) + z)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b + x), $MachinePrecision]}, If[LessEqual[z, -1.45e+72], N[(N[(1.0 - y), $MachinePrecision] * z + t$95$1), $MachinePrecision], If[LessEqual[z, 2.6e+20], N[(N[(1.0 - t), $MachinePrecision] * a + t$95$1), $MachinePrecision], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[((-z) * y + x), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\left(t + y\right) - 2, b, x\right)\\
\mathbf{if}\;z \leq -1.45 \cdot 10^{+72}:\\
\;\;\;\;\mathsf{fma}\left(1 - y, z, t\_1\right)\\
\mathbf{elif}\;z \leq 2.6 \cdot 10^{+20}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, t\_1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(-z, y, x\right) + z\right)\\
\end{array}
\end{array}
if z < -1.45000000000000009e72Initial program 96.0%
Taylor expanded in x around 0
Applied rewrites96.1%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6490.4
Applied rewrites90.4%
if -1.45000000000000009e72 < z < 2.6e20Initial program 98.0%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6495.2
Applied rewrites95.2%
if 2.6e20 < z Initial program 93.8%
Taylor expanded in x around 0
Applied rewrites98.0%
Taylor expanded in b around 0
Applied rewrites91.8%
Final simplification93.6%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -4.2e+75)
(fma (- 1.0 y) z (+ (fma (- y 2.0) b x) a))
(if (<= z 2.6e+20)
(fma (- 1.0 t) a (fma (- (+ t y) 2.0) b x))
(fma (- 1.0 t) a (+ (fma (- z) y x) z)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -4.2e+75) {
tmp = fma((1.0 - y), z, (fma((y - 2.0), b, x) + a));
} else if (z <= 2.6e+20) {
tmp = fma((1.0 - t), a, fma(((t + y) - 2.0), b, x));
} else {
tmp = fma((1.0 - t), a, (fma(-z, y, x) + z));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -4.2e+75) tmp = fma(Float64(1.0 - y), z, Float64(fma(Float64(y - 2.0), b, x) + a)); elseif (z <= 2.6e+20) tmp = fma(Float64(1.0 - t), a, fma(Float64(Float64(t + y) - 2.0), b, x)); else tmp = fma(Float64(1.0 - t), a, Float64(fma(Float64(-z), y, x) + z)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -4.2e+75], N[(N[(1.0 - y), $MachinePrecision] * z + N[(N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision] + a), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2.6e+20], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[((-z) * y + x), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.2 \cdot 10^{+75}:\\
\;\;\;\;\mathsf{fma}\left(1 - y, z, \mathsf{fma}\left(y - 2, b, x\right) + a\right)\\
\mathbf{elif}\;z \leq 2.6 \cdot 10^{+20}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(\left(t + y\right) - 2, b, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(-z, y, x\right) + z\right)\\
\end{array}
\end{array}
if z < -4.19999999999999997e75Initial program 96.0%
Taylor expanded in t around 0
sub-negN/A
+-commutativeN/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
associate-+l+N/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-+.f64N/A
Applied rewrites87.0%
if -4.19999999999999997e75 < z < 2.6e20Initial program 98.0%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6495.2
Applied rewrites95.2%
if 2.6e20 < z Initial program 93.8%
Taylor expanded in x around 0
Applied rewrites98.0%
Taylor expanded in b around 0
Applied rewrites91.8%
Final simplification92.9%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -4.2e+75)
(+ (fma -2.0 b (fma (- b z) y z)) (+ a x))
(if (<= z 2.6e+20)
(fma (- 1.0 t) a (fma (- (+ t y) 2.0) b x))
(fma (- 1.0 t) a (+ (fma (- z) y x) z)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -4.2e+75) {
tmp = fma(-2.0, b, fma((b - z), y, z)) + (a + x);
} else if (z <= 2.6e+20) {
tmp = fma((1.0 - t), a, fma(((t + y) - 2.0), b, x));
} else {
tmp = fma((1.0 - t), a, (fma(-z, y, x) + z));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -4.2e+75) tmp = Float64(fma(-2.0, b, fma(Float64(b - z), y, z)) + Float64(a + x)); elseif (z <= 2.6e+20) tmp = fma(Float64(1.0 - t), a, fma(Float64(Float64(t + y) - 2.0), b, x)); else tmp = fma(Float64(1.0 - t), a, Float64(fma(Float64(-z), y, x) + z)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -4.2e+75], N[(N[(-2.0 * b + N[(N[(b - z), $MachinePrecision] * y + z), $MachinePrecision]), $MachinePrecision] + N[(a + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2.6e+20], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[((-z) * y + x), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.2 \cdot 10^{+75}:\\
\;\;\;\;\mathsf{fma}\left(-2, b, \mathsf{fma}\left(b - z, y, z\right)\right) + \left(a + x\right)\\
\mathbf{elif}\;z \leq 2.6 \cdot 10^{+20}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(\left(t + y\right) - 2, b, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(-z, y, x\right) + z\right)\\
\end{array}
\end{array}
if z < -4.19999999999999997e75Initial program 96.0%
Taylor expanded in x around 0
Applied rewrites96.1%
Taylor expanded in t around 0
Applied rewrites87.0%
if -4.19999999999999997e75 < z < 2.6e20Initial program 98.0%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6495.2
Applied rewrites95.2%
if 2.6e20 < z Initial program 93.8%
Taylor expanded in x around 0
Applied rewrites98.0%
Taylor expanded in b around 0
Applied rewrites91.8%
Final simplification92.9%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -1.9e+39)
(+ (fma -2.0 b (fma (- b z) y z)) (+ a x))
(if (<= z 2e+19)
(+ (fma (- b a) t a) (fma (- y 2.0) b x))
(fma (- 1.0 t) a (+ (fma (- z) y x) z)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -1.9e+39) {
tmp = fma(-2.0, b, fma((b - z), y, z)) + (a + x);
} else if (z <= 2e+19) {
tmp = fma((b - a), t, a) + fma((y - 2.0), b, x);
} else {
tmp = fma((1.0 - t), a, (fma(-z, y, x) + z));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -1.9e+39) tmp = Float64(fma(-2.0, b, fma(Float64(b - z), y, z)) + Float64(a + x)); elseif (z <= 2e+19) tmp = Float64(fma(Float64(b - a), t, a) + fma(Float64(y - 2.0), b, x)); else tmp = fma(Float64(1.0 - t), a, Float64(fma(Float64(-z), y, x) + z)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -1.9e+39], N[(N[(-2.0 * b + N[(N[(b - z), $MachinePrecision] * y + z), $MachinePrecision]), $MachinePrecision] + N[(a + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2e+19], N[(N[(N[(b - a), $MachinePrecision] * t + a), $MachinePrecision] + N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[((-z) * y + x), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.9 \cdot 10^{+39}:\\
\;\;\;\;\mathsf{fma}\left(-2, b, \mathsf{fma}\left(b - z, y, z\right)\right) + \left(a + x\right)\\
\mathbf{elif}\;z \leq 2 \cdot 10^{+19}:\\
\;\;\;\;\mathsf{fma}\left(b - a, t, a\right) + \mathsf{fma}\left(y - 2, b, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(-z, y, x\right) + z\right)\\
\end{array}
\end{array}
if z < -1.8999999999999999e39Initial program 96.3%
Taylor expanded in x around 0
Applied rewrites94.5%
Taylor expanded in t around 0
Applied rewrites86.1%
if -1.8999999999999999e39 < z < 2e19Initial program 98.0%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6495.0
Applied rewrites95.0%
Taylor expanded in t around 0
Applied rewrites95.0%
if 2e19 < z Initial program 93.8%
Taylor expanded in x around 0
Applied rewrites98.0%
Taylor expanded in b around 0
Applied rewrites91.8%
Final simplification92.5%
(FPCore (x y z t a b)
:precision binary64
(if (<= t -8.8e+116)
(fma t b (* (- a) t))
(if (<= t 1.28e+106)
(+ (fma -2.0 b (fma (- b z) y z)) (+ a x))
(fma (- 1.0 t) a (fma (- t 2.0) b x)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -8.8e+116) {
tmp = fma(t, b, (-a * t));
} else if (t <= 1.28e+106) {
tmp = fma(-2.0, b, fma((b - z), y, z)) + (a + x);
} else {
tmp = fma((1.0 - t), a, fma((t - 2.0), b, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -8.8e+116) tmp = fma(t, b, Float64(Float64(-a) * t)); elseif (t <= 1.28e+106) tmp = Float64(fma(-2.0, b, fma(Float64(b - z), y, z)) + Float64(a + x)); else tmp = fma(Float64(1.0 - t), a, fma(Float64(t - 2.0), b, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -8.8e+116], N[(t * b + N[((-a) * t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.28e+106], N[(N[(-2.0 * b + N[(N[(b - z), $MachinePrecision] * y + z), $MachinePrecision]), $MachinePrecision] + N[(a + x), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(t - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -8.8 \cdot 10^{+116}:\\
\;\;\;\;\mathsf{fma}\left(t, b, \left(-a\right) \cdot t\right)\\
\mathbf{elif}\;t \leq 1.28 \cdot 10^{+106}:\\
\;\;\;\;\mathsf{fma}\left(-2, b, \mathsf{fma}\left(b - z, y, z\right)\right) + \left(a + x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(t - 2, b, x\right)\right)\\
\end{array}
\end{array}
if t < -8.799999999999999e116Initial program 95.2%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6474.7
Applied rewrites74.7%
Applied rewrites77.1%
if -8.799999999999999e116 < t < 1.28e106Initial program 98.8%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites94.3%
if 1.28e106 < t Initial program 89.7%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6487.8
Applied rewrites87.8%
Taylor expanded in y around 0
Applied rewrites85.2%
Final simplification90.1%
(FPCore (x y z t a b)
:precision binary64
(if (<= t -3.4e+116)
(fma t b (* (- a) t))
(if (<= t 170000000.0)
(+ (fma -2.0 b (fma (- b z) y z)) a)
(if (<= t 5.2e+148) (+ (fma (- t 2.0) b z) x) (* (- b a) t)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -3.4e+116) {
tmp = fma(t, b, (-a * t));
} else if (t <= 170000000.0) {
tmp = fma(-2.0, b, fma((b - z), y, z)) + a;
} else if (t <= 5.2e+148) {
tmp = fma((t - 2.0), b, z) + x;
} else {
tmp = (b - a) * t;
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -3.4e+116) tmp = fma(t, b, Float64(Float64(-a) * t)); elseif (t <= 170000000.0) tmp = Float64(fma(-2.0, b, fma(Float64(b - z), y, z)) + a); elseif (t <= 5.2e+148) tmp = Float64(fma(Float64(t - 2.0), b, z) + x); else tmp = Float64(Float64(b - a) * t); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -3.4e+116], N[(t * b + N[((-a) * t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 170000000.0], N[(N[(-2.0 * b + N[(N[(b - z), $MachinePrecision] * y + z), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision], If[LessEqual[t, 5.2e+148], N[(N[(N[(t - 2.0), $MachinePrecision] * b + z), $MachinePrecision] + x), $MachinePrecision], N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -3.4 \cdot 10^{+116}:\\
\;\;\;\;\mathsf{fma}\left(t, b, \left(-a\right) \cdot t\right)\\
\mathbf{elif}\;t \leq 170000000:\\
\;\;\;\;\mathsf{fma}\left(-2, b, \mathsf{fma}\left(b - z, y, z\right)\right) + a\\
\mathbf{elif}\;t \leq 5.2 \cdot 10^{+148}:\\
\;\;\;\;\mathsf{fma}\left(t - 2, b, z\right) + x\\
\mathbf{else}:\\
\;\;\;\;\left(b - a\right) \cdot t\\
\end{array}
\end{array}
if t < -3.40000000000000023e116Initial program 95.2%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6474.7
Applied rewrites74.7%
Applied rewrites77.1%
if -3.40000000000000023e116 < t < 1.7e8Initial program 98.7%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites96.2%
Taylor expanded in x around 0
Applied rewrites83.2%
if 1.7e8 < t < 5.2e148Initial program 99.9%
Taylor expanded in x around 0
Applied rewrites99.9%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6491.3
Applied rewrites91.3%
Taylor expanded in y around 0
Applied rewrites74.9%
if 5.2e148 < t Initial program 87.5%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6479.1
Applied rewrites79.1%
(FPCore (x y z t a b)
:precision binary64
(if (<= t -3.4e+116)
(fma t b (* (- a) t))
(if (<= t 3.3e+49)
(+ (fma -2.0 b (fma (- b z) y z)) a)
(fma (- 1.0 t) a (fma (- t 2.0) b x)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -3.4e+116) {
tmp = fma(t, b, (-a * t));
} else if (t <= 3.3e+49) {
tmp = fma(-2.0, b, fma((b - z), y, z)) + a;
} else {
tmp = fma((1.0 - t), a, fma((t - 2.0), b, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -3.4e+116) tmp = fma(t, b, Float64(Float64(-a) * t)); elseif (t <= 3.3e+49) tmp = Float64(fma(-2.0, b, fma(Float64(b - z), y, z)) + a); else tmp = fma(Float64(1.0 - t), a, fma(Float64(t - 2.0), b, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -3.4e+116], N[(t * b + N[((-a) * t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 3.3e+49], N[(N[(-2.0 * b + N[(N[(b - z), $MachinePrecision] * y + z), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(t - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -3.4 \cdot 10^{+116}:\\
\;\;\;\;\mathsf{fma}\left(t, b, \left(-a\right) \cdot t\right)\\
\mathbf{elif}\;t \leq 3.3 \cdot 10^{+49}:\\
\;\;\;\;\mathsf{fma}\left(-2, b, \mathsf{fma}\left(b - z, y, z\right)\right) + a\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(t - 2, b, x\right)\right)\\
\end{array}
\end{array}
if t < -3.40000000000000023e116Initial program 95.2%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6474.7
Applied rewrites74.7%
Applied rewrites77.1%
if -3.40000000000000023e116 < t < 3.2999999999999998e49Initial program 98.8%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites94.6%
Taylor expanded in x around 0
Applied rewrites81.8%
if 3.2999999999999998e49 < t Initial program 91.3%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6485.5
Applied rewrites85.5%
Taylor expanded in y around 0
Applied rewrites83.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma (- 1.0 y) z a)))
(if (<= z -5.8e+29)
t_1
(if (<= z -2.35e-295)
(fma (- y 2.0) b a)
(if (<= z 8.5e+16) (fma (- 1.0 t) a x) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma((1.0 - y), z, a);
double tmp;
if (z <= -5.8e+29) {
tmp = t_1;
} else if (z <= -2.35e-295) {
tmp = fma((y - 2.0), b, a);
} else if (z <= 8.5e+16) {
tmp = fma((1.0 - t), a, x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(Float64(1.0 - y), z, a) tmp = 0.0 if (z <= -5.8e+29) tmp = t_1; elseif (z <= -2.35e-295) tmp = fma(Float64(y - 2.0), b, a); elseif (z <= 8.5e+16) tmp = fma(Float64(1.0 - t), a, x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(1.0 - y), $MachinePrecision] * z + a), $MachinePrecision]}, If[LessEqual[z, -5.8e+29], t$95$1, If[LessEqual[z, -2.35e-295], N[(N[(y - 2.0), $MachinePrecision] * b + a), $MachinePrecision], If[LessEqual[z, 8.5e+16], N[(N[(1.0 - t), $MachinePrecision] * a + x), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(1 - y, z, a\right)\\
\mathbf{if}\;z \leq -5.8 \cdot 10^{+29}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -2.35 \cdot 10^{-295}:\\
\;\;\;\;\mathsf{fma}\left(y - 2, b, a\right)\\
\mathbf{elif}\;z \leq 8.5 \cdot 10^{+16}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -5.7999999999999999e29 or 8.5e16 < z Initial program 95.4%
Taylor expanded in x around 0
Applied rewrites96.3%
Taylor expanded in t around 0
Applied rewrites85.1%
Taylor expanded in x around 0
Applied rewrites80.7%
Taylor expanded in b around 0
Applied rewrites72.4%
if -5.7999999999999999e29 < z < -2.3499999999999999e-295Initial program 98.5%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites70.2%
Taylor expanded in x around 0
Applied rewrites61.3%
Taylor expanded in z around 0
Applied rewrites54.8%
if -2.3499999999999999e-295 < z < 8.5e16Initial program 97.4%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6498.7
Applied rewrites98.7%
Taylor expanded in b around 0
Applied rewrites60.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- b a) t)))
(if (<= t -1.08e+40)
t_1
(if (<= t -1.5e-7)
(* (- y 2.0) b)
(if (<= t 4.6e+61) (+ (+ z x) a) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b - a) * t;
double tmp;
if (t <= -1.08e+40) {
tmp = t_1;
} else if (t <= -1.5e-7) {
tmp = (y - 2.0) * b;
} else if (t <= 4.6e+61) {
tmp = (z + x) + a;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = (b - a) * t
if (t <= (-1.08d+40)) then
tmp = t_1
else if (t <= (-1.5d-7)) then
tmp = (y - 2.0d0) * b
else if (t <= 4.6d+61) then
tmp = (z + x) + a
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 t_1 = (b - a) * t;
double tmp;
if (t <= -1.08e+40) {
tmp = t_1;
} else if (t <= -1.5e-7) {
tmp = (y - 2.0) * b;
} else if (t <= 4.6e+61) {
tmp = (z + x) + a;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (b - a) * t tmp = 0 if t <= -1.08e+40: tmp = t_1 elif t <= -1.5e-7: tmp = (y - 2.0) * b elif t <= 4.6e+61: tmp = (z + x) + a else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(b - a) * t) tmp = 0.0 if (t <= -1.08e+40) tmp = t_1; elseif (t <= -1.5e-7) tmp = Float64(Float64(y - 2.0) * b); elseif (t <= 4.6e+61) tmp = Float64(Float64(z + x) + a); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (b - a) * t; tmp = 0.0; if (t <= -1.08e+40) tmp = t_1; elseif (t <= -1.5e-7) tmp = (y - 2.0) * b; elseif (t <= 4.6e+61) tmp = (z + x) + a; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[t, -1.08e+40], t$95$1, If[LessEqual[t, -1.5e-7], N[(N[(y - 2.0), $MachinePrecision] * b), $MachinePrecision], If[LessEqual[t, 4.6e+61], N[(N[(z + x), $MachinePrecision] + a), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(b - a\right) \cdot t\\
\mathbf{if}\;t \leq -1.08 \cdot 10^{+40}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -1.5 \cdot 10^{-7}:\\
\;\;\;\;\left(y - 2\right) \cdot b\\
\mathbf{elif}\;t \leq 4.6 \cdot 10^{+61}:\\
\;\;\;\;\left(z + x\right) + a\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -1.08000000000000001e40 or 4.5999999999999999e61 < t Initial program 93.7%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6468.5
Applied rewrites68.5%
if -1.08000000000000001e40 < t < -1.4999999999999999e-7Initial program 99.9%
Taylor expanded in x around 0
Applied rewrites99.9%
Taylor expanded in t around 0
Applied rewrites71.0%
Taylor expanded in b around inf
Applied rewrites63.3%
if -1.4999999999999999e-7 < t < 4.5999999999999999e61Initial program 98.6%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites97.7%
Taylor expanded in b around 0
Applied rewrites68.9%
Taylor expanded in y around 0
Applied rewrites52.8%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* (- (+ t y) 2.0) b))) (if (<= b -9e+116) t_1 (if (<= b 7e+92) (fma (- 1.0 y) z (+ a x)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((t + y) - 2.0) * b;
double tmp;
if (b <= -9e+116) {
tmp = t_1;
} else if (b <= 7e+92) {
tmp = fma((1.0 - y), z, (a + x));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(t + y) - 2.0) * b) tmp = 0.0 if (b <= -9e+116) tmp = t_1; elseif (b <= 7e+92) tmp = fma(Float64(1.0 - y), z, Float64(a + x)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[b, -9e+116], t$95$1, If[LessEqual[b, 7e+92], N[(N[(1.0 - y), $MachinePrecision] * z + N[(a + x), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(t + y\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -9 \cdot 10^{+116}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 7 \cdot 10^{+92}:\\
\;\;\;\;\mathsf{fma}\left(1 - y, z, a + x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -9.00000000000000032e116 or 6.99999999999999972e92 < b Initial program 92.1%
Taylor expanded in x around 0
Applied rewrites94.3%
Taylor expanded in t around 0
Applied rewrites68.6%
Taylor expanded in x around 0
Applied rewrites66.0%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
lower-+.f6486.8
Applied rewrites86.8%
if -9.00000000000000032e116 < b < 6.99999999999999972e92Initial program 99.4%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites79.6%
Taylor expanded in b around 0
Applied rewrites71.7%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* (- b z) y))) (if (<= y -240000000.0) t_1 (if (<= y 320000000.0) (+ (+ z x) a) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b - z) * y;
double tmp;
if (y <= -240000000.0) {
tmp = t_1;
} else if (y <= 320000000.0) {
tmp = (z + x) + a;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = (b - z) * y
if (y <= (-240000000.0d0)) then
tmp = t_1
else if (y <= 320000000.0d0) then
tmp = (z + x) + a
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 t_1 = (b - z) * y;
double tmp;
if (y <= -240000000.0) {
tmp = t_1;
} else if (y <= 320000000.0) {
tmp = (z + x) + a;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (b - z) * y tmp = 0 if y <= -240000000.0: tmp = t_1 elif y <= 320000000.0: tmp = (z + x) + a else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(b - z) * y) tmp = 0.0 if (y <= -240000000.0) tmp = t_1; elseif (y <= 320000000.0) tmp = Float64(Float64(z + x) + a); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (b - z) * y; tmp = 0.0; if (y <= -240000000.0) tmp = t_1; elseif (y <= 320000000.0) tmp = (z + x) + a; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - z), $MachinePrecision] * y), $MachinePrecision]}, If[LessEqual[y, -240000000.0], t$95$1, If[LessEqual[y, 320000000.0], N[(N[(z + x), $MachinePrecision] + a), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(b - z\right) \cdot y\\
\mathbf{if}\;y \leq -240000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 320000000:\\
\;\;\;\;\left(z + x\right) + a\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.4e8 or 3.2e8 < y Initial program 95.8%
Taylor expanded in y around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6470.0
Applied rewrites70.0%
if -2.4e8 < y < 3.2e8Initial program 97.7%
Taylor expanded in x around 0
Applied rewrites99.2%
Taylor expanded in t around 0
Applied rewrites67.4%
Taylor expanded in b around 0
Applied rewrites54.3%
Taylor expanded in y around 0
Applied rewrites53.8%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* (- y 2.0) b))) (if (<= b -1.8e+86) t_1 (if (<= b 2.8e+92) (+ (+ z x) a) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (y - 2.0) * b;
double tmp;
if (b <= -1.8e+86) {
tmp = t_1;
} else if (b <= 2.8e+92) {
tmp = (z + x) + a;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = (y - 2.0d0) * b
if (b <= (-1.8d+86)) then
tmp = t_1
else if (b <= 2.8d+92) then
tmp = (z + x) + a
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 t_1 = (y - 2.0) * b;
double tmp;
if (b <= -1.8e+86) {
tmp = t_1;
} else if (b <= 2.8e+92) {
tmp = (z + x) + a;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (y - 2.0) * b tmp = 0 if b <= -1.8e+86: tmp = t_1 elif b <= 2.8e+92: tmp = (z + x) + a else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(y - 2.0) * b) tmp = 0.0 if (b <= -1.8e+86) tmp = t_1; elseif (b <= 2.8e+92) tmp = Float64(Float64(z + x) + a); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (y - 2.0) * b; tmp = 0.0; if (b <= -1.8e+86) tmp = t_1; elseif (b <= 2.8e+92) tmp = (z + x) + a; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(y - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[b, -1.8e+86], t$95$1, If[LessEqual[b, 2.8e+92], N[(N[(z + x), $MachinePrecision] + a), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(y - 2\right) \cdot b\\
\mathbf{if}\;b \leq -1.8 \cdot 10^{+86}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 2.8 \cdot 10^{+92}:\\
\;\;\;\;\left(z + x\right) + a\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.80000000000000003e86 or 2.80000000000000001e92 < b Initial program 92.6%
Taylor expanded in x around 0
Applied rewrites94.7%
Taylor expanded in t around 0
Applied rewrites70.6%
Taylor expanded in b around inf
Applied rewrites57.3%
if -1.80000000000000003e86 < b < 2.80000000000000001e92Initial program 99.3%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites78.8%
Taylor expanded in b around 0
Applied rewrites71.9%
Taylor expanded in y around 0
Applied rewrites51.3%
(FPCore (x y z t a b) :precision binary64 (if (<= y -50000000000.0) (* b y) (if (<= y 4.8e+36) (+ (+ z x) a) (* (- z) y))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -50000000000.0) {
tmp = b * y;
} else if (y <= 4.8e+36) {
tmp = (z + x) + a;
} else {
tmp = -z * y;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (y <= (-50000000000.0d0)) then
tmp = b * y
else if (y <= 4.8d+36) then
tmp = (z + x) + a
else
tmp = -z * y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -50000000000.0) {
tmp = b * y;
} else if (y <= 4.8e+36) {
tmp = (z + x) + a;
} else {
tmp = -z * y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -50000000000.0: tmp = b * y elif y <= 4.8e+36: tmp = (z + x) + a else: tmp = -z * y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -50000000000.0) tmp = Float64(b * y); elseif (y <= 4.8e+36) tmp = Float64(Float64(z + x) + a); else tmp = Float64(Float64(-z) * y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -50000000000.0) tmp = b * y; elseif (y <= 4.8e+36) tmp = (z + x) + a; else tmp = -z * y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -50000000000.0], N[(b * y), $MachinePrecision], If[LessEqual[y, 4.8e+36], N[(N[(z + x), $MachinePrecision] + a), $MachinePrecision], N[((-z) * y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -50000000000:\\
\;\;\;\;b \cdot y\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{+36}:\\
\;\;\;\;\left(z + x\right) + a\\
\mathbf{else}:\\
\;\;\;\;\left(-z\right) \cdot y\\
\end{array}
\end{array}
if y < -5e10Initial program 95.0%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6479.6
Applied rewrites79.6%
Taylor expanded in y around inf
Applied rewrites44.5%
if -5e10 < y < 4.79999999999999985e36Initial program 97.9%
Taylor expanded in x around 0
Applied rewrites99.3%
Taylor expanded in t around 0
Applied rewrites67.0%
Taylor expanded in b around 0
Applied rewrites53.6%
Taylor expanded in y around 0
Applied rewrites52.4%
if 4.79999999999999985e36 < y Initial program 96.1%
Taylor expanded in y around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6480.9
Applied rewrites80.9%
Taylor expanded in z around inf
Applied rewrites51.9%
(FPCore (x y z t a b) :precision binary64 (if (<= b -4.8e+35) (* b y) (if (<= b 7.5e+92) (+ a x) (* b y))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -4.8e+35) {
tmp = b * y;
} else if (b <= 7.5e+92) {
tmp = a + x;
} else {
tmp = b * y;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (b <= (-4.8d+35)) then
tmp = b * y
else if (b <= 7.5d+92) then
tmp = a + x
else
tmp = b * y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -4.8e+35) {
tmp = b * y;
} else if (b <= 7.5e+92) {
tmp = a + x;
} else {
tmp = b * y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -4.8e+35: tmp = b * y elif b <= 7.5e+92: tmp = a + x else: tmp = b * y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -4.8e+35) tmp = Float64(b * y); elseif (b <= 7.5e+92) tmp = Float64(a + x); else tmp = Float64(b * y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -4.8e+35) tmp = b * y; elseif (b <= 7.5e+92) tmp = a + x; else tmp = b * y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -4.8e+35], N[(b * y), $MachinePrecision], If[LessEqual[b, 7.5e+92], N[(a + x), $MachinePrecision], N[(b * y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -4.8 \cdot 10^{+35}:\\
\;\;\;\;b \cdot y\\
\mathbf{elif}\;b \leq 7.5 \cdot 10^{+92}:\\
\;\;\;\;a + x\\
\mathbf{else}:\\
\;\;\;\;b \cdot y\\
\end{array}
\end{array}
if b < -4.80000000000000029e35 or 7.49999999999999946e92 < b Initial program 93.2%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-+.f6488.5
Applied rewrites88.5%
Taylor expanded in y around inf
Applied rewrites38.9%
if -4.80000000000000029e35 < b < 7.49999999999999946e92Initial program 99.3%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites78.8%
Taylor expanded in b around 0
Applied rewrites73.5%
Taylor expanded in z around 0
Applied rewrites35.2%
(FPCore (x y z t a b) :precision binary64 (if (<= b -9.5e+155) (* -2.0 b) (if (<= b 2.05e+93) (+ a x) (* -2.0 b))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -9.5e+155) {
tmp = -2.0 * b;
} else if (b <= 2.05e+93) {
tmp = a + x;
} else {
tmp = -2.0 * b;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (b <= (-9.5d+155)) then
tmp = (-2.0d0) * b
else if (b <= 2.05d+93) then
tmp = a + x
else
tmp = (-2.0d0) * b
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -9.5e+155) {
tmp = -2.0 * b;
} else if (b <= 2.05e+93) {
tmp = a + x;
} else {
tmp = -2.0 * b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -9.5e+155: tmp = -2.0 * b elif b <= 2.05e+93: tmp = a + x else: tmp = -2.0 * b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -9.5e+155) tmp = Float64(-2.0 * b); elseif (b <= 2.05e+93) tmp = Float64(a + x); else tmp = Float64(-2.0 * b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -9.5e+155) tmp = -2.0 * b; elseif (b <= 2.05e+93) tmp = a + x; else tmp = -2.0 * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -9.5e+155], N[(-2.0 * b), $MachinePrecision], If[LessEqual[b, 2.05e+93], N[(a + x), $MachinePrecision], N[(-2.0 * b), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -9.5 \cdot 10^{+155}:\\
\;\;\;\;-2 \cdot b\\
\mathbf{elif}\;b \leq 2.05 \cdot 10^{+93}:\\
\;\;\;\;a + x\\
\mathbf{else}:\\
\;\;\;\;-2 \cdot b\\
\end{array}
\end{array}
if b < -9.5000000000000006e155 or 2.0500000000000001e93 < b Initial program 91.4%
Taylor expanded in x around 0
Applied rewrites93.9%
Taylor expanded in t around 0
Applied rewrites68.4%
Taylor expanded in b around inf
Applied rewrites59.1%
Taylor expanded in y around 0
Applied rewrites20.3%
if -9.5000000000000006e155 < b < 2.0500000000000001e93Initial program 99.4%
Taylor expanded in x around 0
Applied rewrites100.0%
Taylor expanded in t around 0
Applied rewrites79.2%
Taylor expanded in b around 0
Applied rewrites70.0%
Taylor expanded in z around 0
Applied rewrites32.3%
(FPCore (x y z t a b) :precision binary64 (+ a x))
double code(double x, double y, double z, double t, double a, double b) {
return a + x;
}
real(8) function code(x, y, z, t, a, b)
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
code = a + x
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return a + x;
}
def code(x, y, z, t, a, b): return a + x
function code(x, y, z, t, a, b) return Float64(a + x) end
function tmp = code(x, y, z, t, a, b) tmp = a + x; end
code[x_, y_, z_, t_, a_, b_] := N[(a + x), $MachinePrecision]
\begin{array}{l}
\\
a + x
\end{array}
Initial program 96.8%
Taylor expanded in x around 0
Applied rewrites98.0%
Taylor expanded in t around 0
Applied rewrites75.8%
Taylor expanded in b around 0
Applied rewrites53.4%
Taylor expanded in z around 0
Applied rewrites24.2%
herbie shell --seed 2024331
(FPCore (x y z t a b)
:name "Statistics.Distribution.Beta:$centropy from math-functions-0.1.5.2"
:precision binary64
(+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))