
(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 20 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
(let* ((t_1
(+ (+ (+ x (* z (- 1.0 y))) (* a (- 1.0 t))) (* (- (+ y t) 2.0) b))))
(if (<= t_1 INFINITY) t_1 (* t (- b a)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b);
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = t * (b - a);
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b);
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = t * (b - a);
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b) tmp = 0 if t_1 <= math.inf: tmp = t_1 else: tmp = t * (b - a) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(x + Float64(z * Float64(1.0 - y))) + Float64(a * Float64(1.0 - t))) + Float64(Float64(Float64(y + t) - 2.0) * b)) tmp = 0.0 if (t_1 <= Inf) tmp = t_1; else tmp = Float64(t * Float64(b - a)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b); tmp = 0.0; if (t_1 <= Inf) tmp = t_1; else tmp = t * (b - a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, Infinity], t$95$1, N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(x + z \cdot \left(1 - y\right)\right) + a \cdot \left(1 - t\right)\right) + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{if}\;t\_1 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t \cdot \left(b - a\right)\\
\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)) < +inf.0Initial program 100.0%
if +inf.0 < (+.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 0.0%
Taylor expanded in t around inf
lower-*.f64N/A
lower--.f6471.5
Applied rewrites71.5%
Final simplification98.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+ (+ (+ x (* z (- 1.0 y))) (* a (- 1.0 t))) (* (- (+ y t) 2.0) b))))
(if (<= t_1 -2e+301) (* t b) (if (<= t_1 2e+302) (+ z (+ x a)) (* t b)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b);
double tmp;
if (t_1 <= -2e+301) {
tmp = t * b;
} else if (t_1 <= 2e+302) {
tmp = z + (x + a);
} else {
tmp = t * 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) :: t_1
real(8) :: tmp
t_1 = ((x + (z * (1.0d0 - y))) + (a * (1.0d0 - t))) + (((y + t) - 2.0d0) * b)
if (t_1 <= (-2d+301)) then
tmp = t * b
else if (t_1 <= 2d+302) then
tmp = z + (x + a)
else
tmp = t * b
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 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b);
double tmp;
if (t_1 <= -2e+301) {
tmp = t * b;
} else if (t_1 <= 2e+302) {
tmp = z + (x + a);
} else {
tmp = t * b;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b) tmp = 0 if t_1 <= -2e+301: tmp = t * b elif t_1 <= 2e+302: tmp = z + (x + a) else: tmp = t * b return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(x + Float64(z * Float64(1.0 - y))) + Float64(a * Float64(1.0 - t))) + Float64(Float64(Float64(y + t) - 2.0) * b)) tmp = 0.0 if (t_1 <= -2e+301) tmp = Float64(t * b); elseif (t_1 <= 2e+302) tmp = Float64(z + Float64(x + a)); else tmp = Float64(t * b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b); tmp = 0.0; if (t_1 <= -2e+301) tmp = t * b; elseif (t_1 <= 2e+302) tmp = z + (x + a); else tmp = t * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e+301], N[(t * b), $MachinePrecision], If[LessEqual[t$95$1, 2e+302], N[(z + N[(x + a), $MachinePrecision]), $MachinePrecision], N[(t * b), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(x + z \cdot \left(1 - y\right)\right) + a \cdot \left(1 - t\right)\right) + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{+301}:\\
\;\;\;\;t \cdot b\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{+302}:\\
\;\;\;\;z + \left(x + a\right)\\
\mathbf{else}:\\
\;\;\;\;t \cdot b\\
\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)) < -2.00000000000000011e301 or 2.0000000000000002e302 < (+.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 87.0%
Taylor expanded in a around 0
+-commutativeN/A
associate--l+N/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-evalN/A
sub-negN/A
+-commutativeN/A
distribute-rgt-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--.f6474.7
Applied rewrites74.7%
Taylor expanded in t around inf
Applied rewrites33.0%
if -2.00000000000000011e301 < (+.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)) < 2.0000000000000002e302Initial program 100.0%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites77.6%
Taylor expanded in t around 0
Applied rewrites69.6%
Taylor expanded in y around 0
Applied rewrites58.7%
Final simplification47.9%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+ (+ (+ x (* z (- 1.0 y))) (* a (- 1.0 t))) (* (- (+ y t) 2.0) b))))
(if (<= t_1 -1e+290) (* t b) (if (<= t_1 2e+302) (+ x a) (* t b)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b);
double tmp;
if (t_1 <= -1e+290) {
tmp = t * b;
} else if (t_1 <= 2e+302) {
tmp = x + a;
} else {
tmp = t * 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) :: t_1
real(8) :: tmp
t_1 = ((x + (z * (1.0d0 - y))) + (a * (1.0d0 - t))) + (((y + t) - 2.0d0) * b)
if (t_1 <= (-1d+290)) then
tmp = t * b
else if (t_1 <= 2d+302) then
tmp = x + a
else
tmp = t * b
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 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b);
double tmp;
if (t_1 <= -1e+290) {
tmp = t * b;
} else if (t_1 <= 2e+302) {
tmp = x + a;
} else {
tmp = t * b;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b) tmp = 0 if t_1 <= -1e+290: tmp = t * b elif t_1 <= 2e+302: tmp = x + a else: tmp = t * b return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(x + Float64(z * Float64(1.0 - y))) + Float64(a * Float64(1.0 - t))) + Float64(Float64(Float64(y + t) - 2.0) * b)) tmp = 0.0 if (t_1 <= -1e+290) tmp = Float64(t * b); elseif (t_1 <= 2e+302) tmp = Float64(x + a); else tmp = Float64(t * b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (((y + t) - 2.0) * b); tmp = 0.0; if (t_1 <= -1e+290) tmp = t * b; elseif (t_1 <= 2e+302) tmp = x + a; else tmp = t * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -1e+290], N[(t * b), $MachinePrecision], If[LessEqual[t$95$1, 2e+302], N[(x + a), $MachinePrecision], N[(t * b), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(x + z \cdot \left(1 - y\right)\right) + a \cdot \left(1 - t\right)\right) + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{+290}:\\
\;\;\;\;t \cdot b\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{+302}:\\
\;\;\;\;x + a\\
\mathbf{else}:\\
\;\;\;\;t \cdot b\\
\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)) < -1.00000000000000006e290 or 2.0000000000000002e302 < (+.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 87.6%
Taylor expanded in a around 0
+-commutativeN/A
associate--l+N/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-evalN/A
sub-negN/A
+-commutativeN/A
distribute-rgt-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--.f6474.9
Applied rewrites74.9%
Taylor expanded in t around inf
Applied rewrites33.3%
if -1.00000000000000006e290 < (+.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)) < 2.0000000000000002e302Initial program 100.0%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval73.7
Applied rewrites73.7%
Taylor expanded in t around 0
Applied rewrites62.3%
Taylor expanded in b around 0
Applied rewrites45.8%
Final simplification40.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ y (+ t -2.0))) (t_2 (fma z (- 1.0 y) x)))
(if (<= b -2.05e+133)
(fma a (- 1.0 t) (fma b t_1 x))
(if (<= b 1.45e-43) (fma a (- 1.0 t) t_2) (fma b t_1 t_2)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y + (t + -2.0);
double t_2 = fma(z, (1.0 - y), x);
double tmp;
if (b <= -2.05e+133) {
tmp = fma(a, (1.0 - t), fma(b, t_1, x));
} else if (b <= 1.45e-43) {
tmp = fma(a, (1.0 - t), t_2);
} else {
tmp = fma(b, t_1, t_2);
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(y + Float64(t + -2.0)) t_2 = fma(z, Float64(1.0 - y), x) tmp = 0.0 if (b <= -2.05e+133) tmp = fma(a, Float64(1.0 - t), fma(b, t_1, x)); elseif (b <= 1.45e-43) tmp = fma(a, Float64(1.0 - t), t_2); else tmp = fma(b, t_1, t_2); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y + N[(t + -2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(z * N[(1.0 - y), $MachinePrecision] + x), $MachinePrecision]}, If[LessEqual[b, -2.05e+133], N[(a * N[(1.0 - t), $MachinePrecision] + N[(b * t$95$1 + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.45e-43], N[(a * N[(1.0 - t), $MachinePrecision] + t$95$2), $MachinePrecision], N[(b * t$95$1 + t$95$2), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y + \left(t + -2\right)\\
t_2 := \mathsf{fma}\left(z, 1 - y, x\right)\\
\mathbf{if}\;b \leq -2.05 \cdot 10^{+133}:\\
\;\;\;\;\mathsf{fma}\left(a, 1 - t, \mathsf{fma}\left(b, t\_1, x\right)\right)\\
\mathbf{elif}\;b \leq 1.45 \cdot 10^{-43}:\\
\;\;\;\;\mathsf{fma}\left(a, 1 - t, t\_2\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(b, t\_1, t\_2\right)\\
\end{array}
\end{array}
if b < -2.05000000000000002e133Initial program 87.5%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval97.4
Applied rewrites97.4%
if -2.05000000000000002e133 < b < 1.4500000000000001e-43Initial program 98.0%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites90.4%
if 1.4500000000000001e-43 < b Initial program 90.4%
Taylor expanded in a around 0
+-commutativeN/A
associate--l+N/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-evalN/A
sub-negN/A
+-commutativeN/A
distribute-rgt-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--.f6493.0
Applied rewrites93.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma a (- 1.0 t) (fma b (+ y (+ t -2.0)) x))))
(if (<= b -2.05e+133)
t_1
(if (<= b 1.45e-9) (fma a (- 1.0 t) (fma z (- 1.0 y) x)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma(a, (1.0 - t), fma(b, (y + (t + -2.0)), x));
double tmp;
if (b <= -2.05e+133) {
tmp = t_1;
} else if (b <= 1.45e-9) {
tmp = fma(a, (1.0 - t), fma(z, (1.0 - y), x));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(a, Float64(1.0 - t), fma(b, Float64(y + Float64(t + -2.0)), x)) tmp = 0.0 if (b <= -2.05e+133) tmp = t_1; elseif (b <= 1.45e-9) tmp = fma(a, Float64(1.0 - t), fma(z, Float64(1.0 - y), x)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(a * N[(1.0 - t), $MachinePrecision] + N[(b * N[(y + N[(t + -2.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -2.05e+133], t$95$1, If[LessEqual[b, 1.45e-9], N[(a * N[(1.0 - t), $MachinePrecision] + N[(z * N[(1.0 - y), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(a, 1 - t, \mathsf{fma}\left(b, y + \left(t + -2\right), x\right)\right)\\
\mathbf{if}\;b \leq -2.05 \cdot 10^{+133}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.45 \cdot 10^{-9}:\\
\;\;\;\;\mathsf{fma}\left(a, 1 - t, \mathsf{fma}\left(z, 1 - y, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -2.05000000000000002e133 or 1.44999999999999996e-9 < b Initial program 89.1%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval88.7
Applied rewrites88.7%
if -2.05000000000000002e133 < b < 1.44999999999999996e-9Initial program 98.1%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites90.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ t (+ y -2.0))))
(if (<= b -2.05e+133)
(fma b t_1 x)
(if (<= b 2e+44)
(fma a (- 1.0 t) (fma z (- 1.0 y) x))
(fma a (- 1.0 t) (* b t_1))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = t + (y + -2.0);
double tmp;
if (b <= -2.05e+133) {
tmp = fma(b, t_1, x);
} else if (b <= 2e+44) {
tmp = fma(a, (1.0 - t), fma(z, (1.0 - y), x));
} else {
tmp = fma(a, (1.0 - t), (b * t_1));
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(t + Float64(y + -2.0)) tmp = 0.0 if (b <= -2.05e+133) tmp = fma(b, t_1, x); elseif (b <= 2e+44) tmp = fma(a, Float64(1.0 - t), fma(z, Float64(1.0 - y), x)); else tmp = fma(a, Float64(1.0 - t), Float64(b * t_1)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(t + N[(y + -2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -2.05e+133], N[(b * t$95$1 + x), $MachinePrecision], If[LessEqual[b, 2e+44], N[(a * N[(1.0 - t), $MachinePrecision] + N[(z * N[(1.0 - y), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], N[(a * N[(1.0 - t), $MachinePrecision] + N[(b * t$95$1), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t + \left(y + -2\right)\\
\mathbf{if}\;b \leq -2.05 \cdot 10^{+133}:\\
\;\;\;\;\mathsf{fma}\left(b, t\_1, x\right)\\
\mathbf{elif}\;b \leq 2 \cdot 10^{+44}:\\
\;\;\;\;\mathsf{fma}\left(a, 1 - t, \mathsf{fma}\left(z, 1 - y, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(a, 1 - t, b \cdot t\_1\right)\\
\end{array}
\end{array}
if b < -2.05000000000000002e133Initial program 87.5%
Taylor expanded in a around 0
+-commutativeN/A
associate--l+N/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-evalN/A
sub-negN/A
+-commutativeN/A
distribute-rgt-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--.f6487.8
Applied rewrites87.8%
Taylor expanded in z around 0
Applied rewrites93.6%
if -2.05000000000000002e133 < b < 2.0000000000000002e44Initial program 98.2%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites88.4%
if 2.0000000000000002e44 < b Initial program 87.5%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval87.8
Applied rewrites87.8%
Taylor expanded in b around inf
Applied rewrites86.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* t (- b a))))
(if (<= t -4.5e+79)
t_1
(if (<= t 4.5e-212)
(fma b (+ y -2.0) (fma z (- y) x))
(if (<= t 5e+92) (fma z (- 1.0 y) (+ x a)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = t * (b - a);
double tmp;
if (t <= -4.5e+79) {
tmp = t_1;
} else if (t <= 4.5e-212) {
tmp = fma(b, (y + -2.0), fma(z, -y, x));
} else if (t <= 5e+92) {
tmp = fma(z, (1.0 - y), (x + a));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(t * Float64(b - a)) tmp = 0.0 if (t <= -4.5e+79) tmp = t_1; elseif (t <= 4.5e-212) tmp = fma(b, Float64(y + -2.0), fma(z, Float64(-y), x)); elseif (t <= 5e+92) tmp = fma(z, Float64(1.0 - y), Float64(x + a)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -4.5e+79], t$95$1, If[LessEqual[t, 4.5e-212], N[(b * N[(y + -2.0), $MachinePrecision] + N[(z * (-y) + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 5e+92], N[(z * N[(1.0 - y), $MachinePrecision] + N[(x + a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -4.5 \cdot 10^{+79}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 4.5 \cdot 10^{-212}:\\
\;\;\;\;\mathsf{fma}\left(b, y + -2, \mathsf{fma}\left(z, -y, x\right)\right)\\
\mathbf{elif}\;t \leq 5 \cdot 10^{+92}:\\
\;\;\;\;\mathsf{fma}\left(z, 1 - y, x + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -4.49999999999999994e79 or 5.00000000000000022e92 < t Initial program 85.4%
Taylor expanded in t around inf
lower-*.f64N/A
lower--.f6482.5
Applied rewrites82.5%
if -4.49999999999999994e79 < t < 4.4999999999999999e-212Initial program 98.9%
Taylor expanded in a around 0
+-commutativeN/A
associate--l+N/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-evalN/A
sub-negN/A
+-commutativeN/A
distribute-rgt-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--.f6482.4
Applied rewrites82.4%
Taylor expanded in t around 0
Applied rewrites81.5%
Taylor expanded in y around inf
Applied rewrites73.4%
if 4.4999999999999999e-212 < t < 5.00000000000000022e92Initial program 100.0%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites80.4%
Taylor expanded in t around 0
Applied rewrites78.7%
Final simplification78.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* t (- b a))))
(if (<= t -4.5e+79)
t_1
(if (<= t -5.8e-295)
(+ a (fma b (+ y -2.0) x))
(if (<= t 5e+92) (fma z (- 1.0 y) (+ x a)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = t * (b - a);
double tmp;
if (t <= -4.5e+79) {
tmp = t_1;
} else if (t <= -5.8e-295) {
tmp = a + fma(b, (y + -2.0), x);
} else if (t <= 5e+92) {
tmp = fma(z, (1.0 - y), (x + a));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(t * Float64(b - a)) tmp = 0.0 if (t <= -4.5e+79) tmp = t_1; elseif (t <= -5.8e-295) tmp = Float64(a + fma(b, Float64(y + -2.0), x)); elseif (t <= 5e+92) tmp = fma(z, Float64(1.0 - y), Float64(x + a)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -4.5e+79], t$95$1, If[LessEqual[t, -5.8e-295], N[(a + N[(b * N[(y + -2.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 5e+92], N[(z * N[(1.0 - y), $MachinePrecision] + N[(x + a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -4.5 \cdot 10^{+79}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -5.8 \cdot 10^{-295}:\\
\;\;\;\;a + \mathsf{fma}\left(b, y + -2, x\right)\\
\mathbf{elif}\;t \leq 5 \cdot 10^{+92}:\\
\;\;\;\;\mathsf{fma}\left(z, 1 - y, x + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -4.49999999999999994e79 or 5.00000000000000022e92 < t Initial program 85.4%
Taylor expanded in t around inf
lower-*.f64N/A
lower--.f6482.5
Applied rewrites82.5%
if -4.49999999999999994e79 < t < -5.8000000000000003e-295Initial program 98.5%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval76.6
Applied rewrites76.6%
Taylor expanded in t around 0
Applied rewrites72.3%
if -5.8000000000000003e-295 < t < 5.00000000000000022e92Initial program 100.0%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites77.0%
Taylor expanded in t around 0
Applied rewrites75.8%
Final simplification77.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma b (+ t (+ y -2.0)) x)))
(if (<= b -2.05e+133)
t_1
(if (<= b 1.18e+27) (fma a (- 1.0 t) (fma z (- 1.0 y) x)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma(b, (t + (y + -2.0)), x);
double tmp;
if (b <= -2.05e+133) {
tmp = t_1;
} else if (b <= 1.18e+27) {
tmp = fma(a, (1.0 - t), fma(z, (1.0 - y), x));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(b, Float64(t + Float64(y + -2.0)), x) tmp = 0.0 if (b <= -2.05e+133) tmp = t_1; elseif (b <= 1.18e+27) tmp = fma(a, Float64(1.0 - t), fma(z, Float64(1.0 - y), x)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(b * N[(t + N[(y + -2.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]}, If[LessEqual[b, -2.05e+133], t$95$1, If[LessEqual[b, 1.18e+27], N[(a * N[(1.0 - t), $MachinePrecision] + N[(z * N[(1.0 - y), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(b, t + \left(y + -2\right), x\right)\\
\mathbf{if}\;b \leq -2.05 \cdot 10^{+133}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.18 \cdot 10^{+27}:\\
\;\;\;\;\mathsf{fma}\left(a, 1 - t, \mathsf{fma}\left(z, 1 - y, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -2.05000000000000002e133 or 1.18000000000000006e27 < b Initial program 87.9%
Taylor expanded in a around 0
+-commutativeN/A
associate--l+N/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-evalN/A
sub-negN/A
+-commutativeN/A
distribute-rgt-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--.f6492.5
Applied rewrites92.5%
Taylor expanded in z around 0
Applied rewrites86.5%
if -2.05000000000000002e133 < b < 1.18000000000000006e27Initial program 98.2%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites88.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* b (+ y (+ t -2.0)))))
(if (<= b -8e+78)
t_1
(if (<= b 3.5e-305)
(fma a (- 1.0 t) x)
(if (<= b 2e+44) (fma z (- 1.0 y) x) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = b * (y + (t + -2.0));
double tmp;
if (b <= -8e+78) {
tmp = t_1;
} else if (b <= 3.5e-305) {
tmp = fma(a, (1.0 - t), x);
} else if (b <= 2e+44) {
tmp = fma(z, (1.0 - y), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(b * Float64(y + Float64(t + -2.0))) tmp = 0.0 if (b <= -8e+78) tmp = t_1; elseif (b <= 3.5e-305) tmp = fma(a, Float64(1.0 - t), x); elseif (b <= 2e+44) tmp = fma(z, Float64(1.0 - y), x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(b * N[(y + N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -8e+78], t$95$1, If[LessEqual[b, 3.5e-305], N[(a * N[(1.0 - t), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[b, 2e+44], N[(z * N[(1.0 - y), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := b \cdot \left(y + \left(t + -2\right)\right)\\
\mathbf{if}\;b \leq -8 \cdot 10^{+78}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 3.5 \cdot 10^{-305}:\\
\;\;\;\;\mathsf{fma}\left(a, 1 - t, x\right)\\
\mathbf{elif}\;b \leq 2 \cdot 10^{+44}:\\
\;\;\;\;\mathsf{fma}\left(z, 1 - y, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -8.00000000000000007e78 or 2.0000000000000002e44 < b Initial program 87.2%
Taylor expanded in b around inf
lower-*.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval80.8
Applied rewrites80.8%
if -8.00000000000000007e78 < b < 3.4999999999999998e-305Initial program 97.6%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval70.5
Applied rewrites70.5%
Taylor expanded in b around 0
Applied rewrites61.8%
if 3.4999999999999998e-305 < b < 2.0000000000000002e44Initial program 100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate--l+N/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-evalN/A
sub-negN/A
+-commutativeN/A
distribute-rgt-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--.f6472.2
Applied rewrites72.2%
Taylor expanded in b around 0
Applied rewrites61.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* t (- b a))))
(if (<= t -4.5e+79)
t_1
(if (<= t 6.5e-180)
(+ a (fma b (+ y -2.0) x))
(if (<= t 5e+92) (fma z (- 1.0 y) x) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = t * (b - a);
double tmp;
if (t <= -4.5e+79) {
tmp = t_1;
} else if (t <= 6.5e-180) {
tmp = a + fma(b, (y + -2.0), x);
} else if (t <= 5e+92) {
tmp = fma(z, (1.0 - y), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(t * Float64(b - a)) tmp = 0.0 if (t <= -4.5e+79) tmp = t_1; elseif (t <= 6.5e-180) tmp = Float64(a + fma(b, Float64(y + -2.0), x)); elseif (t <= 5e+92) tmp = fma(z, Float64(1.0 - y), x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -4.5e+79], t$95$1, If[LessEqual[t, 6.5e-180], N[(a + N[(b * N[(y + -2.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 5e+92], N[(z * N[(1.0 - y), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -4.5 \cdot 10^{+79}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 6.5 \cdot 10^{-180}:\\
\;\;\;\;a + \mathsf{fma}\left(b, y + -2, x\right)\\
\mathbf{elif}\;t \leq 5 \cdot 10^{+92}:\\
\;\;\;\;\mathsf{fma}\left(z, 1 - y, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -4.49999999999999994e79 or 5.00000000000000022e92 < t Initial program 85.4%
Taylor expanded in t around inf
lower-*.f64N/A
lower--.f6482.5
Applied rewrites82.5%
if -4.49999999999999994e79 < t < 6.50000000000000013e-180Initial program 99.0%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval73.8
Applied rewrites73.8%
Taylor expanded in t around 0
Applied rewrites71.1%
if 6.50000000000000013e-180 < t < 5.00000000000000022e92Initial program 100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate--l+N/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-evalN/A
sub-negN/A
+-commutativeN/A
distribute-rgt-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--.f6484.6
Applied rewrites84.6%
Taylor expanded in b around 0
Applied rewrites65.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* t (- b a))))
(if (<= t -720000000000.0)
t_1
(if (<= t 4.8e-215)
(+ a (fma b -2.0 x))
(if (<= t 1.7e+92) (+ z (+ x a)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = t * (b - a);
double tmp;
if (t <= -720000000000.0) {
tmp = t_1;
} else if (t <= 4.8e-215) {
tmp = a + fma(b, -2.0, x);
} else if (t <= 1.7e+92) {
tmp = z + (x + a);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(t * Float64(b - a)) tmp = 0.0 if (t <= -720000000000.0) tmp = t_1; elseif (t <= 4.8e-215) tmp = Float64(a + fma(b, -2.0, x)); elseif (t <= 1.7e+92) tmp = Float64(z + Float64(x + a)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -720000000000.0], t$95$1, If[LessEqual[t, 4.8e-215], N[(a + N[(b * -2.0 + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.7e+92], N[(z + N[(x + a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -720000000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 4.8 \cdot 10^{-215}:\\
\;\;\;\;a + \mathsf{fma}\left(b, -2, x\right)\\
\mathbf{elif}\;t \leq 1.7 \cdot 10^{+92}:\\
\;\;\;\;z + \left(x + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -7.2e11 or 1.6999999999999999e92 < t Initial program 86.9%
Taylor expanded in t around inf
lower-*.f64N/A
lower--.f6478.1
Applied rewrites78.1%
if -7.2e11 < t < 4.8000000000000002e-215Initial program 98.8%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval73.8
Applied rewrites73.8%
Taylor expanded in t around 0
Applied rewrites73.8%
Taylor expanded in y around 0
Applied rewrites52.6%
if 4.8000000000000002e-215 < t < 1.6999999999999999e92Initial program 100.0%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites79.3%
Taylor expanded in t around 0
Applied rewrites77.7%
Taylor expanded in y around 0
Applied rewrites53.2%
Final simplification62.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* t (- a))))
(if (<= t -1.28e+238)
t_1
(if (<= t -2.6e+81) (* t b) (if (<= t 5e+92) (+ z (+ x a)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = t * -a;
double tmp;
if (t <= -1.28e+238) {
tmp = t_1;
} else if (t <= -2.6e+81) {
tmp = t * b;
} else if (t <= 5e+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 = t * -a
if (t <= (-1.28d+238)) then
tmp = t_1
else if (t <= (-2.6d+81)) then
tmp = t * b
else if (t <= 5d+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 = t * -a;
double tmp;
if (t <= -1.28e+238) {
tmp = t_1;
} else if (t <= -2.6e+81) {
tmp = t * b;
} else if (t <= 5e+92) {
tmp = z + (x + a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = t * -a tmp = 0 if t <= -1.28e+238: tmp = t_1 elif t <= -2.6e+81: tmp = t * b elif t <= 5e+92: tmp = z + (x + a) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(t * Float64(-a)) tmp = 0.0 if (t <= -1.28e+238) tmp = t_1; elseif (t <= -2.6e+81) tmp = Float64(t * b); elseif (t <= 5e+92) tmp = Float64(z + Float64(x + a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = t * -a; tmp = 0.0; if (t <= -1.28e+238) tmp = t_1; elseif (t <= -2.6e+81) tmp = t * b; elseif (t <= 5e+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[(t * (-a)), $MachinePrecision]}, If[LessEqual[t, -1.28e+238], t$95$1, If[LessEqual[t, -2.6e+81], N[(t * b), $MachinePrecision], If[LessEqual[t, 5e+92], N[(z + N[(x + a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(-a\right)\\
\mathbf{if}\;t \leq -1.28 \cdot 10^{+238}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -2.6 \cdot 10^{+81}:\\
\;\;\;\;t \cdot b\\
\mathbf{elif}\;t \leq 5 \cdot 10^{+92}:\\
\;\;\;\;z + \left(x + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -1.28000000000000007e238 or 5.00000000000000022e92 < t Initial program 86.0%
Taylor expanded in t around inf
lower-*.f64N/A
lower--.f6487.9
Applied rewrites87.9%
Taylor expanded in b around 0
Applied rewrites60.6%
if -1.28000000000000007e238 < t < -2.59999999999999992e81Initial program 83.9%
Taylor expanded in a around 0
+-commutativeN/A
associate--l+N/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-evalN/A
sub-negN/A
+-commutativeN/A
distribute-rgt-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--.f6478.0
Applied rewrites78.0%
Taylor expanded in t around inf
Applied rewrites53.2%
if -2.59999999999999992e81 < t < 5.00000000000000022e92Initial program 99.4%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites73.4%
Taylor expanded in t around 0
Applied rewrites71.0%
Taylor expanded in y around 0
Applied rewrites49.1%
Final simplification52.2%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* t (- b a)))) (if (<= t -1.26e+76) t_1 (if (<= t 5e+92) (fma z (- 1.0 y) x) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = t * (b - a);
double tmp;
if (t <= -1.26e+76) {
tmp = t_1;
} else if (t <= 5e+92) {
tmp = fma(z, (1.0 - y), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(t * Float64(b - a)) tmp = 0.0 if (t <= -1.26e+76) tmp = t_1; elseif (t <= 5e+92) tmp = fma(z, Float64(1.0 - y), x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.26e+76], t$95$1, If[LessEqual[t, 5e+92], N[(z * N[(1.0 - y), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -1.26 \cdot 10^{+76}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 5 \cdot 10^{+92}:\\
\;\;\;\;\mathsf{fma}\left(z, 1 - y, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -1.26000000000000007e76 or 5.00000000000000022e92 < t Initial program 85.7%
Taylor expanded in t around inf
lower-*.f64N/A
lower--.f6481.8
Applied rewrites81.8%
if -1.26000000000000007e76 < t < 5.00000000000000022e92Initial program 99.4%
Taylor expanded in a around 0
+-commutativeN/A
associate--l+N/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-evalN/A
sub-negN/A
+-commutativeN/A
distribute-rgt-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--.f6481.6
Applied rewrites81.6%
Taylor expanded in b around 0
Applied rewrites55.8%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* y (- b z)))) (if (<= y -2.05e+80) t_1 (if (<= y 2.5e+125) (fma a (- 1.0 t) x) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y * (b - z);
double tmp;
if (y <= -2.05e+80) {
tmp = t_1;
} else if (y <= 2.5e+125) {
tmp = fma(a, (1.0 - t), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(y * Float64(b - z)) tmp = 0.0 if (y <= -2.05e+80) tmp = t_1; elseif (y <= 2.5e+125) tmp = fma(a, Float64(1.0 - t), x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.05e+80], t$95$1, If[LessEqual[y, 2.5e+125], N[(a * N[(1.0 - t), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(b - z\right)\\
\mathbf{if}\;y \leq -2.05 \cdot 10^{+80}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2.5 \cdot 10^{+125}:\\
\;\;\;\;\mathsf{fma}\left(a, 1 - t, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.05000000000000001e80 or 2.49999999999999981e125 < y Initial program 92.0%
Taylor expanded in y around inf
lower-*.f64N/A
lower--.f6479.1
Applied rewrites79.1%
if -2.05000000000000001e80 < y < 2.49999999999999981e125Initial program 95.8%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval81.4
Applied rewrites81.4%
Taylor expanded in b around 0
Applied rewrites55.4%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* t (- b a)))) (if (<= t -8.8e+77) t_1 (if (<= t 1.7e+92) (+ z (+ x a)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = t * (b - a);
double tmp;
if (t <= -8.8e+77) {
tmp = t_1;
} else if (t <= 1.7e+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 = t * (b - a)
if (t <= (-8.8d+77)) then
tmp = t_1
else if (t <= 1.7d+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 = t * (b - a);
double tmp;
if (t <= -8.8e+77) {
tmp = t_1;
} else if (t <= 1.7e+92) {
tmp = z + (x + a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = t * (b - a) tmp = 0 if t <= -8.8e+77: tmp = t_1 elif t <= 1.7e+92: tmp = z + (x + a) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(t * Float64(b - a)) tmp = 0.0 if (t <= -8.8e+77) tmp = t_1; elseif (t <= 1.7e+92) tmp = Float64(z + Float64(x + a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = t * (b - a); tmp = 0.0; if (t <= -8.8e+77) tmp = t_1; elseif (t <= 1.7e+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[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -8.8e+77], t$95$1, If[LessEqual[t, 1.7e+92], N[(z + N[(x + a), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -8.8 \cdot 10^{+77}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 1.7 \cdot 10^{+92}:\\
\;\;\;\;z + \left(x + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -8.8000000000000002e77 or 1.6999999999999999e92 < t Initial program 85.4%
Taylor expanded in t around inf
lower-*.f64N/A
lower--.f6482.5
Applied rewrites82.5%
if -8.8000000000000002e77 < t < 1.6999999999999999e92Initial program 99.4%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites73.3%
Taylor expanded in t around 0
Applied rewrites71.4%
Taylor expanded in y around 0
Applied rewrites49.4%
Final simplification60.9%
(FPCore (x y z t a b) :precision binary64 (if (<= b -2.4e+133) (* b (+ y -2.0)) (if (<= b 2.25e+45) (+ z (+ x a)) (* b (+ t -2.0)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -2.4e+133) {
tmp = b * (y + -2.0);
} else if (b <= 2.25e+45) {
tmp = z + (x + a);
} else {
tmp = b * (t + -2.0);
}
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 <= (-2.4d+133)) then
tmp = b * (y + (-2.0d0))
else if (b <= 2.25d+45) then
tmp = z + (x + a)
else
tmp = b * (t + (-2.0d0))
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 <= -2.4e+133) {
tmp = b * (y + -2.0);
} else if (b <= 2.25e+45) {
tmp = z + (x + a);
} else {
tmp = b * (t + -2.0);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -2.4e+133: tmp = b * (y + -2.0) elif b <= 2.25e+45: tmp = z + (x + a) else: tmp = b * (t + -2.0) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -2.4e+133) tmp = Float64(b * Float64(y + -2.0)); elseif (b <= 2.25e+45) tmp = Float64(z + Float64(x + a)); else tmp = Float64(b * Float64(t + -2.0)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -2.4e+133) tmp = b * (y + -2.0); elseif (b <= 2.25e+45) tmp = z + (x + a); else tmp = b * (t + -2.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -2.4e+133], N[(b * N[(y + -2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.25e+45], N[(z + N[(x + a), $MachinePrecision]), $MachinePrecision], N[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -2.4 \cdot 10^{+133}:\\
\;\;\;\;b \cdot \left(y + -2\right)\\
\mathbf{elif}\;b \leq 2.25 \cdot 10^{+45}:\\
\;\;\;\;z + \left(x + a\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(t + -2\right)\\
\end{array}
\end{array}
if b < -2.3999999999999999e133Initial program 87.5%
Taylor expanded in b around inf
lower-*.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval85.5
Applied rewrites85.5%
Taylor expanded in t around 0
Applied rewrites65.4%
if -2.3999999999999999e133 < b < 2.2499999999999999e45Initial program 98.2%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites88.4%
Taylor expanded in t around 0
Applied rewrites68.8%
Taylor expanded in y around 0
Applied rewrites48.0%
if 2.2499999999999999e45 < b Initial program 87.5%
Taylor expanded in b around inf
lower-*.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval83.4
Applied rewrites83.4%
Taylor expanded in y around 0
Applied rewrites60.8%
Final simplification52.8%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* b (+ t -2.0)))) (if (<= b -8.4e+133) t_1 (if (<= b 2.25e+45) (+ z (+ x a)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = b * (t + -2.0);
double tmp;
if (b <= -8.4e+133) {
tmp = t_1;
} else if (b <= 2.25e+45) {
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 * (t + (-2.0d0))
if (b <= (-8.4d+133)) then
tmp = t_1
else if (b <= 2.25d+45) 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 * (t + -2.0);
double tmp;
if (b <= -8.4e+133) {
tmp = t_1;
} else if (b <= 2.25e+45) {
tmp = z + (x + a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = b * (t + -2.0) tmp = 0 if b <= -8.4e+133: tmp = t_1 elif b <= 2.25e+45: tmp = z + (x + a) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(b * Float64(t + -2.0)) tmp = 0.0 if (b <= -8.4e+133) tmp = t_1; elseif (b <= 2.25e+45) tmp = Float64(z + Float64(x + a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = b * (t + -2.0); tmp = 0.0; if (b <= -8.4e+133) tmp = t_1; elseif (b <= 2.25e+45) 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[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -8.4e+133], t$95$1, If[LessEqual[b, 2.25e+45], N[(z + N[(x + a), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := b \cdot \left(t + -2\right)\\
\mathbf{if}\;b \leq -8.4 \cdot 10^{+133}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 2.25 \cdot 10^{+45}:\\
\;\;\;\;z + \left(x + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -8.4e133 or 2.2499999999999999e45 < b Initial program 87.5%
Taylor expanded in b around inf
lower-*.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval84.0
Applied rewrites84.0%
Taylor expanded in y around 0
Applied rewrites58.7%
if -8.4e133 < b < 2.2499999999999999e45Initial program 98.2%
Taylor expanded in b around 0
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
distribute-rgt-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
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
Applied rewrites88.4%
Taylor expanded in t around 0
Applied rewrites68.8%
Taylor expanded in y around 0
Applied rewrites48.0%
Final simplification51.7%
(FPCore (x y z t a b) :precision binary64 (if (<= y -5e+80) (* y b) (if (<= y 2.35e+229) (+ x a) (* y b))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -5e+80) {
tmp = y * b;
} else if (y <= 2.35e+229) {
tmp = x + a;
} else {
tmp = y * 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 (y <= (-5d+80)) then
tmp = y * b
else if (y <= 2.35d+229) then
tmp = x + a
else
tmp = y * 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 (y <= -5e+80) {
tmp = y * b;
} else if (y <= 2.35e+229) {
tmp = x + a;
} else {
tmp = y * b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -5e+80: tmp = y * b elif y <= 2.35e+229: tmp = x + a else: tmp = y * b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -5e+80) tmp = Float64(y * b); elseif (y <= 2.35e+229) tmp = Float64(x + a); else tmp = Float64(y * b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -5e+80) tmp = y * b; elseif (y <= 2.35e+229) tmp = x + a; else tmp = y * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -5e+80], N[(y * b), $MachinePrecision], If[LessEqual[y, 2.35e+229], N[(x + a), $MachinePrecision], N[(y * b), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5 \cdot 10^{+80}:\\
\;\;\;\;y \cdot b\\
\mathbf{elif}\;y \leq 2.35 \cdot 10^{+229}:\\
\;\;\;\;x + a\\
\mathbf{else}:\\
\;\;\;\;y \cdot b\\
\end{array}
\end{array}
if y < -4.99999999999999961e80 or 2.35e229 < y Initial program 90.1%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval67.3
Applied rewrites67.3%
Taylor expanded in y around inf
Applied rewrites49.2%
if -4.99999999999999961e80 < y < 2.35e229Initial program 96.2%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval78.5
Applied rewrites78.5%
Taylor expanded in t around 0
Applied rewrites45.4%
Taylor expanded in b around 0
Applied rewrites34.2%
Final simplification38.3%
(FPCore (x y z t a b) :precision binary64 (+ x a))
double code(double x, double y, double z, double t, double a, double b) {
return x + a;
}
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 + a
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x + a;
}
def code(x, y, z, t, a, b): return x + a
function code(x, y, z, t, a, b) return Float64(x + a) end
function tmp = code(x, y, z, t, a, b) tmp = x + a; end
code[x_, y_, z_, t_, a_, b_] := N[(x + a), $MachinePrecision]
\begin{array}{l}
\\
x + a
\end{array}
Initial program 94.5%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-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
+-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
associate-+r-N/A
lower-+.f64N/A
sub-negN/A
lower-+.f64N/A
metadata-eval75.4
Applied rewrites75.4%
Taylor expanded in t around 0
Applied rewrites48.5%
Taylor expanded in b around 0
Applied rewrites27.3%
Final simplification27.3%
herbie shell --seed 2024221
(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)))