
(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 25 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
(+ (+ (* a (- 1.0 t)) (+ x (* z (- 1.0 y)))) (* (- (+ y t) 2.0) b))))
(if (<= t_1 INFINITY) t_1 (* y (- b z)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((a * (1.0 - t)) + (x + (z * (1.0 - y)))) + (((y + t) - 2.0) * b);
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = y * (b - z);
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((a * (1.0 - t)) + (x + (z * (1.0 - y)))) + (((y + t) - 2.0) * b);
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = y * (b - z);
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = ((a * (1.0 - t)) + (x + (z * (1.0 - y)))) + (((y + t) - 2.0) * b) tmp = 0 if t_1 <= math.inf: tmp = t_1 else: tmp = y * (b - z) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(a * Float64(1.0 - t)) + Float64(x + Float64(z * Float64(1.0 - y)))) + Float64(Float64(Float64(y + t) - 2.0) * b)) tmp = 0.0 if (t_1 <= Inf) tmp = t_1; else tmp = Float64(y * Float64(b - z)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((a * (1.0 - t)) + (x + (z * (1.0 - y)))) + (((y + t) - 2.0) * b); tmp = 0.0; if (t_1 <= Inf) tmp = t_1; else tmp = y * (b - z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision] + N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $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[(y * N[(b - z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(a \cdot \left(1 - t\right) + \left(x + z \cdot \left(1 - y\right)\right)\right) + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{if}\;t\_1 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;y \cdot \left(b - z\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 y around inf
*-lowering-*.f64N/A
--lowering--.f6464.3%
Simplified64.3%
Final simplification98.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ z (+ (* y (- b z)) (* b (+ t -2.0))))))
(if (<= b -5.8e+150)
t_1
(if (<= b 1.1e+71)
(+ x (+ (* a (- 1.0 t)) (* z (- 1.0 y))))
(if (<= b 1.6e+154) t_1 (+ x (* (- (+ y t) 2.0) b)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z + ((y * (b - z)) + (b * (t + -2.0)));
double tmp;
if (b <= -5.8e+150) {
tmp = t_1;
} else if (b <= 1.1e+71) {
tmp = x + ((a * (1.0 - t)) + (z * (1.0 - y)));
} else if (b <= 1.6e+154) {
tmp = t_1;
} else {
tmp = x + (((y + t) - 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) :: t_1
real(8) :: tmp
t_1 = z + ((y * (b - z)) + (b * (t + (-2.0d0))))
if (b <= (-5.8d+150)) then
tmp = t_1
else if (b <= 1.1d+71) then
tmp = x + ((a * (1.0d0 - t)) + (z * (1.0d0 - y)))
else if (b <= 1.6d+154) then
tmp = t_1
else
tmp = x + (((y + t) - 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 t_1 = z + ((y * (b - z)) + (b * (t + -2.0)));
double tmp;
if (b <= -5.8e+150) {
tmp = t_1;
} else if (b <= 1.1e+71) {
tmp = x + ((a * (1.0 - t)) + (z * (1.0 - y)));
} else if (b <= 1.6e+154) {
tmp = t_1;
} else {
tmp = x + (((y + t) - 2.0) * b);
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = z + ((y * (b - z)) + (b * (t + -2.0))) tmp = 0 if b <= -5.8e+150: tmp = t_1 elif b <= 1.1e+71: tmp = x + ((a * (1.0 - t)) + (z * (1.0 - y))) elif b <= 1.6e+154: tmp = t_1 else: tmp = x + (((y + t) - 2.0) * b) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(z + Float64(Float64(y * Float64(b - z)) + Float64(b * Float64(t + -2.0)))) tmp = 0.0 if (b <= -5.8e+150) tmp = t_1; elseif (b <= 1.1e+71) tmp = Float64(x + Float64(Float64(a * Float64(1.0 - t)) + Float64(z * Float64(1.0 - y)))); elseif (b <= 1.6e+154) tmp = t_1; else tmp = Float64(x + Float64(Float64(Float64(y + t) - 2.0) * b)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = z + ((y * (b - z)) + (b * (t + -2.0))); tmp = 0.0; if (b <= -5.8e+150) tmp = t_1; elseif (b <= 1.1e+71) tmp = x + ((a * (1.0 - t)) + (z * (1.0 - y))); elseif (b <= 1.6e+154) tmp = t_1; else tmp = x + (((y + t) - 2.0) * b); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(z + N[(N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision] + N[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -5.8e+150], t$95$1, If[LessEqual[b, 1.1e+71], N[(x + N[(N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision] + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.6e+154], t$95$1, N[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z + \left(y \cdot \left(b - z\right) + b \cdot \left(t + -2\right)\right)\\
\mathbf{if}\;b \leq -5.8 \cdot 10^{+150}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.1 \cdot 10^{+71}:\\
\;\;\;\;x + \left(a \cdot \left(1 - t\right) + z \cdot \left(1 - y\right)\right)\\
\mathbf{elif}\;b \leq 1.6 \cdot 10^{+154}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x + \left(\left(y + t\right) - 2\right) \cdot b\\
\end{array}
\end{array}
if b < -5.80000000000000022e150 or 1.09999999999999997e71 < b < 1.6e154Initial program 92.7%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f6485.9%
Simplified85.9%
Taylor expanded in y around 0
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6489.6%
Simplified89.6%
if -5.80000000000000022e150 < b < 1.09999999999999997e71Initial program 97.5%
Taylor expanded in b around 0
sub-negN/A
+-lowering-+.f64N/A
+-commutativeN/A
distribute-neg-inN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
Simplified79.6%
if 1.6e154 < b Initial program 91.9%
Taylor expanded in x around inf
Simplified99.2%
Final simplification84.6%
(FPCore (x y z t a b) :precision binary64 (if (<= t 1e+255) (+ (+ x (* b (+ t -2.0))) (+ (* y (- b z)) (+ z (* a (- 1.0 t))))) (* t (- b a))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= 1e+255) {
tmp = (x + (b * (t + -2.0))) + ((y * (b - z)) + (z + (a * (1.0 - t))));
} else {
tmp = t * (b - a);
}
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 (t <= 1d+255) then
tmp = (x + (b * (t + (-2.0d0)))) + ((y * (b - z)) + (z + (a * (1.0d0 - t))))
else
tmp = t * (b - a)
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 (t <= 1e+255) {
tmp = (x + (b * (t + -2.0))) + ((y * (b - z)) + (z + (a * (1.0 - t))));
} else {
tmp = t * (b - a);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= 1e+255: tmp = (x + (b * (t + -2.0))) + ((y * (b - z)) + (z + (a * (1.0 - t)))) else: tmp = t * (b - a) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= 1e+255) tmp = Float64(Float64(x + Float64(b * Float64(t + -2.0))) + Float64(Float64(y * Float64(b - z)) + Float64(z + Float64(a * Float64(1.0 - t))))); else tmp = Float64(t * Float64(b - a)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= 1e+255) tmp = (x + (b * (t + -2.0))) + ((y * (b - z)) + (z + (a * (1.0 - t)))); else tmp = t * (b - a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, 1e+255], N[(N[(x + N[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision] + N[(z + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 10^{+255}:\\
\;\;\;\;\left(x + b \cdot \left(t + -2\right)\right) + \left(y \cdot \left(b - z\right) + \left(z + a \cdot \left(1 - t\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t \cdot \left(b - a\right)\\
\end{array}
\end{array}
if t < 9.99999999999999988e254Initial program 97.1%
Taylor expanded in y around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
Simplified97.5%
if 9.99999999999999988e254 < t Initial program 69.2%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6492.7%
Simplified92.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ (* (- (+ y t) 2.0) b) (+ x (* z (- 1.0 y))))))
(if (<= y -13200000.0)
t_1
(if (<= y 8.2e+26) (+ (+ (* t (- b a)) (+ z (* b -2.0))) (+ x a)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (((y + t) - 2.0) * b) + (x + (z * (1.0 - y)));
double tmp;
if (y <= -13200000.0) {
tmp = t_1;
} else if (y <= 8.2e+26) {
tmp = ((t * (b - a)) + (z + (b * -2.0))) + (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 + t) - 2.0d0) * b) + (x + (z * (1.0d0 - y)))
if (y <= (-13200000.0d0)) then
tmp = t_1
else if (y <= 8.2d+26) then
tmp = ((t * (b - a)) + (z + (b * (-2.0d0)))) + (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 + t) - 2.0) * b) + (x + (z * (1.0 - y)));
double tmp;
if (y <= -13200000.0) {
tmp = t_1;
} else if (y <= 8.2e+26) {
tmp = ((t * (b - a)) + (z + (b * -2.0))) + (x + a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (((y + t) - 2.0) * b) + (x + (z * (1.0 - y))) tmp = 0 if y <= -13200000.0: tmp = t_1 elif y <= 8.2e+26: tmp = ((t * (b - a)) + (z + (b * -2.0))) + (x + a) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(Float64(y + t) - 2.0) * b) + Float64(x + Float64(z * Float64(1.0 - y)))) tmp = 0.0 if (y <= -13200000.0) tmp = t_1; elseif (y <= 8.2e+26) tmp = Float64(Float64(Float64(t * Float64(b - a)) + Float64(z + Float64(b * -2.0))) + Float64(x + a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (((y + t) - 2.0) * b) + (x + (z * (1.0 - y))); tmp = 0.0; if (y <= -13200000.0) tmp = t_1; elseif (y <= 8.2e+26) tmp = ((t * (b - a)) + (z + (b * -2.0))) + (x + a); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision] + N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -13200000.0], t$95$1, If[LessEqual[y, 8.2e+26], N[(N[(N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision] + N[(z + N[(b * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x + a), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(y + t\right) - 2\right) \cdot b + \left(x + z \cdot \left(1 - y\right)\right)\\
\mathbf{if}\;y \leq -13200000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 8.2 \cdot 10^{+26}:\\
\;\;\;\;\left(t \cdot \left(b - a\right) + \left(z + b \cdot -2\right)\right) + \left(x + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.32e7 or 8.19999999999999967e26 < y Initial program 92.8%
Taylor expanded in a around 0
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f6482.9%
Simplified82.9%
if -1.32e7 < y < 8.19999999999999967e26Initial program 98.4%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified100.0%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6498.9%
Simplified98.9%
Final simplification91.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- (+ y t) 2.0) b)) (t_2 (* a (- 1.0 t))))
(if (<= a -2.3e+121)
(+ t_2 (+ x (* b (+ y -2.0))))
(if (<= a 2.55e+42) (+ t_1 (+ x (* z (- 1.0 y)))) (+ 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) * b;
double t_2 = a * (1.0 - t);
double tmp;
if (a <= -2.3e+121) {
tmp = t_2 + (x + (b * (y + -2.0)));
} else if (a <= 2.55e+42) {
tmp = t_1 + (x + (z * (1.0 - y)));
} else {
tmp = t_1 + t_2;
}
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) :: t_2
real(8) :: tmp
t_1 = ((y + t) - 2.0d0) * b
t_2 = a * (1.0d0 - t)
if (a <= (-2.3d+121)) then
tmp = t_2 + (x + (b * (y + (-2.0d0))))
else if (a <= 2.55d+42) then
tmp = t_1 + (x + (z * (1.0d0 - y)))
else
tmp = t_1 + t_2
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 + t) - 2.0) * b;
double t_2 = a * (1.0 - t);
double tmp;
if (a <= -2.3e+121) {
tmp = t_2 + (x + (b * (y + -2.0)));
} else if (a <= 2.55e+42) {
tmp = t_1 + (x + (z * (1.0 - y)));
} else {
tmp = t_1 + t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = ((y + t) - 2.0) * b t_2 = a * (1.0 - t) tmp = 0 if a <= -2.3e+121: tmp = t_2 + (x + (b * (y + -2.0))) elif a <= 2.55e+42: tmp = t_1 + (x + (z * (1.0 - y))) else: tmp = t_1 + t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(y + t) - 2.0) * b) t_2 = Float64(a * Float64(1.0 - t)) tmp = 0.0 if (a <= -2.3e+121) tmp = Float64(t_2 + Float64(x + Float64(b * Float64(y + -2.0)))); elseif (a <= 2.55e+42) tmp = Float64(t_1 + Float64(x + Float64(z * Float64(1.0 - y)))); else tmp = Float64(t_1 + t_2); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((y + t) - 2.0) * b; t_2 = a * (1.0 - t); tmp = 0.0; if (a <= -2.3e+121) tmp = t_2 + (x + (b * (y + -2.0))); elseif (a <= 2.55e+42) tmp = t_1 + (x + (z * (1.0 - y))); else tmp = t_1 + t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]}, Block[{t$95$2 = N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -2.3e+121], N[(t$95$2 + N[(x + N[(b * N[(y + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.55e+42], N[(t$95$1 + N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 + t$95$2), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(y + t\right) - 2\right) \cdot b\\
t_2 := a \cdot \left(1 - t\right)\\
\mathbf{if}\;a \leq -2.3 \cdot 10^{+121}:\\
\;\;\;\;t\_2 + \left(x + b \cdot \left(y + -2\right)\right)\\
\mathbf{elif}\;a \leq 2.55 \cdot 10^{+42}:\\
\;\;\;\;t\_1 + \left(x + z \cdot \left(1 - y\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1 + t\_2\\
\end{array}
\end{array}
if a < -2.2999999999999999e121Initial program 86.8%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified92.1%
Taylor expanded in a around inf
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f6484.1%
Simplified84.1%
if -2.2999999999999999e121 < a < 2.55e42Initial program 98.8%
Taylor expanded in a around 0
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f6493.1%
Simplified93.1%
if 2.55e42 < a Initial program 91.7%
Taylor expanded in a around inf
sub-negN/A
neg-mul-1N/A
metadata-evalN/A
distribute-lft-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6481.2%
Simplified81.2%
Final simplification89.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* (- (+ y t) 2.0) b))))
(if (<= b -1e+139)
t_1
(if (<= b -1.7e-160)
(- (+ x (* y b)) (* (+ y -1.0) z))
(if (<= b 2.4e-24) (+ (+ x a) (- z (* t a))) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (((y + t) - 2.0) * b);
double tmp;
if (b <= -1e+139) {
tmp = t_1;
} else if (b <= -1.7e-160) {
tmp = (x + (y * b)) - ((y + -1.0) * z);
} else if (b <= 2.4e-24) {
tmp = (x + a) + (z - (t * 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 = x + (((y + t) - 2.0d0) * b)
if (b <= (-1d+139)) then
tmp = t_1
else if (b <= (-1.7d-160)) then
tmp = (x + (y * b)) - ((y + (-1.0d0)) * z)
else if (b <= 2.4d-24) then
tmp = (x + a) + (z - (t * 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 = x + (((y + t) - 2.0) * b);
double tmp;
if (b <= -1e+139) {
tmp = t_1;
} else if (b <= -1.7e-160) {
tmp = (x + (y * b)) - ((y + -1.0) * z);
} else if (b <= 2.4e-24) {
tmp = (x + a) + (z - (t * a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (((y + t) - 2.0) * b) tmp = 0 if b <= -1e+139: tmp = t_1 elif b <= -1.7e-160: tmp = (x + (y * b)) - ((y + -1.0) * z) elif b <= 2.4e-24: tmp = (x + a) + (z - (t * a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(Float64(Float64(y + t) - 2.0) * b)) tmp = 0.0 if (b <= -1e+139) tmp = t_1; elseif (b <= -1.7e-160) tmp = Float64(Float64(x + Float64(y * b)) - Float64(Float64(y + -1.0) * z)); elseif (b <= 2.4e-24) tmp = Float64(Float64(x + a) + Float64(z - Float64(t * a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (((y + t) - 2.0) * b); tmp = 0.0; if (b <= -1e+139) tmp = t_1; elseif (b <= -1.7e-160) tmp = (x + (y * b)) - ((y + -1.0) * z); elseif (b <= 2.4e-24) tmp = (x + a) + (z - (t * a)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1e+139], t$95$1, If[LessEqual[b, -1.7e-160], N[(N[(x + N[(y * b), $MachinePrecision]), $MachinePrecision] - N[(N[(y + -1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.4e-24], N[(N[(x + a), $MachinePrecision] + N[(z - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -1 \cdot 10^{+139}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -1.7 \cdot 10^{-160}:\\
\;\;\;\;\left(x + y \cdot b\right) - \left(y + -1\right) \cdot z\\
\mathbf{elif}\;b \leq 2.4 \cdot 10^{-24}:\\
\;\;\;\;\left(x + a\right) + \left(z - t \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.00000000000000003e139 or 2.3999999999999998e-24 < b Initial program 92.5%
Taylor expanded in x around inf
Simplified80.5%
if -1.00000000000000003e139 < b < -1.70000000000000011e-160Initial program 98.3%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified98.4%
Taylor expanded in z around inf
*-lowering-*.f64N/A
--lowering--.f6470.2%
Simplified70.2%
Taylor expanded in y around inf
*-lowering-*.f6465.6%
Simplified65.6%
if -1.70000000000000011e-160 < b < 2.3999999999999998e-24Initial program 98.6%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified98.6%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6479.1%
Simplified79.1%
Taylor expanded in b around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6476.5%
Simplified76.5%
Final simplification75.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* (- (+ y t) 2.0) b))))
(if (<= b -7.4e+41)
t_1
(if (<= b -3.8e-60)
(+ (* z (- 1.0 y)) (* b (- y 2.0)))
(if (<= b 2.1e-24) (+ (+ x a) (- z (* t a))) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (((y + t) - 2.0) * b);
double tmp;
if (b <= -7.4e+41) {
tmp = t_1;
} else if (b <= -3.8e-60) {
tmp = (z * (1.0 - y)) + (b * (y - 2.0));
} else if (b <= 2.1e-24) {
tmp = (x + a) + (z - (t * 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 = x + (((y + t) - 2.0d0) * b)
if (b <= (-7.4d+41)) then
tmp = t_1
else if (b <= (-3.8d-60)) then
tmp = (z * (1.0d0 - y)) + (b * (y - 2.0d0))
else if (b <= 2.1d-24) then
tmp = (x + a) + (z - (t * 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 = x + (((y + t) - 2.0) * b);
double tmp;
if (b <= -7.4e+41) {
tmp = t_1;
} else if (b <= -3.8e-60) {
tmp = (z * (1.0 - y)) + (b * (y - 2.0));
} else if (b <= 2.1e-24) {
tmp = (x + a) + (z - (t * a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (((y + t) - 2.0) * b) tmp = 0 if b <= -7.4e+41: tmp = t_1 elif b <= -3.8e-60: tmp = (z * (1.0 - y)) + (b * (y - 2.0)) elif b <= 2.1e-24: tmp = (x + a) + (z - (t * a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(Float64(Float64(y + t) - 2.0) * b)) tmp = 0.0 if (b <= -7.4e+41) tmp = t_1; elseif (b <= -3.8e-60) tmp = Float64(Float64(z * Float64(1.0 - y)) + Float64(b * Float64(y - 2.0))); elseif (b <= 2.1e-24) tmp = Float64(Float64(x + a) + Float64(z - Float64(t * a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (((y + t) - 2.0) * b); tmp = 0.0; if (b <= -7.4e+41) tmp = t_1; elseif (b <= -3.8e-60) tmp = (z * (1.0 - y)) + (b * (y - 2.0)); elseif (b <= 2.1e-24) tmp = (x + a) + (z - (t * a)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -7.4e+41], t$95$1, If[LessEqual[b, -3.8e-60], N[(N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision] + N[(b * N[(y - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.1e-24], N[(N[(x + a), $MachinePrecision] + N[(z - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -7.4 \cdot 10^{+41}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -3.8 \cdot 10^{-60}:\\
\;\;\;\;z \cdot \left(1 - y\right) + b \cdot \left(y - 2\right)\\
\mathbf{elif}\;b \leq 2.1 \cdot 10^{-24}:\\
\;\;\;\;\left(x + a\right) + \left(z - t \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -7.39999999999999962e41 or 2.0999999999999999e-24 < b Initial program 92.8%
Taylor expanded in x around inf
Simplified78.8%
if -7.39999999999999962e41 < b < -3.79999999999999994e-60Initial program 100.0%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f6482.7%
Simplified82.7%
Taylor expanded in y around inf
Simplified77.1%
if -3.79999999999999994e-60 < b < 2.0999999999999999e-24Initial program 99.0%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified99.0%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6475.3%
Simplified75.3%
Taylor expanded in b around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6470.5%
Simplified70.5%
Final simplification75.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* a (- 1.0 t))))
(if (<= b -5.8e+150)
(+ z (+ (* y (- b z)) (* b (+ t -2.0))))
(if (<= b 2.8e+26)
(+ x (+ t_1 (* z (- 1.0 y))))
(+ (* (- (+ y t) 2.0) b) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = a * (1.0 - t);
double tmp;
if (b <= -5.8e+150) {
tmp = z + ((y * (b - z)) + (b * (t + -2.0)));
} else if (b <= 2.8e+26) {
tmp = x + (t_1 + (z * (1.0 - y)));
} else {
tmp = (((y + t) - 2.0) * b) + 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 = a * (1.0d0 - t)
if (b <= (-5.8d+150)) then
tmp = z + ((y * (b - z)) + (b * (t + (-2.0d0))))
else if (b <= 2.8d+26) then
tmp = x + (t_1 + (z * (1.0d0 - y)))
else
tmp = (((y + t) - 2.0d0) * b) + 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 = a * (1.0 - t);
double tmp;
if (b <= -5.8e+150) {
tmp = z + ((y * (b - z)) + (b * (t + -2.0)));
} else if (b <= 2.8e+26) {
tmp = x + (t_1 + (z * (1.0 - y)));
} else {
tmp = (((y + t) - 2.0) * b) + t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = a * (1.0 - t) tmp = 0 if b <= -5.8e+150: tmp = z + ((y * (b - z)) + (b * (t + -2.0))) elif b <= 2.8e+26: tmp = x + (t_1 + (z * (1.0 - y))) else: tmp = (((y + t) - 2.0) * b) + t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(a * Float64(1.0 - t)) tmp = 0.0 if (b <= -5.8e+150) tmp = Float64(z + Float64(Float64(y * Float64(b - z)) + Float64(b * Float64(t + -2.0)))); elseif (b <= 2.8e+26) tmp = Float64(x + Float64(t_1 + Float64(z * Float64(1.0 - y)))); else tmp = Float64(Float64(Float64(Float64(y + t) - 2.0) * b) + t_1); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = a * (1.0 - t); tmp = 0.0; if (b <= -5.8e+150) tmp = z + ((y * (b - z)) + (b * (t + -2.0))); elseif (b <= 2.8e+26) tmp = x + (t_1 + (z * (1.0 - y))); else tmp = (((y + t) - 2.0) * b) + t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -5.8e+150], N[(z + N[(N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision] + N[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.8e+26], N[(x + N[(t$95$1 + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision] + t$95$1), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \left(1 - t\right)\\
\mathbf{if}\;b \leq -5.8 \cdot 10^{+150}:\\
\;\;\;\;z + \left(y \cdot \left(b - z\right) + b \cdot \left(t + -2\right)\right)\\
\mathbf{elif}\;b \leq 2.8 \cdot 10^{+26}:\\
\;\;\;\;x + \left(t\_1 + z \cdot \left(1 - y\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(y + t\right) - 2\right) \cdot b + t\_1\\
\end{array}
\end{array}
if b < -5.80000000000000022e150Initial program 96.8%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f6491.0%
Simplified91.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6491.0%
Simplified91.0%
if -5.80000000000000022e150 < b < 2.8e26Initial program 98.0%
Taylor expanded in b around 0
sub-negN/A
+-lowering-+.f64N/A
+-commutativeN/A
distribute-neg-inN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
Simplified81.6%
if 2.8e26 < b Initial program 90.4%
Taylor expanded in a around inf
sub-negN/A
neg-mul-1N/A
metadata-evalN/A
distribute-lft-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6483.4%
Simplified83.4%
Final simplification83.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* (- (+ y t) 2.0) b))))
(if (<= b -2.05e-25)
t_1
(if (<= b 9.5e+68) (+ x (+ (* a (- 1.0 t)) (* z (- 1.0 y)))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (((y + t) - 2.0) * b);
double tmp;
if (b <= -2.05e-25) {
tmp = t_1;
} else if (b <= 9.5e+68) {
tmp = x + ((a * (1.0 - t)) + (z * (1.0 - y)));
} 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 = x + (((y + t) - 2.0d0) * b)
if (b <= (-2.05d-25)) then
tmp = t_1
else if (b <= 9.5d+68) then
tmp = x + ((a * (1.0d0 - t)) + (z * (1.0d0 - y)))
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 = x + (((y + t) - 2.0) * b);
double tmp;
if (b <= -2.05e-25) {
tmp = t_1;
} else if (b <= 9.5e+68) {
tmp = x + ((a * (1.0 - t)) + (z * (1.0 - y)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (((y + t) - 2.0) * b) tmp = 0 if b <= -2.05e-25: tmp = t_1 elif b <= 9.5e+68: tmp = x + ((a * (1.0 - t)) + (z * (1.0 - y))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(Float64(Float64(y + t) - 2.0) * b)) tmp = 0.0 if (b <= -2.05e-25) tmp = t_1; elseif (b <= 9.5e+68) tmp = Float64(x + Float64(Float64(a * Float64(1.0 - t)) + Float64(z * Float64(1.0 - y)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (((y + t) - 2.0) * b); tmp = 0.0; if (b <= -2.05e-25) tmp = t_1; elseif (b <= 9.5e+68) tmp = x + ((a * (1.0 - t)) + (z * (1.0 - y))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -2.05e-25], t$95$1, If[LessEqual[b, 9.5e+68], N[(x + N[(N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision] + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -2.05 \cdot 10^{-25}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 9.5 \cdot 10^{+68}:\\
\;\;\;\;x + \left(a \cdot \left(1 - t\right) + z \cdot \left(1 - y\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -2.04999999999999994e-25 or 9.50000000000000069e68 < b Initial program 92.9%
Taylor expanded in x around inf
Simplified81.0%
if -2.04999999999999994e-25 < b < 9.50000000000000069e68Initial program 98.4%
Taylor expanded in b around 0
sub-negN/A
+-lowering-+.f64N/A
+-commutativeN/A
distribute-neg-inN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
Simplified83.9%
Final simplification82.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* y (- b z))))
(if (<= y -1900000000.0)
t_1
(if (<= y -1.45e-69)
(+ z (+ x a))
(if (<= y 600000000.0) (+ z (* b (+ t -2.0))) 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 <= -1900000000.0) {
tmp = t_1;
} else if (y <= -1.45e-69) {
tmp = z + (x + a);
} else if (y <= 600000000.0) {
tmp = z + (b * (t + -2.0));
} 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 * (b - z)
if (y <= (-1900000000.0d0)) then
tmp = t_1
else if (y <= (-1.45d-69)) then
tmp = z + (x + a)
else if (y <= 600000000.0d0) then
tmp = z + (b * (t + (-2.0d0)))
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 * (b - z);
double tmp;
if (y <= -1900000000.0) {
tmp = t_1;
} else if (y <= -1.45e-69) {
tmp = z + (x + a);
} else if (y <= 600000000.0) {
tmp = z + (b * (t + -2.0));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y * (b - z) tmp = 0 if y <= -1900000000.0: tmp = t_1 elif y <= -1.45e-69: tmp = z + (x + a) elif y <= 600000000.0: tmp = z + (b * (t + -2.0)) 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 <= -1900000000.0) tmp = t_1; elseif (y <= -1.45e-69) tmp = Float64(z + Float64(x + a)); elseif (y <= 600000000.0) tmp = Float64(z + Float64(b * Float64(t + -2.0))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y * (b - z); tmp = 0.0; if (y <= -1900000000.0) tmp = t_1; elseif (y <= -1.45e-69) tmp = z + (x + a); elseif (y <= 600000000.0) tmp = z + (b * (t + -2.0)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1900000000.0], t$95$1, If[LessEqual[y, -1.45e-69], N[(z + N[(x + a), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 600000000.0], N[(z + N[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(b - z\right)\\
\mathbf{if}\;y \leq -1900000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.45 \cdot 10^{-69}:\\
\;\;\;\;z + \left(x + a\right)\\
\mathbf{elif}\;y \leq 600000000:\\
\;\;\;\;z + b \cdot \left(t + -2\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.9e9 or 6e8 < y Initial program 93.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6467.1%
Simplified67.1%
if -1.9e9 < y < -1.4499999999999999e-69Initial program 99.9%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified99.9%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6498.3%
Simplified98.3%
Taylor expanded in z around inf
Simplified65.3%
if -1.4499999999999999e-69 < y < 6e8Initial program 98.0%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f6455.6%
Simplified55.6%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6455.7%
Simplified55.7%
Final simplification62.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* y (- b z))))
(if (<= y -1700000000.0)
t_1
(if (<= y 1.55e-229)
(+ z (+ x a))
(if (<= y 5.8e+74) (* t (- b a)) 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 <= -1700000000.0) {
tmp = t_1;
} else if (y <= 1.55e-229) {
tmp = z + (x + a);
} else if (y <= 5.8e+74) {
tmp = t * (b - 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 * (b - z)
if (y <= (-1700000000.0d0)) then
tmp = t_1
else if (y <= 1.55d-229) then
tmp = z + (x + a)
else if (y <= 5.8d+74) then
tmp = t * (b - 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 * (b - z);
double tmp;
if (y <= -1700000000.0) {
tmp = t_1;
} else if (y <= 1.55e-229) {
tmp = z + (x + a);
} else if (y <= 5.8e+74) {
tmp = t * (b - a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y * (b - z) tmp = 0 if y <= -1700000000.0: tmp = t_1 elif y <= 1.55e-229: tmp = z + (x + a) elif y <= 5.8e+74: tmp = t * (b - a) 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 <= -1700000000.0) tmp = t_1; elseif (y <= 1.55e-229) tmp = Float64(z + Float64(x + a)); elseif (y <= 5.8e+74) tmp = Float64(t * Float64(b - a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y * (b - z); tmp = 0.0; if (y <= -1700000000.0) tmp = t_1; elseif (y <= 1.55e-229) tmp = z + (x + a); elseif (y <= 5.8e+74) tmp = t * (b - a); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1700000000.0], t$95$1, If[LessEqual[y, 1.55e-229], N[(z + N[(x + a), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5.8e+74], N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(b - z\right)\\
\mathbf{if}\;y \leq -1700000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.55 \cdot 10^{-229}:\\
\;\;\;\;z + \left(x + a\right)\\
\mathbf{elif}\;y \leq 5.8 \cdot 10^{+74}:\\
\;\;\;\;t \cdot \left(b - a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.7e9 or 5.8000000000000005e74 < y Initial program 92.3%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6470.0%
Simplified70.0%
if -1.7e9 < y < 1.55e-229Initial program 98.7%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified100.0%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6499.5%
Simplified99.5%
Taylor expanded in z around inf
Simplified54.6%
if 1.55e-229 < y < 5.8000000000000005e74Initial program 98.2%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6452.2%
Simplified52.2%
Final simplification61.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* y (- b z))))
(if (<= y -2050000000.0)
t_1
(if (<= y 3.2e-242)
(+ x (* t b))
(if (<= y 2.1e+76) (* t (- b a)) 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 <= -2050000000.0) {
tmp = t_1;
} else if (y <= 3.2e-242) {
tmp = x + (t * b);
} else if (y <= 2.1e+76) {
tmp = t * (b - 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 * (b - z)
if (y <= (-2050000000.0d0)) then
tmp = t_1
else if (y <= 3.2d-242) then
tmp = x + (t * b)
else if (y <= 2.1d+76) then
tmp = t * (b - 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 * (b - z);
double tmp;
if (y <= -2050000000.0) {
tmp = t_1;
} else if (y <= 3.2e-242) {
tmp = x + (t * b);
} else if (y <= 2.1e+76) {
tmp = t * (b - a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y * (b - z) tmp = 0 if y <= -2050000000.0: tmp = t_1 elif y <= 3.2e-242: tmp = x + (t * b) elif y <= 2.1e+76: tmp = t * (b - a) 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 <= -2050000000.0) tmp = t_1; elseif (y <= 3.2e-242) tmp = Float64(x + Float64(t * b)); elseif (y <= 2.1e+76) tmp = Float64(t * Float64(b - a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y * (b - z); tmp = 0.0; if (y <= -2050000000.0) tmp = t_1; elseif (y <= 3.2e-242) tmp = x + (t * b); elseif (y <= 2.1e+76) tmp = t * (b - a); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2050000000.0], t$95$1, If[LessEqual[y, 3.2e-242], N[(x + N[(t * b), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.1e+76], N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(b - z\right)\\
\mathbf{if}\;y \leq -2050000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.2 \cdot 10^{-242}:\\
\;\;\;\;x + t \cdot b\\
\mathbf{elif}\;y \leq 2.1 \cdot 10^{+76}:\\
\;\;\;\;t \cdot \left(b - a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.05e9 or 2.10000000000000007e76 < y Initial program 92.3%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6470.0%
Simplified70.0%
if -2.05e9 < y < 3.19999999999999999e-242Initial program 98.7%
Taylor expanded in x around inf
Simplified55.1%
Taylor expanded in t around inf
Simplified41.6%
if 3.19999999999999999e-242 < y < 2.10000000000000007e76Initial program 98.3%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6451.3%
Simplified51.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* b (+ t -2.0))))
(if (<= b -1.35e-72)
t_1
(if (<= b 1.3e+84) (* a (- 1.0 t)) (if (<= b 3.95e+105) (* y b) 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 <= -1.35e-72) {
tmp = t_1;
} else if (b <= 1.3e+84) {
tmp = a * (1.0 - t);
} else if (b <= 3.95e+105) {
tmp = y * b;
} 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 <= (-1.35d-72)) then
tmp = t_1
else if (b <= 1.3d+84) then
tmp = a * (1.0d0 - t)
else if (b <= 3.95d+105) then
tmp = y * b
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 <= -1.35e-72) {
tmp = t_1;
} else if (b <= 1.3e+84) {
tmp = a * (1.0 - t);
} else if (b <= 3.95e+105) {
tmp = y * b;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = b * (t + -2.0) tmp = 0 if b <= -1.35e-72: tmp = t_1 elif b <= 1.3e+84: tmp = a * (1.0 - t) elif b <= 3.95e+105: tmp = y * b 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 <= -1.35e-72) tmp = t_1; elseif (b <= 1.3e+84) tmp = Float64(a * Float64(1.0 - t)); elseif (b <= 3.95e+105) tmp = Float64(y * b); 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 <= -1.35e-72) tmp = t_1; elseif (b <= 1.3e+84) tmp = a * (1.0 - t); elseif (b <= 3.95e+105) tmp = y * b; 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, -1.35e-72], t$95$1, If[LessEqual[b, 1.3e+84], N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 3.95e+105], N[(y * b), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := b \cdot \left(t + -2\right)\\
\mathbf{if}\;b \leq -1.35 \cdot 10^{-72}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.3 \cdot 10^{+84}:\\
\;\;\;\;a \cdot \left(1 - t\right)\\
\mathbf{elif}\;b \leq 3.95 \cdot 10^{+105}:\\
\;\;\;\;y \cdot b\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.35e-72 or 3.95e105 < b Initial program 93.7%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified92.2%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6469.9%
Simplified69.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6443.9%
Simplified43.9%
if -1.35e-72 < b < 1.3000000000000001e84Initial program 97.5%
Taylor expanded in a around inf
sub-negN/A
neg-mul-1N/A
metadata-evalN/A
distribute-lft-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6439.9%
Simplified39.9%
if 1.3000000000000001e84 < b < 3.95e105Initial program 100.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f64100.0%
Simplified100.0%
Taylor expanded in b around inf
Simplified100.0%
Final simplification43.0%
(FPCore (x y z t a b) :precision binary64 (if (<= b -3e+204) (* t b) (if (<= b -4.4e+38) x (if (<= b 3.8e+85) (* a (- 1.0 t)) (* y b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -3e+204) {
tmp = t * b;
} else if (b <= -4.4e+38) {
tmp = x;
} else if (b <= 3.8e+85) {
tmp = a * (1.0 - t);
} 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 (b <= (-3d+204)) then
tmp = t * b
else if (b <= (-4.4d+38)) then
tmp = x
else if (b <= 3.8d+85) then
tmp = a * (1.0d0 - t)
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 (b <= -3e+204) {
tmp = t * b;
} else if (b <= -4.4e+38) {
tmp = x;
} else if (b <= 3.8e+85) {
tmp = a * (1.0 - t);
} else {
tmp = y * b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -3e+204: tmp = t * b elif b <= -4.4e+38: tmp = x elif b <= 3.8e+85: tmp = a * (1.0 - t) else: tmp = y * b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -3e+204) tmp = Float64(t * b); elseif (b <= -4.4e+38) tmp = x; elseif (b <= 3.8e+85) tmp = Float64(a * Float64(1.0 - t)); else tmp = Float64(y * b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -3e+204) tmp = t * b; elseif (b <= -4.4e+38) tmp = x; elseif (b <= 3.8e+85) tmp = a * (1.0 - t); else tmp = y * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -3e+204], N[(t * b), $MachinePrecision], If[LessEqual[b, -4.4e+38], x, If[LessEqual[b, 3.8e+85], N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision], N[(y * b), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3 \cdot 10^{+204}:\\
\;\;\;\;t \cdot b\\
\mathbf{elif}\;b \leq -4.4 \cdot 10^{+38}:\\
\;\;\;\;x\\
\mathbf{elif}\;b \leq 3.8 \cdot 10^{+85}:\\
\;\;\;\;a \cdot \left(1 - t\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot b\\
\end{array}
\end{array}
if b < -2.99999999999999983e204Initial program 95.6%
Taylor expanded in x around inf
Simplified91.5%
Taylor expanded in t around inf
*-lowering-*.f6453.2%
Simplified53.2%
if -2.99999999999999983e204 < b < -4.40000000000000013e38Initial program 94.2%
Taylor expanded in x around inf
Simplified30.4%
if -4.40000000000000013e38 < b < 3.79999999999999992e85Initial program 97.9%
Taylor expanded in a around inf
sub-negN/A
neg-mul-1N/A
metadata-evalN/A
distribute-lft-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6437.6%
Simplified37.6%
if 3.79999999999999992e85 < b Initial program 91.1%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6451.1%
Simplified51.1%
Taylor expanded in b around inf
Simplified45.5%
Final simplification39.8%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (+ x (* (- (+ y t) 2.0) b)))) (if (<= b -8.2e-64) t_1 (if (<= b 2.1e-24) (+ (+ x a) (- z (* t a))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (((y + t) - 2.0) * b);
double tmp;
if (b <= -8.2e-64) {
tmp = t_1;
} else if (b <= 2.1e-24) {
tmp = (x + a) + (z - (t * 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 = x + (((y + t) - 2.0d0) * b)
if (b <= (-8.2d-64)) then
tmp = t_1
else if (b <= 2.1d-24) then
tmp = (x + a) + (z - (t * 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 = x + (((y + t) - 2.0) * b);
double tmp;
if (b <= -8.2e-64) {
tmp = t_1;
} else if (b <= 2.1e-24) {
tmp = (x + a) + (z - (t * a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (((y + t) - 2.0) * b) tmp = 0 if b <= -8.2e-64: tmp = t_1 elif b <= 2.1e-24: tmp = (x + a) + (z - (t * a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(Float64(Float64(y + t) - 2.0) * b)) tmp = 0.0 if (b <= -8.2e-64) tmp = t_1; elseif (b <= 2.1e-24) tmp = Float64(Float64(x + a) + Float64(z - Float64(t * a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (((y + t) - 2.0) * b); tmp = 0.0; if (b <= -8.2e-64) tmp = t_1; elseif (b <= 2.1e-24) tmp = (x + a) + (z - (t * a)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -8.2e-64], t$95$1, If[LessEqual[b, 2.1e-24], N[(N[(x + a), $MachinePrecision] + N[(z - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -8.2 \cdot 10^{-64}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 2.1 \cdot 10^{-24}:\\
\;\;\;\;\left(x + a\right) + \left(z - t \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -8.2000000000000001e-64 or 2.0999999999999999e-24 < b Initial program 93.6%
Taylor expanded in x around inf
Simplified75.1%
if -8.2000000000000001e-64 < b < 2.0999999999999999e-24Initial program 99.0%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified99.0%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6475.0%
Simplified75.0%
Taylor expanded in b around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6471.1%
Simplified71.1%
Final simplification73.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* z (- 1.0 y)))))
(if (<= z -1.65e+132)
t_1
(if (<= z 3.4e+48) (+ x (* (- (+ y t) 2.0) b)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (z * (1.0 - y));
double tmp;
if (z <= -1.65e+132) {
tmp = t_1;
} else if (z <= 3.4e+48) {
tmp = x + (((y + t) - 2.0) * b);
} 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 = x + (z * (1.0d0 - y))
if (z <= (-1.65d+132)) then
tmp = t_1
else if (z <= 3.4d+48) then
tmp = x + (((y + t) - 2.0d0) * b)
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 = x + (z * (1.0 - y));
double tmp;
if (z <= -1.65e+132) {
tmp = t_1;
} else if (z <= 3.4e+48) {
tmp = x + (((y + t) - 2.0) * b);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (z * (1.0 - y)) tmp = 0 if z <= -1.65e+132: tmp = t_1 elif z <= 3.4e+48: tmp = x + (((y + t) - 2.0) * b) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(z * Float64(1.0 - y))) tmp = 0.0 if (z <= -1.65e+132) tmp = t_1; elseif (z <= 3.4e+48) tmp = Float64(x + Float64(Float64(Float64(y + t) - 2.0) * b)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (z * (1.0 - y)); tmp = 0.0; if (z <= -1.65e+132) tmp = t_1; elseif (z <= 3.4e+48) tmp = x + (((y + t) - 2.0) * b); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.65e+132], t$95$1, If[LessEqual[z, 3.4e+48], N[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + z \cdot \left(1 - y\right)\\
\mathbf{if}\;z \leq -1.65 \cdot 10^{+132}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 3.4 \cdot 10^{+48}:\\
\;\;\;\;x + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.65000000000000015e132 or 3.4000000000000003e48 < z Initial program 88.4%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified91.0%
Taylor expanded in z around inf
*-lowering-*.f64N/A
--lowering--.f6472.4%
Simplified72.4%
Taylor expanded in b around 0
*-lft-identityN/A
fma-defineN/A
remove-double-negN/A
distribute-rgt-neg-outN/A
mul-1-negN/A
fmm-undefN/A
*-lft-identityN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
associate--r-N/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6470.3%
Simplified70.3%
if -1.65000000000000015e132 < z < 3.4000000000000003e48Initial program 98.9%
Taylor expanded in x around inf
Simplified68.9%
Final simplification69.3%
(FPCore (x y z t a b) :precision binary64 (if (<= y -2700000000.0) (* y b) (if (<= y -7e-70) x (if (<= y 1150000000000.0) (* t b) (* y b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -2700000000.0) {
tmp = y * b;
} else if (y <= -7e-70) {
tmp = x;
} else if (y <= 1150000000000.0) {
tmp = t * b;
} 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 <= (-2700000000.0d0)) then
tmp = y * b
else if (y <= (-7d-70)) then
tmp = x
else if (y <= 1150000000000.0d0) then
tmp = t * b
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 <= -2700000000.0) {
tmp = y * b;
} else if (y <= -7e-70) {
tmp = x;
} else if (y <= 1150000000000.0) {
tmp = t * b;
} else {
tmp = y * b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -2700000000.0: tmp = y * b elif y <= -7e-70: tmp = x elif y <= 1150000000000.0: tmp = t * b else: tmp = y * b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -2700000000.0) tmp = Float64(y * b); elseif (y <= -7e-70) tmp = x; elseif (y <= 1150000000000.0) tmp = Float64(t * b); else tmp = Float64(y * b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -2700000000.0) tmp = y * b; elseif (y <= -7e-70) tmp = x; elseif (y <= 1150000000000.0) tmp = t * b; else tmp = y * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -2700000000.0], N[(y * b), $MachinePrecision], If[LessEqual[y, -7e-70], x, If[LessEqual[y, 1150000000000.0], N[(t * b), $MachinePrecision], N[(y * b), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2700000000:\\
\;\;\;\;y \cdot b\\
\mathbf{elif}\;y \leq -7 \cdot 10^{-70}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 1150000000000:\\
\;\;\;\;t \cdot b\\
\mathbf{else}:\\
\;\;\;\;y \cdot b\\
\end{array}
\end{array}
if y < -2.7e9 or 1.15e12 < y Initial program 93.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6467.1%
Simplified67.1%
Taylor expanded in b around inf
Simplified40.8%
if -2.7e9 < y < -6.99999999999999949e-70Initial program 99.9%
Taylor expanded in x around inf
Simplified35.3%
if -6.99999999999999949e-70 < y < 1.15e12Initial program 98.0%
Taylor expanded in x around inf
Simplified56.3%
Taylor expanded in t around inf
*-lowering-*.f6424.7%
Simplified24.7%
Final simplification33.8%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* t (- b a)))) (if (<= t -8.5e+72) t_1 (if (<= t 2e+110) (+ x (* b (+ y -2.0))) 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.5e+72) {
tmp = t_1;
} else if (t <= 2e+110) {
tmp = x + (b * (y + -2.0));
} 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.5d+72)) then
tmp = t_1
else if (t <= 2d+110) then
tmp = x + (b * (y + (-2.0d0)))
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.5e+72) {
tmp = t_1;
} else if (t <= 2e+110) {
tmp = x + (b * (y + -2.0));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = t * (b - a) tmp = 0 if t <= -8.5e+72: tmp = t_1 elif t <= 2e+110: tmp = x + (b * (y + -2.0)) 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.5e+72) tmp = t_1; elseif (t <= 2e+110) tmp = Float64(x + Float64(b * Float64(y + -2.0))); 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.5e+72) tmp = t_1; elseif (t <= 2e+110) tmp = x + (b * (y + -2.0)); 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.5e+72], t$95$1, If[LessEqual[t, 2e+110], N[(x + N[(b * N[(y + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -8.5 \cdot 10^{+72}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2 \cdot 10^{+110}:\\
\;\;\;\;x + b \cdot \left(y + -2\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -8.5000000000000004e72 or 2e110 < t Initial program 92.0%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6480.6%
Simplified80.6%
if -8.5000000000000004e72 < t < 2e110Initial program 97.6%
Taylor expanded in x around inf
Simplified61.0%
Taylor expanded in t around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6456.8%
Simplified56.8%
Final simplification65.0%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* t (- b a)))) (if (<= t -6e+72) t_1 (if (<= t 2e+110) (+ x (* y b)) 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 <= -6e+72) {
tmp = t_1;
} else if (t <= 2e+110) {
tmp = x + (y * b);
} 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 <= (-6d+72)) then
tmp = t_1
else if (t <= 2d+110) then
tmp = x + (y * b)
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 <= -6e+72) {
tmp = t_1;
} else if (t <= 2e+110) {
tmp = x + (y * b);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = t * (b - a) tmp = 0 if t <= -6e+72: tmp = t_1 elif t <= 2e+110: tmp = x + (y * b) 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 <= -6e+72) tmp = t_1; elseif (t <= 2e+110) tmp = Float64(x + Float64(y * b)); 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 <= -6e+72) tmp = t_1; elseif (t <= 2e+110) tmp = x + (y * b); 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, -6e+72], t$95$1, If[LessEqual[t, 2e+110], N[(x + N[(y * b), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -6 \cdot 10^{+72}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2 \cdot 10^{+110}:\\
\;\;\;\;x + y \cdot b\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -6.00000000000000006e72 or 2e110 < t Initial program 92.0%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6480.6%
Simplified80.6%
if -6.00000000000000006e72 < t < 2e110Initial program 97.6%
Taylor expanded in x around inf
Simplified61.0%
Taylor expanded in y around inf
Simplified45.1%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* t (- b a)))) (if (<= t -1.1e+73) t_1 (if (<= t 2.55e+17) (* y (- b z)) 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.1e+73) {
tmp = t_1;
} else if (t <= 2.55e+17) {
tmp = y * (b - z);
} 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 <= (-1.1d+73)) then
tmp = t_1
else if (t <= 2.55d+17) then
tmp = y * (b - z)
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 <= -1.1e+73) {
tmp = t_1;
} else if (t <= 2.55e+17) {
tmp = y * (b - z);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = t * (b - a) tmp = 0 if t <= -1.1e+73: tmp = t_1 elif t <= 2.55e+17: tmp = y * (b - z) 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.1e+73) tmp = t_1; elseif (t <= 2.55e+17) tmp = Float64(y * Float64(b - z)); 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 <= -1.1e+73) tmp = t_1; elseif (t <= 2.55e+17) tmp = y * (b - z); 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, -1.1e+73], t$95$1, If[LessEqual[t, 2.55e+17], N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -1.1 \cdot 10^{+73}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2.55 \cdot 10^{+17}:\\
\;\;\;\;y \cdot \left(b - z\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -1.1e73 or 2.55e17 < t Initial program 92.8%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6472.6%
Simplified72.6%
if -1.1e73 < t < 2.55e17Initial program 97.9%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6441.5%
Simplified41.5%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* t (- b a)))) (if (<= t -900000000.0) t_1 (if (<= t 1.75e+15) (* b (+ y -2.0)) 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 <= -900000000.0) {
tmp = t_1;
} else if (t <= 1.75e+15) {
tmp = b * (y + -2.0);
} 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 <= (-900000000.0d0)) then
tmp = t_1
else if (t <= 1.75d+15) then
tmp = b * (y + (-2.0d0))
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 <= -900000000.0) {
tmp = t_1;
} else if (t <= 1.75e+15) {
tmp = b * (y + -2.0);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = t * (b - a) tmp = 0 if t <= -900000000.0: tmp = t_1 elif t <= 1.75e+15: tmp = b * (y + -2.0) 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 <= -900000000.0) tmp = t_1; elseif (t <= 1.75e+15) tmp = Float64(b * Float64(y + -2.0)); 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 <= -900000000.0) tmp = t_1; elseif (t <= 1.75e+15) tmp = b * (y + -2.0); 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, -900000000.0], t$95$1, If[LessEqual[t, 1.75e+15], N[(b * N[(y + -2.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -900000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 1.75 \cdot 10^{+15}:\\
\;\;\;\;b \cdot \left(y + -2\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -9e8 or 1.75e15 < t Initial program 93.5%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6468.7%
Simplified68.7%
if -9e8 < t < 1.75e15Initial program 97.7%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f6465.3%
Simplified65.3%
Taylor expanded in y around inf
Simplified65.3%
Taylor expanded in z around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6440.5%
Simplified40.5%
Final simplification54.1%
(FPCore (x y z t a b) :precision binary64 (if (<= b -4.2e-73) (* b (+ t -2.0)) (if (<= b 1.4e+69) (* a (- 1.0 t)) (* b (+ y -2.0)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -4.2e-73) {
tmp = b * (t + -2.0);
} else if (b <= 1.4e+69) {
tmp = a * (1.0 - t);
} else {
tmp = b * (y + -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 <= (-4.2d-73)) then
tmp = b * (t + (-2.0d0))
else if (b <= 1.4d+69) then
tmp = a * (1.0d0 - t)
else
tmp = b * (y + (-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 <= -4.2e-73) {
tmp = b * (t + -2.0);
} else if (b <= 1.4e+69) {
tmp = a * (1.0 - t);
} else {
tmp = b * (y + -2.0);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -4.2e-73: tmp = b * (t + -2.0) elif b <= 1.4e+69: tmp = a * (1.0 - t) else: tmp = b * (y + -2.0) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -4.2e-73) tmp = Float64(b * Float64(t + -2.0)); elseif (b <= 1.4e+69) tmp = Float64(a * Float64(1.0 - t)); else tmp = Float64(b * Float64(y + -2.0)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -4.2e-73) tmp = b * (t + -2.0); elseif (b <= 1.4e+69) tmp = a * (1.0 - t); else tmp = b * (y + -2.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -4.2e-73], N[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.4e+69], N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision], N[(b * N[(y + -2.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -4.2 \cdot 10^{-73}:\\
\;\;\;\;b \cdot \left(t + -2\right)\\
\mathbf{elif}\;b \leq 1.4 \cdot 10^{+69}:\\
\;\;\;\;a \cdot \left(1 - t\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(y + -2\right)\\
\end{array}
\end{array}
if b < -4.1999999999999997e-73Initial program 96.1%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified97.4%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6473.1%
Simplified73.1%
Taylor expanded in b around inf
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6438.3%
Simplified38.3%
if -4.1999999999999997e-73 < b < 1.39999999999999991e69Initial program 98.3%
Taylor expanded in a around inf
sub-negN/A
neg-mul-1N/A
metadata-evalN/A
distribute-lft-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6440.3%
Simplified40.3%
if 1.39999999999999991e69 < b Initial program 90.0%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
--lowering--.f6483.4%
Simplified83.4%
Taylor expanded in y around inf
Simplified59.3%
Taylor expanded in z around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6457.7%
Simplified57.7%
Final simplification43.7%
(FPCore (x y z t a b) :precision binary64 (if (<= t -5.6e+72) (* t b) (if (<= t 1.02e+33) x (* t b))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -5.6e+72) {
tmp = t * b;
} else if (t <= 1.02e+33) {
tmp = x;
} 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) :: tmp
if (t <= (-5.6d+72)) then
tmp = t * b
else if (t <= 1.02d+33) then
tmp = x
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 tmp;
if (t <= -5.6e+72) {
tmp = t * b;
} else if (t <= 1.02e+33) {
tmp = x;
} else {
tmp = t * b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= -5.6e+72: tmp = t * b elif t <= 1.02e+33: tmp = x else: tmp = t * b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -5.6e+72) tmp = Float64(t * b); elseif (t <= 1.02e+33) tmp = x; else tmp = Float64(t * b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= -5.6e+72) tmp = t * b; elseif (t <= 1.02e+33) tmp = x; else tmp = t * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -5.6e+72], N[(t * b), $MachinePrecision], If[LessEqual[t, 1.02e+33], x, N[(t * b), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -5.6 \cdot 10^{+72}:\\
\;\;\;\;t \cdot b\\
\mathbf{elif}\;t \leq 1.02 \cdot 10^{+33}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;t \cdot b\\
\end{array}
\end{array}
if t < -5.5999999999999998e72 or 1.02000000000000001e33 < t Initial program 92.8%
Taylor expanded in x around inf
Simplified55.1%
Taylor expanded in t around inf
*-lowering-*.f6439.7%
Simplified39.7%
if -5.5999999999999998e72 < t < 1.02000000000000001e33Initial program 97.9%
Taylor expanded in x around inf
Simplified22.4%
Final simplification29.9%
(FPCore (x y z t a b) :precision binary64 (if (<= z -1.55e+118) z (if (<= z 1.05e+143) x z)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -1.55e+118) {
tmp = z;
} else if (z <= 1.05e+143) {
tmp = x;
} else {
tmp = z;
}
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 (z <= (-1.55d+118)) then
tmp = z
else if (z <= 1.05d+143) then
tmp = x
else
tmp = z
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 (z <= -1.55e+118) {
tmp = z;
} else if (z <= 1.05e+143) {
tmp = x;
} else {
tmp = z;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if z <= -1.55e+118: tmp = z elif z <= 1.05e+143: tmp = x else: tmp = z return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -1.55e+118) tmp = z; elseif (z <= 1.05e+143) tmp = x; else tmp = z; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (z <= -1.55e+118) tmp = z; elseif (z <= 1.05e+143) tmp = x; else tmp = z; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -1.55e+118], z, If[LessEqual[z, 1.05e+143], x, z]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.55 \cdot 10^{+118}:\\
\;\;\;\;z\\
\mathbf{elif}\;z \leq 1.05 \cdot 10^{+143}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;z\\
\end{array}
\end{array}
if z < -1.54999999999999993e118 or 1.04999999999999994e143 < z Initial program 87.5%
Taylor expanded in t around 0
associate-+r+N/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
mul-1-negN/A
Simplified88.8%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6460.8%
Simplified60.8%
Taylor expanded in z around inf
Simplified28.4%
if -1.54999999999999993e118 < z < 1.04999999999999994e143Initial program 98.9%
Taylor expanded in x around inf
Simplified19.3%
(FPCore (x y z t a b) :precision binary64 x)
double code(double x, double y, double z, double t, double a, double b) {
return 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 = x
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x;
}
def code(x, y, z, t, a, b): return x
function code(x, y, z, t, a, b) return x end
function tmp = code(x, y, z, t, a, b) tmp = x; end
code[x_, y_, z_, t_, a_, b_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 95.7%
Taylor expanded in x around inf
Simplified15.9%
herbie shell --seed 2024160
(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)))