
(FPCore (x y z t a b) :precision binary64 (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))
double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a)) + (((y + t) - 2.0d0) * b)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
def code(x, y, z, t, a, b): return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b)) end
function tmp = code(x, y, z, t, a, b) tmp = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 27 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b) :precision binary64 (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))
double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a)) + (((y + t) - 2.0d0) * b)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
def code(x, y, z, t, a, b): return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b)) end
function tmp = code(x, y, z, t, a, b) tmp = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b
\end{array}
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* a (- 1.0 t))) (t_2 (+ x (* z (- 1.0 y)))))
(if (<= (+ (+ t_2 t_1) (* (- (+ y t) 2.0) b)) INFINITY)
(+ (+ t_2 (* b (+ (+ y t) -2.0))) t_1)
(+ a (* y (- b z))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = a * (1.0 - t);
double t_2 = x + (z * (1.0 - y));
double tmp;
if (((t_2 + t_1) + (((y + t) - 2.0) * b)) <= ((double) INFINITY)) {
tmp = (t_2 + (b * ((y + t) + -2.0))) + t_1;
} else {
tmp = a + (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);
double t_2 = x + (z * (1.0 - y));
double tmp;
if (((t_2 + t_1) + (((y + t) - 2.0) * b)) <= Double.POSITIVE_INFINITY) {
tmp = (t_2 + (b * ((y + t) + -2.0))) + t_1;
} else {
tmp = a + (y * (b - z));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = a * (1.0 - t) t_2 = x + (z * (1.0 - y)) tmp = 0 if ((t_2 + t_1) + (((y + t) - 2.0) * b)) <= math.inf: tmp = (t_2 + (b * ((y + t) + -2.0))) + t_1 else: tmp = a + (y * (b - z)) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(a * Float64(1.0 - t)) t_2 = Float64(x + Float64(z * Float64(1.0 - y))) tmp = 0.0 if (Float64(Float64(t_2 + t_1) + Float64(Float64(Float64(y + t) - 2.0) * b)) <= Inf) tmp = Float64(Float64(t_2 + Float64(b * Float64(Float64(y + t) + -2.0))) + t_1); else tmp = Float64(a + 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); t_2 = x + (z * (1.0 - y)); tmp = 0.0; if (((t_2 + t_1) + (((y + t) - 2.0) * b)) <= Inf) tmp = (t_2 + (b * ((y + t) + -2.0))) + t_1; else tmp = a + (y * (b - z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(t$95$2 + t$95$1), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(t$95$2 + N[(b * N[(N[(y + t), $MachinePrecision] + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$1), $MachinePrecision], N[(a + N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \left(1 - t\right)\\
t_2 := x + z \cdot \left(1 - y\right)\\
\mathbf{if}\;\left(t\_2 + t\_1\right) + \left(\left(y + t\right) - 2\right) \cdot b \leq \infty:\\
\;\;\;\;\left(t\_2 + b \cdot \left(\left(y + t\right) + -2\right)\right) + t\_1\\
\mathbf{else}:\\
\;\;\;\;a + 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%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval100.0%
Applied egg-rr100.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%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval0.0%
Applied egg-rr0.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6466.7%
Simplified66.7%
Taylor expanded in t around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6468.6%
Simplified68.6%
Final simplification98.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 INFINITY) t_1 (+ a (* y (- b z))))))
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 = a + (y * (b - z));
}
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 = a + (y * (b - z));
}
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 = a + (y * (b - z)) 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(a + Float64(y * Float64(b - z))); 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 = a + (y * (b - z)); 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[(a + N[(y * N[(b - z), $MachinePrecision]), $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}:\\
\;\;\;\;a + 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%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval0.0%
Applied egg-rr0.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6466.7%
Simplified66.7%
Taylor expanded in t around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6468.6%
Simplified68.6%
Final simplification98.9%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* z (- 1.0 y))) (t_2 (- (+ (* b (+ y (+ t -2.0))) t_1) (* t a))))
(if (<= b -3.3e+113)
(+ x (* (- (+ y t) 2.0) b))
(if (<= b -1.5e-59)
t_2
(if (<= b 2.15e-19) (+ x (+ t_1 (* a (- 1.0 t)))) t_2)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z * (1.0 - y);
double t_2 = ((b * (y + (t + -2.0))) + t_1) - (t * a);
double tmp;
if (b <= -3.3e+113) {
tmp = x + (((y + t) - 2.0) * b);
} else if (b <= -1.5e-59) {
tmp = t_2;
} else if (b <= 2.15e-19) {
tmp = x + (t_1 + (a * (1.0 - t)));
} else {
tmp = 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 = z * (1.0d0 - y)
t_2 = ((b * (y + (t + (-2.0d0)))) + t_1) - (t * a)
if (b <= (-3.3d+113)) then
tmp = x + (((y + t) - 2.0d0) * b)
else if (b <= (-1.5d-59)) then
tmp = t_2
else if (b <= 2.15d-19) then
tmp = x + (t_1 + (a * (1.0d0 - t)))
else
tmp = 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 = z * (1.0 - y);
double t_2 = ((b * (y + (t + -2.0))) + t_1) - (t * a);
double tmp;
if (b <= -3.3e+113) {
tmp = x + (((y + t) - 2.0) * b);
} else if (b <= -1.5e-59) {
tmp = t_2;
} else if (b <= 2.15e-19) {
tmp = x + (t_1 + (a * (1.0 - t)));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = z * (1.0 - y) t_2 = ((b * (y + (t + -2.0))) + t_1) - (t * a) tmp = 0 if b <= -3.3e+113: tmp = x + (((y + t) - 2.0) * b) elif b <= -1.5e-59: tmp = t_2 elif b <= 2.15e-19: tmp = x + (t_1 + (a * (1.0 - t))) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(z * Float64(1.0 - y)) t_2 = Float64(Float64(Float64(b * Float64(y + Float64(t + -2.0))) + t_1) - Float64(t * a)) tmp = 0.0 if (b <= -3.3e+113) tmp = Float64(x + Float64(Float64(Float64(y + t) - 2.0) * b)); elseif (b <= -1.5e-59) tmp = t_2; elseif (b <= 2.15e-19) tmp = Float64(x + Float64(t_1 + Float64(a * Float64(1.0 - t)))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = z * (1.0 - y); t_2 = ((b * (y + (t + -2.0))) + t_1) - (t * a); tmp = 0.0; if (b <= -3.3e+113) tmp = x + (((y + t) - 2.0) * b); elseif (b <= -1.5e-59) tmp = t_2; elseif (b <= 2.15e-19) tmp = x + (t_1 + (a * (1.0 - t))); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(b * N[(y + N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$1), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -3.3e+113], N[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -1.5e-59], t$95$2, If[LessEqual[b, 2.15e-19], N[(x + N[(t$95$1 + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \left(1 - y\right)\\
t_2 := \left(b \cdot \left(y + \left(t + -2\right)\right) + t\_1\right) - t \cdot a\\
\mathbf{if}\;b \leq -3.3 \cdot 10^{+113}:\\
\;\;\;\;x + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{elif}\;b \leq -1.5 \cdot 10^{-59}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;b \leq 2.15 \cdot 10^{-19}:\\
\;\;\;\;x + \left(t\_1 + a \cdot \left(1 - t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if b < -3.3000000000000003e113Initial program 88.1%
Taylor expanded in x around inf
Simplified88.2%
if -3.3000000000000003e113 < b < -1.5e-59 or 2.15e-19 < b Initial program 95.9%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval95.9%
Applied egg-rr95.9%
Taylor expanded in x around 0
--lowering--.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+r-N/A
+-commutativeN/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6493.9%
Simplified93.9%
Taylor expanded in t around inf
Simplified89.7%
if -1.5e-59 < b < 2.15e-19Initial program 100.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
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/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
Simplified96.4%
Final simplification92.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* z (- 1.0 y))))
(if (<= z -1.38e+107)
t_1
(if (<= z -0.0275)
(+ x (* t b))
(if (<= z 1.25e-48)
(* b (+ t (+ y -2.0)))
(if (<= z 9.5e+51)
(+ x z)
(if (<= z 4.8e+173) (* a (- 1.0 t)) t_1)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z * (1.0 - y);
double tmp;
if (z <= -1.38e+107) {
tmp = t_1;
} else if (z <= -0.0275) {
tmp = x + (t * b);
} else if (z <= 1.25e-48) {
tmp = b * (t + (y + -2.0));
} else if (z <= 9.5e+51) {
tmp = x + z;
} else if (z <= 4.8e+173) {
tmp = a * (1.0 - t);
} 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 = z * (1.0d0 - y)
if (z <= (-1.38d+107)) then
tmp = t_1
else if (z <= (-0.0275d0)) then
tmp = x + (t * b)
else if (z <= 1.25d-48) then
tmp = b * (t + (y + (-2.0d0)))
else if (z <= 9.5d+51) then
tmp = x + z
else if (z <= 4.8d+173) then
tmp = a * (1.0d0 - t)
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 = z * (1.0 - y);
double tmp;
if (z <= -1.38e+107) {
tmp = t_1;
} else if (z <= -0.0275) {
tmp = x + (t * b);
} else if (z <= 1.25e-48) {
tmp = b * (t + (y + -2.0));
} else if (z <= 9.5e+51) {
tmp = x + z;
} else if (z <= 4.8e+173) {
tmp = a * (1.0 - t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = z * (1.0 - y) tmp = 0 if z <= -1.38e+107: tmp = t_1 elif z <= -0.0275: tmp = x + (t * b) elif z <= 1.25e-48: tmp = b * (t + (y + -2.0)) elif z <= 9.5e+51: tmp = x + z elif z <= 4.8e+173: tmp = a * (1.0 - t) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(z * Float64(1.0 - y)) tmp = 0.0 if (z <= -1.38e+107) tmp = t_1; elseif (z <= -0.0275) tmp = Float64(x + Float64(t * b)); elseif (z <= 1.25e-48) tmp = Float64(b * Float64(t + Float64(y + -2.0))); elseif (z <= 9.5e+51) tmp = Float64(x + z); elseif (z <= 4.8e+173) tmp = Float64(a * Float64(1.0 - t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = z * (1.0 - y); tmp = 0.0; if (z <= -1.38e+107) tmp = t_1; elseif (z <= -0.0275) tmp = x + (t * b); elseif (z <= 1.25e-48) tmp = b * (t + (y + -2.0)); elseif (z <= 9.5e+51) tmp = x + z; elseif (z <= 4.8e+173) tmp = a * (1.0 - t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.38e+107], t$95$1, If[LessEqual[z, -0.0275], N[(x + N[(t * b), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 1.25e-48], N[(b * N[(t + N[(y + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 9.5e+51], N[(x + z), $MachinePrecision], If[LessEqual[z, 4.8e+173], N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \left(1 - y\right)\\
\mathbf{if}\;z \leq -1.38 \cdot 10^{+107}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -0.0275:\\
\;\;\;\;x + t \cdot b\\
\mathbf{elif}\;z \leq 1.25 \cdot 10^{-48}:\\
\;\;\;\;b \cdot \left(t + \left(y + -2\right)\right)\\
\mathbf{elif}\;z \leq 9.5 \cdot 10^{+51}:\\
\;\;\;\;x + z\\
\mathbf{elif}\;z \leq 4.8 \cdot 10^{+173}:\\
\;\;\;\;a \cdot \left(1 - t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.38e107 or 4.7999999999999998e173 < z Initial program 89.9%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6472.7%
Simplified72.7%
if -1.38e107 < z < -0.0275000000000000001Initial program 99.9%
Taylor expanded in x around inf
Simplified54.5%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6447.8%
Simplified47.8%
if -0.0275000000000000001 < z < 1.25e-48Initial program 99.2%
Taylor expanded in b around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval55.7%
Simplified55.7%
if 1.25e-48 < z < 9.4999999999999999e51Initial program 94.1%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified100.0%
Taylor expanded in b around 0
associate-+r+N/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
*-rgt-identityN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f6488.6%
Simplified88.6%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
cancel-sign-sub-invN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6469.5%
Simplified69.5%
Taylor expanded in y around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f6463.5%
Simplified63.5%
if 9.4999999999999999e51 < z < 4.7999999999999998e173Initial program 100.0%
Taylor expanded in a around inf
sub-negN/A
metadata-evalN/A
neg-mul-1N/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--.f6455.1%
Simplified55.1%
Final simplification60.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* z (- 1.0 y))))
(if (<= a -3.8e+212)
(+ a (+ (* b (+ y -2.0)) (* t (- b a))))
(if (<= a 1e-16)
(+ (* (- (+ y t) 2.0) b) (- (+ x t_1) (* t a)))
(+ (+ (* b (+ y (+ t -2.0))) t_1) (* a (- 1.0 t)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z * (1.0 - y);
double tmp;
if (a <= -3.8e+212) {
tmp = a + ((b * (y + -2.0)) + (t * (b - a)));
} else if (a <= 1e-16) {
tmp = (((y + t) - 2.0) * b) + ((x + t_1) - (t * a));
} else {
tmp = ((b * (y + (t + -2.0))) + t_1) + (a * (1.0 - t));
}
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 * (1.0d0 - y)
if (a <= (-3.8d+212)) then
tmp = a + ((b * (y + (-2.0d0))) + (t * (b - a)))
else if (a <= 1d-16) then
tmp = (((y + t) - 2.0d0) * b) + ((x + t_1) - (t * a))
else
tmp = ((b * (y + (t + (-2.0d0)))) + t_1) + (a * (1.0d0 - t))
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 * (1.0 - y);
double tmp;
if (a <= -3.8e+212) {
tmp = a + ((b * (y + -2.0)) + (t * (b - a)));
} else if (a <= 1e-16) {
tmp = (((y + t) - 2.0) * b) + ((x + t_1) - (t * a));
} else {
tmp = ((b * (y + (t + -2.0))) + t_1) + (a * (1.0 - t));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = z * (1.0 - y) tmp = 0 if a <= -3.8e+212: tmp = a + ((b * (y + -2.0)) + (t * (b - a))) elif a <= 1e-16: tmp = (((y + t) - 2.0) * b) + ((x + t_1) - (t * a)) else: tmp = ((b * (y + (t + -2.0))) + t_1) + (a * (1.0 - t)) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(z * Float64(1.0 - y)) tmp = 0.0 if (a <= -3.8e+212) tmp = Float64(a + Float64(Float64(b * Float64(y + -2.0)) + Float64(t * Float64(b - a)))); elseif (a <= 1e-16) tmp = Float64(Float64(Float64(Float64(y + t) - 2.0) * b) + Float64(Float64(x + t_1) - Float64(t * a))); else tmp = Float64(Float64(Float64(b * Float64(y + Float64(t + -2.0))) + t_1) + Float64(a * Float64(1.0 - t))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = z * (1.0 - y); tmp = 0.0; if (a <= -3.8e+212) tmp = a + ((b * (y + -2.0)) + (t * (b - a))); elseif (a <= 1e-16) tmp = (((y + t) - 2.0) * b) + ((x + t_1) - (t * a)); else tmp = ((b * (y + (t + -2.0))) + t_1) + (a * (1.0 - t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -3.8e+212], N[(a + N[(N[(b * N[(y + -2.0), $MachinePrecision]), $MachinePrecision] + N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1e-16], N[(N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision] + N[(N[(x + t$95$1), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b * N[(y + N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$1), $MachinePrecision] + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \left(1 - y\right)\\
\mathbf{if}\;a \leq -3.8 \cdot 10^{+212}:\\
\;\;\;\;a + \left(b \cdot \left(y + -2\right) + t \cdot \left(b - a\right)\right)\\
\mathbf{elif}\;a \leq 10^{-16}:\\
\;\;\;\;\left(\left(y + t\right) - 2\right) \cdot b + \left(\left(x + t\_1\right) - t \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;\left(b \cdot \left(y + \left(t + -2\right)\right) + t\_1\right) + a \cdot \left(1 - t\right)\\
\end{array}
\end{array}
if a < -3.79999999999999988e212Initial program 90.0%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified95.0%
Taylor expanded in a around inf
Simplified95.3%
if -3.79999999999999988e212 < a < 9.9999999999999998e-17Initial program 97.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6493.8%
Simplified93.8%
if 9.9999999999999998e-17 < a Initial program 97.0%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval97.0%
Applied egg-rr97.0%
Taylor expanded in x around 0
--lowering--.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+r-N/A
+-commutativeN/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6491.7%
Simplified91.7%
Final simplification93.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* z (- 1.0 y))))
(if (<= a -2.25e+210)
(+ a (+ (* b (+ y -2.0)) (* t (- b a))))
(if (<= a 1.7e+69)
(+ (* (- (+ y t) 2.0) b) (- (+ x t_1) (* t a)))
(+ x (+ t_1 (* a (- 1.0 t))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z * (1.0 - y);
double tmp;
if (a <= -2.25e+210) {
tmp = a + ((b * (y + -2.0)) + (t * (b - a)));
} else if (a <= 1.7e+69) {
tmp = (((y + t) - 2.0) * b) + ((x + t_1) - (t * a));
} else {
tmp = x + (t_1 + (a * (1.0 - t)));
}
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 * (1.0d0 - y)
if (a <= (-2.25d+210)) then
tmp = a + ((b * (y + (-2.0d0))) + (t * (b - a)))
else if (a <= 1.7d+69) then
tmp = (((y + t) - 2.0d0) * b) + ((x + t_1) - (t * a))
else
tmp = x + (t_1 + (a * (1.0d0 - t)))
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 * (1.0 - y);
double tmp;
if (a <= -2.25e+210) {
tmp = a + ((b * (y + -2.0)) + (t * (b - a)));
} else if (a <= 1.7e+69) {
tmp = (((y + t) - 2.0) * b) + ((x + t_1) - (t * a));
} else {
tmp = x + (t_1 + (a * (1.0 - t)));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = z * (1.0 - y) tmp = 0 if a <= -2.25e+210: tmp = a + ((b * (y + -2.0)) + (t * (b - a))) elif a <= 1.7e+69: tmp = (((y + t) - 2.0) * b) + ((x + t_1) - (t * a)) else: tmp = x + (t_1 + (a * (1.0 - t))) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(z * Float64(1.0 - y)) tmp = 0.0 if (a <= -2.25e+210) tmp = Float64(a + Float64(Float64(b * Float64(y + -2.0)) + Float64(t * Float64(b - a)))); elseif (a <= 1.7e+69) tmp = Float64(Float64(Float64(Float64(y + t) - 2.0) * b) + Float64(Float64(x + t_1) - Float64(t * a))); else tmp = Float64(x + Float64(t_1 + Float64(a * Float64(1.0 - t)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = z * (1.0 - y); tmp = 0.0; if (a <= -2.25e+210) tmp = a + ((b * (y + -2.0)) + (t * (b - a))); elseif (a <= 1.7e+69) tmp = (((y + t) - 2.0) * b) + ((x + t_1) - (t * a)); else tmp = x + (t_1 + (a * (1.0 - t))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -2.25e+210], N[(a + N[(N[(b * N[(y + -2.0), $MachinePrecision]), $MachinePrecision] + N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.7e+69], N[(N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision] + N[(N[(x + t$95$1), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(t$95$1 + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \left(1 - y\right)\\
\mathbf{if}\;a \leq -2.25 \cdot 10^{+210}:\\
\;\;\;\;a + \left(b \cdot \left(y + -2\right) + t \cdot \left(b - a\right)\right)\\
\mathbf{elif}\;a \leq 1.7 \cdot 10^{+69}:\\
\;\;\;\;\left(\left(y + t\right) - 2\right) \cdot b + \left(\left(x + t\_1\right) - t \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;x + \left(t\_1 + a \cdot \left(1 - t\right)\right)\\
\end{array}
\end{array}
if a < -2.25000000000000002e210Initial program 90.0%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified95.0%
Taylor expanded in a around inf
Simplified95.3%
if -2.25000000000000002e210 < a < 1.69999999999999993e69Initial program 97.3%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6493.9%
Simplified93.9%
if 1.69999999999999993e69 < a Initial program 96.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
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/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
Simplified88.7%
Final simplification93.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -3.1e+114)
(+ x (* (- (+ y t) 2.0) b))
(if (<= b -1.5e-59)
(+ z (+ (* y (- b z)) (* b (+ t -2.0))))
(if (<= b 6e+71)
(+ x (+ (* z (- 1.0 y)) (* a (- 1.0 t))))
(+ a (+ (* b (+ y -2.0)) (* t (- b a))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -3.1e+114) {
tmp = x + (((y + t) - 2.0) * b);
} else if (b <= -1.5e-59) {
tmp = z + ((y * (b - z)) + (b * (t + -2.0)));
} else if (b <= 6e+71) {
tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t)));
} else {
tmp = a + ((b * (y + -2.0)) + (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 (b <= (-3.1d+114)) then
tmp = x + (((y + t) - 2.0d0) * b)
else if (b <= (-1.5d-59)) then
tmp = z + ((y * (b - z)) + (b * (t + (-2.0d0))))
else if (b <= 6d+71) then
tmp = x + ((z * (1.0d0 - y)) + (a * (1.0d0 - t)))
else
tmp = a + ((b * (y + (-2.0d0))) + (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 (b <= -3.1e+114) {
tmp = x + (((y + t) - 2.0) * b);
} else if (b <= -1.5e-59) {
tmp = z + ((y * (b - z)) + (b * (t + -2.0)));
} else if (b <= 6e+71) {
tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t)));
} else {
tmp = a + ((b * (y + -2.0)) + (t * (b - a)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -3.1e+114: tmp = x + (((y + t) - 2.0) * b) elif b <= -1.5e-59: tmp = z + ((y * (b - z)) + (b * (t + -2.0))) elif b <= 6e+71: tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t))) else: tmp = a + ((b * (y + -2.0)) + (t * (b - a))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -3.1e+114) tmp = Float64(x + Float64(Float64(Float64(y + t) - 2.0) * b)); elseif (b <= -1.5e-59) tmp = Float64(z + Float64(Float64(y * Float64(b - z)) + Float64(b * Float64(t + -2.0)))); elseif (b <= 6e+71) tmp = Float64(x + Float64(Float64(z * Float64(1.0 - y)) + Float64(a * Float64(1.0 - t)))); else tmp = Float64(a + Float64(Float64(b * Float64(y + -2.0)) + Float64(t * Float64(b - a)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -3.1e+114) tmp = x + (((y + t) - 2.0) * b); elseif (b <= -1.5e-59) tmp = z + ((y * (b - z)) + (b * (t + -2.0))); elseif (b <= 6e+71) tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t))); else tmp = a + ((b * (y + -2.0)) + (t * (b - a))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -3.1e+114], N[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -1.5e-59], N[(z + N[(N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision] + N[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 6e+71], N[(x + N[(N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision] + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a + N[(N[(b * N[(y + -2.0), $MachinePrecision]), $MachinePrecision] + N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.1 \cdot 10^{+114}:\\
\;\;\;\;x + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{elif}\;b \leq -1.5 \cdot 10^{-59}:\\
\;\;\;\;z + \left(y \cdot \left(b - z\right) + b \cdot \left(t + -2\right)\right)\\
\mathbf{elif}\;b \leq 6 \cdot 10^{+71}:\\
\;\;\;\;x + \left(z \cdot \left(1 - y\right) + a \cdot \left(1 - t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;a + \left(b \cdot \left(y + -2\right) + t \cdot \left(b - a\right)\right)\\
\end{array}
\end{array}
if b < -3.1e114Initial program 90.2%
Taylor expanded in x around inf
Simplified90.4%
if -3.1e114 < b < -1.5e-59Initial program 97.1%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6486.2%
Simplified86.2%
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.1%
Simplified89.1%
if -1.5e-59 < b < 6.00000000000000025e71Initial program 100.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
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/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
Simplified93.8%
if 6.00000000000000025e71 < b Initial program 91.6%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified93.7%
Taylor expanded in a around inf
Simplified85.4%
Final simplification91.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* a (- 1.0 t))) (t_2 (* (- (+ y t) 2.0) b)))
(if (<= b -2.4e+114)
(+ x t_2)
(if (<= b -1.5e-59)
(+ z (+ (* y (- b z)) (* b (+ t -2.0))))
(if (<= b 4.5e+71) (+ x (+ (* z (- 1.0 y)) t_1)) (+ t_2 t_1))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = a * (1.0 - t);
double t_2 = ((y + t) - 2.0) * b;
double tmp;
if (b <= -2.4e+114) {
tmp = x + t_2;
} else if (b <= -1.5e-59) {
tmp = z + ((y * (b - z)) + (b * (t + -2.0)));
} else if (b <= 4.5e+71) {
tmp = x + ((z * (1.0 - y)) + t_1);
} else {
tmp = t_2 + 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) :: t_2
real(8) :: tmp
t_1 = a * (1.0d0 - t)
t_2 = ((y + t) - 2.0d0) * b
if (b <= (-2.4d+114)) then
tmp = x + t_2
else if (b <= (-1.5d-59)) then
tmp = z + ((y * (b - z)) + (b * (t + (-2.0d0))))
else if (b <= 4.5d+71) then
tmp = x + ((z * (1.0d0 - y)) + t_1)
else
tmp = t_2 + 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 t_2 = ((y + t) - 2.0) * b;
double tmp;
if (b <= -2.4e+114) {
tmp = x + t_2;
} else if (b <= -1.5e-59) {
tmp = z + ((y * (b - z)) + (b * (t + -2.0)));
} else if (b <= 4.5e+71) {
tmp = x + ((z * (1.0 - y)) + t_1);
} else {
tmp = t_2 + t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = a * (1.0 - t) t_2 = ((y + t) - 2.0) * b tmp = 0 if b <= -2.4e+114: tmp = x + t_2 elif b <= -1.5e-59: tmp = z + ((y * (b - z)) + (b * (t + -2.0))) elif b <= 4.5e+71: tmp = x + ((z * (1.0 - y)) + t_1) else: tmp = t_2 + t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(a * Float64(1.0 - t)) t_2 = Float64(Float64(Float64(y + t) - 2.0) * b) tmp = 0.0 if (b <= -2.4e+114) tmp = Float64(x + t_2); elseif (b <= -1.5e-59) tmp = Float64(z + Float64(Float64(y * Float64(b - z)) + Float64(b * Float64(t + -2.0)))); elseif (b <= 4.5e+71) tmp = Float64(x + Float64(Float64(z * Float64(1.0 - y)) + t_1)); else tmp = Float64(t_2 + t_1); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = a * (1.0 - t); t_2 = ((y + t) - 2.0) * b; tmp = 0.0; if (b <= -2.4e+114) tmp = x + t_2; elseif (b <= -1.5e-59) tmp = z + ((y * (b - z)) + (b * (t + -2.0))); elseif (b <= 4.5e+71) tmp = x + ((z * (1.0 - y)) + t_1); else tmp = t_2 + 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]}, Block[{t$95$2 = N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[b, -2.4e+114], N[(x + t$95$2), $MachinePrecision], If[LessEqual[b, -1.5e-59], N[(z + N[(N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision] + N[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 4.5e+71], N[(x + N[(N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision] + t$95$1), $MachinePrecision]), $MachinePrecision], N[(t$95$2 + t$95$1), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \left(1 - t\right)\\
t_2 := \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -2.4 \cdot 10^{+114}:\\
\;\;\;\;x + t\_2\\
\mathbf{elif}\;b \leq -1.5 \cdot 10^{-59}:\\
\;\;\;\;z + \left(y \cdot \left(b - z\right) + b \cdot \left(t + -2\right)\right)\\
\mathbf{elif}\;b \leq 4.5 \cdot 10^{+71}:\\
\;\;\;\;x + \left(z \cdot \left(1 - y\right) + t\_1\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2 + t\_1\\
\end{array}
\end{array}
if b < -2.4e114Initial program 90.2%
Taylor expanded in x around inf
Simplified90.4%
if -2.4e114 < b < -1.5e-59Initial program 97.1%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6486.2%
Simplified86.2%
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.1%
Simplified89.1%
if -1.5e-59 < b < 4.50000000000000043e71Initial program 100.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
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/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
Simplified93.8%
if 4.50000000000000043e71 < b Initial program 91.6%
Taylor expanded in a around inf
sub-negN/A
metadata-evalN/A
neg-mul-1N/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--.f6485.4%
Simplified85.4%
Final simplification91.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* (- (+ y t) 2.0) b))))
(if (<= b -3.1e+114)
t_1
(if (<= b -1.5e-59)
(+ z (+ (* y (- b z)) (* b (+ t -2.0))))
(if (<= b 1.3e+194) (+ x (+ (* z (- 1.0 y)) (* a (- 1.0 t)))) 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 <= -3.1e+114) {
tmp = t_1;
} else if (b <= -1.5e-59) {
tmp = z + ((y * (b - z)) + (b * (t + -2.0)));
} else if (b <= 1.3e+194) {
tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t)));
} 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 <= (-3.1d+114)) then
tmp = t_1
else if (b <= (-1.5d-59)) then
tmp = z + ((y * (b - z)) + (b * (t + (-2.0d0))))
else if (b <= 1.3d+194) then
tmp = x + ((z * (1.0d0 - y)) + (a * (1.0d0 - t)))
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 <= -3.1e+114) {
tmp = t_1;
} else if (b <= -1.5e-59) {
tmp = z + ((y * (b - z)) + (b * (t + -2.0)));
} else if (b <= 1.3e+194) {
tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t)));
} 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 <= -3.1e+114: tmp = t_1 elif b <= -1.5e-59: tmp = z + ((y * (b - z)) + (b * (t + -2.0))) elif b <= 1.3e+194: tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t))) 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 <= -3.1e+114) tmp = t_1; elseif (b <= -1.5e-59) tmp = Float64(z + Float64(Float64(y * Float64(b - z)) + Float64(b * Float64(t + -2.0)))); elseif (b <= 1.3e+194) tmp = Float64(x + Float64(Float64(z * Float64(1.0 - y)) + Float64(a * Float64(1.0 - t)))); 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 <= -3.1e+114) tmp = t_1; elseif (b <= -1.5e-59) tmp = z + ((y * (b - z)) + (b * (t + -2.0))); elseif (b <= 1.3e+194) tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t))); 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, -3.1e+114], t$95$1, If[LessEqual[b, -1.5e-59], N[(z + N[(N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision] + N[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.3e+194], N[(x + N[(N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision] + N[(a * N[(1.0 - t), $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 -3.1 \cdot 10^{+114}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -1.5 \cdot 10^{-59}:\\
\;\;\;\;z + \left(y \cdot \left(b - z\right) + b \cdot \left(t + -2\right)\right)\\
\mathbf{elif}\;b \leq 1.3 \cdot 10^{+194}:\\
\;\;\;\;x + \left(z \cdot \left(1 - y\right) + a \cdot \left(1 - t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -3.1e114 or 1.2999999999999999e194 < b Initial program 91.2%
Taylor expanded in x around inf
Simplified87.5%
if -3.1e114 < b < -1.5e-59Initial program 97.1%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6486.2%
Simplified86.2%
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.1%
Simplified89.1%
if -1.5e-59 < b < 1.2999999999999999e194Initial program 98.7%
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
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/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
Simplified89.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* b (+ t (+ y -2.0)))))
(if (<= b -1.5e-59)
t_1
(if (<= b -1.35e-254)
(+ x (* a (- 1.0 t)))
(if (<= b 6.4e+17)
(+ x (* z (- 1.0 y)))
(if (<= b 1.15e+135) (* t (- b a)) t_1))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = b * (t + (y + -2.0));
double tmp;
if (b <= -1.5e-59) {
tmp = t_1;
} else if (b <= -1.35e-254) {
tmp = x + (a * (1.0 - t));
} else if (b <= 6.4e+17) {
tmp = x + (z * (1.0 - y));
} else if (b <= 1.15e+135) {
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 = b * (t + (y + (-2.0d0)))
if (b <= (-1.5d-59)) then
tmp = t_1
else if (b <= (-1.35d-254)) then
tmp = x + (a * (1.0d0 - t))
else if (b <= 6.4d+17) then
tmp = x + (z * (1.0d0 - y))
else if (b <= 1.15d+135) 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 = b * (t + (y + -2.0));
double tmp;
if (b <= -1.5e-59) {
tmp = t_1;
} else if (b <= -1.35e-254) {
tmp = x + (a * (1.0 - t));
} else if (b <= 6.4e+17) {
tmp = x + (z * (1.0 - y));
} else if (b <= 1.15e+135) {
tmp = t * (b - a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = b * (t + (y + -2.0)) tmp = 0 if b <= -1.5e-59: tmp = t_1 elif b <= -1.35e-254: tmp = x + (a * (1.0 - t)) elif b <= 6.4e+17: tmp = x + (z * (1.0 - y)) elif b <= 1.15e+135: tmp = t * (b - a) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(b * Float64(t + Float64(y + -2.0))) tmp = 0.0 if (b <= -1.5e-59) tmp = t_1; elseif (b <= -1.35e-254) tmp = Float64(x + Float64(a * Float64(1.0 - t))); elseif (b <= 6.4e+17) tmp = Float64(x + Float64(z * Float64(1.0 - y))); elseif (b <= 1.15e+135) 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 = b * (t + (y + -2.0)); tmp = 0.0; if (b <= -1.5e-59) tmp = t_1; elseif (b <= -1.35e-254) tmp = x + (a * (1.0 - t)); elseif (b <= 6.4e+17) tmp = x + (z * (1.0 - y)); elseif (b <= 1.15e+135) 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[(b * N[(t + N[(y + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.5e-59], t$95$1, If[LessEqual[b, -1.35e-254], N[(x + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 6.4e+17], N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.15e+135], N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := b \cdot \left(t + \left(y + -2\right)\right)\\
\mathbf{if}\;b \leq -1.5 \cdot 10^{-59}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -1.35 \cdot 10^{-254}:\\
\;\;\;\;x + a \cdot \left(1 - t\right)\\
\mathbf{elif}\;b \leq 6.4 \cdot 10^{+17}:\\
\;\;\;\;x + z \cdot \left(1 - y\right)\\
\mathbf{elif}\;b \leq 1.15 \cdot 10^{+135}:\\
\;\;\;\;t \cdot \left(b - a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.5e-59 or 1.1500000000000001e135 < b Initial program 92.1%
Taylor expanded in b around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval68.9%
Simplified68.9%
if -1.5e-59 < b < -1.35000000000000003e-254Initial program 100.0%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval100.0%
Applied egg-rr100.0%
Taylor expanded in x around inf
Simplified67.5%
if -1.35000000000000003e-254 < b < 6.4e17Initial program 100.0%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified100.0%
Taylor expanded in b around 0
associate-+r+N/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
*-rgt-identityN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f6495.1%
Simplified95.1%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
cancel-sign-sub-invN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6464.6%
Simplified64.6%
if 6.4e17 < b < 1.1500000000000001e135Initial program 99.9%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6466.7%
Simplified66.7%
Final simplification67.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* (- (+ y t) 2.0) b))))
(if (<= b -3.2e-35)
t_1
(if (<= b 3.8e-48)
(- (+ x a) (* z (+ y -1.0)))
(if (<= b 3.5e+196) (+ (* y (- b z)) (* a (- 1.0 t))) 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 <= -3.2e-35) {
tmp = t_1;
} else if (b <= 3.8e-48) {
tmp = (x + a) - (z * (y + -1.0));
} else if (b <= 3.5e+196) {
tmp = (y * (b - z)) + (a * (1.0 - t));
} 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 <= (-3.2d-35)) then
tmp = t_1
else if (b <= 3.8d-48) then
tmp = (x + a) - (z * (y + (-1.0d0)))
else if (b <= 3.5d+196) then
tmp = (y * (b - z)) + (a * (1.0d0 - t))
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 <= -3.2e-35) {
tmp = t_1;
} else if (b <= 3.8e-48) {
tmp = (x + a) - (z * (y + -1.0));
} else if (b <= 3.5e+196) {
tmp = (y * (b - z)) + (a * (1.0 - t));
} 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 <= -3.2e-35: tmp = t_1 elif b <= 3.8e-48: tmp = (x + a) - (z * (y + -1.0)) elif b <= 3.5e+196: tmp = (y * (b - z)) + (a * (1.0 - t)) 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 <= -3.2e-35) tmp = t_1; elseif (b <= 3.8e-48) tmp = Float64(Float64(x + a) - Float64(z * Float64(y + -1.0))); elseif (b <= 3.5e+196) tmp = Float64(Float64(y * Float64(b - z)) + Float64(a * Float64(1.0 - t))); 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 <= -3.2e-35) tmp = t_1; elseif (b <= 3.8e-48) tmp = (x + a) - (z * (y + -1.0)); elseif (b <= 3.5e+196) tmp = (y * (b - z)) + (a * (1.0 - t)); 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, -3.2e-35], t$95$1, If[LessEqual[b, 3.8e-48], N[(N[(x + a), $MachinePrecision] - N[(z * N[(y + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 3.5e+196], N[(N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision] + N[(a * N[(1.0 - t), $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 -3.2 \cdot 10^{-35}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 3.8 \cdot 10^{-48}:\\
\;\;\;\;\left(x + a\right) - z \cdot \left(y + -1\right)\\
\mathbf{elif}\;b \leq 3.5 \cdot 10^{+196}:\\
\;\;\;\;y \cdot \left(b - z\right) + a \cdot \left(1 - t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -3.1999999999999998e-35 or 3.4999999999999998e196 < b Initial program 92.8%
Taylor expanded in x around inf
Simplified78.8%
if -3.1999999999999998e-35 < b < 3.80000000000000002e-48Initial program 100.0%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified100.0%
Taylor expanded in b around 0
associate-+r+N/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
*-rgt-identityN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f6495.6%
Simplified95.6%
Taylor expanded in t around 0
+-commutativeN/A
+-lowering-+.f6477.3%
Simplified77.3%
if 3.80000000000000002e-48 < b < 3.4999999999999998e196Initial program 95.2%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval95.2%
Applied egg-rr95.2%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6464.7%
Simplified64.7%
Final simplification75.7%
(FPCore (x y z t a b) :precision binary64 (if (<= y -6.5e+134) (+ (* y (- b z)) (* a (- 1.0 t))) (+ (+ (* b (+ y -2.0)) (* t (- b a))) (- (+ x a) (* z (+ y -1.0))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -6.5e+134) {
tmp = (y * (b - z)) + (a * (1.0 - t));
} else {
tmp = ((b * (y + -2.0)) + (t * (b - a))) + ((x + a) - (z * (y + -1.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 (y <= (-6.5d+134)) then
tmp = (y * (b - z)) + (a * (1.0d0 - t))
else
tmp = ((b * (y + (-2.0d0))) + (t * (b - a))) + ((x + a) - (z * (y + (-1.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 (y <= -6.5e+134) {
tmp = (y * (b - z)) + (a * (1.0 - t));
} else {
tmp = ((b * (y + -2.0)) + (t * (b - a))) + ((x + a) - (z * (y + -1.0)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -6.5e+134: tmp = (y * (b - z)) + (a * (1.0 - t)) else: tmp = ((b * (y + -2.0)) + (t * (b - a))) + ((x + a) - (z * (y + -1.0))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -6.5e+134) tmp = Float64(Float64(y * Float64(b - z)) + Float64(a * Float64(1.0 - t))); else tmp = Float64(Float64(Float64(b * Float64(y + -2.0)) + Float64(t * Float64(b - a))) + Float64(Float64(x + a) - Float64(z * Float64(y + -1.0)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -6.5e+134) tmp = (y * (b - z)) + (a * (1.0 - t)); else tmp = ((b * (y + -2.0)) + (t * (b - a))) + ((x + a) - (z * (y + -1.0))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -6.5e+134], N[(N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision] + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b * N[(y + -2.0), $MachinePrecision]), $MachinePrecision] + N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(x + a), $MachinePrecision] - N[(z * N[(y + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6.5 \cdot 10^{+134}:\\
\;\;\;\;y \cdot \left(b - z\right) + a \cdot \left(1 - t\right)\\
\mathbf{else}:\\
\;\;\;\;\left(b \cdot \left(y + -2\right) + t \cdot \left(b - a\right)\right) + \left(\left(x + a\right) - z \cdot \left(y + -1\right)\right)\\
\end{array}
\end{array}
if y < -6.5e134Initial program 83.3%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval83.3%
Applied egg-rr83.3%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6494.4%
Simplified94.4%
if -6.5e134 < y Initial program 98.6%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified99.5%
Final simplification98.8%
(FPCore (x y z t a b) :precision binary64 (if (<= b -1.95e+263) (+ x (* (- (+ y t) 2.0) b)) (+ (+ x (* b (+ t -2.0))) (+ (* y (- b z)) (+ z (* a (- 1.0 t)))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.95e+263) {
tmp = x + (((y + t) - 2.0) * b);
} else {
tmp = (x + (b * (t + -2.0))) + ((y * (b - z)) + (z + (a * (1.0 - t))));
}
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 <= (-1.95d+263)) then
tmp = x + (((y + t) - 2.0d0) * b)
else
tmp = (x + (b * (t + (-2.0d0)))) + ((y * (b - z)) + (z + (a * (1.0d0 - t))))
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 <= -1.95e+263) {
tmp = x + (((y + t) - 2.0) * b);
} else {
tmp = (x + (b * (t + -2.0))) + ((y * (b - z)) + (z + (a * (1.0 - t))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -1.95e+263: tmp = x + (((y + t) - 2.0) * b) else: tmp = (x + (b * (t + -2.0))) + ((y * (b - z)) + (z + (a * (1.0 - t)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -1.95e+263) tmp = Float64(x + Float64(Float64(Float64(y + t) - 2.0) * b)); else tmp = Float64(Float64(x + Float64(b * Float64(t + -2.0))) + Float64(Float64(y * Float64(b - z)) + Float64(z + Float64(a * Float64(1.0 - t))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -1.95e+263) tmp = x + (((y + t) - 2.0) * b); else tmp = (x + (b * (t + -2.0))) + ((y * (b - z)) + (z + (a * (1.0 - t)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -1.95e+263], N[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision], 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]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.95 \cdot 10^{+263}:\\
\;\;\;\;x + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{else}:\\
\;\;\;\;\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)\\
\end{array}
\end{array}
if b < -1.95000000000000014e263Initial program 76.9%
Taylor expanded in x around inf
Simplified100.0%
if -1.95000000000000014e263 < b Initial program 97.5%
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
Simplified98.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* (- (+ y t) 2.0) b))))
(if (<= b -3.3e-35)
t_1
(if (<= b 1.85e+18)
(- (+ x a) (* z (+ y -1.0)))
(if (<= b 1.15e+135) (* t (- b 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 <= -3.3e-35) {
tmp = t_1;
} else if (b <= 1.85e+18) {
tmp = (x + a) - (z * (y + -1.0));
} else if (b <= 1.15e+135) {
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 = x + (((y + t) - 2.0d0) * b)
if (b <= (-3.3d-35)) then
tmp = t_1
else if (b <= 1.85d+18) then
tmp = (x + a) - (z * (y + (-1.0d0)))
else if (b <= 1.15d+135) 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 = x + (((y + t) - 2.0) * b);
double tmp;
if (b <= -3.3e-35) {
tmp = t_1;
} else if (b <= 1.85e+18) {
tmp = (x + a) - (z * (y + -1.0));
} else if (b <= 1.15e+135) {
tmp = t * (b - 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 <= -3.3e-35: tmp = t_1 elif b <= 1.85e+18: tmp = (x + a) - (z * (y + -1.0)) elif b <= 1.15e+135: tmp = t * (b - 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 <= -3.3e-35) tmp = t_1; elseif (b <= 1.85e+18) tmp = Float64(Float64(x + a) - Float64(z * Float64(y + -1.0))); elseif (b <= 1.15e+135) 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 = x + (((y + t) - 2.0) * b); tmp = 0.0; if (b <= -3.3e-35) tmp = t_1; elseif (b <= 1.85e+18) tmp = (x + a) - (z * (y + -1.0)); elseif (b <= 1.15e+135) 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[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -3.3e-35], t$95$1, If[LessEqual[b, 1.85e+18], N[(N[(x + a), $MachinePrecision] - N[(z * N[(y + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.15e+135], N[(t * N[(b - a), $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 -3.3 \cdot 10^{-35}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.85 \cdot 10^{+18}:\\
\;\;\;\;\left(x + a\right) - z \cdot \left(y + -1\right)\\
\mathbf{elif}\;b \leq 1.15 \cdot 10^{+135}:\\
\;\;\;\;t \cdot \left(b - a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -3.3e-35 or 1.1500000000000001e135 < b Initial program 91.6%
Taylor expanded in x around inf
Simplified76.2%
if -3.3e-35 < b < 1.85e18Initial program 100.0%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified100.0%
Taylor expanded in b around 0
associate-+r+N/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
*-rgt-identityN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f6493.8%
Simplified93.8%
Taylor expanded in t around 0
+-commutativeN/A
+-lowering-+.f6476.4%
Simplified76.4%
if 1.85e18 < b < 1.1500000000000001e135Initial program 99.9%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6466.7%
Simplified66.7%
Final simplification75.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* (- (+ y t) 2.0) b))))
(if (<= b -3.4e-35)
t_1
(if (<= b 1.3e+194) (+ x (+ (* z (- 1.0 y)) (* a (- 1.0 t)))) 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 <= -3.4e-35) {
tmp = t_1;
} else if (b <= 1.3e+194) {
tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t)));
} 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 <= (-3.4d-35)) then
tmp = t_1
else if (b <= 1.3d+194) then
tmp = x + ((z * (1.0d0 - y)) + (a * (1.0d0 - t)))
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 <= -3.4e-35) {
tmp = t_1;
} else if (b <= 1.3e+194) {
tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t)));
} 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 <= -3.4e-35: tmp = t_1 elif b <= 1.3e+194: tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t))) 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 <= -3.4e-35) tmp = t_1; elseif (b <= 1.3e+194) tmp = Float64(x + Float64(Float64(z * Float64(1.0 - y)) + Float64(a * Float64(1.0 - t)))); 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 <= -3.4e-35) tmp = t_1; elseif (b <= 1.3e+194) tmp = x + ((z * (1.0 - y)) + (a * (1.0 - t))); 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, -3.4e-35], t$95$1, If[LessEqual[b, 1.3e+194], N[(x + N[(N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision] + N[(a * N[(1.0 - t), $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 -3.4 \cdot 10^{-35}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.3 \cdot 10^{+194}:\\
\;\;\;\;x + \left(z \cdot \left(1 - y\right) + a \cdot \left(1 - t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -3.4000000000000003e-35 or 1.2999999999999999e194 < b Initial program 92.8%
Taylor expanded in x around inf
Simplified78.8%
if -3.4000000000000003e-35 < b < 1.2999999999999999e194Initial program 98.7%
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
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/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
Simplified88.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* z (- 1.0 y)))))
(if (<= z -1.38e+107)
t_1
(if (<= z 4.4e-55)
(+ x (* (- (+ y t) 2.0) b))
(if (<= z 4.8e+173) (+ x (* a (- 1.0 t))) 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.38e+107) {
tmp = t_1;
} else if (z <= 4.4e-55) {
tmp = x + (((y + t) - 2.0) * b);
} else if (z <= 4.8e+173) {
tmp = x + (a * (1.0 - t));
} 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.38d+107)) then
tmp = t_1
else if (z <= 4.4d-55) then
tmp = x + (((y + t) - 2.0d0) * b)
else if (z <= 4.8d+173) then
tmp = x + (a * (1.0d0 - t))
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.38e+107) {
tmp = t_1;
} else if (z <= 4.4e-55) {
tmp = x + (((y + t) - 2.0) * b);
} else if (z <= 4.8e+173) {
tmp = x + (a * (1.0 - t));
} 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.38e+107: tmp = t_1 elif z <= 4.4e-55: tmp = x + (((y + t) - 2.0) * b) elif z <= 4.8e+173: tmp = x + (a * (1.0 - t)) 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.38e+107) tmp = t_1; elseif (z <= 4.4e-55) tmp = Float64(x + Float64(Float64(Float64(y + t) - 2.0) * b)); elseif (z <= 4.8e+173) tmp = Float64(x + Float64(a * Float64(1.0 - t))); 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.38e+107) tmp = t_1; elseif (z <= 4.4e-55) tmp = x + (((y + t) - 2.0) * b); elseif (z <= 4.8e+173) tmp = x + (a * (1.0 - t)); 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.38e+107], t$95$1, If[LessEqual[z, 4.4e-55], N[(x + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 4.8e+173], N[(x + N[(a * N[(1.0 - t), $MachinePrecision]), $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.38 \cdot 10^{+107}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 4.4 \cdot 10^{-55}:\\
\;\;\;\;x + \left(\left(y + t\right) - 2\right) \cdot b\\
\mathbf{elif}\;z \leq 4.8 \cdot 10^{+173}:\\
\;\;\;\;x + a \cdot \left(1 - t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.38e107 or 4.7999999999999998e173 < z Initial program 89.9%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified91.4%
Taylor expanded in b around 0
associate-+r+N/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
*-rgt-identityN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f6485.3%
Simplified85.3%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
cancel-sign-sub-invN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6476.1%
Simplified76.1%
if -1.38e107 < z < 4.3999999999999999e-55Initial program 99.3%
Taylor expanded in x around inf
Simplified67.2%
if 4.3999999999999999e-55 < z < 4.7999999999999998e173Initial program 97.7%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval97.7%
Applied egg-rr97.7%
Taylor expanded in x around inf
Simplified61.7%
Final simplification68.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* b (+ t (+ y -2.0)))))
(if (<= b -1.65e-37)
t_1
(if (<= b 1.32e+14)
(+ x (* z (- 1.0 y)))
(if (<= b 1.15e+135) (* t (- b a)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = b * (t + (y + -2.0));
double tmp;
if (b <= -1.65e-37) {
tmp = t_1;
} else if (b <= 1.32e+14) {
tmp = x + (z * (1.0 - y));
} else if (b <= 1.15e+135) {
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 = b * (t + (y + (-2.0d0)))
if (b <= (-1.65d-37)) then
tmp = t_1
else if (b <= 1.32d+14) then
tmp = x + (z * (1.0d0 - y))
else if (b <= 1.15d+135) 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 = b * (t + (y + -2.0));
double tmp;
if (b <= -1.65e-37) {
tmp = t_1;
} else if (b <= 1.32e+14) {
tmp = x + (z * (1.0 - y));
} else if (b <= 1.15e+135) {
tmp = t * (b - a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = b * (t + (y + -2.0)) tmp = 0 if b <= -1.65e-37: tmp = t_1 elif b <= 1.32e+14: tmp = x + (z * (1.0 - y)) elif b <= 1.15e+135: tmp = t * (b - a) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(b * Float64(t + Float64(y + -2.0))) tmp = 0.0 if (b <= -1.65e-37) tmp = t_1; elseif (b <= 1.32e+14) tmp = Float64(x + Float64(z * Float64(1.0 - y))); elseif (b <= 1.15e+135) 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 = b * (t + (y + -2.0)); tmp = 0.0; if (b <= -1.65e-37) tmp = t_1; elseif (b <= 1.32e+14) tmp = x + (z * (1.0 - y)); elseif (b <= 1.15e+135) 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[(b * N[(t + N[(y + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.65e-37], t$95$1, If[LessEqual[b, 1.32e+14], N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.15e+135], N[(t * N[(b - a), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := b \cdot \left(t + \left(y + -2\right)\right)\\
\mathbf{if}\;b \leq -1.65 \cdot 10^{-37}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.32 \cdot 10^{+14}:\\
\;\;\;\;x + z \cdot \left(1 - y\right)\\
\mathbf{elif}\;b \leq 1.15 \cdot 10^{+135}:\\
\;\;\;\;t \cdot \left(b - a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.64999999999999991e-37 or 1.1500000000000001e135 < b Initial program 91.8%
Taylor expanded in b around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval69.7%
Simplified69.7%
if -1.64999999999999991e-37 < b < 1.32e14Initial program 100.0%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified100.0%
Taylor expanded in b around 0
associate-+r+N/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
*-rgt-identityN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f6494.4%
Simplified94.4%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
cancel-sign-sub-invN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6461.0%
Simplified61.0%
if 1.32e14 < b < 1.1500000000000001e135Initial program 99.9%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6466.7%
Simplified66.7%
Final simplification65.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* y (- b z))))
(if (<= y -33.0)
t_1
(if (<= y 2.9e-174)
(+ x z)
(if (<= y 245000000000.0) (* a (- 1.0 t)) 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 <= -33.0) {
tmp = t_1;
} else if (y <= 2.9e-174) {
tmp = x + z;
} else if (y <= 245000000000.0) {
tmp = a * (1.0 - t);
} 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 <= (-33.0d0)) then
tmp = t_1
else if (y <= 2.9d-174) then
tmp = x + z
else if (y <= 245000000000.0d0) then
tmp = a * (1.0d0 - t)
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 <= -33.0) {
tmp = t_1;
} else if (y <= 2.9e-174) {
tmp = x + z;
} else if (y <= 245000000000.0) {
tmp = a * (1.0 - t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y * (b - z) tmp = 0 if y <= -33.0: tmp = t_1 elif y <= 2.9e-174: tmp = x + z elif y <= 245000000000.0: tmp = a * (1.0 - t) 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 <= -33.0) tmp = t_1; elseif (y <= 2.9e-174) tmp = Float64(x + z); elseif (y <= 245000000000.0) tmp = Float64(a * Float64(1.0 - t)); 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 <= -33.0) tmp = t_1; elseif (y <= 2.9e-174) tmp = x + z; elseif (y <= 245000000000.0) tmp = a * (1.0 - t); 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, -33.0], t$95$1, If[LessEqual[y, 2.9e-174], N[(x + z), $MachinePrecision], If[LessEqual[y, 245000000000.0], N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(b - z\right)\\
\mathbf{if}\;y \leq -33:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2.9 \cdot 10^{-174}:\\
\;\;\;\;x + z\\
\mathbf{elif}\;y \leq 245000000000:\\
\;\;\;\;a \cdot \left(1 - t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -33 or 2.45e11 < y Initial program 95.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6467.5%
Simplified67.5%
if -33 < y < 2.9000000000000001e-174Initial program 96.8%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified100.0%
Taylor expanded in b around 0
associate-+r+N/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
*-rgt-identityN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f6467.8%
Simplified67.8%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
cancel-sign-sub-invN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6444.8%
Simplified44.8%
Taylor expanded in y around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f6443.5%
Simplified43.5%
if 2.9000000000000001e-174 < y < 2.45e11Initial program 100.0%
Taylor expanded in a around inf
sub-negN/A
metadata-evalN/A
neg-mul-1N/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--.f6443.2%
Simplified43.2%
Final simplification54.8%
(FPCore (x y z t a b) :precision binary64 (if (<= b -1.75e-60) (* y b) (if (<= b 3.6e+17) (+ x z) (if (<= b 1.6e+246) (* t b) (* y b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.75e-60) {
tmp = y * b;
} else if (b <= 3.6e+17) {
tmp = x + z;
} else if (b <= 1.6e+246) {
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 (b <= (-1.75d-60)) then
tmp = y * b
else if (b <= 3.6d+17) then
tmp = x + z
else if (b <= 1.6d+246) 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 (b <= -1.75e-60) {
tmp = y * b;
} else if (b <= 3.6e+17) {
tmp = x + z;
} else if (b <= 1.6e+246) {
tmp = t * b;
} else {
tmp = y * b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -1.75e-60: tmp = y * b elif b <= 3.6e+17: tmp = x + z elif b <= 1.6e+246: tmp = t * b else: tmp = y * b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -1.75e-60) tmp = Float64(y * b); elseif (b <= 3.6e+17) tmp = Float64(x + z); elseif (b <= 1.6e+246) 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 (b <= -1.75e-60) tmp = y * b; elseif (b <= 3.6e+17) tmp = x + z; elseif (b <= 1.6e+246) tmp = t * b; else tmp = y * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -1.75e-60], N[(y * b), $MachinePrecision], If[LessEqual[b, 3.6e+17], N[(x + z), $MachinePrecision], If[LessEqual[b, 1.6e+246], N[(t * b), $MachinePrecision], N[(y * b), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.75 \cdot 10^{-60}:\\
\;\;\;\;y \cdot b\\
\mathbf{elif}\;b \leq 3.6 \cdot 10^{+17}:\\
\;\;\;\;x + z\\
\mathbf{elif}\;b \leq 1.6 \cdot 10^{+246}:\\
\;\;\;\;t \cdot b\\
\mathbf{else}:\\
\;\;\;\;y \cdot b\\
\end{array}
\end{array}
if b < -1.74999999999999988e-60 or 1.60000000000000007e246 < b Initial program 92.2%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6449.7%
Simplified49.7%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f6438.1%
Simplified38.1%
if -1.74999999999999988e-60 < b < 3.6e17Initial program 100.0%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified100.0%
Taylor expanded in b around 0
associate-+r+N/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
*-rgt-identityN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f6495.1%
Simplified95.1%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
cancel-sign-sub-invN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6460.8%
Simplified60.8%
Taylor expanded in y around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f6440.7%
Simplified40.7%
if 3.6e17 < b < 1.60000000000000007e246Initial program 95.2%
Taylor expanded in x around inf
Simplified52.2%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6433.3%
Simplified33.3%
Final simplification38.5%
(FPCore (x y z t a b) :precision binary64 (if (<= b -9.2e-41) (* t b) (if (<= b -2.3e-143) z (if (<= b 2.5) x (* t b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -9.2e-41) {
tmp = t * b;
} else if (b <= -2.3e-143) {
tmp = z;
} else if (b <= 2.5) {
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 (b <= (-9.2d-41)) then
tmp = t * b
else if (b <= (-2.3d-143)) then
tmp = z
else if (b <= 2.5d0) 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 (b <= -9.2e-41) {
tmp = t * b;
} else if (b <= -2.3e-143) {
tmp = z;
} else if (b <= 2.5) {
tmp = x;
} else {
tmp = t * b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -9.2e-41: tmp = t * b elif b <= -2.3e-143: tmp = z elif b <= 2.5: tmp = x else: tmp = t * b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -9.2e-41) tmp = Float64(t * b); elseif (b <= -2.3e-143) tmp = z; elseif (b <= 2.5) 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 (b <= -9.2e-41) tmp = t * b; elseif (b <= -2.3e-143) tmp = z; elseif (b <= 2.5) tmp = x; else tmp = t * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -9.2e-41], N[(t * b), $MachinePrecision], If[LessEqual[b, -2.3e-143], z, If[LessEqual[b, 2.5], x, N[(t * b), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -9.2 \cdot 10^{-41}:\\
\;\;\;\;t \cdot b\\
\mathbf{elif}\;b \leq -2.3 \cdot 10^{-143}:\\
\;\;\;\;z\\
\mathbf{elif}\;b \leq 2.5:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;t \cdot b\\
\end{array}
\end{array}
if b < -9.20000000000000041e-41 or 2.5 < b Initial program 93.2%
Taylor expanded in x around inf
Simplified67.7%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6429.7%
Simplified29.7%
if -9.20000000000000041e-41 < b < -2.30000000000000011e-143Initial program 100.0%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6447.2%
Simplified47.2%
Taylor expanded in y around 0
Simplified30.9%
if -2.30000000000000011e-143 < b < 2.5Initial program 100.0%
Taylor expanded in x around inf
Simplified27.4%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (+ a (* y (- b z))))) (if (<= y -3.2e+14) t_1 (if (<= y 8500.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 = a + (y * (b - z));
double tmp;
if (y <= -3.2e+14) {
tmp = t_1;
} else if (y <= 8500.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 = a + (y * (b - z))
if (y <= (-3.2d+14)) then
tmp = t_1
else if (y <= 8500.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 = a + (y * (b - z));
double tmp;
if (y <= -3.2e+14) {
tmp = t_1;
} else if (y <= 8500.0) {
tmp = z + (b * (t + -2.0));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = a + (y * (b - z)) tmp = 0 if y <= -3.2e+14: tmp = t_1 elif y <= 8500.0: tmp = z + (b * (t + -2.0)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(a + Float64(y * Float64(b - z))) tmp = 0.0 if (y <= -3.2e+14) tmp = t_1; elseif (y <= 8500.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 = a + (y * (b - z)); tmp = 0.0; if (y <= -3.2e+14) tmp = t_1; elseif (y <= 8500.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[(a + N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.2e+14], t$95$1, If[LessEqual[y, 8500.0], N[(z + N[(b * N[(t + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a + y \cdot \left(b - z\right)\\
\mathbf{if}\;y \leq -3.2 \cdot 10^{+14}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 8500:\\
\;\;\;\;z + b \cdot \left(t + -2\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -3.2e14 or 8500 < y Initial program 94.9%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval94.9%
Applied egg-rr94.9%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6486.1%
Simplified86.1%
Taylor expanded in t around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6473.8%
Simplified73.8%
if -3.2e14 < y < 8500Initial program 97.8%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6454.0%
Simplified54.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6452.4%
Simplified52.4%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* t (- b a)))) (if (<= t -7e+27) t_1 (if (<= t 72000000.0) (+ a (* 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 <= -7e+27) {
tmp = t_1;
} else if (t <= 72000000.0) {
tmp = a + (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 <= (-7d+27)) then
tmp = t_1
else if (t <= 72000000.0d0) then
tmp = a + (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 <= -7e+27) {
tmp = t_1;
} else if (t <= 72000000.0) {
tmp = a + (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 <= -7e+27: tmp = t_1 elif t <= 72000000.0: tmp = a + (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 <= -7e+27) tmp = t_1; elseif (t <= 72000000.0) tmp = Float64(a + 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 <= -7e+27) tmp = t_1; elseif (t <= 72000000.0) tmp = a + (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, -7e+27], t$95$1, If[LessEqual[t, 72000000.0], N[(a + N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -7 \cdot 10^{+27}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 72000000:\\
\;\;\;\;a + y \cdot \left(b - z\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -7.0000000000000004e27 or 7.2e7 < t Initial program 92.8%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6466.8%
Simplified66.8%
if -7.0000000000000004e27 < t < 7.2e7Initial program 99.3%
+-commutativeN/A
associate-+r-N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.3%
Applied egg-rr99.3%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6458.3%
Simplified58.3%
Taylor expanded in t around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f6457.9%
Simplified57.9%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* t (- b a)))) (if (<= t -6.2e-6) t_1 (if (<= t 45000000.0) (+ x 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 <= -6.2e-6) {
tmp = t_1;
} else if (t <= 45000000.0) {
tmp = x + 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 <= (-6.2d-6)) then
tmp = t_1
else if (t <= 45000000.0d0) then
tmp = x + 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 <= -6.2e-6) {
tmp = t_1;
} else if (t <= 45000000.0) {
tmp = x + z;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = t * (b - a) tmp = 0 if t <= -6.2e-6: tmp = t_1 elif t <= 45000000.0: tmp = x + 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 <= -6.2e-6) tmp = t_1; elseif (t <= 45000000.0) tmp = Float64(x + 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 <= -6.2e-6) tmp = t_1; elseif (t <= 45000000.0) tmp = x + 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, -6.2e-6], t$95$1, If[LessEqual[t, 45000000.0], N[(x + z), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(b - a\right)\\
\mathbf{if}\;t \leq -6.2 \cdot 10^{-6}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 45000000:\\
\;\;\;\;x + z\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -6.1999999999999999e-6 or 4.5e7 < t Initial program 93.0%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f6464.7%
Simplified64.7%
if -6.1999999999999999e-6 < t < 4.5e7Initial program 99.2%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified99.2%
Taylor expanded in b around 0
associate-+r+N/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
*-rgt-identityN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f6470.7%
Simplified70.7%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
cancel-sign-sub-invN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6452.9%
Simplified52.9%
Taylor expanded in y around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f6431.5%
Simplified31.5%
Final simplification46.4%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* a (- 1.0 t)))) (if (<= a -1.4e+132) t_1 (if (<= a 1.75e-15) (+ x z) 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 (a <= -1.4e+132) {
tmp = t_1;
} else if (a <= 1.75e-15) {
tmp = x + 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 = a * (1.0d0 - t)
if (a <= (-1.4d+132)) then
tmp = t_1
else if (a <= 1.75d-15) then
tmp = x + 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 = a * (1.0 - t);
double tmp;
if (a <= -1.4e+132) {
tmp = t_1;
} else if (a <= 1.75e-15) {
tmp = x + z;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = a * (1.0 - t) tmp = 0 if a <= -1.4e+132: tmp = t_1 elif a <= 1.75e-15: tmp = x + z else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(a * Float64(1.0 - t)) tmp = 0.0 if (a <= -1.4e+132) tmp = t_1; elseif (a <= 1.75e-15) tmp = Float64(x + z); else tmp = 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 (a <= -1.4e+132) tmp = t_1; elseif (a <= 1.75e-15) tmp = x + z; else tmp = 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[a, -1.4e+132], t$95$1, If[LessEqual[a, 1.75e-15], N[(x + z), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \left(1 - t\right)\\
\mathbf{if}\;a \leq -1.4 \cdot 10^{+132}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 1.75 \cdot 10^{-15}:\\
\;\;\;\;x + z\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.4e132 or 1.75e-15 < a Initial program 95.3%
Taylor expanded in a around inf
sub-negN/A
metadata-evalN/A
neg-mul-1N/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--.f6452.1%
Simplified52.1%
if -1.4e132 < a < 1.75e-15Initial program 97.3%
Taylor expanded in t around 0
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate--r+N/A
sub-negN/A
+-lowering-+.f64N/A
sub-negN/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
Simplified96.6%
Taylor expanded in b around 0
associate-+r+N/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
*-rgt-identityN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f6457.4%
Simplified57.4%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
cancel-sign-sub-invN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6449.3%
Simplified49.3%
Taylor expanded in y around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f6433.0%
Simplified33.0%
Final simplification41.0%
(FPCore (x y z t a b) :precision binary64 (if (<= y -3.5e+40) (* y b) (if (<= y 1.05e+96) (* t b) (* y b))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -3.5e+40) {
tmp = y * b;
} else if (y <= 1.05e+96) {
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 <= (-3.5d+40)) then
tmp = y * b
else if (y <= 1.05d+96) 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 <= -3.5e+40) {
tmp = y * b;
} else if (y <= 1.05e+96) {
tmp = t * b;
} else {
tmp = y * b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -3.5e+40: tmp = y * b elif y <= 1.05e+96: tmp = t * b else: tmp = y * b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -3.5e+40) tmp = Float64(y * b); elseif (y <= 1.05e+96) 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 <= -3.5e+40) tmp = y * b; elseif (y <= 1.05e+96) tmp = t * b; else tmp = y * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -3.5e+40], N[(y * b), $MachinePrecision], If[LessEqual[y, 1.05e+96], N[(t * b), $MachinePrecision], N[(y * b), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -3.5 \cdot 10^{+40}:\\
\;\;\;\;y \cdot b\\
\mathbf{elif}\;y \leq 1.05 \cdot 10^{+96}:\\
\;\;\;\;t \cdot b\\
\mathbf{else}:\\
\;\;\;\;y \cdot b\\
\end{array}
\end{array}
if y < -3.4999999999999999e40 or 1.0500000000000001e96 < y Initial program 93.6%
Taylor expanded in y around inf
*-lowering-*.f64N/A
--lowering--.f6474.1%
Simplified74.1%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f6441.9%
Simplified41.9%
if -3.4999999999999999e40 < y < 1.0500000000000001e96Initial program 98.1%
Taylor expanded in x around inf
Simplified49.9%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6422.0%
Simplified22.0%
(FPCore (x y z t a b) :precision binary64 (if (<= z -5.2e+78) z (if (<= z 3.8e+191) x z)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -5.2e+78) {
tmp = z;
} else if (z <= 3.8e+191) {
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 <= (-5.2d+78)) then
tmp = z
else if (z <= 3.8d+191) 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 <= -5.2e+78) {
tmp = z;
} else if (z <= 3.8e+191) {
tmp = x;
} else {
tmp = z;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if z <= -5.2e+78: tmp = z elif z <= 3.8e+191: tmp = x else: tmp = z return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -5.2e+78) tmp = z; elseif (z <= 3.8e+191) tmp = x; else tmp = z; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (z <= -5.2e+78) tmp = z; elseif (z <= 3.8e+191) tmp = x; else tmp = z; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -5.2e+78], z, If[LessEqual[z, 3.8e+191], x, z]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -5.2 \cdot 10^{+78}:\\
\;\;\;\;z\\
\mathbf{elif}\;z \leq 3.8 \cdot 10^{+191}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;z\\
\end{array}
\end{array}
if z < -5.2e78 or 3.7999999999999998e191 < z Initial program 91.2%
Taylor expanded in z around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
neg-mul-1N/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6469.5%
Simplified69.5%
Taylor expanded in y around 0
Simplified29.5%
if -5.2e78 < z < 3.7999999999999998e191Initial program 98.4%
Taylor expanded in x around inf
Simplified18.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 96.4%
Taylor expanded in x around inf
Simplified14.5%
herbie shell --seed 2024152
(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)))