
(FPCore (x y z t a b) :precision binary64 (/ (- (+ (* (+ x y) z) (* (+ t y) a)) (* y b)) (+ (+ x t) y)))
double code(double x, double y, double z, double t, double a, double b) {
return ((((x + y) * z) + ((t + y) * a)) - (y * b)) / ((x + t) + y);
}
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) * z) + ((t + y) * a)) - (y * b)) / ((x + t) + y)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return ((((x + y) * z) + ((t + y) * a)) - (y * b)) / ((x + t) + y);
}
def code(x, y, z, t, a, b): return ((((x + y) * z) + ((t + y) * a)) - (y * b)) / ((x + t) + y)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(Float64(Float64(x + y) * z) + Float64(Float64(t + y) * a)) - Float64(y * b)) / Float64(Float64(x + t) + y)) end
function tmp = code(x, y, z, t, a, b) tmp = ((((x + y) * z) + ((t + y) * a)) - (y * b)) / ((x + t) + y); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(N[(N[(x + y), $MachinePrecision] * z), $MachinePrecision] + N[(N[(t + y), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / N[(N[(x + t), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(x + y\right) \cdot z + \left(t + y\right) \cdot a\right) - y \cdot b}{\left(x + t\right) + y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 24 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b) :precision binary64 (/ (- (+ (* (+ x y) z) (* (+ t y) a)) (* y b)) (+ (+ x t) y)))
double code(double x, double y, double z, double t, double a, double b) {
return ((((x + y) * z) + ((t + y) * a)) - (y * b)) / ((x + t) + y);
}
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) * z) + ((t + y) * a)) - (y * b)) / ((x + t) + y)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return ((((x + y) * z) + ((t + y) * a)) - (y * b)) / ((x + t) + y);
}
def code(x, y, z, t, a, b): return ((((x + y) * z) + ((t + y) * a)) - (y * b)) / ((x + t) + y)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(Float64(Float64(x + y) * z) + Float64(Float64(t + y) * a)) - Float64(y * b)) / Float64(Float64(x + t) + y)) end
function tmp = code(x, y, z, t, a, b) tmp = ((((x + y) * z) + ((t + y) * a)) - (y * b)) / ((x + t) + y); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(N[(N[(x + y), $MachinePrecision] * z), $MachinePrecision] + N[(N[(t + y), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / N[(N[(x + t), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(x + y\right) \cdot z + \left(t + y\right) \cdot a\right) - y \cdot b}{\left(x + t\right) + y}
\end{array}
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ y (+ x t)))
(t_2 (* (+ y t) a))
(t_3 (/ (- (+ t_2 (* z (+ x y))) (* y b)) t_1))
(t_4 (+ t (+ x y))))
(if (<= t_3 (- INFINITY))
(*
a
(+
(/ t t_1)
(- (fma z (/ (+ x y) (* a t_1)) (/ y t_1)) (* (/ b t_1) (/ y a)))))
(if (<= t_3 1e+259)
(- (+ (* z (+ (/ x t_4) (/ y t_4))) (/ t_2 t_4)) (/ (* y b) t_4))
(- (+ z a) b)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y + (x + t);
double t_2 = (y + t) * a;
double t_3 = ((t_2 + (z * (x + y))) - (y * b)) / t_1;
double t_4 = t + (x + y);
double tmp;
if (t_3 <= -((double) INFINITY)) {
tmp = a * ((t / t_1) + (fma(z, ((x + y) / (a * t_1)), (y / t_1)) - ((b / t_1) * (y / a))));
} else if (t_3 <= 1e+259) {
tmp = ((z * ((x / t_4) + (y / t_4))) + (t_2 / t_4)) - ((y * b) / t_4);
} else {
tmp = (z + a) - b;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(y + Float64(x + t)) t_2 = Float64(Float64(y + t) * a) t_3 = Float64(Float64(Float64(t_2 + Float64(z * Float64(x + y))) - Float64(y * b)) / t_1) t_4 = Float64(t + Float64(x + y)) tmp = 0.0 if (t_3 <= Float64(-Inf)) tmp = Float64(a * Float64(Float64(t / t_1) + Float64(fma(z, Float64(Float64(x + y) / Float64(a * t_1)), Float64(y / t_1)) - Float64(Float64(b / t_1) * Float64(y / a))))); elseif (t_3 <= 1e+259) tmp = Float64(Float64(Float64(z * Float64(Float64(x / t_4) + Float64(y / t_4))) + Float64(t_2 / t_4)) - Float64(Float64(y * b) / t_4)); else tmp = Float64(Float64(z + a) - b); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(y + t), $MachinePrecision] * a), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(t$95$2 + N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision]}, Block[{t$95$4 = N[(t + N[(x + y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, (-Infinity)], N[(a * N[(N[(t / t$95$1), $MachinePrecision] + N[(N[(z * N[(N[(x + y), $MachinePrecision] / N[(a * t$95$1), $MachinePrecision]), $MachinePrecision] + N[(y / t$95$1), $MachinePrecision]), $MachinePrecision] - N[(N[(b / t$95$1), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 1e+259], N[(N[(N[(z * N[(N[(x / t$95$4), $MachinePrecision] + N[(y / t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$2 / t$95$4), $MachinePrecision]), $MachinePrecision] - N[(N[(y * b), $MachinePrecision] / t$95$4), $MachinePrecision]), $MachinePrecision], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y + \left(x + t\right)\\
t_2 := \left(y + t\right) \cdot a\\
t_3 := \frac{\left(t\_2 + z \cdot \left(x + y\right)\right) - y \cdot b}{t\_1}\\
t_4 := t + \left(x + y\right)\\
\mathbf{if}\;t\_3 \leq -\infty:\\
\;\;\;\;a \cdot \left(\frac{t}{t\_1} + \left(\mathsf{fma}\left(z, \frac{x + y}{a \cdot t\_1}, \frac{y}{t\_1}\right) - \frac{b}{t\_1} \cdot \frac{y}{a}\right)\right)\\
\mathbf{elif}\;t\_3 \leq 10^{+259}:\\
\;\;\;\;\left(z \cdot \left(\frac{x}{t\_4} + \frac{y}{t\_4}\right) + \frac{t\_2}{t\_4}\right) - \frac{y \cdot b}{t\_4}\\
\mathbf{else}:\\
\;\;\;\;\left(z + a\right) - b\\
\end{array}
\end{array}
if (/.f64 (-.f64 (+.f64 (*.f64 (+.f64 x y) z) (*.f64 (+.f64 t y) a)) (*.f64 y b)) (+.f64 (+.f64 x t) y)) < -inf.0Initial program 5.9%
Taylor expanded in a around inf 40.5%
associate--l+40.5%
+-commutative40.5%
+-commutative40.5%
associate-+l+40.5%
+-commutative40.5%
Simplified66.3%
if -inf.0 < (/.f64 (-.f64 (+.f64 (*.f64 (+.f64 x y) z) (*.f64 (+.f64 t y) a)) (*.f64 y b)) (+.f64 (+.f64 x t) y)) < 9.999999999999999e258Initial program 99.1%
Taylor expanded in z around 0 99.1%
if 9.999999999999999e258 < (/.f64 (-.f64 (+.f64 (*.f64 (+.f64 x y) z) (*.f64 (+.f64 t y) a)) (*.f64 y b)) (+.f64 (+.f64 x t) y)) Initial program 9.1%
Taylor expanded in y around inf 71.9%
Final simplification87.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ t (+ x y)))
(t_2 (* (+ y t) a))
(t_3 (/ (- (+ t_2 (* z (+ x y))) (* y b)) (+ y (+ x t)))))
(if (or (<= t_3 -1e+194) (not (<= t_3 1e+259)))
(- (+ z a) b)
(- (+ (* z (+ (/ x t_1) (/ y t_1))) (/ t_2 t_1)) (/ (* y b) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = t + (x + y);
double t_2 = (y + t) * a;
double t_3 = ((t_2 + (z * (x + y))) - (y * b)) / (y + (x + t));
double tmp;
if ((t_3 <= -1e+194) || !(t_3 <= 1e+259)) {
tmp = (z + a) - b;
} else {
tmp = ((z * ((x / t_1) + (y / t_1))) + (t_2 / t_1)) - ((y * b) / t_1);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = t + (x + y)
t_2 = (y + t) * a
t_3 = ((t_2 + (z * (x + y))) - (y * b)) / (y + (x + t))
if ((t_3 <= (-1d+194)) .or. (.not. (t_3 <= 1d+259))) then
tmp = (z + a) - b
else
tmp = ((z * ((x / t_1) + (y / t_1))) + (t_2 / t_1)) - ((y * b) / t_1)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = t + (x + y);
double t_2 = (y + t) * a;
double t_3 = ((t_2 + (z * (x + y))) - (y * b)) / (y + (x + t));
double tmp;
if ((t_3 <= -1e+194) || !(t_3 <= 1e+259)) {
tmp = (z + a) - b;
} else {
tmp = ((z * ((x / t_1) + (y / t_1))) + (t_2 / t_1)) - ((y * b) / t_1);
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = t + (x + y) t_2 = (y + t) * a t_3 = ((t_2 + (z * (x + y))) - (y * b)) / (y + (x + t)) tmp = 0 if (t_3 <= -1e+194) or not (t_3 <= 1e+259): tmp = (z + a) - b else: tmp = ((z * ((x / t_1) + (y / t_1))) + (t_2 / t_1)) - ((y * b) / t_1) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(t + Float64(x + y)) t_2 = Float64(Float64(y + t) * a) t_3 = Float64(Float64(Float64(t_2 + Float64(z * Float64(x + y))) - Float64(y * b)) / Float64(y + Float64(x + t))) tmp = 0.0 if ((t_3 <= -1e+194) || !(t_3 <= 1e+259)) tmp = Float64(Float64(z + a) - b); else tmp = Float64(Float64(Float64(z * Float64(Float64(x / t_1) + Float64(y / t_1))) + Float64(t_2 / t_1)) - Float64(Float64(y * b) / t_1)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = t + (x + y); t_2 = (y + t) * a; t_3 = ((t_2 + (z * (x + y))) - (y * b)) / (y + (x + t)); tmp = 0.0; if ((t_3 <= -1e+194) || ~((t_3 <= 1e+259))) tmp = (z + a) - b; else tmp = ((z * ((x / t_1) + (y / t_1))) + (t_2 / t_1)) - ((y * b) / t_1); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(t + N[(x + y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(y + t), $MachinePrecision] * a), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(t$95$2 + N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t$95$3, -1e+194], N[Not[LessEqual[t$95$3, 1e+259]], $MachinePrecision]], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision], N[(N[(N[(z * N[(N[(x / t$95$1), $MachinePrecision] + N[(y / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$2 / t$95$1), $MachinePrecision]), $MachinePrecision] - N[(N[(y * b), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t + \left(x + y\right)\\
t_2 := \left(y + t\right) \cdot a\\
t_3 := \frac{\left(t\_2 + z \cdot \left(x + y\right)\right) - y \cdot b}{y + \left(x + t\right)}\\
\mathbf{if}\;t\_3 \leq -1 \cdot 10^{+194} \lor \neg \left(t\_3 \leq 10^{+259}\right):\\
\;\;\;\;\left(z + a\right) - b\\
\mathbf{else}:\\
\;\;\;\;\left(z \cdot \left(\frac{x}{t\_1} + \frac{y}{t\_1}\right) + \frac{t\_2}{t\_1}\right) - \frac{y \cdot b}{t\_1}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (+.f64 (*.f64 (+.f64 x y) z) (*.f64 (+.f64 t y) a)) (*.f64 y b)) (+.f64 (+.f64 x t) y)) < -9.99999999999999945e193 or 9.999999999999999e258 < (/.f64 (-.f64 (+.f64 (*.f64 (+.f64 x y) z) (*.f64 (+.f64 t y) a)) (*.f64 y b)) (+.f64 (+.f64 x t) y)) Initial program 15.1%
Taylor expanded in y around inf 71.9%
if -9.99999999999999945e193 < (/.f64 (-.f64 (+.f64 (*.f64 (+.f64 x y) z) (*.f64 (+.f64 t y) a)) (*.f64 y b)) (+.f64 (+.f64 x t) y)) < 9.999999999999999e258Initial program 99.0%
Taylor expanded in z around 0 99.0%
Final simplification86.8%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (/ (- (+ (* (+ y t) a) (* z (+ x y))) (* y b)) (+ y (+ x t))))) (if (or (<= t_1 -1e+194) (not (<= t_1 1e+259))) (- (+ z a) b) t_1)))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((((y + t) * a) + (z * (x + y))) - (y * b)) / (y + (x + t));
double tmp;
if ((t_1 <= -1e+194) || !(t_1 <= 1e+259)) {
tmp = (z + a) - b;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = ((((y + t) * a) + (z * (x + y))) - (y * b)) / (y + (x + t))
if ((t_1 <= (-1d+194)) .or. (.not. (t_1 <= 1d+259))) then
tmp = (z + a) - b
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((((y + t) * a) + (z * (x + y))) - (y * b)) / (y + (x + t));
double tmp;
if ((t_1 <= -1e+194) || !(t_1 <= 1e+259)) {
tmp = (z + a) - b;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = ((((y + t) * a) + (z * (x + y))) - (y * b)) / (y + (x + t)) tmp = 0 if (t_1 <= -1e+194) or not (t_1 <= 1e+259): tmp = (z + a) - b else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(Float64(Float64(y + t) * a) + Float64(z * Float64(x + y))) - Float64(y * b)) / Float64(y + Float64(x + t))) tmp = 0.0 if ((t_1 <= -1e+194) || !(t_1 <= 1e+259)) tmp = Float64(Float64(z + a) - b); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((((y + t) * a) + (z * (x + y))) - (y * b)) / (y + (x + t)); tmp = 0.0; if ((t_1 <= -1e+194) || ~((t_1 <= 1e+259))) tmp = (z + a) - b; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(N[(N[(y + t), $MachinePrecision] * a), $MachinePrecision] + N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t$95$1, -1e+194], N[Not[LessEqual[t$95$1, 1e+259]], $MachinePrecision]], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision], t$95$1]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\left(\left(y + t\right) \cdot a + z \cdot \left(x + y\right)\right) - y \cdot b}{y + \left(x + t\right)}\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{+194} \lor \neg \left(t\_1 \leq 10^{+259}\right):\\
\;\;\;\;\left(z + a\right) - b\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 (-.f64 (+.f64 (*.f64 (+.f64 x y) z) (*.f64 (+.f64 t y) a)) (*.f64 y b)) (+.f64 (+.f64 x t) y)) < -9.99999999999999945e193 or 9.999999999999999e258 < (/.f64 (-.f64 (+.f64 (*.f64 (+.f64 x y) z) (*.f64 (+.f64 t y) a)) (*.f64 y b)) (+.f64 (+.f64 x t) y)) Initial program 15.1%
Taylor expanded in y around inf 71.9%
if -9.99999999999999945e193 < (/.f64 (-.f64 (+.f64 (*.f64 (+.f64 x y) z) (*.f64 (+.f64 t y) a)) (*.f64 y b)) (+.f64 (+.f64 x t) y)) < 9.999999999999999e258Initial program 99.0%
Final simplification86.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ y (+ x t))) (t_2 (- (+ z a) b)) (t_3 (* z (+ x y))))
(if (<= y -1.06e+154)
t_2
(if (<= y -1.2e+76)
(+ z (- (* a (/ y x)) (* b (/ y x))))
(if (<= y -3.4e+14)
(/ (- (+ (* y a) t_3) (* y b)) (+ x y))
(if (<= y -2.65e-36)
(* a (/ (+ y t) t_1))
(if (<= y -4.6e-134)
(/ 1.0 (/ (+ t (+ x y)) (- t_3 (* y b))))
(if (<= y 1.75e-60)
(/ (+ (* t a) (* x z)) (+ x t))
(if (<= y 6e+69)
(* a (+ 1.0 (- (/ z a) (* (/ b t_1) (/ y a)))))
t_2)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y + (x + t);
double t_2 = (z + a) - b;
double t_3 = z * (x + y);
double tmp;
if (y <= -1.06e+154) {
tmp = t_2;
} else if (y <= -1.2e+76) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -3.4e+14) {
tmp = (((y * a) + t_3) - (y * b)) / (x + y);
} else if (y <= -2.65e-36) {
tmp = a * ((y + t) / t_1);
} else if (y <= -4.6e-134) {
tmp = 1.0 / ((t + (x + y)) / (t_3 - (y * b)));
} else if (y <= 1.75e-60) {
tmp = ((t * a) + (x * z)) / (x + t);
} else if (y <= 6e+69) {
tmp = a * (1.0 + ((z / a) - ((b / t_1) * (y / a))));
} 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) :: t_3
real(8) :: tmp
t_1 = y + (x + t)
t_2 = (z + a) - b
t_3 = z * (x + y)
if (y <= (-1.06d+154)) then
tmp = t_2
else if (y <= (-1.2d+76)) then
tmp = z + ((a * (y / x)) - (b * (y / x)))
else if (y <= (-3.4d+14)) then
tmp = (((y * a) + t_3) - (y * b)) / (x + y)
else if (y <= (-2.65d-36)) then
tmp = a * ((y + t) / t_1)
else if (y <= (-4.6d-134)) then
tmp = 1.0d0 / ((t + (x + y)) / (t_3 - (y * b)))
else if (y <= 1.75d-60) then
tmp = ((t * a) + (x * z)) / (x + t)
else if (y <= 6d+69) then
tmp = a * (1.0d0 + ((z / a) - ((b / t_1) * (y / a))))
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 = y + (x + t);
double t_2 = (z + a) - b;
double t_3 = z * (x + y);
double tmp;
if (y <= -1.06e+154) {
tmp = t_2;
} else if (y <= -1.2e+76) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -3.4e+14) {
tmp = (((y * a) + t_3) - (y * b)) / (x + y);
} else if (y <= -2.65e-36) {
tmp = a * ((y + t) / t_1);
} else if (y <= -4.6e-134) {
tmp = 1.0 / ((t + (x + y)) / (t_3 - (y * b)));
} else if (y <= 1.75e-60) {
tmp = ((t * a) + (x * z)) / (x + t);
} else if (y <= 6e+69) {
tmp = a * (1.0 + ((z / a) - ((b / t_1) * (y / a))));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y + (x + t) t_2 = (z + a) - b t_3 = z * (x + y) tmp = 0 if y <= -1.06e+154: tmp = t_2 elif y <= -1.2e+76: tmp = z + ((a * (y / x)) - (b * (y / x))) elif y <= -3.4e+14: tmp = (((y * a) + t_3) - (y * b)) / (x + y) elif y <= -2.65e-36: tmp = a * ((y + t) / t_1) elif y <= -4.6e-134: tmp = 1.0 / ((t + (x + y)) / (t_3 - (y * b))) elif y <= 1.75e-60: tmp = ((t * a) + (x * z)) / (x + t) elif y <= 6e+69: tmp = a * (1.0 + ((z / a) - ((b / t_1) * (y / a)))) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(y + Float64(x + t)) t_2 = Float64(Float64(z + a) - b) t_3 = Float64(z * Float64(x + y)) tmp = 0.0 if (y <= -1.06e+154) tmp = t_2; elseif (y <= -1.2e+76) tmp = Float64(z + Float64(Float64(a * Float64(y / x)) - Float64(b * Float64(y / x)))); elseif (y <= -3.4e+14) tmp = Float64(Float64(Float64(Float64(y * a) + t_3) - Float64(y * b)) / Float64(x + y)); elseif (y <= -2.65e-36) tmp = Float64(a * Float64(Float64(y + t) / t_1)); elseif (y <= -4.6e-134) tmp = Float64(1.0 / Float64(Float64(t + Float64(x + y)) / Float64(t_3 - Float64(y * b)))); elseif (y <= 1.75e-60) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); elseif (y <= 6e+69) tmp = Float64(a * Float64(1.0 + Float64(Float64(z / a) - Float64(Float64(b / t_1) * Float64(y / a))))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y + (x + t); t_2 = (z + a) - b; t_3 = z * (x + y); tmp = 0.0; if (y <= -1.06e+154) tmp = t_2; elseif (y <= -1.2e+76) tmp = z + ((a * (y / x)) - (b * (y / x))); elseif (y <= -3.4e+14) tmp = (((y * a) + t_3) - (y * b)) / (x + y); elseif (y <= -2.65e-36) tmp = a * ((y + t) / t_1); elseif (y <= -4.6e-134) tmp = 1.0 / ((t + (x + y)) / (t_3 - (y * b))); elseif (y <= 1.75e-60) tmp = ((t * a) + (x * z)) / (x + t); elseif (y <= 6e+69) tmp = a * (1.0 + ((z / a) - ((b / t_1) * (y / a)))); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, Block[{t$95$3 = N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.06e+154], t$95$2, If[LessEqual[y, -1.2e+76], N[(z + N[(N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision] - N[(b * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -3.4e+14], N[(N[(N[(N[(y * a), $MachinePrecision] + t$95$3), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / N[(x + y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -2.65e-36], N[(a * N[(N[(y + t), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -4.6e-134], N[(1.0 / N[(N[(t + N[(x + y), $MachinePrecision]), $MachinePrecision] / N[(t$95$3 - N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.75e-60], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 6e+69], N[(a * N[(1.0 + N[(N[(z / a), $MachinePrecision] - N[(N[(b / t$95$1), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y + \left(x + t\right)\\
t_2 := \left(z + a\right) - b\\
t_3 := z \cdot \left(x + y\right)\\
\mathbf{if}\;y \leq -1.06 \cdot 10^{+154}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.2 \cdot 10^{+76}:\\
\;\;\;\;z + \left(a \cdot \frac{y}{x} - b \cdot \frac{y}{x}\right)\\
\mathbf{elif}\;y \leq -3.4 \cdot 10^{+14}:\\
\;\;\;\;\frac{\left(y \cdot a + t\_3\right) - y \cdot b}{x + y}\\
\mathbf{elif}\;y \leq -2.65 \cdot 10^{-36}:\\
\;\;\;\;a \cdot \frac{y + t}{t\_1}\\
\mathbf{elif}\;y \leq -4.6 \cdot 10^{-134}:\\
\;\;\;\;\frac{1}{\frac{t + \left(x + y\right)}{t\_3 - y \cdot b}}\\
\mathbf{elif}\;y \leq 1.75 \cdot 10^{-60}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{elif}\;y \leq 6 \cdot 10^{+69}:\\
\;\;\;\;a \cdot \left(1 + \left(\frac{z}{a} - \frac{b}{t\_1} \cdot \frac{y}{a}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y < -1.06e154 or 5.99999999999999967e69 < y Initial program 29.2%
Taylor expanded in y around inf 82.3%
if -1.06e154 < y < -1.2e76Initial program 41.1%
Taylor expanded in x around inf 32.9%
associate--l+40.4%
associate-/l*47.7%
associate-/l*47.7%
associate-/l*58.4%
associate-/l*58.4%
Simplified58.4%
Taylor expanded in t around 0 48.3%
associate-/l*55.7%
associate-*r/66.4%
Simplified66.4%
if -1.2e76 < y < -3.4e14Initial program 76.1%
Taylor expanded in t around 0 67.7%
if -3.4e14 < y < -2.6499999999999999e-36Initial program 65.3%
Taylor expanded in a around inf 47.2%
associate-/l*73.2%
+-commutative73.2%
+-commutative73.2%
associate-+l+73.2%
+-commutative73.2%
Simplified73.2%
if -2.6499999999999999e-36 < y < -4.6000000000000001e-134Initial program 86.3%
Taylor expanded in a around 0 76.6%
+-commutative76.6%
*-commutative76.6%
Simplified76.6%
clear-num76.9%
inv-pow76.9%
+-commutative76.9%
+-commutative76.9%
+-commutative76.9%
*-commutative76.9%
Applied egg-rr76.9%
unpow-176.9%
+-commutative76.9%
+-commutative76.9%
associate-+r+76.9%
+-commutative76.9%
Simplified76.9%
if -4.6000000000000001e-134 < y < 1.74999999999999988e-60Initial program 83.7%
Taylor expanded in y around 0 70.4%
if 1.74999999999999988e-60 < y < 5.99999999999999967e69Initial program 74.2%
Taylor expanded in a around inf 81.6%
associate--l+81.6%
+-commutative81.6%
+-commutative81.6%
associate-+l+81.6%
+-commutative81.6%
Simplified88.5%
Taylor expanded in x around inf 76.1%
Taylor expanded in t around inf 77.2%
Final simplification75.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (+ y t))) (t_2 (- (+ z a) b)) (t_3 (* z (+ x y))))
(if (<= y -1.06e+154)
t_2
(if (<= y -1.28e+76)
(+ z (- (* a (/ y x)) (* b (/ y x))))
(if (<= y -3.3e+28)
(/ (- (+ (* y a) t_3) (* y b)) (+ x y))
(if (<= y -8.6e-36)
(* a (- (/ t t_1) (* b (/ (/ y a) t_1))))
(if (<= y -4.8e-134)
(/ 1.0 (/ (+ t (+ x y)) (- t_3 (* y b))))
(if (<= y 1.55e-61)
(/ (+ (* t a) (* x z)) (+ x t))
(if (<= y 6.6e+69)
(* a (+ 1.0 (- (/ z a) (* (/ b (+ y (+ x t))) (/ y a)))))
t_2)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y + t);
double t_2 = (z + a) - b;
double t_3 = z * (x + y);
double tmp;
if (y <= -1.06e+154) {
tmp = t_2;
} else if (y <= -1.28e+76) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -3.3e+28) {
tmp = (((y * a) + t_3) - (y * b)) / (x + y);
} else if (y <= -8.6e-36) {
tmp = a * ((t / t_1) - (b * ((y / a) / t_1)));
} else if (y <= -4.8e-134) {
tmp = 1.0 / ((t + (x + y)) / (t_3 - (y * b)));
} else if (y <= 1.55e-61) {
tmp = ((t * a) + (x * z)) / (x + t);
} else if (y <= 6.6e+69) {
tmp = a * (1.0 + ((z / a) - ((b / (y + (x + t))) * (y / a))));
} 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) :: t_3
real(8) :: tmp
t_1 = x + (y + t)
t_2 = (z + a) - b
t_3 = z * (x + y)
if (y <= (-1.06d+154)) then
tmp = t_2
else if (y <= (-1.28d+76)) then
tmp = z + ((a * (y / x)) - (b * (y / x)))
else if (y <= (-3.3d+28)) then
tmp = (((y * a) + t_3) - (y * b)) / (x + y)
else if (y <= (-8.6d-36)) then
tmp = a * ((t / t_1) - (b * ((y / a) / t_1)))
else if (y <= (-4.8d-134)) then
tmp = 1.0d0 / ((t + (x + y)) / (t_3 - (y * b)))
else if (y <= 1.55d-61) then
tmp = ((t * a) + (x * z)) / (x + t)
else if (y <= 6.6d+69) then
tmp = a * (1.0d0 + ((z / a) - ((b / (y + (x + t))) * (y / a))))
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 = x + (y + t);
double t_2 = (z + a) - b;
double t_3 = z * (x + y);
double tmp;
if (y <= -1.06e+154) {
tmp = t_2;
} else if (y <= -1.28e+76) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -3.3e+28) {
tmp = (((y * a) + t_3) - (y * b)) / (x + y);
} else if (y <= -8.6e-36) {
tmp = a * ((t / t_1) - (b * ((y / a) / t_1)));
} else if (y <= -4.8e-134) {
tmp = 1.0 / ((t + (x + y)) / (t_3 - (y * b)));
} else if (y <= 1.55e-61) {
tmp = ((t * a) + (x * z)) / (x + t);
} else if (y <= 6.6e+69) {
tmp = a * (1.0 + ((z / a) - ((b / (y + (x + t))) * (y / a))));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (y + t) t_2 = (z + a) - b t_3 = z * (x + y) tmp = 0 if y <= -1.06e+154: tmp = t_2 elif y <= -1.28e+76: tmp = z + ((a * (y / x)) - (b * (y / x))) elif y <= -3.3e+28: tmp = (((y * a) + t_3) - (y * b)) / (x + y) elif y <= -8.6e-36: tmp = a * ((t / t_1) - (b * ((y / a) / t_1))) elif y <= -4.8e-134: tmp = 1.0 / ((t + (x + y)) / (t_3 - (y * b))) elif y <= 1.55e-61: tmp = ((t * a) + (x * z)) / (x + t) elif y <= 6.6e+69: tmp = a * (1.0 + ((z / a) - ((b / (y + (x + t))) * (y / a)))) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(y + t)) t_2 = Float64(Float64(z + a) - b) t_3 = Float64(z * Float64(x + y)) tmp = 0.0 if (y <= -1.06e+154) tmp = t_2; elseif (y <= -1.28e+76) tmp = Float64(z + Float64(Float64(a * Float64(y / x)) - Float64(b * Float64(y / x)))); elseif (y <= -3.3e+28) tmp = Float64(Float64(Float64(Float64(y * a) + t_3) - Float64(y * b)) / Float64(x + y)); elseif (y <= -8.6e-36) tmp = Float64(a * Float64(Float64(t / t_1) - Float64(b * Float64(Float64(y / a) / t_1)))); elseif (y <= -4.8e-134) tmp = Float64(1.0 / Float64(Float64(t + Float64(x + y)) / Float64(t_3 - Float64(y * b)))); elseif (y <= 1.55e-61) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); elseif (y <= 6.6e+69) tmp = Float64(a * Float64(1.0 + Float64(Float64(z / a) - Float64(Float64(b / Float64(y + Float64(x + t))) * Float64(y / a))))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (y + t); t_2 = (z + a) - b; t_3 = z * (x + y); tmp = 0.0; if (y <= -1.06e+154) tmp = t_2; elseif (y <= -1.28e+76) tmp = z + ((a * (y / x)) - (b * (y / x))); elseif (y <= -3.3e+28) tmp = (((y * a) + t_3) - (y * b)) / (x + y); elseif (y <= -8.6e-36) tmp = a * ((t / t_1) - (b * ((y / a) / t_1))); elseif (y <= -4.8e-134) tmp = 1.0 / ((t + (x + y)) / (t_3 - (y * b))); elseif (y <= 1.55e-61) tmp = ((t * a) + (x * z)) / (x + t); elseif (y <= 6.6e+69) tmp = a * (1.0 + ((z / a) - ((b / (y + (x + t))) * (y / a)))); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(y + t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, Block[{t$95$3 = N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.06e+154], t$95$2, If[LessEqual[y, -1.28e+76], N[(z + N[(N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision] - N[(b * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -3.3e+28], N[(N[(N[(N[(y * a), $MachinePrecision] + t$95$3), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / N[(x + y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -8.6e-36], N[(a * N[(N[(t / t$95$1), $MachinePrecision] - N[(b * N[(N[(y / a), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -4.8e-134], N[(1.0 / N[(N[(t + N[(x + y), $MachinePrecision]), $MachinePrecision] / N[(t$95$3 - N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.55e-61], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 6.6e+69], N[(a * N[(1.0 + N[(N[(z / a), $MachinePrecision] - N[(N[(b / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(y + t\right)\\
t_2 := \left(z + a\right) - b\\
t_3 := z \cdot \left(x + y\right)\\
\mathbf{if}\;y \leq -1.06 \cdot 10^{+154}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.28 \cdot 10^{+76}:\\
\;\;\;\;z + \left(a \cdot \frac{y}{x} - b \cdot \frac{y}{x}\right)\\
\mathbf{elif}\;y \leq -3.3 \cdot 10^{+28}:\\
\;\;\;\;\frac{\left(y \cdot a + t\_3\right) - y \cdot b}{x + y}\\
\mathbf{elif}\;y \leq -8.6 \cdot 10^{-36}:\\
\;\;\;\;a \cdot \left(\frac{t}{t\_1} - b \cdot \frac{\frac{y}{a}}{t\_1}\right)\\
\mathbf{elif}\;y \leq -4.8 \cdot 10^{-134}:\\
\;\;\;\;\frac{1}{\frac{t + \left(x + y\right)}{t\_3 - y \cdot b}}\\
\mathbf{elif}\;y \leq 1.55 \cdot 10^{-61}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{elif}\;y \leq 6.6 \cdot 10^{+69}:\\
\;\;\;\;a \cdot \left(1 + \left(\frac{z}{a} - \frac{b}{y + \left(x + t\right)} \cdot \frac{y}{a}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y < -1.06e154 or 6.5999999999999997e69 < y Initial program 29.2%
Taylor expanded in y around inf 82.3%
if -1.06e154 < y < -1.27999999999999994e76Initial program 41.1%
Taylor expanded in x around inf 32.9%
associate--l+40.4%
associate-/l*47.7%
associate-/l*47.7%
associate-/l*58.4%
associate-/l*58.4%
Simplified58.4%
Taylor expanded in t around 0 48.3%
associate-/l*55.7%
associate-*r/66.4%
Simplified66.4%
if -1.27999999999999994e76 < y < -3.3e28Initial program 86.3%
Taylor expanded in t around 0 86.3%
if -3.3e28 < y < -8.6000000000000004e-36Initial program 64.1%
Taylor expanded in a around inf 82.3%
associate--l+82.3%
+-commutative82.3%
+-commutative82.3%
associate-+l+82.3%
+-commutative82.3%
Simplified88.0%
Taylor expanded in x around inf 57.7%
Taylor expanded in z around 0 64.1%
+-commutative64.1%
associate-+r+64.1%
+-commutative64.1%
associate-/l*70.0%
associate-/r*70.0%
+-commutative70.0%
associate-+r+70.0%
+-commutative70.0%
Simplified70.0%
if -8.6000000000000004e-36 < y < -4.80000000000000019e-134Initial program 86.3%
Taylor expanded in a around 0 76.6%
+-commutative76.6%
*-commutative76.6%
Simplified76.6%
clear-num76.9%
inv-pow76.9%
+-commutative76.9%
+-commutative76.9%
+-commutative76.9%
*-commutative76.9%
Applied egg-rr76.9%
unpow-176.9%
+-commutative76.9%
+-commutative76.9%
associate-+r+76.9%
+-commutative76.9%
Simplified76.9%
if -4.80000000000000019e-134 < y < 1.54999999999999997e-61Initial program 83.7%
Taylor expanded in y around 0 70.4%
if 1.54999999999999997e-61 < y < 6.5999999999999997e69Initial program 74.2%
Taylor expanded in a around inf 81.6%
associate--l+81.6%
+-commutative81.6%
+-commutative81.6%
associate-+l+81.6%
+-commutative81.6%
Simplified88.5%
Taylor expanded in x around inf 76.1%
Taylor expanded in t around inf 77.2%
Final simplification75.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (- (+ z a) b)) (t_2 (+ y (+ x t))))
(if (<= y -1.06e+154)
t_1
(if (<= y -1.25e+76)
(+ z (- (* a (/ y x)) (* b (/ y x))))
(if (<= y -2.7e+14)
(/ (* y t_1) t_2)
(if (<= y -2.9e-40)
(* a (/ (+ y t) t_2))
(if (<= y -4.6e-134)
(/ 1.0 (/ (+ t (+ x y)) (- (* z (+ x y)) (* y b))))
(if (<= y 1.18e-29) (/ (+ (* t a) (* x z)) (+ x t)) t_1))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (z + a) - b;
double t_2 = y + (x + t);
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -1.25e+76) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -2.7e+14) {
tmp = (y * t_1) / t_2;
} else if (y <= -2.9e-40) {
tmp = a * ((y + t) / t_2);
} else if (y <= -4.6e-134) {
tmp = 1.0 / ((t + (x + y)) / ((z * (x + y)) - (y * b)));
} else if (y <= 1.18e-29) {
tmp = ((t * a) + (x * z)) / (x + 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) :: t_2
real(8) :: tmp
t_1 = (z + a) - b
t_2 = y + (x + t)
if (y <= (-1.06d+154)) then
tmp = t_1
else if (y <= (-1.25d+76)) then
tmp = z + ((a * (y / x)) - (b * (y / x)))
else if (y <= (-2.7d+14)) then
tmp = (y * t_1) / t_2
else if (y <= (-2.9d-40)) then
tmp = a * ((y + t) / t_2)
else if (y <= (-4.6d-134)) then
tmp = 1.0d0 / ((t + (x + y)) / ((z * (x + y)) - (y * b)))
else if (y <= 1.18d-29) then
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b;
double t_2 = y + (x + t);
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -1.25e+76) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -2.7e+14) {
tmp = (y * t_1) / t_2;
} else if (y <= -2.9e-40) {
tmp = a * ((y + t) / t_2);
} else if (y <= -4.6e-134) {
tmp = 1.0 / ((t + (x + y)) / ((z * (x + y)) - (y * b)));
} else if (y <= 1.18e-29) {
tmp = ((t * a) + (x * z)) / (x + t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (z + a) - b t_2 = y + (x + t) tmp = 0 if y <= -1.06e+154: tmp = t_1 elif y <= -1.25e+76: tmp = z + ((a * (y / x)) - (b * (y / x))) elif y <= -2.7e+14: tmp = (y * t_1) / t_2 elif y <= -2.9e-40: tmp = a * ((y + t) / t_2) elif y <= -4.6e-134: tmp = 1.0 / ((t + (x + y)) / ((z * (x + y)) - (y * b))) elif y <= 1.18e-29: tmp = ((t * a) + (x * z)) / (x + t) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(z + a) - b) t_2 = Float64(y + Float64(x + t)) tmp = 0.0 if (y <= -1.06e+154) tmp = t_1; elseif (y <= -1.25e+76) tmp = Float64(z + Float64(Float64(a * Float64(y / x)) - Float64(b * Float64(y / x)))); elseif (y <= -2.7e+14) tmp = Float64(Float64(y * t_1) / t_2); elseif (y <= -2.9e-40) tmp = Float64(a * Float64(Float64(y + t) / t_2)); elseif (y <= -4.6e-134) tmp = Float64(1.0 / Float64(Float64(t + Float64(x + y)) / Float64(Float64(z * Float64(x + y)) - Float64(y * b)))); elseif (y <= 1.18e-29) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (z + a) - b; t_2 = y + (x + t); tmp = 0.0; if (y <= -1.06e+154) tmp = t_1; elseif (y <= -1.25e+76) tmp = z + ((a * (y / x)) - (b * (y / x))); elseif (y <= -2.7e+14) tmp = (y * t_1) / t_2; elseif (y <= -2.9e-40) tmp = a * ((y + t) / t_2); elseif (y <= -4.6e-134) tmp = 1.0 / ((t + (x + y)) / ((z * (x + y)) - (y * b))); elseif (y <= 1.18e-29) tmp = ((t * a) + (x * z)) / (x + t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, Block[{t$95$2 = N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.06e+154], t$95$1, If[LessEqual[y, -1.25e+76], N[(z + N[(N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision] - N[(b * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -2.7e+14], N[(N[(y * t$95$1), $MachinePrecision] / t$95$2), $MachinePrecision], If[LessEqual[y, -2.9e-40], N[(a * N[(N[(y + t), $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -4.6e-134], N[(1.0 / N[(N[(t + N[(x + y), $MachinePrecision]), $MachinePrecision] / N[(N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.18e-29], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(z + a\right) - b\\
t_2 := y + \left(x + t\right)\\
\mathbf{if}\;y \leq -1.06 \cdot 10^{+154}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.25 \cdot 10^{+76}:\\
\;\;\;\;z + \left(a \cdot \frac{y}{x} - b \cdot \frac{y}{x}\right)\\
\mathbf{elif}\;y \leq -2.7 \cdot 10^{+14}:\\
\;\;\;\;\frac{y \cdot t\_1}{t\_2}\\
\mathbf{elif}\;y \leq -2.9 \cdot 10^{-40}:\\
\;\;\;\;a \cdot \frac{y + t}{t\_2}\\
\mathbf{elif}\;y \leq -4.6 \cdot 10^{-134}:\\
\;\;\;\;\frac{1}{\frac{t + \left(x + y\right)}{z \cdot \left(x + y\right) - y \cdot b}}\\
\mathbf{elif}\;y \leq 1.18 \cdot 10^{-29}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.06e154 or 1.17999999999999996e-29 < y Initial program 37.7%
Taylor expanded in y around inf 76.7%
if -1.06e154 < y < -1.24999999999999998e76Initial program 41.1%
Taylor expanded in x around inf 32.9%
associate--l+40.4%
associate-/l*47.7%
associate-/l*47.7%
associate-/l*58.4%
associate-/l*58.4%
Simplified58.4%
Taylor expanded in t around 0 48.3%
associate-/l*55.7%
associate-*r/66.4%
Simplified66.4%
if -1.24999999999999998e76 < y < -2.7e14Initial program 76.1%
Taylor expanded in y around inf 67.7%
if -2.7e14 < y < -2.8999999999999999e-40Initial program 65.3%
Taylor expanded in a around inf 47.2%
associate-/l*73.2%
+-commutative73.2%
+-commutative73.2%
associate-+l+73.2%
+-commutative73.2%
Simplified73.2%
if -2.8999999999999999e-40 < y < -4.6000000000000001e-134Initial program 86.3%
Taylor expanded in a around 0 76.6%
+-commutative76.6%
*-commutative76.6%
Simplified76.6%
clear-num76.9%
inv-pow76.9%
+-commutative76.9%
+-commutative76.9%
+-commutative76.9%
*-commutative76.9%
Applied egg-rr76.9%
unpow-176.9%
+-commutative76.9%
+-commutative76.9%
associate-+r+76.9%
+-commutative76.9%
Simplified76.9%
if -4.6000000000000001e-134 < y < 1.17999999999999996e-29Initial program 83.5%
Taylor expanded in y around 0 69.9%
Final simplification73.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (- (+ z a) b)) (t_2 (* z (+ x y))))
(if (<= y -1.06e+154)
t_1
(if (<= y -6.1e+75)
(+ z (- (* a (/ y x)) (* b (/ y x))))
(if (<= y -8.2e+15)
(/ (- (+ (* y a) t_2) (* y b)) (+ x y))
(if (<= y -8.2e-40)
(* a (/ (+ y t) (+ y (+ x t))))
(if (<= y -4.6e-134)
(/ 1.0 (/ (+ t (+ x y)) (- t_2 (* y b))))
(if (<= y 1.15e-29) (/ (+ (* t a) (* x z)) (+ x t)) t_1))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (z + a) - b;
double t_2 = z * (x + y);
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -6.1e+75) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -8.2e+15) {
tmp = (((y * a) + t_2) - (y * b)) / (x + y);
} else if (y <= -8.2e-40) {
tmp = a * ((y + t) / (y + (x + t)));
} else if (y <= -4.6e-134) {
tmp = 1.0 / ((t + (x + y)) / (t_2 - (y * b)));
} else if (y <= 1.15e-29) {
tmp = ((t * a) + (x * z)) / (x + 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) :: t_2
real(8) :: tmp
t_1 = (z + a) - b
t_2 = z * (x + y)
if (y <= (-1.06d+154)) then
tmp = t_1
else if (y <= (-6.1d+75)) then
tmp = z + ((a * (y / x)) - (b * (y / x)))
else if (y <= (-8.2d+15)) then
tmp = (((y * a) + t_2) - (y * b)) / (x + y)
else if (y <= (-8.2d-40)) then
tmp = a * ((y + t) / (y + (x + t)))
else if (y <= (-4.6d-134)) then
tmp = 1.0d0 / ((t + (x + y)) / (t_2 - (y * b)))
else if (y <= 1.15d-29) then
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b;
double t_2 = z * (x + y);
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -6.1e+75) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -8.2e+15) {
tmp = (((y * a) + t_2) - (y * b)) / (x + y);
} else if (y <= -8.2e-40) {
tmp = a * ((y + t) / (y + (x + t)));
} else if (y <= -4.6e-134) {
tmp = 1.0 / ((t + (x + y)) / (t_2 - (y * b)));
} else if (y <= 1.15e-29) {
tmp = ((t * a) + (x * z)) / (x + t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (z + a) - b t_2 = z * (x + y) tmp = 0 if y <= -1.06e+154: tmp = t_1 elif y <= -6.1e+75: tmp = z + ((a * (y / x)) - (b * (y / x))) elif y <= -8.2e+15: tmp = (((y * a) + t_2) - (y * b)) / (x + y) elif y <= -8.2e-40: tmp = a * ((y + t) / (y + (x + t))) elif y <= -4.6e-134: tmp = 1.0 / ((t + (x + y)) / (t_2 - (y * b))) elif y <= 1.15e-29: tmp = ((t * a) + (x * z)) / (x + t) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(z + a) - b) t_2 = Float64(z * Float64(x + y)) tmp = 0.0 if (y <= -1.06e+154) tmp = t_1; elseif (y <= -6.1e+75) tmp = Float64(z + Float64(Float64(a * Float64(y / x)) - Float64(b * Float64(y / x)))); elseif (y <= -8.2e+15) tmp = Float64(Float64(Float64(Float64(y * a) + t_2) - Float64(y * b)) / Float64(x + y)); elseif (y <= -8.2e-40) tmp = Float64(a * Float64(Float64(y + t) / Float64(y + Float64(x + t)))); elseif (y <= -4.6e-134) tmp = Float64(1.0 / Float64(Float64(t + Float64(x + y)) / Float64(t_2 - Float64(y * b)))); elseif (y <= 1.15e-29) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (z + a) - b; t_2 = z * (x + y); tmp = 0.0; if (y <= -1.06e+154) tmp = t_1; elseif (y <= -6.1e+75) tmp = z + ((a * (y / x)) - (b * (y / x))); elseif (y <= -8.2e+15) tmp = (((y * a) + t_2) - (y * b)) / (x + y); elseif (y <= -8.2e-40) tmp = a * ((y + t) / (y + (x + t))); elseif (y <= -4.6e-134) tmp = 1.0 / ((t + (x + y)) / (t_2 - (y * b))); elseif (y <= 1.15e-29) tmp = ((t * a) + (x * z)) / (x + t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, Block[{t$95$2 = N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.06e+154], t$95$1, If[LessEqual[y, -6.1e+75], N[(z + N[(N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision] - N[(b * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -8.2e+15], N[(N[(N[(N[(y * a), $MachinePrecision] + t$95$2), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / N[(x + y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -8.2e-40], N[(a * N[(N[(y + t), $MachinePrecision] / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -4.6e-134], N[(1.0 / N[(N[(t + N[(x + y), $MachinePrecision]), $MachinePrecision] / N[(t$95$2 - N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.15e-29], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(z + a\right) - b\\
t_2 := z \cdot \left(x + y\right)\\
\mathbf{if}\;y \leq -1.06 \cdot 10^{+154}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -6.1 \cdot 10^{+75}:\\
\;\;\;\;z + \left(a \cdot \frac{y}{x} - b \cdot \frac{y}{x}\right)\\
\mathbf{elif}\;y \leq -8.2 \cdot 10^{+15}:\\
\;\;\;\;\frac{\left(y \cdot a + t\_2\right) - y \cdot b}{x + y}\\
\mathbf{elif}\;y \leq -8.2 \cdot 10^{-40}:\\
\;\;\;\;a \cdot \frac{y + t}{y + \left(x + t\right)}\\
\mathbf{elif}\;y \leq -4.6 \cdot 10^{-134}:\\
\;\;\;\;\frac{1}{\frac{t + \left(x + y\right)}{t\_2 - y \cdot b}}\\
\mathbf{elif}\;y \leq 1.15 \cdot 10^{-29}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.06e154 or 1.14999999999999996e-29 < y Initial program 37.7%
Taylor expanded in y around inf 76.7%
if -1.06e154 < y < -6.10000000000000009e75Initial program 41.1%
Taylor expanded in x around inf 32.9%
associate--l+40.4%
associate-/l*47.7%
associate-/l*47.7%
associate-/l*58.4%
associate-/l*58.4%
Simplified58.4%
Taylor expanded in t around 0 48.3%
associate-/l*55.7%
associate-*r/66.4%
Simplified66.4%
if -6.10000000000000009e75 < y < -8.2e15Initial program 76.1%
Taylor expanded in t around 0 67.7%
if -8.2e15 < y < -8.19999999999999926e-40Initial program 65.3%
Taylor expanded in a around inf 47.2%
associate-/l*73.2%
+-commutative73.2%
+-commutative73.2%
associate-+l+73.2%
+-commutative73.2%
Simplified73.2%
if -8.19999999999999926e-40 < y < -4.6000000000000001e-134Initial program 86.3%
Taylor expanded in a around 0 76.6%
+-commutative76.6%
*-commutative76.6%
Simplified76.6%
clear-num76.9%
inv-pow76.9%
+-commutative76.9%
+-commutative76.9%
+-commutative76.9%
*-commutative76.9%
Applied egg-rr76.9%
unpow-176.9%
+-commutative76.9%
+-commutative76.9%
associate-+r+76.9%
+-commutative76.9%
Simplified76.9%
if -4.6000000000000001e-134 < y < 1.14999999999999996e-29Initial program 83.5%
Taylor expanded in y around 0 69.9%
Final simplification73.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ y (+ x t))) (t_2 (- (+ z a) b)))
(if (<= y -1.1e+154)
t_2
(if (<= y -1.28e+76)
(+ z (* y (- (/ a x) (/ b x))))
(if (<= y -1.9e+54)
t_2
(if (<= y -1.5e-59)
(* a (/ (+ y t) t_1))
(if (<= y -6.2e-119)
(* z (/ (+ x y) t_1))
(if (<= y 8.6e-30) (/ (+ (* t a) (* x z)) (+ x t)) t_2))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y + (x + t);
double t_2 = (z + a) - b;
double tmp;
if (y <= -1.1e+154) {
tmp = t_2;
} else if (y <= -1.28e+76) {
tmp = z + (y * ((a / x) - (b / x)));
} else if (y <= -1.9e+54) {
tmp = t_2;
} else if (y <= -1.5e-59) {
tmp = a * ((y + t) / t_1);
} else if (y <= -6.2e-119) {
tmp = z * ((x + y) / t_1);
} else if (y <= 8.6e-30) {
tmp = ((t * a) + (x * z)) / (x + 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 = y + (x + t)
t_2 = (z + a) - b
if (y <= (-1.1d+154)) then
tmp = t_2
else if (y <= (-1.28d+76)) then
tmp = z + (y * ((a / x) - (b / x)))
else if (y <= (-1.9d+54)) then
tmp = t_2
else if (y <= (-1.5d-59)) then
tmp = a * ((y + t) / t_1)
else if (y <= (-6.2d-119)) then
tmp = z * ((x + y) / t_1)
else if (y <= 8.6d-30) then
tmp = ((t * a) + (x * z)) / (x + 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 = y + (x + t);
double t_2 = (z + a) - b;
double tmp;
if (y <= -1.1e+154) {
tmp = t_2;
} else if (y <= -1.28e+76) {
tmp = z + (y * ((a / x) - (b / x)));
} else if (y <= -1.9e+54) {
tmp = t_2;
} else if (y <= -1.5e-59) {
tmp = a * ((y + t) / t_1);
} else if (y <= -6.2e-119) {
tmp = z * ((x + y) / t_1);
} else if (y <= 8.6e-30) {
tmp = ((t * a) + (x * z)) / (x + t);
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y + (x + t) t_2 = (z + a) - b tmp = 0 if y <= -1.1e+154: tmp = t_2 elif y <= -1.28e+76: tmp = z + (y * ((a / x) - (b / x))) elif y <= -1.9e+54: tmp = t_2 elif y <= -1.5e-59: tmp = a * ((y + t) / t_1) elif y <= -6.2e-119: tmp = z * ((x + y) / t_1) elif y <= 8.6e-30: tmp = ((t * a) + (x * z)) / (x + t) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(y + Float64(x + t)) t_2 = Float64(Float64(z + a) - b) tmp = 0.0 if (y <= -1.1e+154) tmp = t_2; elseif (y <= -1.28e+76) tmp = Float64(z + Float64(y * Float64(Float64(a / x) - Float64(b / x)))); elseif (y <= -1.9e+54) tmp = t_2; elseif (y <= -1.5e-59) tmp = Float64(a * Float64(Float64(y + t) / t_1)); elseif (y <= -6.2e-119) tmp = Float64(z * Float64(Float64(x + y) / t_1)); elseif (y <= 8.6e-30) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y + (x + t); t_2 = (z + a) - b; tmp = 0.0; if (y <= -1.1e+154) tmp = t_2; elseif (y <= -1.28e+76) tmp = z + (y * ((a / x) - (b / x))); elseif (y <= -1.9e+54) tmp = t_2; elseif (y <= -1.5e-59) tmp = a * ((y + t) / t_1); elseif (y <= -6.2e-119) tmp = z * ((x + y) / t_1); elseif (y <= 8.6e-30) tmp = ((t * a) + (x * z)) / (x + t); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, If[LessEqual[y, -1.1e+154], t$95$2, If[LessEqual[y, -1.28e+76], N[(z + N[(y * N[(N[(a / x), $MachinePrecision] - N[(b / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -1.9e+54], t$95$2, If[LessEqual[y, -1.5e-59], N[(a * N[(N[(y + t), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -6.2e-119], N[(z * N[(N[(x + y), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8.6e-30], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y + \left(x + t\right)\\
t_2 := \left(z + a\right) - b\\
\mathbf{if}\;y \leq -1.1 \cdot 10^{+154}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.28 \cdot 10^{+76}:\\
\;\;\;\;z + y \cdot \left(\frac{a}{x} - \frac{b}{x}\right)\\
\mathbf{elif}\;y \leq -1.9 \cdot 10^{+54}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.5 \cdot 10^{-59}:\\
\;\;\;\;a \cdot \frac{y + t}{t\_1}\\
\mathbf{elif}\;y \leq -6.2 \cdot 10^{-119}:\\
\;\;\;\;z \cdot \frac{x + y}{t\_1}\\
\mathbf{elif}\;y \leq 8.6 \cdot 10^{-30}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y < -1.1000000000000001e154 or -1.27999999999999994e76 < y < -1.9000000000000001e54 or 8.59999999999999932e-30 < y Initial program 39.1%
Taylor expanded in y around inf 77.5%
if -1.1000000000000001e154 < y < -1.27999999999999994e76Initial program 41.1%
Taylor expanded in x around inf 32.9%
associate--l+40.4%
associate-/l*47.7%
associate-/l*47.7%
associate-/l*58.4%
associate-/l*58.4%
Simplified58.4%
Taylor expanded in y around inf 66.2%
if -1.9000000000000001e54 < y < -1.5e-59Initial program 72.0%
Taylor expanded in a around inf 31.3%
associate-/l*51.1%
+-commutative51.1%
+-commutative51.1%
associate-+l+51.1%
+-commutative51.1%
Simplified51.1%
if -1.5e-59 < y < -6.19999999999999956e-119Initial program 85.5%
Taylor expanded in z around inf 55.2%
associate-/l*62.2%
+-commutative62.2%
+-commutative62.2%
+-commutative62.2%
associate-+l+62.2%
+-commutative62.2%
Simplified62.2%
if -6.19999999999999956e-119 < y < 8.59999999999999932e-30Initial program 84.0%
Taylor expanded in y around 0 68.8%
Final simplification70.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (- (+ z a) b)))
(if (<= y -1.06e+154)
t_1
(if (<= y -1.28e+76)
(+ z (* y (- (/ a x) (/ b x))))
(if (<= y -1.9e+54)
t_1
(if (<= y -2.35e-42)
(* a (/ (+ y t) (+ y (+ x t))))
(if (<= y -8.5e-88)
(/ (- (* y z) (* y b)) (+ y t))
(if (<= y 1.4e-29) (/ (+ (* t a) (* x z)) (+ x t)) t_1))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (z + a) - b;
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -1.28e+76) {
tmp = z + (y * ((a / x) - (b / x)));
} else if (y <= -1.9e+54) {
tmp = t_1;
} else if (y <= -2.35e-42) {
tmp = a * ((y + t) / (y + (x + t)));
} else if (y <= -8.5e-88) {
tmp = ((y * z) - (y * b)) / (y + t);
} else if (y <= 1.4e-29) {
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b
if (y <= (-1.06d+154)) then
tmp = t_1
else if (y <= (-1.28d+76)) then
tmp = z + (y * ((a / x) - (b / x)))
else if (y <= (-1.9d+54)) then
tmp = t_1
else if (y <= (-2.35d-42)) then
tmp = a * ((y + t) / (y + (x + t)))
else if (y <= (-8.5d-88)) then
tmp = ((y * z) - (y * b)) / (y + t)
else if (y <= 1.4d-29) then
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b;
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -1.28e+76) {
tmp = z + (y * ((a / x) - (b / x)));
} else if (y <= -1.9e+54) {
tmp = t_1;
} else if (y <= -2.35e-42) {
tmp = a * ((y + t) / (y + (x + t)));
} else if (y <= -8.5e-88) {
tmp = ((y * z) - (y * b)) / (y + t);
} else if (y <= 1.4e-29) {
tmp = ((t * a) + (x * z)) / (x + t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (z + a) - b tmp = 0 if y <= -1.06e+154: tmp = t_1 elif y <= -1.28e+76: tmp = z + (y * ((a / x) - (b / x))) elif y <= -1.9e+54: tmp = t_1 elif y <= -2.35e-42: tmp = a * ((y + t) / (y + (x + t))) elif y <= -8.5e-88: tmp = ((y * z) - (y * b)) / (y + t) elif y <= 1.4e-29: tmp = ((t * a) + (x * z)) / (x + t) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(z + a) - b) tmp = 0.0 if (y <= -1.06e+154) tmp = t_1; elseif (y <= -1.28e+76) tmp = Float64(z + Float64(y * Float64(Float64(a / x) - Float64(b / x)))); elseif (y <= -1.9e+54) tmp = t_1; elseif (y <= -2.35e-42) tmp = Float64(a * Float64(Float64(y + t) / Float64(y + Float64(x + t)))); elseif (y <= -8.5e-88) tmp = Float64(Float64(Float64(y * z) - Float64(y * b)) / Float64(y + t)); elseif (y <= 1.4e-29) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (z + a) - b; tmp = 0.0; if (y <= -1.06e+154) tmp = t_1; elseif (y <= -1.28e+76) tmp = z + (y * ((a / x) - (b / x))); elseif (y <= -1.9e+54) tmp = t_1; elseif (y <= -2.35e-42) tmp = a * ((y + t) / (y + (x + t))); elseif (y <= -8.5e-88) tmp = ((y * z) - (y * b)) / (y + t); elseif (y <= 1.4e-29) tmp = ((t * a) + (x * z)) / (x + t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, If[LessEqual[y, -1.06e+154], t$95$1, If[LessEqual[y, -1.28e+76], N[(z + N[(y * N[(N[(a / x), $MachinePrecision] - N[(b / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -1.9e+54], t$95$1, If[LessEqual[y, -2.35e-42], N[(a * N[(N[(y + t), $MachinePrecision] / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -8.5e-88], N[(N[(N[(y * z), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / N[(y + t), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.4e-29], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(z + a\right) - b\\
\mathbf{if}\;y \leq -1.06 \cdot 10^{+154}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.28 \cdot 10^{+76}:\\
\;\;\;\;z + y \cdot \left(\frac{a}{x} - \frac{b}{x}\right)\\
\mathbf{elif}\;y \leq -1.9 \cdot 10^{+54}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -2.35 \cdot 10^{-42}:\\
\;\;\;\;a \cdot \frac{y + t}{y + \left(x + t\right)}\\
\mathbf{elif}\;y \leq -8.5 \cdot 10^{-88}:\\
\;\;\;\;\frac{y \cdot z - y \cdot b}{y + t}\\
\mathbf{elif}\;y \leq 1.4 \cdot 10^{-29}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.06e154 or -1.27999999999999994e76 < y < -1.9000000000000001e54 or 1.4000000000000001e-29 < y Initial program 39.1%
Taylor expanded in y around inf 77.5%
if -1.06e154 < y < -1.27999999999999994e76Initial program 41.1%
Taylor expanded in x around inf 32.9%
associate--l+40.4%
associate-/l*47.7%
associate-/l*47.7%
associate-/l*58.4%
associate-/l*58.4%
Simplified58.4%
Taylor expanded in y around inf 66.2%
if -1.9000000000000001e54 < y < -2.35e-42Initial program 72.5%
Taylor expanded in a around inf 35.3%
associate-/l*53.4%
+-commutative53.4%
+-commutative53.4%
associate-+l+53.4%
+-commutative53.4%
Simplified53.4%
if -2.35e-42 < y < -8.4999999999999996e-88Initial program 81.0%
Taylor expanded in a around 0 70.9%
+-commutative70.9%
*-commutative70.9%
Simplified70.9%
Taylor expanded in x around 0 51.0%
*-commutative51.0%
Simplified51.0%
if -8.4999999999999996e-88 < y < 1.4000000000000001e-29Initial program 84.0%
Taylor expanded in y around 0 68.9%
Final simplification70.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (- (+ z a) b)))
(if (<= y -1.06e+154)
t_1
(if (<= y -2.85e+75)
(+ z (- (* a (/ y x)) (* b (/ y x))))
(if (<= y -2.45e+54)
t_1
(if (<= y -2.1e-42)
(* a (/ (+ y t) (+ y (+ x t))))
(if (<= y -3.25e-87)
(/ (- (* y z) (* y b)) (+ y t))
(if (<= y 9.2e-30) (/ (+ (* t a) (* x z)) (+ x t)) t_1))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (z + a) - b;
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -2.85e+75) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -2.45e+54) {
tmp = t_1;
} else if (y <= -2.1e-42) {
tmp = a * ((y + t) / (y + (x + t)));
} else if (y <= -3.25e-87) {
tmp = ((y * z) - (y * b)) / (y + t);
} else if (y <= 9.2e-30) {
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b
if (y <= (-1.06d+154)) then
tmp = t_1
else if (y <= (-2.85d+75)) then
tmp = z + ((a * (y / x)) - (b * (y / x)))
else if (y <= (-2.45d+54)) then
tmp = t_1
else if (y <= (-2.1d-42)) then
tmp = a * ((y + t) / (y + (x + t)))
else if (y <= (-3.25d-87)) then
tmp = ((y * z) - (y * b)) / (y + t)
else if (y <= 9.2d-30) then
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b;
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -2.85e+75) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -2.45e+54) {
tmp = t_1;
} else if (y <= -2.1e-42) {
tmp = a * ((y + t) / (y + (x + t)));
} else if (y <= -3.25e-87) {
tmp = ((y * z) - (y * b)) / (y + t);
} else if (y <= 9.2e-30) {
tmp = ((t * a) + (x * z)) / (x + t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (z + a) - b tmp = 0 if y <= -1.06e+154: tmp = t_1 elif y <= -2.85e+75: tmp = z + ((a * (y / x)) - (b * (y / x))) elif y <= -2.45e+54: tmp = t_1 elif y <= -2.1e-42: tmp = a * ((y + t) / (y + (x + t))) elif y <= -3.25e-87: tmp = ((y * z) - (y * b)) / (y + t) elif y <= 9.2e-30: tmp = ((t * a) + (x * z)) / (x + t) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(z + a) - b) tmp = 0.0 if (y <= -1.06e+154) tmp = t_1; elseif (y <= -2.85e+75) tmp = Float64(z + Float64(Float64(a * Float64(y / x)) - Float64(b * Float64(y / x)))); elseif (y <= -2.45e+54) tmp = t_1; elseif (y <= -2.1e-42) tmp = Float64(a * Float64(Float64(y + t) / Float64(y + Float64(x + t)))); elseif (y <= -3.25e-87) tmp = Float64(Float64(Float64(y * z) - Float64(y * b)) / Float64(y + t)); elseif (y <= 9.2e-30) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (z + a) - b; tmp = 0.0; if (y <= -1.06e+154) tmp = t_1; elseif (y <= -2.85e+75) tmp = z + ((a * (y / x)) - (b * (y / x))); elseif (y <= -2.45e+54) tmp = t_1; elseif (y <= -2.1e-42) tmp = a * ((y + t) / (y + (x + t))); elseif (y <= -3.25e-87) tmp = ((y * z) - (y * b)) / (y + t); elseif (y <= 9.2e-30) tmp = ((t * a) + (x * z)) / (x + t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, If[LessEqual[y, -1.06e+154], t$95$1, If[LessEqual[y, -2.85e+75], N[(z + N[(N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision] - N[(b * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -2.45e+54], t$95$1, If[LessEqual[y, -2.1e-42], N[(a * N[(N[(y + t), $MachinePrecision] / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -3.25e-87], N[(N[(N[(y * z), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / N[(y + t), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 9.2e-30], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(z + a\right) - b\\
\mathbf{if}\;y \leq -1.06 \cdot 10^{+154}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -2.85 \cdot 10^{+75}:\\
\;\;\;\;z + \left(a \cdot \frac{y}{x} - b \cdot \frac{y}{x}\right)\\
\mathbf{elif}\;y \leq -2.45 \cdot 10^{+54}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -2.1 \cdot 10^{-42}:\\
\;\;\;\;a \cdot \frac{y + t}{y + \left(x + t\right)}\\
\mathbf{elif}\;y \leq -3.25 \cdot 10^{-87}:\\
\;\;\;\;\frac{y \cdot z - y \cdot b}{y + t}\\
\mathbf{elif}\;y \leq 9.2 \cdot 10^{-30}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.06e154 or -2.8500000000000002e75 < y < -2.45e54 or 9.19999999999999937e-30 < y Initial program 39.1%
Taylor expanded in y around inf 77.5%
if -1.06e154 < y < -2.8500000000000002e75Initial program 41.1%
Taylor expanded in x around inf 32.9%
associate--l+40.4%
associate-/l*47.7%
associate-/l*47.7%
associate-/l*58.4%
associate-/l*58.4%
Simplified58.4%
Taylor expanded in t around 0 48.3%
associate-/l*55.7%
associate-*r/66.4%
Simplified66.4%
if -2.45e54 < y < -2.10000000000000006e-42Initial program 72.5%
Taylor expanded in a around inf 35.3%
associate-/l*53.4%
+-commutative53.4%
+-commutative53.4%
associate-+l+53.4%
+-commutative53.4%
Simplified53.4%
if -2.10000000000000006e-42 < y < -3.2500000000000001e-87Initial program 81.0%
Taylor expanded in a around 0 70.9%
+-commutative70.9%
*-commutative70.9%
Simplified70.9%
Taylor expanded in x around 0 51.0%
*-commutative51.0%
Simplified51.0%
if -3.2500000000000001e-87 < y < 9.19999999999999937e-30Initial program 84.0%
Taylor expanded in y around 0 68.9%
Final simplification70.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (- (+ z a) b)) (t_2 (+ y (+ x t))))
(if (<= y -1.06e+154)
t_1
(if (<= y -1.12e+76)
(+ z (- (* a (/ y x)) (* b (/ y x))))
(if (<= y -2.45e+14)
(/ (* y t_1) t_2)
(if (<= y -8e-36)
(* a (/ (+ y t) t_2))
(if (<= y -4.6e-134)
(/ (- (* z (+ x y)) (* y b)) t_2)
(if (<= y 1.18e-29) (/ (+ (* t a) (* x z)) (+ x t)) t_1))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (z + a) - b;
double t_2 = y + (x + t);
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -1.12e+76) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -2.45e+14) {
tmp = (y * t_1) / t_2;
} else if (y <= -8e-36) {
tmp = a * ((y + t) / t_2);
} else if (y <= -4.6e-134) {
tmp = ((z * (x + y)) - (y * b)) / t_2;
} else if (y <= 1.18e-29) {
tmp = ((t * a) + (x * z)) / (x + 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) :: t_2
real(8) :: tmp
t_1 = (z + a) - b
t_2 = y + (x + t)
if (y <= (-1.06d+154)) then
tmp = t_1
else if (y <= (-1.12d+76)) then
tmp = z + ((a * (y / x)) - (b * (y / x)))
else if (y <= (-2.45d+14)) then
tmp = (y * t_1) / t_2
else if (y <= (-8d-36)) then
tmp = a * ((y + t) / t_2)
else if (y <= (-4.6d-134)) then
tmp = ((z * (x + y)) - (y * b)) / t_2
else if (y <= 1.18d-29) then
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b;
double t_2 = y + (x + t);
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -1.12e+76) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -2.45e+14) {
tmp = (y * t_1) / t_2;
} else if (y <= -8e-36) {
tmp = a * ((y + t) / t_2);
} else if (y <= -4.6e-134) {
tmp = ((z * (x + y)) - (y * b)) / t_2;
} else if (y <= 1.18e-29) {
tmp = ((t * a) + (x * z)) / (x + t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (z + a) - b t_2 = y + (x + t) tmp = 0 if y <= -1.06e+154: tmp = t_1 elif y <= -1.12e+76: tmp = z + ((a * (y / x)) - (b * (y / x))) elif y <= -2.45e+14: tmp = (y * t_1) / t_2 elif y <= -8e-36: tmp = a * ((y + t) / t_2) elif y <= -4.6e-134: tmp = ((z * (x + y)) - (y * b)) / t_2 elif y <= 1.18e-29: tmp = ((t * a) + (x * z)) / (x + t) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(z + a) - b) t_2 = Float64(y + Float64(x + t)) tmp = 0.0 if (y <= -1.06e+154) tmp = t_1; elseif (y <= -1.12e+76) tmp = Float64(z + Float64(Float64(a * Float64(y / x)) - Float64(b * Float64(y / x)))); elseif (y <= -2.45e+14) tmp = Float64(Float64(y * t_1) / t_2); elseif (y <= -8e-36) tmp = Float64(a * Float64(Float64(y + t) / t_2)); elseif (y <= -4.6e-134) tmp = Float64(Float64(Float64(z * Float64(x + y)) - Float64(y * b)) / t_2); elseif (y <= 1.18e-29) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (z + a) - b; t_2 = y + (x + t); tmp = 0.0; if (y <= -1.06e+154) tmp = t_1; elseif (y <= -1.12e+76) tmp = z + ((a * (y / x)) - (b * (y / x))); elseif (y <= -2.45e+14) tmp = (y * t_1) / t_2; elseif (y <= -8e-36) tmp = a * ((y + t) / t_2); elseif (y <= -4.6e-134) tmp = ((z * (x + y)) - (y * b)) / t_2; elseif (y <= 1.18e-29) tmp = ((t * a) + (x * z)) / (x + t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, Block[{t$95$2 = N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.06e+154], t$95$1, If[LessEqual[y, -1.12e+76], N[(z + N[(N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision] - N[(b * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -2.45e+14], N[(N[(y * t$95$1), $MachinePrecision] / t$95$2), $MachinePrecision], If[LessEqual[y, -8e-36], N[(a * N[(N[(y + t), $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -4.6e-134], N[(N[(N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision], If[LessEqual[y, 1.18e-29], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(z + a\right) - b\\
t_2 := y + \left(x + t\right)\\
\mathbf{if}\;y \leq -1.06 \cdot 10^{+154}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.12 \cdot 10^{+76}:\\
\;\;\;\;z + \left(a \cdot \frac{y}{x} - b \cdot \frac{y}{x}\right)\\
\mathbf{elif}\;y \leq -2.45 \cdot 10^{+14}:\\
\;\;\;\;\frac{y \cdot t\_1}{t\_2}\\
\mathbf{elif}\;y \leq -8 \cdot 10^{-36}:\\
\;\;\;\;a \cdot \frac{y + t}{t\_2}\\
\mathbf{elif}\;y \leq -4.6 \cdot 10^{-134}:\\
\;\;\;\;\frac{z \cdot \left(x + y\right) - y \cdot b}{t\_2}\\
\mathbf{elif}\;y \leq 1.18 \cdot 10^{-29}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.06e154 or 1.17999999999999996e-29 < y Initial program 37.7%
Taylor expanded in y around inf 76.7%
if -1.06e154 < y < -1.12000000000000005e76Initial program 41.1%
Taylor expanded in x around inf 32.9%
associate--l+40.4%
associate-/l*47.7%
associate-/l*47.7%
associate-/l*58.4%
associate-/l*58.4%
Simplified58.4%
Taylor expanded in t around 0 48.3%
associate-/l*55.7%
associate-*r/66.4%
Simplified66.4%
if -1.12000000000000005e76 < y < -2.45e14Initial program 76.1%
Taylor expanded in y around inf 67.7%
if -2.45e14 < y < -7.9999999999999995e-36Initial program 65.3%
Taylor expanded in a around inf 47.2%
associate-/l*73.2%
+-commutative73.2%
+-commutative73.2%
associate-+l+73.2%
+-commutative73.2%
Simplified73.2%
if -7.9999999999999995e-36 < y < -4.6000000000000001e-134Initial program 86.3%
Taylor expanded in a around 0 76.6%
+-commutative76.6%
*-commutative76.6%
Simplified76.6%
if -4.6000000000000001e-134 < y < 1.17999999999999996e-29Initial program 83.5%
Taylor expanded in y around 0 69.9%
Final simplification73.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ y (+ x y))))
(if (<= y -7.2e-26)
(* b (- (+ (* (/ a b) t_1) (/ z b)) t_1))
(if (<= y -4.8e-134)
(/ 1.0 (/ (+ t (+ x y)) (- (* z (+ x y)) (* y b))))
(if (<= y 1.1e-60)
(/ (+ (* t a) (* x z)) (+ x t))
(if (<= y 2.15e+68)
(* a (+ 1.0 (- (/ z a) (* (/ b (+ y (+ x t))) (/ y a)))))
(- (+ z a) b)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y / (x + y);
double tmp;
if (y <= -7.2e-26) {
tmp = b * ((((a / b) * t_1) + (z / b)) - t_1);
} else if (y <= -4.8e-134) {
tmp = 1.0 / ((t + (x + y)) / ((z * (x + y)) - (y * b)));
} else if (y <= 1.1e-60) {
tmp = ((t * a) + (x * z)) / (x + t);
} else if (y <= 2.15e+68) {
tmp = a * (1.0 + ((z / a) - ((b / (y + (x + t))) * (y / a))));
} else {
tmp = (z + a) - b;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = y / (x + y)
if (y <= (-7.2d-26)) then
tmp = b * ((((a / b) * t_1) + (z / b)) - t_1)
else if (y <= (-4.8d-134)) then
tmp = 1.0d0 / ((t + (x + y)) / ((z * (x + y)) - (y * b)))
else if (y <= 1.1d-60) then
tmp = ((t * a) + (x * z)) / (x + t)
else if (y <= 2.15d+68) then
tmp = a * (1.0d0 + ((z / a) - ((b / (y + (x + t))) * (y / a))))
else
tmp = (z + a) - b
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y / (x + y);
double tmp;
if (y <= -7.2e-26) {
tmp = b * ((((a / b) * t_1) + (z / b)) - t_1);
} else if (y <= -4.8e-134) {
tmp = 1.0 / ((t + (x + y)) / ((z * (x + y)) - (y * b)));
} else if (y <= 1.1e-60) {
tmp = ((t * a) + (x * z)) / (x + t);
} else if (y <= 2.15e+68) {
tmp = a * (1.0 + ((z / a) - ((b / (y + (x + t))) * (y / a))));
} else {
tmp = (z + a) - b;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y / (x + y) tmp = 0 if y <= -7.2e-26: tmp = b * ((((a / b) * t_1) + (z / b)) - t_1) elif y <= -4.8e-134: tmp = 1.0 / ((t + (x + y)) / ((z * (x + y)) - (y * b))) elif y <= 1.1e-60: tmp = ((t * a) + (x * z)) / (x + t) elif y <= 2.15e+68: tmp = a * (1.0 + ((z / a) - ((b / (y + (x + t))) * (y / a)))) else: tmp = (z + a) - b return tmp
function code(x, y, z, t, a, b) t_1 = Float64(y / Float64(x + y)) tmp = 0.0 if (y <= -7.2e-26) tmp = Float64(b * Float64(Float64(Float64(Float64(a / b) * t_1) + Float64(z / b)) - t_1)); elseif (y <= -4.8e-134) tmp = Float64(1.0 / Float64(Float64(t + Float64(x + y)) / Float64(Float64(z * Float64(x + y)) - Float64(y * b)))); elseif (y <= 1.1e-60) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); elseif (y <= 2.15e+68) tmp = Float64(a * Float64(1.0 + Float64(Float64(z / a) - Float64(Float64(b / Float64(y + Float64(x + t))) * Float64(y / a))))); else tmp = Float64(Float64(z + a) - b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y / (x + y); tmp = 0.0; if (y <= -7.2e-26) tmp = b * ((((a / b) * t_1) + (z / b)) - t_1); elseif (y <= -4.8e-134) tmp = 1.0 / ((t + (x + y)) / ((z * (x + y)) - (y * b))); elseif (y <= 1.1e-60) tmp = ((t * a) + (x * z)) / (x + t); elseif (y <= 2.15e+68) tmp = a * (1.0 + ((z / a) - ((b / (y + (x + t))) * (y / a)))); else tmp = (z + a) - b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y / N[(x + y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -7.2e-26], N[(b * N[(N[(N[(N[(a / b), $MachinePrecision] * t$95$1), $MachinePrecision] + N[(z / b), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -4.8e-134], N[(1.0 / N[(N[(t + N[(x + y), $MachinePrecision]), $MachinePrecision] / N[(N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.1e-60], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.15e+68], N[(a * N[(1.0 + N[(N[(z / a), $MachinePrecision] - N[(N[(b / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y}{x + y}\\
\mathbf{if}\;y \leq -7.2 \cdot 10^{-26}:\\
\;\;\;\;b \cdot \left(\left(\frac{a}{b} \cdot t\_1 + \frac{z}{b}\right) - t\_1\right)\\
\mathbf{elif}\;y \leq -4.8 \cdot 10^{-134}:\\
\;\;\;\;\frac{1}{\frac{t + \left(x + y\right)}{z \cdot \left(x + y\right) - y \cdot b}}\\
\mathbf{elif}\;y \leq 1.1 \cdot 10^{-60}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{elif}\;y \leq 2.15 \cdot 10^{+68}:\\
\;\;\;\;a \cdot \left(1 + \left(\frac{z}{a} - \frac{b}{y + \left(x + t\right)} \cdot \frac{y}{a}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(z + a\right) - b\\
\end{array}
\end{array}
if y < -7.2000000000000003e-26Initial program 39.7%
Taylor expanded in b around inf 51.5%
Taylor expanded in t around 0 56.8%
+-commutative56.8%
mul-1-neg56.8%
unsub-neg56.8%
+-commutative56.8%
times-frac69.8%
+-commutative69.8%
+-commutative69.8%
Simplified69.8%
if -7.2000000000000003e-26 < y < -4.80000000000000019e-134Initial program 84.0%
Taylor expanded in a around 0 71.3%
+-commutative71.3%
*-commutative71.3%
Simplified71.3%
clear-num71.4%
inv-pow71.4%
+-commutative71.4%
+-commutative71.4%
+-commutative71.4%
*-commutative71.4%
Applied egg-rr71.4%
unpow-171.4%
+-commutative71.4%
+-commutative71.4%
associate-+r+71.4%
+-commutative71.4%
Simplified71.4%
if -4.80000000000000019e-134 < y < 1.0999999999999999e-60Initial program 83.7%
Taylor expanded in y around 0 70.4%
if 1.0999999999999999e-60 < y < 2.1500000000000001e68Initial program 74.2%
Taylor expanded in a around inf 81.6%
associate--l+81.6%
+-commutative81.6%
+-commutative81.6%
associate-+l+81.6%
+-commutative81.6%
Simplified88.5%
Taylor expanded in x around inf 76.1%
Taylor expanded in t around inf 77.2%
if 2.1500000000000001e68 < y Initial program 35.4%
Taylor expanded in y around inf 82.3%
Final simplification73.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (- (+ z a) b)))
(if (<= y -1.06e+154)
t_1
(if (<= y -2.65e+75)
(+ z (- (* a (/ y x)) (* b (/ y x))))
(if (<= y -2.2e+69)
t_1
(if (<= y -1.38e-86)
(/ (* y t_1) (+ y (+ x t)))
(if (<= y 8.6e-30) (/ (+ (* t a) (* x z)) (+ x t)) t_1)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (z + a) - b;
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -2.65e+75) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -2.2e+69) {
tmp = t_1;
} else if (y <= -1.38e-86) {
tmp = (y * t_1) / (y + (x + t));
} else if (y <= 8.6e-30) {
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b
if (y <= (-1.06d+154)) then
tmp = t_1
else if (y <= (-2.65d+75)) then
tmp = z + ((a * (y / x)) - (b * (y / x)))
else if (y <= (-2.2d+69)) then
tmp = t_1
else if (y <= (-1.38d-86)) then
tmp = (y * t_1) / (y + (x + t))
else if (y <= 8.6d-30) then
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b;
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -2.65e+75) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -2.2e+69) {
tmp = t_1;
} else if (y <= -1.38e-86) {
tmp = (y * t_1) / (y + (x + t));
} else if (y <= 8.6e-30) {
tmp = ((t * a) + (x * z)) / (x + t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (z + a) - b tmp = 0 if y <= -1.06e+154: tmp = t_1 elif y <= -2.65e+75: tmp = z + ((a * (y / x)) - (b * (y / x))) elif y <= -2.2e+69: tmp = t_1 elif y <= -1.38e-86: tmp = (y * t_1) / (y + (x + t)) elif y <= 8.6e-30: tmp = ((t * a) + (x * z)) / (x + t) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(z + a) - b) tmp = 0.0 if (y <= -1.06e+154) tmp = t_1; elseif (y <= -2.65e+75) tmp = Float64(z + Float64(Float64(a * Float64(y / x)) - Float64(b * Float64(y / x)))); elseif (y <= -2.2e+69) tmp = t_1; elseif (y <= -1.38e-86) tmp = Float64(Float64(y * t_1) / Float64(y + Float64(x + t))); elseif (y <= 8.6e-30) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (z + a) - b; tmp = 0.0; if (y <= -1.06e+154) tmp = t_1; elseif (y <= -2.65e+75) tmp = z + ((a * (y / x)) - (b * (y / x))); elseif (y <= -2.2e+69) tmp = t_1; elseif (y <= -1.38e-86) tmp = (y * t_1) / (y + (x + t)); elseif (y <= 8.6e-30) tmp = ((t * a) + (x * z)) / (x + t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, If[LessEqual[y, -1.06e+154], t$95$1, If[LessEqual[y, -2.65e+75], N[(z + N[(N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision] - N[(b * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -2.2e+69], t$95$1, If[LessEqual[y, -1.38e-86], N[(N[(y * t$95$1), $MachinePrecision] / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8.6e-30], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(z + a\right) - b\\
\mathbf{if}\;y \leq -1.06 \cdot 10^{+154}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -2.65 \cdot 10^{+75}:\\
\;\;\;\;z + \left(a \cdot \frac{y}{x} - b \cdot \frac{y}{x}\right)\\
\mathbf{elif}\;y \leq -2.2 \cdot 10^{+69}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.38 \cdot 10^{-86}:\\
\;\;\;\;\frac{y \cdot t\_1}{y + \left(x + t\right)}\\
\mathbf{elif}\;y \leq 8.6 \cdot 10^{-30}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.06e154 or -2.6499999999999999e75 < y < -2.2000000000000002e69 or 8.59999999999999932e-30 < y Initial program 38.6%
Taylor expanded in y around inf 77.3%
if -1.06e154 < y < -2.6499999999999999e75Initial program 41.1%
Taylor expanded in x around inf 32.9%
associate--l+40.4%
associate-/l*47.7%
associate-/l*47.7%
associate-/l*58.4%
associate-/l*58.4%
Simplified58.4%
Taylor expanded in t around 0 48.3%
associate-/l*55.7%
associate-*r/66.4%
Simplified66.4%
if -2.2000000000000002e69 < y < -1.38000000000000008e-86Initial program 76.0%
Taylor expanded in y around inf 55.9%
if -1.38000000000000008e-86 < y < 8.59999999999999932e-30Initial program 84.0%
Taylor expanded in y around 0 68.9%
Final simplification70.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (- (+ z a) b)))
(if (<= y -1.06e+154)
t_1
(if (<= y -6.5e+74)
(+ z (- (* a (/ y x)) (* b (/ y x))))
(if (<= y -9e+55)
t_1
(if (<= y -3.7e-87)
(/ (- (* (+ y t) a) (* y b)) (+ y (+ x t)))
(if (<= y 9.2e-30) (/ (+ (* t a) (* x z)) (+ x t)) t_1)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (z + a) - b;
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -6.5e+74) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -9e+55) {
tmp = t_1;
} else if (y <= -3.7e-87) {
tmp = (((y + t) * a) - (y * b)) / (y + (x + t));
} else if (y <= 9.2e-30) {
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b
if (y <= (-1.06d+154)) then
tmp = t_1
else if (y <= (-6.5d+74)) then
tmp = z + ((a * (y / x)) - (b * (y / x)))
else if (y <= (-9d+55)) then
tmp = t_1
else if (y <= (-3.7d-87)) then
tmp = (((y + t) * a) - (y * b)) / (y + (x + t))
else if (y <= 9.2d-30) then
tmp = ((t * a) + (x * z)) / (x + 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 + a) - b;
double tmp;
if (y <= -1.06e+154) {
tmp = t_1;
} else if (y <= -6.5e+74) {
tmp = z + ((a * (y / x)) - (b * (y / x)));
} else if (y <= -9e+55) {
tmp = t_1;
} else if (y <= -3.7e-87) {
tmp = (((y + t) * a) - (y * b)) / (y + (x + t));
} else if (y <= 9.2e-30) {
tmp = ((t * a) + (x * z)) / (x + t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (z + a) - b tmp = 0 if y <= -1.06e+154: tmp = t_1 elif y <= -6.5e+74: tmp = z + ((a * (y / x)) - (b * (y / x))) elif y <= -9e+55: tmp = t_1 elif y <= -3.7e-87: tmp = (((y + t) * a) - (y * b)) / (y + (x + t)) elif y <= 9.2e-30: tmp = ((t * a) + (x * z)) / (x + t) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(z + a) - b) tmp = 0.0 if (y <= -1.06e+154) tmp = t_1; elseif (y <= -6.5e+74) tmp = Float64(z + Float64(Float64(a * Float64(y / x)) - Float64(b * Float64(y / x)))); elseif (y <= -9e+55) tmp = t_1; elseif (y <= -3.7e-87) tmp = Float64(Float64(Float64(Float64(y + t) * a) - Float64(y * b)) / Float64(y + Float64(x + t))); elseif (y <= 9.2e-30) tmp = Float64(Float64(Float64(t * a) + Float64(x * z)) / Float64(x + t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (z + a) - b; tmp = 0.0; if (y <= -1.06e+154) tmp = t_1; elseif (y <= -6.5e+74) tmp = z + ((a * (y / x)) - (b * (y / x))); elseif (y <= -9e+55) tmp = t_1; elseif (y <= -3.7e-87) tmp = (((y + t) * a) - (y * b)) / (y + (x + t)); elseif (y <= 9.2e-30) tmp = ((t * a) + (x * z)) / (x + t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, If[LessEqual[y, -1.06e+154], t$95$1, If[LessEqual[y, -6.5e+74], N[(z + N[(N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision] - N[(b * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -9e+55], t$95$1, If[LessEqual[y, -3.7e-87], N[(N[(N[(N[(y + t), $MachinePrecision] * a), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision] / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 9.2e-30], N[(N[(N[(t * a), $MachinePrecision] + N[(x * z), $MachinePrecision]), $MachinePrecision] / N[(x + t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(z + a\right) - b\\
\mathbf{if}\;y \leq -1.06 \cdot 10^{+154}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -6.5 \cdot 10^{+74}:\\
\;\;\;\;z + \left(a \cdot \frac{y}{x} - b \cdot \frac{y}{x}\right)\\
\mathbf{elif}\;y \leq -9 \cdot 10^{+55}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -3.7 \cdot 10^{-87}:\\
\;\;\;\;\frac{\left(y + t\right) \cdot a - y \cdot b}{y + \left(x + t\right)}\\
\mathbf{elif}\;y \leq 9.2 \cdot 10^{-30}:\\
\;\;\;\;\frac{t \cdot a + x \cdot z}{x + t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.06e154 or -6.49999999999999962e74 < y < -8.99999999999999996e55 or 9.19999999999999937e-30 < y Initial program 39.1%
Taylor expanded in y around inf 77.5%
if -1.06e154 < y < -6.49999999999999962e74Initial program 41.1%
Taylor expanded in x around inf 32.9%
associate--l+40.4%
associate-/l*47.7%
associate-/l*47.7%
associate-/l*58.4%
associate-/l*58.4%
Simplified58.4%
Taylor expanded in t around 0 48.3%
associate-/l*55.7%
associate-*r/66.4%
Simplified66.4%
if -8.99999999999999996e55 < y < -3.7000000000000002e-87Initial program 75.3%
Taylor expanded in z around 0 56.6%
*-commutative56.6%
Simplified56.6%
if -3.7000000000000002e-87 < y < 9.19999999999999937e-30Initial program 84.0%
Taylor expanded in y around 0 68.9%
Final simplification71.0%
(FPCore (x y z t a b) :precision binary64 (if (<= x -1.35e+73) (- z (* y (/ b x))) (if (<= x 2.8e+136) (- (+ z a) b) (* z (/ (+ x y) (+ y (+ x t)))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (x <= -1.35e+73) {
tmp = z - (y * (b / x));
} else if (x <= 2.8e+136) {
tmp = (z + a) - b;
} else {
tmp = z * ((x + y) / (y + (x + 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 (x <= (-1.35d+73)) then
tmp = z - (y * (b / x))
else if (x <= 2.8d+136) then
tmp = (z + a) - b
else
tmp = z * ((x + y) / (y + (x + 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 (x <= -1.35e+73) {
tmp = z - (y * (b / x));
} else if (x <= 2.8e+136) {
tmp = (z + a) - b;
} else {
tmp = z * ((x + y) / (y + (x + t)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if x <= -1.35e+73: tmp = z - (y * (b / x)) elif x <= 2.8e+136: tmp = (z + a) - b else: tmp = z * ((x + y) / (y + (x + t))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (x <= -1.35e+73) tmp = Float64(z - Float64(y * Float64(b / x))); elseif (x <= 2.8e+136) tmp = Float64(Float64(z + a) - b); else tmp = Float64(z * Float64(Float64(x + y) / Float64(y + Float64(x + t)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (x <= -1.35e+73) tmp = z - (y * (b / x)); elseif (x <= 2.8e+136) tmp = (z + a) - b; else tmp = z * ((x + y) / (y + (x + t))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[x, -1.35e+73], N[(z - N[(y * N[(b / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 2.8e+136], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision], N[(z * N[(N[(x + y), $MachinePrecision] / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.35 \cdot 10^{+73}:\\
\;\;\;\;z - y \cdot \frac{b}{x}\\
\mathbf{elif}\;x \leq 2.8 \cdot 10^{+136}:\\
\;\;\;\;\left(z + a\right) - b\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{x + y}{y + \left(x + t\right)}\\
\end{array}
\end{array}
if x < -1.35e73Initial program 55.8%
Taylor expanded in x around inf 58.0%
associate--l+58.1%
associate-/l*57.9%
associate-/l*58.2%
associate-/l*59.1%
associate-/l*71.5%
Simplified71.5%
Taylor expanded in y around inf 70.4%
Taylor expanded in a around 0 64.9%
neg-mul-164.9%
distribute-neg-frac64.9%
Simplified64.9%
if -1.35e73 < x < 2.8000000000000002e136Initial program 65.2%
Taylor expanded in y around inf 63.0%
if 2.8000000000000002e136 < x Initial program 47.2%
Taylor expanded in z around inf 33.2%
associate-/l*68.6%
+-commutative68.6%
+-commutative68.6%
+-commutative68.6%
associate-+l+68.6%
+-commutative68.6%
Simplified68.6%
Final simplification64.0%
(FPCore (x y z t a b) :precision binary64 (if (<= x -1.05e+72) (+ z (* y (- (/ a x) (/ b x)))) (if (<= x 3.5e+135) (- (+ z a) b) (* z (/ (+ x y) (+ y (+ x t)))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (x <= -1.05e+72) {
tmp = z + (y * ((a / x) - (b / x)));
} else if (x <= 3.5e+135) {
tmp = (z + a) - b;
} else {
tmp = z * ((x + y) / (y + (x + 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 (x <= (-1.05d+72)) then
tmp = z + (y * ((a / x) - (b / x)))
else if (x <= 3.5d+135) then
tmp = (z + a) - b
else
tmp = z * ((x + y) / (y + (x + 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 (x <= -1.05e+72) {
tmp = z + (y * ((a / x) - (b / x)));
} else if (x <= 3.5e+135) {
tmp = (z + a) - b;
} else {
tmp = z * ((x + y) / (y + (x + t)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if x <= -1.05e+72: tmp = z + (y * ((a / x) - (b / x))) elif x <= 3.5e+135: tmp = (z + a) - b else: tmp = z * ((x + y) / (y + (x + t))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (x <= -1.05e+72) tmp = Float64(z + Float64(y * Float64(Float64(a / x) - Float64(b / x)))); elseif (x <= 3.5e+135) tmp = Float64(Float64(z + a) - b); else tmp = Float64(z * Float64(Float64(x + y) / Float64(y + Float64(x + t)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (x <= -1.05e+72) tmp = z + (y * ((a / x) - (b / x))); elseif (x <= 3.5e+135) tmp = (z + a) - b; else tmp = z * ((x + y) / (y + (x + t))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[x, -1.05e+72], N[(z + N[(y * N[(N[(a / x), $MachinePrecision] - N[(b / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 3.5e+135], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision], N[(z * N[(N[(x + y), $MachinePrecision] / N[(y + N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.05 \cdot 10^{+72}:\\
\;\;\;\;z + y \cdot \left(\frac{a}{x} - \frac{b}{x}\right)\\
\mathbf{elif}\;x \leq 3.5 \cdot 10^{+135}:\\
\;\;\;\;\left(z + a\right) - b\\
\mathbf{else}:\\
\;\;\;\;z \cdot \frac{x + y}{y + \left(x + t\right)}\\
\end{array}
\end{array}
if x < -1.0500000000000001e72Initial program 55.8%
Taylor expanded in x around inf 58.0%
associate--l+58.1%
associate-/l*57.9%
associate-/l*58.2%
associate-/l*59.1%
associate-/l*71.5%
Simplified71.5%
Taylor expanded in y around inf 70.4%
if -1.0500000000000001e72 < x < 3.5000000000000003e135Initial program 65.2%
Taylor expanded in y around inf 63.0%
if 3.5000000000000003e135 < x Initial program 47.2%
Taylor expanded in z around inf 33.2%
associate-/l*68.6%
+-commutative68.6%
+-commutative68.6%
+-commutative68.6%
associate-+l+68.6%
+-commutative68.6%
Simplified68.6%
Final simplification65.1%
(FPCore (x y z t a b) :precision binary64 (if (or (<= x -2.2e+162) (not (<= x 7.5e+139))) (+ z (* y (/ a x))) (- (+ z a) b)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((x <= -2.2e+162) || !(x <= 7.5e+139)) {
tmp = z + (y * (a / x));
} else {
tmp = (z + a) - 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 ((x <= (-2.2d+162)) .or. (.not. (x <= 7.5d+139))) then
tmp = z + (y * (a / x))
else
tmp = (z + a) - 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 ((x <= -2.2e+162) || !(x <= 7.5e+139)) {
tmp = z + (y * (a / x));
} else {
tmp = (z + a) - b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (x <= -2.2e+162) or not (x <= 7.5e+139): tmp = z + (y * (a / x)) else: tmp = (z + a) - b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((x <= -2.2e+162) || !(x <= 7.5e+139)) tmp = Float64(z + Float64(y * Float64(a / x))); else tmp = Float64(Float64(z + a) - b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((x <= -2.2e+162) || ~((x <= 7.5e+139))) tmp = z + (y * (a / x)); else tmp = (z + a) - b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[x, -2.2e+162], N[Not[LessEqual[x, 7.5e+139]], $MachinePrecision]], N[(z + N[(y * N[(a / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2.2 \cdot 10^{+162} \lor \neg \left(x \leq 7.5 \cdot 10^{+139}\right):\\
\;\;\;\;z + y \cdot \frac{a}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(z + a\right) - b\\
\end{array}
\end{array}
if x < -2.2000000000000002e162 or 7.49999999999999992e139 < x Initial program 56.0%
Taylor expanded in x around inf 60.2%
associate--l+60.2%
associate-/l*62.1%
associate-/l*62.4%
associate-/l*62.6%
associate-/l*76.9%
Simplified76.9%
Taylor expanded in y around inf 72.7%
Taylor expanded in a around inf 59.2%
if -2.2000000000000002e162 < x < 7.49999999999999992e139Initial program 63.2%
Taylor expanded in y around inf 61.8%
Final simplification61.1%
(FPCore (x y z t a b) :precision binary64 (if (or (<= x -1.5e+71) (not (<= x 2.4e+125))) (- z (* b (/ y x))) (- (+ z a) b)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((x <= -1.5e+71) || !(x <= 2.4e+125)) {
tmp = z - (b * (y / x));
} else {
tmp = (z + a) - 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 ((x <= (-1.5d+71)) .or. (.not. (x <= 2.4d+125))) then
tmp = z - (b * (y / x))
else
tmp = (z + a) - 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 ((x <= -1.5e+71) || !(x <= 2.4e+125)) {
tmp = z - (b * (y / x));
} else {
tmp = (z + a) - b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (x <= -1.5e+71) or not (x <= 2.4e+125): tmp = z - (b * (y / x)) else: tmp = (z + a) - b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((x <= -1.5e+71) || !(x <= 2.4e+125)) tmp = Float64(z - Float64(b * Float64(y / x))); else tmp = Float64(Float64(z + a) - b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((x <= -1.5e+71) || ~((x <= 2.4e+125))) tmp = z - (b * (y / x)); else tmp = (z + a) - b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[x, -1.5e+71], N[Not[LessEqual[x, 2.4e+125]], $MachinePrecision]], N[(z - N[(b * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.5 \cdot 10^{+71} \lor \neg \left(x \leq 2.4 \cdot 10^{+125}\right):\\
\;\;\;\;z - b \cdot \frac{y}{x}\\
\mathbf{else}:\\
\;\;\;\;\left(z + a\right) - b\\
\end{array}
\end{array}
if x < -1.50000000000000006e71 or 2.4e125 < x Initial program 52.9%
Taylor expanded in x around inf 56.1%
associate--l+56.1%
associate-/l*57.7%
associate-/l*58.0%
associate-/l*58.7%
associate-/l*71.7%
Simplified71.7%
Taylor expanded in y around inf 69.7%
Taylor expanded in a around 0 55.8%
mul-1-neg55.8%
unsub-neg55.8%
associate-/l*62.6%
Simplified62.6%
if -1.50000000000000006e71 < x < 2.4e125Initial program 65.2%
Taylor expanded in y around inf 63.0%
Final simplification62.9%
(FPCore (x y z t a b) :precision binary64 (if (<= x -7.5e+162) z (if (<= x 7.5e+134) (- (+ z a) b) (* x (/ z (+ x t))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (x <= -7.5e+162) {
tmp = z;
} else if (x <= 7.5e+134) {
tmp = (z + a) - b;
} else {
tmp = x * (z / (x + 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 (x <= (-7.5d+162)) then
tmp = z
else if (x <= 7.5d+134) then
tmp = (z + a) - b
else
tmp = x * (z / (x + 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 (x <= -7.5e+162) {
tmp = z;
} else if (x <= 7.5e+134) {
tmp = (z + a) - b;
} else {
tmp = x * (z / (x + t));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if x <= -7.5e+162: tmp = z elif x <= 7.5e+134: tmp = (z + a) - b else: tmp = x * (z / (x + t)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (x <= -7.5e+162) tmp = z; elseif (x <= 7.5e+134) tmp = Float64(Float64(z + a) - b); else tmp = Float64(x * Float64(z / Float64(x + t))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (x <= -7.5e+162) tmp = z; elseif (x <= 7.5e+134) tmp = (z + a) - b; else tmp = x * (z / (x + t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[x, -7.5e+162], z, If[LessEqual[x, 7.5e+134], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision], N[(x * N[(z / N[(x + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -7.5 \cdot 10^{+162}:\\
\;\;\;\;z\\
\mathbf{elif}\;x \leq 7.5 \cdot 10^{+134}:\\
\;\;\;\;\left(z + a\right) - b\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{z}{x + t}\\
\end{array}
\end{array}
if x < -7.50000000000000033e162Initial program 62.0%
Taylor expanded in x around inf 51.0%
if -7.50000000000000033e162 < x < 7.5000000000000001e134Initial program 63.2%
Taylor expanded in y around inf 61.8%
if 7.5000000000000001e134 < x Initial program 47.2%
Taylor expanded in a around 0 36.5%
+-commutative36.5%
*-commutative36.5%
Simplified36.5%
Taylor expanded in y around 0 33.2%
+-commutative33.2%
Simplified33.2%
associate-/l*54.9%
+-commutative54.9%
Applied egg-rr54.9%
Final simplification59.4%
(FPCore (x y z t a b) :precision binary64 (if (<= x -1.3e+73) (- z (* y (/ b x))) (if (<= x 1.3e+140) (- (+ z a) b) (- z (* b (/ y x))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (x <= -1.3e+73) {
tmp = z - (y * (b / x));
} else if (x <= 1.3e+140) {
tmp = (z + a) - b;
} else {
tmp = z - (b * (y / x));
}
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 (x <= (-1.3d+73)) then
tmp = z - (y * (b / x))
else if (x <= 1.3d+140) then
tmp = (z + a) - b
else
tmp = z - (b * (y / x))
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 (x <= -1.3e+73) {
tmp = z - (y * (b / x));
} else if (x <= 1.3e+140) {
tmp = (z + a) - b;
} else {
tmp = z - (b * (y / x));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if x <= -1.3e+73: tmp = z - (y * (b / x)) elif x <= 1.3e+140: tmp = (z + a) - b else: tmp = z - (b * (y / x)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (x <= -1.3e+73) tmp = Float64(z - Float64(y * Float64(b / x))); elseif (x <= 1.3e+140) tmp = Float64(Float64(z + a) - b); else tmp = Float64(z - Float64(b * Float64(y / x))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (x <= -1.3e+73) tmp = z - (y * (b / x)); elseif (x <= 1.3e+140) tmp = (z + a) - b; else tmp = z - (b * (y / x)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[x, -1.3e+73], N[(z - N[(y * N[(b / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.3e+140], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision], N[(z - N[(b * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.3 \cdot 10^{+73}:\\
\;\;\;\;z - y \cdot \frac{b}{x}\\
\mathbf{elif}\;x \leq 1.3 \cdot 10^{+140}:\\
\;\;\;\;\left(z + a\right) - b\\
\mathbf{else}:\\
\;\;\;\;z - b \cdot \frac{y}{x}\\
\end{array}
\end{array}
if x < -1.3e73Initial program 55.8%
Taylor expanded in x around inf 58.0%
associate--l+58.1%
associate-/l*57.9%
associate-/l*58.2%
associate-/l*59.1%
associate-/l*71.5%
Simplified71.5%
Taylor expanded in y around inf 70.4%
Taylor expanded in a around 0 64.9%
neg-mul-164.9%
distribute-neg-frac64.9%
Simplified64.9%
if -1.3e73 < x < 1.3000000000000001e140Initial program 65.2%
Taylor expanded in y around inf 63.0%
if 1.3000000000000001e140 < x Initial program 47.2%
Taylor expanded in x around inf 52.2%
associate--l+52.2%
associate-/l*57.2%
associate-/l*57.7%
associate-/l*58.0%
associate-/l*72.2%
Simplified72.2%
Taylor expanded in y around inf 68.3%
Taylor expanded in a around 0 53.3%
mul-1-neg53.3%
unsub-neg53.3%
associate-/l*61.5%
Simplified61.5%
Final simplification63.2%
(FPCore (x y z t a b) :precision binary64 (if (<= x -5.9e+162) z (if (<= x 5e+135) (- (+ z a) b) z)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (x <= -5.9e+162) {
tmp = z;
} else if (x <= 5e+135) {
tmp = (z + a) - b;
} 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 (x <= (-5.9d+162)) then
tmp = z
else if (x <= 5d+135) then
tmp = (z + a) - b
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 (x <= -5.9e+162) {
tmp = z;
} else if (x <= 5e+135) {
tmp = (z + a) - b;
} else {
tmp = z;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if x <= -5.9e+162: tmp = z elif x <= 5e+135: tmp = (z + a) - b else: tmp = z return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (x <= -5.9e+162) tmp = z; elseif (x <= 5e+135) tmp = Float64(Float64(z + a) - b); else tmp = z; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (x <= -5.9e+162) tmp = z; elseif (x <= 5e+135) tmp = (z + a) - b; else tmp = z; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[x, -5.9e+162], z, If[LessEqual[x, 5e+135], N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision], z]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5.9 \cdot 10^{+162}:\\
\;\;\;\;z\\
\mathbf{elif}\;x \leq 5 \cdot 10^{+135}:\\
\;\;\;\;\left(z + a\right) - b\\
\mathbf{else}:\\
\;\;\;\;z\\
\end{array}
\end{array}
if x < -5.90000000000000027e162 or 5.00000000000000029e135 < x Initial program 56.0%
Taylor expanded in x around inf 52.3%
if -5.90000000000000027e162 < x < 5.00000000000000029e135Initial program 63.2%
Taylor expanded in y around inf 61.8%
Final simplification59.3%
(FPCore (x y z t a b) :precision binary64 (if (<= a -5.4e-48) a (if (<= a 1.4e+20) (- z b) a)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (a <= -5.4e-48) {
tmp = a;
} else if (a <= 1.4e+20) {
tmp = z - b;
} else {
tmp = 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 (a <= (-5.4d-48)) then
tmp = a
else if (a <= 1.4d+20) then
tmp = z - b
else
tmp = 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 (a <= -5.4e-48) {
tmp = a;
} else if (a <= 1.4e+20) {
tmp = z - b;
} else {
tmp = a;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if a <= -5.4e-48: tmp = a elif a <= 1.4e+20: tmp = z - b else: tmp = a return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (a <= -5.4e-48) tmp = a; elseif (a <= 1.4e+20) tmp = Float64(z - b); else tmp = a; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (a <= -5.4e-48) tmp = a; elseif (a <= 1.4e+20) tmp = z - b; else tmp = a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[a, -5.4e-48], a, If[LessEqual[a, 1.4e+20], N[(z - b), $MachinePrecision], a]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -5.4 \cdot 10^{-48}:\\
\;\;\;\;a\\
\mathbf{elif}\;a \leq 1.4 \cdot 10^{+20}:\\
\;\;\;\;z - b\\
\mathbf{else}:\\
\;\;\;\;a\\
\end{array}
\end{array}
if a < -5.40000000000000023e-48 or 1.4e20 < a Initial program 53.7%
Taylor expanded in t around inf 52.3%
if -5.40000000000000023e-48 < a < 1.4e20Initial program 68.2%
Taylor expanded in a around 0 59.2%
+-commutative59.2%
*-commutative59.2%
Simplified59.2%
Taylor expanded in y around inf 50.1%
Final simplification51.1%
(FPCore (x y z t a b) :precision binary64 (if (<= a -2.55e-48) a (if (<= a 7500000000.0) z a)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (a <= -2.55e-48) {
tmp = a;
} else if (a <= 7500000000.0) {
tmp = z;
} else {
tmp = 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 (a <= (-2.55d-48)) then
tmp = a
else if (a <= 7500000000.0d0) then
tmp = z
else
tmp = 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 (a <= -2.55e-48) {
tmp = a;
} else if (a <= 7500000000.0) {
tmp = z;
} else {
tmp = a;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if a <= -2.55e-48: tmp = a elif a <= 7500000000.0: tmp = z else: tmp = a return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (a <= -2.55e-48) tmp = a; elseif (a <= 7500000000.0) tmp = z; else tmp = a; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (a <= -2.55e-48) tmp = a; elseif (a <= 7500000000.0) tmp = z; else tmp = a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[a, -2.55e-48], a, If[LessEqual[a, 7500000000.0], z, a]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2.55 \cdot 10^{-48}:\\
\;\;\;\;a\\
\mathbf{elif}\;a \leq 7500000000:\\
\;\;\;\;z\\
\mathbf{else}:\\
\;\;\;\;a\\
\end{array}
\end{array}
if a < -2.55000000000000006e-48 or 7.5e9 < a Initial program 53.3%
Taylor expanded in t around inf 51.9%
if -2.55000000000000006e-48 < a < 7.5e9Initial program 68.7%
Taylor expanded in x around inf 43.9%
Final simplification47.7%
(FPCore (x y z t a b) :precision binary64 a)
double code(double x, double y, double z, double t, double a, double b) {
return a;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = a
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return a;
}
def code(x, y, z, t, a, b): return a
function code(x, y, z, t, a, b) return a end
function tmp = code(x, y, z, t, a, b) tmp = a; end
code[x_, y_, z_, t_, a_, b_] := a
\begin{array}{l}
\\
a
\end{array}
Initial program 61.3%
Taylor expanded in t around inf 29.9%
Final simplification29.9%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ (+ x t) y))
(t_2 (- (+ (* (+ x y) z) (* (+ t y) a)) (* y b)))
(t_3 (/ t_2 t_1))
(t_4 (- (+ z a) b)))
(if (< t_3 -3.5813117084150564e+153)
t_4
(if (< t_3 1.2285964308315609e+82) (/ 1.0 (/ t_1 t_2)) t_4))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x + t) + y;
double t_2 = (((x + y) * z) + ((t + y) * a)) - (y * b);
double t_3 = t_2 / t_1;
double t_4 = (z + a) - b;
double tmp;
if (t_3 < -3.5813117084150564e+153) {
tmp = t_4;
} else if (t_3 < 1.2285964308315609e+82) {
tmp = 1.0 / (t_1 / t_2);
} else {
tmp = t_4;
}
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) :: t_3
real(8) :: t_4
real(8) :: tmp
t_1 = (x + t) + y
t_2 = (((x + y) * z) + ((t + y) * a)) - (y * b)
t_3 = t_2 / t_1
t_4 = (z + a) - b
if (t_3 < (-3.5813117084150564d+153)) then
tmp = t_4
else if (t_3 < 1.2285964308315609d+82) then
tmp = 1.0d0 / (t_1 / t_2)
else
tmp = t_4
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 + t) + y;
double t_2 = (((x + y) * z) + ((t + y) * a)) - (y * b);
double t_3 = t_2 / t_1;
double t_4 = (z + a) - b;
double tmp;
if (t_3 < -3.5813117084150564e+153) {
tmp = t_4;
} else if (t_3 < 1.2285964308315609e+82) {
tmp = 1.0 / (t_1 / t_2);
} else {
tmp = t_4;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x + t) + y t_2 = (((x + y) * z) + ((t + y) * a)) - (y * b) t_3 = t_2 / t_1 t_4 = (z + a) - b tmp = 0 if t_3 < -3.5813117084150564e+153: tmp = t_4 elif t_3 < 1.2285964308315609e+82: tmp = 1.0 / (t_1 / t_2) else: tmp = t_4 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x + t) + y) t_2 = Float64(Float64(Float64(Float64(x + y) * z) + Float64(Float64(t + y) * a)) - Float64(y * b)) t_3 = Float64(t_2 / t_1) t_4 = Float64(Float64(z + a) - b) tmp = 0.0 if (t_3 < -3.5813117084150564e+153) tmp = t_4; elseif (t_3 < 1.2285964308315609e+82) tmp = Float64(1.0 / Float64(t_1 / t_2)); else tmp = t_4; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x + t) + y; t_2 = (((x + y) * z) + ((t + y) * a)) - (y * b); t_3 = t_2 / t_1; t_4 = (z + a) - b; tmp = 0.0; if (t_3 < -3.5813117084150564e+153) tmp = t_4; elseif (t_3 < 1.2285964308315609e+82) tmp = 1.0 / (t_1 / t_2); else tmp = t_4; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x + t), $MachinePrecision] + y), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(N[(x + y), $MachinePrecision] * z), $MachinePrecision] + N[(N[(t + y), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] - N[(y * b), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(t$95$2 / t$95$1), $MachinePrecision]}, Block[{t$95$4 = N[(N[(z + a), $MachinePrecision] - b), $MachinePrecision]}, If[Less[t$95$3, -3.5813117084150564e+153], t$95$4, If[Less[t$95$3, 1.2285964308315609e+82], N[(1.0 / N[(t$95$1 / t$95$2), $MachinePrecision]), $MachinePrecision], t$95$4]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x + t\right) + y\\
t_2 := \left(\left(x + y\right) \cdot z + \left(t + y\right) \cdot a\right) - y \cdot b\\
t_3 := \frac{t\_2}{t\_1}\\
t_4 := \left(z + a\right) - b\\
\mathbf{if}\;t\_3 < -3.5813117084150564 \cdot 10^{+153}:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;t\_3 < 1.2285964308315609 \cdot 10^{+82}:\\
\;\;\;\;\frac{1}{\frac{t\_1}{t\_2}}\\
\mathbf{else}:\\
\;\;\;\;t\_4\\
\end{array}
\end{array}
herbie shell --seed 2024077
(FPCore (x y z t a b)
:name "AI.Clustering.Hierarchical.Internal:ward from clustering-0.2.1"
:precision binary64
:alt
(if (< (/ (- (+ (* (+ x y) z) (* (+ t y) a)) (* y b)) (+ (+ x t) y)) -3.5813117084150564e+153) (- (+ z a) b) (if (< (/ (- (+ (* (+ x y) z) (* (+ t y) a)) (* y b)) (+ (+ x t) y)) 1.2285964308315609e+82) (/ 1.0 (/ (+ (+ x t) y) (- (+ (* (+ x y) z) (* (+ t y) a)) (* y b)))) (- (+ z a) b)))
(/ (- (+ (* (+ x y) z) (* (+ t y) a)) (* y b)) (+ (+ x t) y)))