
(FPCore (x y z t a b c i j k) :precision binary64 (- (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) (* (* j 27.0) k)))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = (((((((x * 18.0d0) * y) * z) * t) - ((a * 4.0d0) * t)) + (b * c)) - ((x * 4.0d0) * i)) - ((j * 27.0d0) * k)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
def code(x, y, z, t, a, b, c, i, j, k): return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k)
function code(x, y, z, t, a, b, c, i, j, k) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 18.0) * y) * z) * t) - Float64(Float64(a * 4.0) * t)) + Float64(b * c)) - Float64(Float64(x * 4.0) * i)) - Float64(Float64(j * 27.0) * k)) end
function tmp = code(x, y, z, t, a, b, c, i, j, k) tmp = (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k); end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(N[(N[(N[(N[(N[(N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] - N[(N[(a * 4.0), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(x * 4.0), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 25 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c i j k) :precision binary64 (- (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) (* (* j 27.0) k)))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = (((((((x * 18.0d0) * y) * z) * t) - ((a * 4.0d0) * t)) + (b * c)) - ((x * 4.0d0) * i)) - ((j * 27.0d0) * k)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
def code(x, y, z, t, a, b, c, i, j, k): return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k)
function code(x, y, z, t, a, b, c, i, j, k) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 18.0) * y) * z) * t) - Float64(Float64(a * 4.0) * t)) + Float64(b * c)) - Float64(Float64(x * 4.0) * i)) - Float64(Float64(j * 27.0) * k)) end
function tmp = code(x, y, z, t, a, b, c, i, j, k) tmp = (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k); end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(N[(N[(N[(N[(N[(N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] - N[(N[(a * 4.0), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(x * 4.0), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k
\end{array}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= y -1.2e+104)
(*
y
(-
(fma 18.0 (* (* t x) z) (* b (/ c y)))
(fma 4.0 (* a (/ t y)) (fma 4.0 (* i (/ x y)) (* (* j (/ k y)) 27.0)))))
(-
(-
(+ (- (* t (* z (* y (* 18.0 x)))) (* t (* 4.0 a))) (* b c))
(* i (* x 4.0)))
(* k (* j 27.0)))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (y <= -1.2e+104) {
tmp = y * (fma(18.0, ((t * x) * z), (b * (c / y))) - fma(4.0, (a * (t / y)), fma(4.0, (i * (x / y)), ((j * (k / y)) * 27.0))));
} else {
tmp = ((((t * (z * (y * (18.0 * x)))) - (t * (4.0 * a))) + (b * c)) - (i * (x * 4.0))) - (k * (j * 27.0));
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (y <= -1.2e+104) tmp = Float64(y * Float64(fma(18.0, Float64(Float64(t * x) * z), Float64(b * Float64(c / y))) - fma(4.0, Float64(a * Float64(t / y)), fma(4.0, Float64(i * Float64(x / y)), Float64(Float64(j * Float64(k / y)) * 27.0))))); else tmp = Float64(Float64(Float64(Float64(Float64(t * Float64(z * Float64(y * Float64(18.0 * x)))) - Float64(t * Float64(4.0 * a))) + Float64(b * c)) - Float64(i * Float64(x * 4.0))) - Float64(k * Float64(j * 27.0))); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[y, -1.2e+104], N[(y * N[(N[(18.0 * N[(N[(t * x), $MachinePrecision] * z), $MachinePrecision] + N[(b * N[(c / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(4.0 * N[(a * N[(t / y), $MachinePrecision]), $MachinePrecision] + N[(4.0 * N[(i * N[(x / y), $MachinePrecision]), $MachinePrecision] + N[(N[(j * N[(k / y), $MachinePrecision]), $MachinePrecision] * 27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(t * N[(z * N[(y * N[(18.0 * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t * N[(4.0 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(i * N[(x * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(k * N[(j * 27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.2 \cdot 10^{+104}:\\
\;\;\;\;y \cdot \left(\mathsf{fma}\left(18, \left(t \cdot x\right) \cdot z, b \cdot \frac{c}{y}\right) - \mathsf{fma}\left(4, a \cdot \frac{t}{y}, \mathsf{fma}\left(4, i \cdot \frac{x}{y}, \left(j \cdot \frac{k}{y}\right) \cdot 27\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\left(t \cdot \left(z \cdot \left(y \cdot \left(18 \cdot x\right)\right)\right) - t \cdot \left(4 \cdot a\right)\right) + b \cdot c\right) - i \cdot \left(x \cdot 4\right)\right) - k \cdot \left(j \cdot 27\right)\\
\end{array}
\end{array}
if y < -1.2e104Initial program 68.5%
Simplified74.9%
Taylor expanded in y around inf 89.3%
sub-neg89.3%
+-commutative89.3%
metadata-eval89.3%
cancel-sign-sub-inv89.3%
sub-neg89.3%
associate--r+89.3%
Simplified95.6%
if -1.2e104 < y Initial program 87.6%
Final simplification89.1%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1
(-
(-
(+ (- (* t (* z (* y (* 18.0 x)))) (* t (* 4.0 a))) (* b c))
(* i (* x 4.0)))
(* k (* j 27.0)))))
(if (<= t_1 INFINITY) t_1 (* x (+ (* (* 18.0 t) (* y z)) (* i -4.0))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = ((((t * (z * (y * (18.0 * x)))) - (t * (4.0 * a))) + (b * c)) - (i * (x * 4.0))) - (k * (j * 27.0));
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = x * (((18.0 * t) * (y * z)) + (i * -4.0));
}
return tmp;
}
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = ((((t * (z * (y * (18.0 * x)))) - (t * (4.0 * a))) + (b * c)) - (i * (x * 4.0))) - (k * (j * 27.0));
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = x * (((18.0 * t) * (y * z)) + (i * -4.0));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = ((((t * (z * (y * (18.0 * x)))) - (t * (4.0 * a))) + (b * c)) - (i * (x * 4.0))) - (k * (j * 27.0)) tmp = 0 if t_1 <= math.inf: tmp = t_1 else: tmp = x * (((18.0 * t) * (y * z)) + (i * -4.0)) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(Float64(Float64(Float64(t * Float64(z * Float64(y * Float64(18.0 * x)))) - Float64(t * Float64(4.0 * a))) + Float64(b * c)) - Float64(i * Float64(x * 4.0))) - Float64(k * Float64(j * 27.0))) tmp = 0.0 if (t_1 <= Inf) tmp = t_1; else tmp = Float64(x * Float64(Float64(Float64(18.0 * t) * Float64(y * z)) + Float64(i * -4.0))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = ((((t * (z * (y * (18.0 * x)))) - (t * (4.0 * a))) + (b * c)) - (i * (x * 4.0))) - (k * (j * 27.0));
tmp = 0.0;
if (t_1 <= Inf)
tmp = t_1;
else
tmp = x * (((18.0 * t) * (y * z)) + (i * -4.0));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(N[(N[(N[(t * N[(z * N[(y * N[(18.0 * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t * N[(4.0 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(i * N[(x * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(k * N[(j * 27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, Infinity], t$95$1, N[(x * N[(N[(N[(18.0 * t), $MachinePrecision] * N[(y * z), $MachinePrecision]), $MachinePrecision] + N[(i * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := \left(\left(\left(t \cdot \left(z \cdot \left(y \cdot \left(18 \cdot x\right)\right)\right) - t \cdot \left(4 \cdot a\right)\right) + b \cdot c\right) - i \cdot \left(x \cdot 4\right)\right) - k \cdot \left(j \cdot 27\right)\\
\mathbf{if}\;t\_1 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\left(18 \cdot t\right) \cdot \left(y \cdot z\right) + i \cdot -4\right)\\
\end{array}
\end{array}
if (-.f64 (-.f64 (+.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 x #s(literal 18 binary64)) y) z) t) (*.f64 (*.f64 a #s(literal 4 binary64)) t)) (*.f64 b c)) (*.f64 (*.f64 x #s(literal 4 binary64)) i)) (*.f64 (*.f64 j #s(literal 27 binary64)) k)) < +inf.0Initial program 97.4%
if +inf.0 < (-.f64 (-.f64 (+.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 x #s(literal 18 binary64)) y) z) t) (*.f64 (*.f64 a #s(literal 4 binary64)) t)) (*.f64 b c)) (*.f64 (*.f64 x #s(literal 4 binary64)) i)) (*.f64 (*.f64 j #s(literal 27 binary64)) k)) Initial program 0.0%
Simplified28.6%
Taylor expanded in a around 0 25.9%
Taylor expanded in b around 0 28.6%
Taylor expanded in x around inf 65.9%
cancel-sign-sub-inv65.9%
associate-*r*65.9%
*-commutative65.9%
metadata-eval65.9%
*-commutative65.9%
Simplified65.9%
Final simplification93.1%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (+ (* b c) (* -4.0 (* t a))))
(t_2 (* 4.0 (* x i)))
(t_3 (* k (* j 27.0))))
(if (<= t_3 -1e-15)
(- (* b c) (* k (+ (* j 27.0) (* 4.0 (/ (* x i) k)))))
(if (<= t_3 200000000000.0)
(- t_1 t_2)
(if (<= t_3 5e+48)
(- (* 18.0 (* t (* x (* y z)))) (+ t_2 (* 27.0 (* j k))))
(- t_1 t_3))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (b * c) + (-4.0 * (t * a));
double t_2 = 4.0 * (x * i);
double t_3 = k * (j * 27.0);
double tmp;
if (t_3 <= -1e-15) {
tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k))));
} else if (t_3 <= 200000000000.0) {
tmp = t_1 - t_2;
} else if (t_3 <= 5e+48) {
tmp = (18.0 * (t * (x * (y * z)))) - (t_2 + (27.0 * (j * k)));
} else {
tmp = t_1 - t_3;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = (b * c) + ((-4.0d0) * (t * a))
t_2 = 4.0d0 * (x * i)
t_3 = k * (j * 27.0d0)
if (t_3 <= (-1d-15)) then
tmp = (b * c) - (k * ((j * 27.0d0) + (4.0d0 * ((x * i) / k))))
else if (t_3 <= 200000000000.0d0) then
tmp = t_1 - t_2
else if (t_3 <= 5d+48) then
tmp = (18.0d0 * (t * (x * (y * z)))) - (t_2 + (27.0d0 * (j * k)))
else
tmp = t_1 - t_3
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (b * c) + (-4.0 * (t * a));
double t_2 = 4.0 * (x * i);
double t_3 = k * (j * 27.0);
double tmp;
if (t_3 <= -1e-15) {
tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k))));
} else if (t_3 <= 200000000000.0) {
tmp = t_1 - t_2;
} else if (t_3 <= 5e+48) {
tmp = (18.0 * (t * (x * (y * z)))) - (t_2 + (27.0 * (j * k)));
} else {
tmp = t_1 - t_3;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = (b * c) + (-4.0 * (t * a)) t_2 = 4.0 * (x * i) t_3 = k * (j * 27.0) tmp = 0 if t_3 <= -1e-15: tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k)))) elif t_3 <= 200000000000.0: tmp = t_1 - t_2 elif t_3 <= 5e+48: tmp = (18.0 * (t * (x * (y * z)))) - (t_2 + (27.0 * (j * k))) else: tmp = t_1 - t_3 return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))) t_2 = Float64(4.0 * Float64(x * i)) t_3 = Float64(k * Float64(j * 27.0)) tmp = 0.0 if (t_3 <= -1e-15) tmp = Float64(Float64(b * c) - Float64(k * Float64(Float64(j * 27.0) + Float64(4.0 * Float64(Float64(x * i) / k))))); elseif (t_3 <= 200000000000.0) tmp = Float64(t_1 - t_2); elseif (t_3 <= 5e+48) tmp = Float64(Float64(18.0 * Float64(t * Float64(x * Float64(y * z)))) - Float64(t_2 + Float64(27.0 * Float64(j * k)))); else tmp = Float64(t_1 - t_3); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = (b * c) + (-4.0 * (t * a));
t_2 = 4.0 * (x * i);
t_3 = k * (j * 27.0);
tmp = 0.0;
if (t_3 <= -1e-15)
tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k))));
elseif (t_3 <= 200000000000.0)
tmp = t_1 - t_2;
elseif (t_3 <= 5e+48)
tmp = (18.0 * (t * (x * (y * z)))) - (t_2 + (27.0 * (j * k)));
else
tmp = t_1 - t_3;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(k * N[(j * 27.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, -1e-15], N[(N[(b * c), $MachinePrecision] - N[(k * N[(N[(j * 27.0), $MachinePrecision] + N[(4.0 * N[(N[(x * i), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 200000000000.0], N[(t$95$1 - t$95$2), $MachinePrecision], If[LessEqual[t$95$3, 5e+48], N[(N[(18.0 * N[(t * N[(x * N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t$95$2 + N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 - t$95$3), $MachinePrecision]]]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := b \cdot c + -4 \cdot \left(t \cdot a\right)\\
t_2 := 4 \cdot \left(x \cdot i\right)\\
t_3 := k \cdot \left(j \cdot 27\right)\\
\mathbf{if}\;t\_3 \leq -1 \cdot 10^{-15}:\\
\;\;\;\;b \cdot c - k \cdot \left(j \cdot 27 + 4 \cdot \frac{x \cdot i}{k}\right)\\
\mathbf{elif}\;t\_3 \leq 200000000000:\\
\;\;\;\;t\_1 - t\_2\\
\mathbf{elif}\;t\_3 \leq 5 \cdot 10^{+48}:\\
\;\;\;\;18 \cdot \left(t \cdot \left(x \cdot \left(y \cdot z\right)\right)\right) - \left(t\_2 + 27 \cdot \left(j \cdot k\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1 - t\_3\\
\end{array}
\end{array}
if (*.f64 (*.f64 j #s(literal 27 binary64)) k) < -1.0000000000000001e-15Initial program 74.2%
Simplified76.0%
Taylor expanded in t around 0 70.5%
Taylor expanded in k around inf 72.2%
if -1.0000000000000001e-15 < (*.f64 (*.f64 j #s(literal 27 binary64)) k) < 2e11Initial program 87.1%
distribute-rgt-out--90.0%
associate-*r*88.0%
*-commutative88.0%
associate-*l*88.0%
associate-*r*88.0%
Applied egg-rr88.0%
Taylor expanded in y around 0 80.3%
Taylor expanded in j around 0 78.9%
if 2e11 < (*.f64 (*.f64 j #s(literal 27 binary64)) k) < 4.99999999999999973e48Initial program 87.7%
Simplified99.4%
Taylor expanded in a around 0 99.6%
Taylor expanded in b around 0 99.6%
if 4.99999999999999973e48 < (*.f64 (*.f64 j #s(literal 27 binary64)) k) Initial program 87.6%
distribute-rgt-out--89.7%
associate-*r*89.7%
*-commutative89.7%
associate-*l*89.7%
associate-*r*89.7%
Applied egg-rr89.7%
Taylor expanded in x around 0 75.0%
Final simplification77.2%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* j (* k -27.0)))
(t_2 (* t (- (* 4.0 (- a)) (* (* y x) (* z -18.0))))))
(if (<= t -4.3e-5)
t_2
(if (<= t -5.5e-86)
(+ (* 18.0 (* (* t x) (* y z))) t_1)
(if (<= t 2.8e-286)
(- (* b c) (* x (* 4.0 i)))
(if (<= t 1.6e-214)
(- (* b c) (* 27.0 (* j k)))
(if (<= t 1e+41) (+ t_1 (* -4.0 (* x i))) t_2)))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = j * (k * -27.0);
double t_2 = t * ((4.0 * -a) - ((y * x) * (z * -18.0)));
double tmp;
if (t <= -4.3e-5) {
tmp = t_2;
} else if (t <= -5.5e-86) {
tmp = (18.0 * ((t * x) * (y * z))) + t_1;
} else if (t <= 2.8e-286) {
tmp = (b * c) - (x * (4.0 * i));
} else if (t <= 1.6e-214) {
tmp = (b * c) - (27.0 * (j * k));
} else if (t <= 1e+41) {
tmp = t_1 + (-4.0 * (x * i));
} else {
tmp = t_2;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = j * (k * (-27.0d0))
t_2 = t * ((4.0d0 * -a) - ((y * x) * (z * (-18.0d0))))
if (t <= (-4.3d-5)) then
tmp = t_2
else if (t <= (-5.5d-86)) then
tmp = (18.0d0 * ((t * x) * (y * z))) + t_1
else if (t <= 2.8d-286) then
tmp = (b * c) - (x * (4.0d0 * i))
else if (t <= 1.6d-214) then
tmp = (b * c) - (27.0d0 * (j * k))
else if (t <= 1d+41) then
tmp = t_1 + ((-4.0d0) * (x * i))
else
tmp = t_2
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = j * (k * -27.0);
double t_2 = t * ((4.0 * -a) - ((y * x) * (z * -18.0)));
double tmp;
if (t <= -4.3e-5) {
tmp = t_2;
} else if (t <= -5.5e-86) {
tmp = (18.0 * ((t * x) * (y * z))) + t_1;
} else if (t <= 2.8e-286) {
tmp = (b * c) - (x * (4.0 * i));
} else if (t <= 1.6e-214) {
tmp = (b * c) - (27.0 * (j * k));
} else if (t <= 1e+41) {
tmp = t_1 + (-4.0 * (x * i));
} else {
tmp = t_2;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = j * (k * -27.0) t_2 = t * ((4.0 * -a) - ((y * x) * (z * -18.0))) tmp = 0 if t <= -4.3e-5: tmp = t_2 elif t <= -5.5e-86: tmp = (18.0 * ((t * x) * (y * z))) + t_1 elif t <= 2.8e-286: tmp = (b * c) - (x * (4.0 * i)) elif t <= 1.6e-214: tmp = (b * c) - (27.0 * (j * k)) elif t <= 1e+41: tmp = t_1 + (-4.0 * (x * i)) else: tmp = t_2 return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(j * Float64(k * -27.0)) t_2 = Float64(t * Float64(Float64(4.0 * Float64(-a)) - Float64(Float64(y * x) * Float64(z * -18.0)))) tmp = 0.0 if (t <= -4.3e-5) tmp = t_2; elseif (t <= -5.5e-86) tmp = Float64(Float64(18.0 * Float64(Float64(t * x) * Float64(y * z))) + t_1); elseif (t <= 2.8e-286) tmp = Float64(Float64(b * c) - Float64(x * Float64(4.0 * i))); elseif (t <= 1.6e-214) tmp = Float64(Float64(b * c) - Float64(27.0 * Float64(j * k))); elseif (t <= 1e+41) tmp = Float64(t_1 + Float64(-4.0 * Float64(x * i))); else tmp = t_2; end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = j * (k * -27.0);
t_2 = t * ((4.0 * -a) - ((y * x) * (z * -18.0)));
tmp = 0.0;
if (t <= -4.3e-5)
tmp = t_2;
elseif (t <= -5.5e-86)
tmp = (18.0 * ((t * x) * (y * z))) + t_1;
elseif (t <= 2.8e-286)
tmp = (b * c) - (x * (4.0 * i));
elseif (t <= 1.6e-214)
tmp = (b * c) - (27.0 * (j * k));
elseif (t <= 1e+41)
tmp = t_1 + (-4.0 * (x * i));
else
tmp = t_2;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(j * N[(k * -27.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t * N[(N[(4.0 * (-a)), $MachinePrecision] - N[(N[(y * x), $MachinePrecision] * N[(z * -18.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -4.3e-5], t$95$2, If[LessEqual[t, -5.5e-86], N[(N[(18.0 * N[(N[(t * x), $MachinePrecision] * N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$1), $MachinePrecision], If[LessEqual[t, 2.8e-286], N[(N[(b * c), $MachinePrecision] - N[(x * N[(4.0 * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.6e-214], N[(N[(b * c), $MachinePrecision] - N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1e+41], N[(t$95$1 + N[(-4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := j \cdot \left(k \cdot -27\right)\\
t_2 := t \cdot \left(4 \cdot \left(-a\right) - \left(y \cdot x\right) \cdot \left(z \cdot -18\right)\right)\\
\mathbf{if}\;t \leq -4.3 \cdot 10^{-5}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t \leq -5.5 \cdot 10^{-86}:\\
\;\;\;\;18 \cdot \left(\left(t \cdot x\right) \cdot \left(y \cdot z\right)\right) + t\_1\\
\mathbf{elif}\;t \leq 2.8 \cdot 10^{-286}:\\
\;\;\;\;b \cdot c - x \cdot \left(4 \cdot i\right)\\
\mathbf{elif}\;t \leq 1.6 \cdot 10^{-214}:\\
\;\;\;\;b \cdot c - 27 \cdot \left(j \cdot k\right)\\
\mathbf{elif}\;t \leq 10^{+41}:\\
\;\;\;\;t\_1 + -4 \cdot \left(x \cdot i\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if t < -4.3000000000000002e-5 or 1.00000000000000001e41 < t Initial program 84.5%
Taylor expanded in t around -inf 65.5%
associate-*r*65.5%
neg-mul-165.5%
cancel-sign-sub-inv65.5%
metadata-eval65.5%
*-commutative65.5%
*-commutative65.5%
Simplified65.5%
Taylor expanded in x around 0 65.5%
*-commutative65.5%
associate-*r*65.5%
associate-*l*65.5%
associate-*r*66.3%
Simplified66.3%
if -4.3000000000000002e-5 < t < -5.5e-86Initial program 93.4%
Simplified86.3%
Taylor expanded in y around inf 57.8%
associate-*r*64.2%
Simplified64.2%
if -5.5e-86 < t < 2.8e-286Initial program 75.5%
Simplified71.9%
Taylor expanded in t around 0 81.5%
Taylor expanded in i around inf 79.8%
associate-*r*79.8%
*-commutative79.8%
Simplified79.8%
if 2.8e-286 < t < 1.60000000000000007e-214Initial program 87.2%
Simplified81.2%
Taylor expanded in t around 0 87.8%
Taylor expanded in i around 0 94.0%
if 1.60000000000000007e-214 < t < 1.00000000000000001e41Initial program 88.1%
Simplified89.7%
Taylor expanded in i around inf 57.7%
Final simplification68.8%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= t -4e-88)
(-
(+ (* b c) (* t (- (* (* 18.0 x) (* y z)) (* 4.0 a))))
(+ (* x (* 4.0 i)) (* j (* k 27.0))))
(if (<= t 1.95e-72)
(- (+ (* b c) (* -4.0 (* t a))) (+ (* 4.0 (* x i)) (* 27.0 (* j k))))
(-
(-
(+ (* b c) (* t (- (* x (* z (* y 18.0))) (* 4.0 a))))
(* i (* x 4.0)))
(* k (* j 27.0))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (t <= -4e-88) {
tmp = ((b * c) + (t * (((18.0 * x) * (y * z)) - (4.0 * a)))) - ((x * (4.0 * i)) + (j * (k * 27.0)));
} else if (t <= 1.95e-72) {
tmp = ((b * c) + (-4.0 * (t * a))) - ((4.0 * (x * i)) + (27.0 * (j * k)));
} else {
tmp = (((b * c) + (t * ((x * (z * (y * 18.0))) - (4.0 * a)))) - (i * (x * 4.0))) - (k * (j * 27.0));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if (t <= (-4d-88)) then
tmp = ((b * c) + (t * (((18.0d0 * x) * (y * z)) - (4.0d0 * a)))) - ((x * (4.0d0 * i)) + (j * (k * 27.0d0)))
else if (t <= 1.95d-72) then
tmp = ((b * c) + ((-4.0d0) * (t * a))) - ((4.0d0 * (x * i)) + (27.0d0 * (j * k)))
else
tmp = (((b * c) + (t * ((x * (z * (y * 18.0d0))) - (4.0d0 * a)))) - (i * (x * 4.0d0))) - (k * (j * 27.0d0))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (t <= -4e-88) {
tmp = ((b * c) + (t * (((18.0 * x) * (y * z)) - (4.0 * a)))) - ((x * (4.0 * i)) + (j * (k * 27.0)));
} else if (t <= 1.95e-72) {
tmp = ((b * c) + (-4.0 * (t * a))) - ((4.0 * (x * i)) + (27.0 * (j * k)));
} else {
tmp = (((b * c) + (t * ((x * (z * (y * 18.0))) - (4.0 * a)))) - (i * (x * 4.0))) - (k * (j * 27.0));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if t <= -4e-88: tmp = ((b * c) + (t * (((18.0 * x) * (y * z)) - (4.0 * a)))) - ((x * (4.0 * i)) + (j * (k * 27.0))) elif t <= 1.95e-72: tmp = ((b * c) + (-4.0 * (t * a))) - ((4.0 * (x * i)) + (27.0 * (j * k))) else: tmp = (((b * c) + (t * ((x * (z * (y * 18.0))) - (4.0 * a)))) - (i * (x * 4.0))) - (k * (j * 27.0)) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (t <= -4e-88) tmp = Float64(Float64(Float64(b * c) + Float64(t * Float64(Float64(Float64(18.0 * x) * Float64(y * z)) - Float64(4.0 * a)))) - Float64(Float64(x * Float64(4.0 * i)) + Float64(j * Float64(k * 27.0)))); elseif (t <= 1.95e-72) tmp = Float64(Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))) - Float64(Float64(4.0 * Float64(x * i)) + Float64(27.0 * Float64(j * k)))); else tmp = Float64(Float64(Float64(Float64(b * c) + Float64(t * Float64(Float64(x * Float64(z * Float64(y * 18.0))) - Float64(4.0 * a)))) - Float64(i * Float64(x * 4.0))) - Float64(k * Float64(j * 27.0))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if (t <= -4e-88)
tmp = ((b * c) + (t * (((18.0 * x) * (y * z)) - (4.0 * a)))) - ((x * (4.0 * i)) + (j * (k * 27.0)));
elseif (t <= 1.95e-72)
tmp = ((b * c) + (-4.0 * (t * a))) - ((4.0 * (x * i)) + (27.0 * (j * k)));
else
tmp = (((b * c) + (t * ((x * (z * (y * 18.0))) - (4.0 * a)))) - (i * (x * 4.0))) - (k * (j * 27.0));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[t, -4e-88], N[(N[(N[(b * c), $MachinePrecision] + N[(t * N[(N[(N[(18.0 * x), $MachinePrecision] * N[(y * z), $MachinePrecision]), $MachinePrecision] - N[(4.0 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(x * N[(4.0 * i), $MachinePrecision]), $MachinePrecision] + N[(j * N[(k * 27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.95e-72], N[(N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision] + N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(b * c), $MachinePrecision] + N[(t * N[(N[(x * N[(z * N[(y * 18.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(4.0 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(i * N[(x * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(k * N[(j * 27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;t \leq -4 \cdot 10^{-88}:\\
\;\;\;\;\left(b \cdot c + t \cdot \left(\left(18 \cdot x\right) \cdot \left(y \cdot z\right) - 4 \cdot a\right)\right) - \left(x \cdot \left(4 \cdot i\right) + j \cdot \left(k \cdot 27\right)\right)\\
\mathbf{elif}\;t \leq 1.95 \cdot 10^{-72}:\\
\;\;\;\;\left(b \cdot c + -4 \cdot \left(t \cdot a\right)\right) - \left(4 \cdot \left(x \cdot i\right) + 27 \cdot \left(j \cdot k\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(b \cdot c + t \cdot \left(x \cdot \left(z \cdot \left(y \cdot 18\right)\right) - 4 \cdot a\right)\right) - i \cdot \left(x \cdot 4\right)\right) - k \cdot \left(j \cdot 27\right)\\
\end{array}
\end{array}
if t < -3.99999999999999974e-88Initial program 88.5%
Simplified92.2%
if -3.99999999999999974e-88 < t < 1.95e-72Initial program 80.9%
Simplified78.0%
Taylor expanded in y around 0 86.2%
if 1.95e-72 < t Initial program 83.2%
distribute-rgt-out--87.4%
associate-*r*88.8%
*-commutative88.8%
associate-*l*88.8%
associate-*r*88.8%
Applied egg-rr88.8%
Final simplification89.0%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= (* k (* j 27.0)) -1e+292)
(- (* b c) (* k (+ (* j 27.0) (* 4.0 (/ (* x i) k)))))
(-
(+ (* b c) (* t (- (* (* 18.0 x) (* y z)) (* 4.0 a))))
(+ (* x (* 4.0 i)) (* j (* k 27.0))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((k * (j * 27.0)) <= -1e+292) {
tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k))));
} else {
tmp = ((b * c) + (t * (((18.0 * x) * (y * z)) - (4.0 * a)))) - ((x * (4.0 * i)) + (j * (k * 27.0)));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if ((k * (j * 27.0d0)) <= (-1d+292)) then
tmp = (b * c) - (k * ((j * 27.0d0) + (4.0d0 * ((x * i) / k))))
else
tmp = ((b * c) + (t * (((18.0d0 * x) * (y * z)) - (4.0d0 * a)))) - ((x * (4.0d0 * i)) + (j * (k * 27.0d0)))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((k * (j * 27.0)) <= -1e+292) {
tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k))));
} else {
tmp = ((b * c) + (t * (((18.0 * x) * (y * z)) - (4.0 * a)))) - ((x * (4.0 * i)) + (j * (k * 27.0)));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if (k * (j * 27.0)) <= -1e+292: tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k)))) else: tmp = ((b * c) + (t * (((18.0 * x) * (y * z)) - (4.0 * a)))) - ((x * (4.0 * i)) + (j * (k * 27.0))) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (Float64(k * Float64(j * 27.0)) <= -1e+292) tmp = Float64(Float64(b * c) - Float64(k * Float64(Float64(j * 27.0) + Float64(4.0 * Float64(Float64(x * i) / k))))); else tmp = Float64(Float64(Float64(b * c) + Float64(t * Float64(Float64(Float64(18.0 * x) * Float64(y * z)) - Float64(4.0 * a)))) - Float64(Float64(x * Float64(4.0 * i)) + Float64(j * Float64(k * 27.0)))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if ((k * (j * 27.0)) <= -1e+292)
tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k))));
else
tmp = ((b * c) + (t * (((18.0 * x) * (y * z)) - (4.0 * a)))) - ((x * (4.0 * i)) + (j * (k * 27.0)));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[N[(k * N[(j * 27.0), $MachinePrecision]), $MachinePrecision], -1e+292], N[(N[(b * c), $MachinePrecision] - N[(k * N[(N[(j * 27.0), $MachinePrecision] + N[(4.0 * N[(N[(x * i), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b * c), $MachinePrecision] + N[(t * N[(N[(N[(18.0 * x), $MachinePrecision] * N[(y * z), $MachinePrecision]), $MachinePrecision] - N[(4.0 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(x * N[(4.0 * i), $MachinePrecision]), $MachinePrecision] + N[(j * N[(k * 27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;k \cdot \left(j \cdot 27\right) \leq -1 \cdot 10^{+292}:\\
\;\;\;\;b \cdot c - k \cdot \left(j \cdot 27 + 4 \cdot \frac{x \cdot i}{k}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(b \cdot c + t \cdot \left(\left(18 \cdot x\right) \cdot \left(y \cdot z\right) - 4 \cdot a\right)\right) - \left(x \cdot \left(4 \cdot i\right) + j \cdot \left(k \cdot 27\right)\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 j #s(literal 27 binary64)) k) < -1e292Initial program 47.1%
Simplified41.8%
Taylor expanded in t around 0 59.4%
Taylor expanded in k around inf 71.2%
if -1e292 < (*.f64 (*.f64 j #s(literal 27 binary64)) k) Initial program 86.7%
Simplified89.3%
Final simplification88.1%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* t (- (* 4.0 (- a)) (* (* y x) (* z -18.0))))))
(if (<= t -2.9e-87)
t_1
(if (<= t 1.7e-283)
(- (* b c) (* x (* 4.0 i)))
(if (<= t 5.4e-222)
(- (* b c) (* 27.0 (* j k)))
(if (<= t 1e+41) (+ (* j (* k -27.0)) (* -4.0 (* x i))) t_1))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = t * ((4.0 * -a) - ((y * x) * (z * -18.0)));
double tmp;
if (t <= -2.9e-87) {
tmp = t_1;
} else if (t <= 1.7e-283) {
tmp = (b * c) - (x * (4.0 * i));
} else if (t <= 5.4e-222) {
tmp = (b * c) - (27.0 * (j * k));
} else if (t <= 1e+41) {
tmp = (j * (k * -27.0)) + (-4.0 * (x * i));
} else {
tmp = t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = t * ((4.0d0 * -a) - ((y * x) * (z * (-18.0d0))))
if (t <= (-2.9d-87)) then
tmp = t_1
else if (t <= 1.7d-283) then
tmp = (b * c) - (x * (4.0d0 * i))
else if (t <= 5.4d-222) then
tmp = (b * c) - (27.0d0 * (j * k))
else if (t <= 1d+41) then
tmp = (j * (k * (-27.0d0))) + ((-4.0d0) * (x * i))
else
tmp = t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = t * ((4.0 * -a) - ((y * x) * (z * -18.0)));
double tmp;
if (t <= -2.9e-87) {
tmp = t_1;
} else if (t <= 1.7e-283) {
tmp = (b * c) - (x * (4.0 * i));
} else if (t <= 5.4e-222) {
tmp = (b * c) - (27.0 * (j * k));
} else if (t <= 1e+41) {
tmp = (j * (k * -27.0)) + (-4.0 * (x * i));
} else {
tmp = t_1;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = t * ((4.0 * -a) - ((y * x) * (z * -18.0))) tmp = 0 if t <= -2.9e-87: tmp = t_1 elif t <= 1.7e-283: tmp = (b * c) - (x * (4.0 * i)) elif t <= 5.4e-222: tmp = (b * c) - (27.0 * (j * k)) elif t <= 1e+41: tmp = (j * (k * -27.0)) + (-4.0 * (x * i)) else: tmp = t_1 return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(t * Float64(Float64(4.0 * Float64(-a)) - Float64(Float64(y * x) * Float64(z * -18.0)))) tmp = 0.0 if (t <= -2.9e-87) tmp = t_1; elseif (t <= 1.7e-283) tmp = Float64(Float64(b * c) - Float64(x * Float64(4.0 * i))); elseif (t <= 5.4e-222) tmp = Float64(Float64(b * c) - Float64(27.0 * Float64(j * k))); elseif (t <= 1e+41) tmp = Float64(Float64(j * Float64(k * -27.0)) + Float64(-4.0 * Float64(x * i))); else tmp = t_1; end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = t * ((4.0 * -a) - ((y * x) * (z * -18.0)));
tmp = 0.0;
if (t <= -2.9e-87)
tmp = t_1;
elseif (t <= 1.7e-283)
tmp = (b * c) - (x * (4.0 * i));
elseif (t <= 5.4e-222)
tmp = (b * c) - (27.0 * (j * k));
elseif (t <= 1e+41)
tmp = (j * (k * -27.0)) + (-4.0 * (x * i));
else
tmp = t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(t * N[(N[(4.0 * (-a)), $MachinePrecision] - N[(N[(y * x), $MachinePrecision] * N[(z * -18.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2.9e-87], t$95$1, If[LessEqual[t, 1.7e-283], N[(N[(b * c), $MachinePrecision] - N[(x * N[(4.0 * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 5.4e-222], N[(N[(b * c), $MachinePrecision] - N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1e+41], N[(N[(j * N[(k * -27.0), $MachinePrecision]), $MachinePrecision] + N[(-4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := t \cdot \left(4 \cdot \left(-a\right) - \left(y \cdot x\right) \cdot \left(z \cdot -18\right)\right)\\
\mathbf{if}\;t \leq -2.9 \cdot 10^{-87}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 1.7 \cdot 10^{-283}:\\
\;\;\;\;b \cdot c - x \cdot \left(4 \cdot i\right)\\
\mathbf{elif}\;t \leq 5.4 \cdot 10^{-222}:\\
\;\;\;\;b \cdot c - 27 \cdot \left(j \cdot k\right)\\
\mathbf{elif}\;t \leq 10^{+41}:\\
\;\;\;\;j \cdot \left(k \cdot -27\right) + -4 \cdot \left(x \cdot i\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -2.8999999999999999e-87 or 1.00000000000000001e41 < t Initial program 85.4%
Taylor expanded in t around -inf 61.6%
associate-*r*61.6%
neg-mul-161.6%
cancel-sign-sub-inv61.6%
metadata-eval61.6%
*-commutative61.6%
*-commutative61.6%
Simplified61.6%
Taylor expanded in x around 0 61.6%
*-commutative61.6%
associate-*r*61.6%
associate-*l*61.6%
associate-*r*62.3%
Simplified62.3%
if -2.8999999999999999e-87 < t < 1.6999999999999999e-283Initial program 75.5%
Simplified71.9%
Taylor expanded in t around 0 81.5%
Taylor expanded in i around inf 79.8%
associate-*r*79.8%
*-commutative79.8%
Simplified79.8%
if 1.6999999999999999e-283 < t < 5.4e-222Initial program 86.5%
Simplified86.5%
Taylor expanded in t around 0 93.2%
Taylor expanded in i around 0 99.9%
if 5.4e-222 < t < 1.00000000000000001e41Initial program 88.3%
Simplified88.2%
Taylor expanded in i around inf 56.9%
Final simplification66.9%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* x (+ (* (* 18.0 t) (* y z)) (* i -4.0)))))
(if (<= x -9.2e+82)
t_1
(if (<= x -1.85e+34)
(+ (* b c) (* j (* k -27.0)))
(if (<= x -6e-225)
(+ (* b c) (* -4.0 (* t a)))
(if (<= x 2.8e+48) (- (* b c) (* 27.0 (* j k))) t_1))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = x * (((18.0 * t) * (y * z)) + (i * -4.0));
double tmp;
if (x <= -9.2e+82) {
tmp = t_1;
} else if (x <= -1.85e+34) {
tmp = (b * c) + (j * (k * -27.0));
} else if (x <= -6e-225) {
tmp = (b * c) + (-4.0 * (t * a));
} else if (x <= 2.8e+48) {
tmp = (b * c) - (27.0 * (j * k));
} else {
tmp = t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = x * (((18.0d0 * t) * (y * z)) + (i * (-4.0d0)))
if (x <= (-9.2d+82)) then
tmp = t_1
else if (x <= (-1.85d+34)) then
tmp = (b * c) + (j * (k * (-27.0d0)))
else if (x <= (-6d-225)) then
tmp = (b * c) + ((-4.0d0) * (t * a))
else if (x <= 2.8d+48) then
tmp = (b * c) - (27.0d0 * (j * k))
else
tmp = t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = x * (((18.0 * t) * (y * z)) + (i * -4.0));
double tmp;
if (x <= -9.2e+82) {
tmp = t_1;
} else if (x <= -1.85e+34) {
tmp = (b * c) + (j * (k * -27.0));
} else if (x <= -6e-225) {
tmp = (b * c) + (-4.0 * (t * a));
} else if (x <= 2.8e+48) {
tmp = (b * c) - (27.0 * (j * k));
} else {
tmp = t_1;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = x * (((18.0 * t) * (y * z)) + (i * -4.0)) tmp = 0 if x <= -9.2e+82: tmp = t_1 elif x <= -1.85e+34: tmp = (b * c) + (j * (k * -27.0)) elif x <= -6e-225: tmp = (b * c) + (-4.0 * (t * a)) elif x <= 2.8e+48: tmp = (b * c) - (27.0 * (j * k)) else: tmp = t_1 return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(x * Float64(Float64(Float64(18.0 * t) * Float64(y * z)) + Float64(i * -4.0))) tmp = 0.0 if (x <= -9.2e+82) tmp = t_1; elseif (x <= -1.85e+34) tmp = Float64(Float64(b * c) + Float64(j * Float64(k * -27.0))); elseif (x <= -6e-225) tmp = Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))); elseif (x <= 2.8e+48) tmp = Float64(Float64(b * c) - Float64(27.0 * Float64(j * k))); else tmp = t_1; end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = x * (((18.0 * t) * (y * z)) + (i * -4.0));
tmp = 0.0;
if (x <= -9.2e+82)
tmp = t_1;
elseif (x <= -1.85e+34)
tmp = (b * c) + (j * (k * -27.0));
elseif (x <= -6e-225)
tmp = (b * c) + (-4.0 * (t * a));
elseif (x <= 2.8e+48)
tmp = (b * c) - (27.0 * (j * k));
else
tmp = t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(x * N[(N[(N[(18.0 * t), $MachinePrecision] * N[(y * z), $MachinePrecision]), $MachinePrecision] + N[(i * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -9.2e+82], t$95$1, If[LessEqual[x, -1.85e+34], N[(N[(b * c), $MachinePrecision] + N[(j * N[(k * -27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -6e-225], N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 2.8e+48], N[(N[(b * c), $MachinePrecision] - N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := x \cdot \left(\left(18 \cdot t\right) \cdot \left(y \cdot z\right) + i \cdot -4\right)\\
\mathbf{if}\;x \leq -9.2 \cdot 10^{+82}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq -1.85 \cdot 10^{+34}:\\
\;\;\;\;b \cdot c + j \cdot \left(k \cdot -27\right)\\
\mathbf{elif}\;x \leq -6 \cdot 10^{-225}:\\
\;\;\;\;b \cdot c + -4 \cdot \left(t \cdot a\right)\\
\mathbf{elif}\;x \leq 2.8 \cdot 10^{+48}:\\
\;\;\;\;b \cdot c - 27 \cdot \left(j \cdot k\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -9.19999999999999953e82 or 2.80000000000000012e48 < x Initial program 73.0%
Simplified81.4%
Taylor expanded in a around 0 72.4%
Taylor expanded in b around 0 63.1%
Taylor expanded in x around inf 73.2%
cancel-sign-sub-inv73.2%
associate-*r*73.2%
*-commutative73.2%
metadata-eval73.2%
*-commutative73.2%
Simplified73.2%
if -9.19999999999999953e82 < x < -1.85000000000000004e34Initial program 72.3%
Simplified81.7%
Taylor expanded in b around inf 68.3%
if -1.85000000000000004e34 < x < -6.00000000000000035e-225Initial program 92.7%
distribute-rgt-out--94.5%
associate-*r*92.8%
*-commutative92.8%
associate-*l*92.8%
associate-*r*92.8%
Applied egg-rr92.8%
Taylor expanded in y around 0 80.6%
Taylor expanded in j around 0 71.9%
Taylor expanded in i around 0 63.0%
if -6.00000000000000035e-225 < x < 2.80000000000000012e48Initial program 94.0%
Simplified88.5%
Taylor expanded in t around 0 67.9%
Taylor expanded in i around 0 61.0%
Final simplification66.8%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* 4.0 (* x i))))
(if (<= z -8.2e+179)
(-
(- (+ (* b c) (* (* 18.0 t) (* z (* y x)))) (* i (* x 4.0)))
(* k (* j 27.0)))
(if (<= z 3.9e+139)
(- (+ (* b c) (* -4.0 (* t a))) (+ t_1 (* 27.0 (* j k))))
(- (+ (* b c) (* t (- (* 18.0 (* x (* y z))) (* 4.0 a)))) t_1)))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = 4.0 * (x * i);
double tmp;
if (z <= -8.2e+179) {
tmp = (((b * c) + ((18.0 * t) * (z * (y * x)))) - (i * (x * 4.0))) - (k * (j * 27.0));
} else if (z <= 3.9e+139) {
tmp = ((b * c) + (-4.0 * (t * a))) - (t_1 + (27.0 * (j * k)));
} else {
tmp = ((b * c) + (t * ((18.0 * (x * (y * z))) - (4.0 * a)))) - t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = 4.0d0 * (x * i)
if (z <= (-8.2d+179)) then
tmp = (((b * c) + ((18.0d0 * t) * (z * (y * x)))) - (i * (x * 4.0d0))) - (k * (j * 27.0d0))
else if (z <= 3.9d+139) then
tmp = ((b * c) + ((-4.0d0) * (t * a))) - (t_1 + (27.0d0 * (j * k)))
else
tmp = ((b * c) + (t * ((18.0d0 * (x * (y * z))) - (4.0d0 * a)))) - t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = 4.0 * (x * i);
double tmp;
if (z <= -8.2e+179) {
tmp = (((b * c) + ((18.0 * t) * (z * (y * x)))) - (i * (x * 4.0))) - (k * (j * 27.0));
} else if (z <= 3.9e+139) {
tmp = ((b * c) + (-4.0 * (t * a))) - (t_1 + (27.0 * (j * k)));
} else {
tmp = ((b * c) + (t * ((18.0 * (x * (y * z))) - (4.0 * a)))) - t_1;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = 4.0 * (x * i) tmp = 0 if z <= -8.2e+179: tmp = (((b * c) + ((18.0 * t) * (z * (y * x)))) - (i * (x * 4.0))) - (k * (j * 27.0)) elif z <= 3.9e+139: tmp = ((b * c) + (-4.0 * (t * a))) - (t_1 + (27.0 * (j * k))) else: tmp = ((b * c) + (t * ((18.0 * (x * (y * z))) - (4.0 * a)))) - t_1 return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(4.0 * Float64(x * i)) tmp = 0.0 if (z <= -8.2e+179) tmp = Float64(Float64(Float64(Float64(b * c) + Float64(Float64(18.0 * t) * Float64(z * Float64(y * x)))) - Float64(i * Float64(x * 4.0))) - Float64(k * Float64(j * 27.0))); elseif (z <= 3.9e+139) tmp = Float64(Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))) - Float64(t_1 + Float64(27.0 * Float64(j * k)))); else tmp = Float64(Float64(Float64(b * c) + Float64(t * Float64(Float64(18.0 * Float64(x * Float64(y * z))) - Float64(4.0 * a)))) - t_1); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = 4.0 * (x * i);
tmp = 0.0;
if (z <= -8.2e+179)
tmp = (((b * c) + ((18.0 * t) * (z * (y * x)))) - (i * (x * 4.0))) - (k * (j * 27.0));
elseif (z <= 3.9e+139)
tmp = ((b * c) + (-4.0 * (t * a))) - (t_1 + (27.0 * (j * k)));
else
tmp = ((b * c) + (t * ((18.0 * (x * (y * z))) - (4.0 * a)))) - t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -8.2e+179], N[(N[(N[(N[(b * c), $MachinePrecision] + N[(N[(18.0 * t), $MachinePrecision] * N[(z * N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(i * N[(x * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(k * N[(j * 27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.9e+139], N[(N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t$95$1 + N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b * c), $MachinePrecision] + N[(t * N[(N[(18.0 * N[(x * N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(4.0 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := 4 \cdot \left(x \cdot i\right)\\
\mathbf{if}\;z \leq -8.2 \cdot 10^{+179}:\\
\;\;\;\;\left(\left(b \cdot c + \left(18 \cdot t\right) \cdot \left(z \cdot \left(y \cdot x\right)\right)\right) - i \cdot \left(x \cdot 4\right)\right) - k \cdot \left(j \cdot 27\right)\\
\mathbf{elif}\;z \leq 3.9 \cdot 10^{+139}:\\
\;\;\;\;\left(b \cdot c + -4 \cdot \left(t \cdot a\right)\right) - \left(t\_1 + 27 \cdot \left(j \cdot k\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(b \cdot c + t \cdot \left(18 \cdot \left(x \cdot \left(y \cdot z\right)\right) - 4 \cdot a\right)\right) - t\_1\\
\end{array}
\end{array}
if z < -8.20000000000000021e179Initial program 72.6%
distribute-rgt-out--83.7%
associate-*r*75.9%
*-commutative75.9%
associate-*l*75.9%
associate-*r*75.9%
Applied egg-rr75.9%
Taylor expanded in x around inf 65.8%
associate-*r*65.8%
associate-*r*73.6%
Simplified73.6%
if -8.20000000000000021e179 < z < 3.90000000000000006e139Initial program 88.4%
Simplified90.7%
Taylor expanded in y around 0 89.4%
if 3.90000000000000006e139 < z Initial program 74.9%
Simplified75.0%
Taylor expanded in j around 0 68.2%
Final simplification83.9%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (+ (* b c) (* -4.0 (* t a)))) (t_2 (* k (* j 27.0))))
(if (<= t_2 -1e-15)
(- (* b c) (* k (+ (* j 27.0) (* 4.0 (/ (* x i) k)))))
(if (<= t_2 2e-9) (- t_1 (* 4.0 (* x i))) (- t_1 (* 27.0 (* j k)))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (b * c) + (-4.0 * (t * a));
double t_2 = k * (j * 27.0);
double tmp;
if (t_2 <= -1e-15) {
tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k))));
} else if (t_2 <= 2e-9) {
tmp = t_1 - (4.0 * (x * i));
} else {
tmp = t_1 - (27.0 * (j * k));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (b * c) + ((-4.0d0) * (t * a))
t_2 = k * (j * 27.0d0)
if (t_2 <= (-1d-15)) then
tmp = (b * c) - (k * ((j * 27.0d0) + (4.0d0 * ((x * i) / k))))
else if (t_2 <= 2d-9) then
tmp = t_1 - (4.0d0 * (x * i))
else
tmp = t_1 - (27.0d0 * (j * k))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (b * c) + (-4.0 * (t * a));
double t_2 = k * (j * 27.0);
double tmp;
if (t_2 <= -1e-15) {
tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k))));
} else if (t_2 <= 2e-9) {
tmp = t_1 - (4.0 * (x * i));
} else {
tmp = t_1 - (27.0 * (j * k));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = (b * c) + (-4.0 * (t * a)) t_2 = k * (j * 27.0) tmp = 0 if t_2 <= -1e-15: tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k)))) elif t_2 <= 2e-9: tmp = t_1 - (4.0 * (x * i)) else: tmp = t_1 - (27.0 * (j * k)) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))) t_2 = Float64(k * Float64(j * 27.0)) tmp = 0.0 if (t_2 <= -1e-15) tmp = Float64(Float64(b * c) - Float64(k * Float64(Float64(j * 27.0) + Float64(4.0 * Float64(Float64(x * i) / k))))); elseif (t_2 <= 2e-9) tmp = Float64(t_1 - Float64(4.0 * Float64(x * i))); else tmp = Float64(t_1 - Float64(27.0 * Float64(j * k))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = (b * c) + (-4.0 * (t * a));
t_2 = k * (j * 27.0);
tmp = 0.0;
if (t_2 <= -1e-15)
tmp = (b * c) - (k * ((j * 27.0) + (4.0 * ((x * i) / k))));
elseif (t_2 <= 2e-9)
tmp = t_1 - (4.0 * (x * i));
else
tmp = t_1 - (27.0 * (j * k));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(k * N[(j * 27.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -1e-15], N[(N[(b * c), $MachinePrecision] - N[(k * N[(N[(j * 27.0), $MachinePrecision] + N[(4.0 * N[(N[(x * i), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 2e-9], N[(t$95$1 - N[(4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 - N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := b \cdot c + -4 \cdot \left(t \cdot a\right)\\
t_2 := k \cdot \left(j \cdot 27\right)\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{-15}:\\
\;\;\;\;b \cdot c - k \cdot \left(j \cdot 27 + 4 \cdot \frac{x \cdot i}{k}\right)\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{-9}:\\
\;\;\;\;t\_1 - 4 \cdot \left(x \cdot i\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1 - 27 \cdot \left(j \cdot k\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 j #s(literal 27 binary64)) k) < -1.0000000000000001e-15Initial program 74.2%
Simplified76.0%
Taylor expanded in t around 0 70.5%
Taylor expanded in k around inf 72.2%
if -1.0000000000000001e-15 < (*.f64 (*.f64 j #s(literal 27 binary64)) k) < 2.00000000000000012e-9Initial program 86.6%
distribute-rgt-out--89.6%
associate-*r*87.6%
*-commutative87.6%
associate-*l*87.6%
associate-*r*87.6%
Applied egg-rr87.6%
Taylor expanded in y around 0 79.6%
Taylor expanded in j around 0 78.9%
if 2.00000000000000012e-9 < (*.f64 (*.f64 j #s(literal 27 binary64)) k) Initial program 88.6%
Simplified93.4%
Taylor expanded in x around 0 72.6%
Final simplification75.8%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (+ (* b c) (* -4.0 (* t a)))))
(if (<= (* b c) -5e+85)
t_1
(if (<= (* b c) 1e-321)
(* -4.0 (+ (* x i) (* t a)))
(if (<= (* b c) 7.8e+161)
(+ (* j (* k -27.0)) (* -4.0 (* x i)))
t_1)))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (b * c) + (-4.0 * (t * a));
double tmp;
if ((b * c) <= -5e+85) {
tmp = t_1;
} else if ((b * c) <= 1e-321) {
tmp = -4.0 * ((x * i) + (t * a));
} else if ((b * c) <= 7.8e+161) {
tmp = (j * (k * -27.0)) + (-4.0 * (x * i));
} else {
tmp = t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = (b * c) + ((-4.0d0) * (t * a))
if ((b * c) <= (-5d+85)) then
tmp = t_1
else if ((b * c) <= 1d-321) then
tmp = (-4.0d0) * ((x * i) + (t * a))
else if ((b * c) <= 7.8d+161) then
tmp = (j * (k * (-27.0d0))) + ((-4.0d0) * (x * i))
else
tmp = t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (b * c) + (-4.0 * (t * a));
double tmp;
if ((b * c) <= -5e+85) {
tmp = t_1;
} else if ((b * c) <= 1e-321) {
tmp = -4.0 * ((x * i) + (t * a));
} else if ((b * c) <= 7.8e+161) {
tmp = (j * (k * -27.0)) + (-4.0 * (x * i));
} else {
tmp = t_1;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = (b * c) + (-4.0 * (t * a)) tmp = 0 if (b * c) <= -5e+85: tmp = t_1 elif (b * c) <= 1e-321: tmp = -4.0 * ((x * i) + (t * a)) elif (b * c) <= 7.8e+161: tmp = (j * (k * -27.0)) + (-4.0 * (x * i)) else: tmp = t_1 return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))) tmp = 0.0 if (Float64(b * c) <= -5e+85) tmp = t_1; elseif (Float64(b * c) <= 1e-321) tmp = Float64(-4.0 * Float64(Float64(x * i) + Float64(t * a))); elseif (Float64(b * c) <= 7.8e+161) tmp = Float64(Float64(j * Float64(k * -27.0)) + Float64(-4.0 * Float64(x * i))); else tmp = t_1; end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = (b * c) + (-4.0 * (t * a));
tmp = 0.0;
if ((b * c) <= -5e+85)
tmp = t_1;
elseif ((b * c) <= 1e-321)
tmp = -4.0 * ((x * i) + (t * a));
elseif ((b * c) <= 7.8e+161)
tmp = (j * (k * -27.0)) + (-4.0 * (x * i));
else
tmp = t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(b * c), $MachinePrecision], -5e+85], t$95$1, If[LessEqual[N[(b * c), $MachinePrecision], 1e-321], N[(-4.0 * N[(N[(x * i), $MachinePrecision] + N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 7.8e+161], N[(N[(j * N[(k * -27.0), $MachinePrecision]), $MachinePrecision] + N[(-4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := b \cdot c + -4 \cdot \left(t \cdot a\right)\\
\mathbf{if}\;b \cdot c \leq -5 \cdot 10^{+85}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \cdot c \leq 10^{-321}:\\
\;\;\;\;-4 \cdot \left(x \cdot i + t \cdot a\right)\\
\mathbf{elif}\;b \cdot c \leq 7.8 \cdot 10^{+161}:\\
\;\;\;\;j \cdot \left(k \cdot -27\right) + -4 \cdot \left(x \cdot i\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (*.f64 b c) < -5.0000000000000001e85 or 7.8000000000000004e161 < (*.f64 b c) Initial program 77.4%
distribute-rgt-out--80.6%
associate-*r*79.5%
*-commutative79.5%
associate-*l*79.5%
associate-*r*79.5%
Applied egg-rr79.5%
Taylor expanded in y around 0 79.4%
Taylor expanded in j around 0 79.1%
Taylor expanded in i around 0 71.3%
if -5.0000000000000001e85 < (*.f64 b c) < 9.98013e-322Initial program 86.1%
distribute-rgt-out--88.2%
associate-*r*86.3%
*-commutative86.3%
associate-*l*86.3%
associate-*r*86.3%
Applied egg-rr86.3%
Taylor expanded in y around 0 76.2%
Taylor expanded in j around 0 59.5%
Taylor expanded in b around 0 58.2%
cancel-sign-sub-inv58.2%
metadata-eval58.2%
distribute-lft-out58.2%
*-commutative58.2%
*-commutative58.2%
Simplified58.2%
if 9.98013e-322 < (*.f64 b c) < 7.8000000000000004e161Initial program 90.8%
Simplified94.1%
Taylor expanded in i around inf 62.6%
Final simplification64.3%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* 4.0 (* x i))))
(if (or (<= (* 4.0 a) -5e-29) (not (<= (* 4.0 a) 1e-29)))
(- (+ (* b c) (* -4.0 (* t a))) t_1)
(- (* b c) (+ t_1 (* 27.0 (* j k)))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = 4.0 * (x * i);
double tmp;
if (((4.0 * a) <= -5e-29) || !((4.0 * a) <= 1e-29)) {
tmp = ((b * c) + (-4.0 * (t * a))) - t_1;
} else {
tmp = (b * c) - (t_1 + (27.0 * (j * k)));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = 4.0d0 * (x * i)
if (((4.0d0 * a) <= (-5d-29)) .or. (.not. ((4.0d0 * a) <= 1d-29))) then
tmp = ((b * c) + ((-4.0d0) * (t * a))) - t_1
else
tmp = (b * c) - (t_1 + (27.0d0 * (j * k)))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = 4.0 * (x * i);
double tmp;
if (((4.0 * a) <= -5e-29) || !((4.0 * a) <= 1e-29)) {
tmp = ((b * c) + (-4.0 * (t * a))) - t_1;
} else {
tmp = (b * c) - (t_1 + (27.0 * (j * k)));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = 4.0 * (x * i) tmp = 0 if ((4.0 * a) <= -5e-29) or not ((4.0 * a) <= 1e-29): tmp = ((b * c) + (-4.0 * (t * a))) - t_1 else: tmp = (b * c) - (t_1 + (27.0 * (j * k))) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(4.0 * Float64(x * i)) tmp = 0.0 if ((Float64(4.0 * a) <= -5e-29) || !(Float64(4.0 * a) <= 1e-29)) tmp = Float64(Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))) - t_1); else tmp = Float64(Float64(b * c) - Float64(t_1 + Float64(27.0 * Float64(j * k)))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = 4.0 * (x * i);
tmp = 0.0;
if (((4.0 * a) <= -5e-29) || ~(((4.0 * a) <= 1e-29)))
tmp = ((b * c) + (-4.0 * (t * a))) - t_1;
else
tmp = (b * c) - (t_1 + (27.0 * (j * k)));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[N[(4.0 * a), $MachinePrecision], -5e-29], N[Not[LessEqual[N[(4.0 * a), $MachinePrecision], 1e-29]], $MachinePrecision]], N[(N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision], N[(N[(b * c), $MachinePrecision] - N[(t$95$1 + N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := 4 \cdot \left(x \cdot i\right)\\
\mathbf{if}\;4 \cdot a \leq -5 \cdot 10^{-29} \lor \neg \left(4 \cdot a \leq 10^{-29}\right):\\
\;\;\;\;\left(b \cdot c + -4 \cdot \left(t \cdot a\right)\right) - t\_1\\
\mathbf{else}:\\
\;\;\;\;b \cdot c - \left(t\_1 + 27 \cdot \left(j \cdot k\right)\right)\\
\end{array}
\end{array}
if (*.f64 a #s(literal 4 binary64)) < -4.99999999999999986e-29 or 9.99999999999999943e-30 < (*.f64 a #s(literal 4 binary64)) Initial program 82.4%
distribute-rgt-out--86.9%
associate-*r*86.9%
*-commutative86.9%
associate-*l*86.9%
associate-*r*86.9%
Applied egg-rr86.9%
Taylor expanded in y around 0 78.0%
Taylor expanded in j around 0 72.2%
if -4.99999999999999986e-29 < (*.f64 a #s(literal 4 binary64)) < 9.99999999999999943e-30Initial program 85.9%
Simplified84.5%
Taylor expanded in t around 0 77.3%
Final simplification74.6%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* 4.0 (* x i))))
(if (<= z 4.2e+140)
(- (+ (* b c) (* -4.0 (* t a))) (+ t_1 (* 27.0 (* j k))))
(- (+ (* b c) (* t (- (* 18.0 (* x (* y z))) (* 4.0 a)))) t_1))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = 4.0 * (x * i);
double tmp;
if (z <= 4.2e+140) {
tmp = ((b * c) + (-4.0 * (t * a))) - (t_1 + (27.0 * (j * k)));
} else {
tmp = ((b * c) + (t * ((18.0 * (x * (y * z))) - (4.0 * a)))) - t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = 4.0d0 * (x * i)
if (z <= 4.2d+140) then
tmp = ((b * c) + ((-4.0d0) * (t * a))) - (t_1 + (27.0d0 * (j * k)))
else
tmp = ((b * c) + (t * ((18.0d0 * (x * (y * z))) - (4.0d0 * a)))) - t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = 4.0 * (x * i);
double tmp;
if (z <= 4.2e+140) {
tmp = ((b * c) + (-4.0 * (t * a))) - (t_1 + (27.0 * (j * k)));
} else {
tmp = ((b * c) + (t * ((18.0 * (x * (y * z))) - (4.0 * a)))) - t_1;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = 4.0 * (x * i) tmp = 0 if z <= 4.2e+140: tmp = ((b * c) + (-4.0 * (t * a))) - (t_1 + (27.0 * (j * k))) else: tmp = ((b * c) + (t * ((18.0 * (x * (y * z))) - (4.0 * a)))) - t_1 return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(4.0 * Float64(x * i)) tmp = 0.0 if (z <= 4.2e+140) tmp = Float64(Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))) - Float64(t_1 + Float64(27.0 * Float64(j * k)))); else tmp = Float64(Float64(Float64(b * c) + Float64(t * Float64(Float64(18.0 * Float64(x * Float64(y * z))) - Float64(4.0 * a)))) - t_1); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = 4.0 * (x * i);
tmp = 0.0;
if (z <= 4.2e+140)
tmp = ((b * c) + (-4.0 * (t * a))) - (t_1 + (27.0 * (j * k)));
else
tmp = ((b * c) + (t * ((18.0 * (x * (y * z))) - (4.0 * a)))) - t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, 4.2e+140], N[(N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t$95$1 + N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b * c), $MachinePrecision] + N[(t * N[(N[(18.0 * N[(x * N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(4.0 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := 4 \cdot \left(x \cdot i\right)\\
\mathbf{if}\;z \leq 4.2 \cdot 10^{+140}:\\
\;\;\;\;\left(b \cdot c + -4 \cdot \left(t \cdot a\right)\right) - \left(t\_1 + 27 \cdot \left(j \cdot k\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(b \cdot c + t \cdot \left(18 \cdot \left(x \cdot \left(y \cdot z\right)\right) - 4 \cdot a\right)\right) - t\_1\\
\end{array}
\end{array}
if z < 4.2000000000000004e140Initial program 85.8%
Simplified88.2%
Taylor expanded in y around 0 83.1%
if 4.2000000000000004e140 < z Initial program 74.9%
Simplified75.0%
Taylor expanded in j around 0 68.2%
Final simplification80.8%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= (* b c) -5.2e+139)
(* b c)
(if (<= (* b c) -2.6e-272)
(* -4.0 (* t a))
(if (<= (* b c) 2.6e+21) (* (* j k) -27.0) (* b c)))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((b * c) <= -5.2e+139) {
tmp = b * c;
} else if ((b * c) <= -2.6e-272) {
tmp = -4.0 * (t * a);
} else if ((b * c) <= 2.6e+21) {
tmp = (j * k) * -27.0;
} else {
tmp = b * c;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if ((b * c) <= (-5.2d+139)) then
tmp = b * c
else if ((b * c) <= (-2.6d-272)) then
tmp = (-4.0d0) * (t * a)
else if ((b * c) <= 2.6d+21) then
tmp = (j * k) * (-27.0d0)
else
tmp = b * c
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((b * c) <= -5.2e+139) {
tmp = b * c;
} else if ((b * c) <= -2.6e-272) {
tmp = -4.0 * (t * a);
} else if ((b * c) <= 2.6e+21) {
tmp = (j * k) * -27.0;
} else {
tmp = b * c;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if (b * c) <= -5.2e+139: tmp = b * c elif (b * c) <= -2.6e-272: tmp = -4.0 * (t * a) elif (b * c) <= 2.6e+21: tmp = (j * k) * -27.0 else: tmp = b * c return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (Float64(b * c) <= -5.2e+139) tmp = Float64(b * c); elseif (Float64(b * c) <= -2.6e-272) tmp = Float64(-4.0 * Float64(t * a)); elseif (Float64(b * c) <= 2.6e+21) tmp = Float64(Float64(j * k) * -27.0); else tmp = Float64(b * c); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if ((b * c) <= -5.2e+139)
tmp = b * c;
elseif ((b * c) <= -2.6e-272)
tmp = -4.0 * (t * a);
elseif ((b * c) <= 2.6e+21)
tmp = (j * k) * -27.0;
else
tmp = b * c;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[N[(b * c), $MachinePrecision], -5.2e+139], N[(b * c), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], -2.6e-272], N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 2.6e+21], N[(N[(j * k), $MachinePrecision] * -27.0), $MachinePrecision], N[(b * c), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -5.2 \cdot 10^{+139}:\\
\;\;\;\;b \cdot c\\
\mathbf{elif}\;b \cdot c \leq -2.6 \cdot 10^{-272}:\\
\;\;\;\;-4 \cdot \left(t \cdot a\right)\\
\mathbf{elif}\;b \cdot c \leq 2.6 \cdot 10^{+21}:\\
\;\;\;\;\left(j \cdot k\right) \cdot -27\\
\mathbf{else}:\\
\;\;\;\;b \cdot c\\
\end{array}
\end{array}
if (*.f64 b c) < -5.20000000000000044e139 or 2.6e21 < (*.f64 b c) Initial program 79.4%
distribute-rgt-out--82.2%
associate-*r*82.2%
*-commutative82.2%
associate-*l*82.2%
associate-*r*82.2%
Applied egg-rr82.2%
Taylor expanded in b around inf 54.5%
if -5.20000000000000044e139 < (*.f64 b c) < -2.59999999999999992e-272Initial program 87.2%
distribute-rgt-out--90.1%
associate-*r*90.2%
*-commutative90.2%
associate-*l*90.2%
associate-*r*90.2%
Applied egg-rr90.2%
Taylor expanded in y around 0 79.4%
Taylor expanded in a around inf 31.9%
*-commutative31.9%
Simplified31.9%
if -2.59999999999999992e-272 < (*.f64 b c) < 2.6e21Initial program 87.5%
Simplified86.6%
Taylor expanded in j around inf 35.1%
Final simplification42.1%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (+ (* 4.0 (* x i)) (* 27.0 (* j k)))))
(if (<= y -8.8e+256)
(- (* 18.0 (* t (* x (* y z)))) t_1)
(- (+ (* b c) (* -4.0 (* t a))) t_1))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (4.0 * (x * i)) + (27.0 * (j * k));
double tmp;
if (y <= -8.8e+256) {
tmp = (18.0 * (t * (x * (y * z)))) - t_1;
} else {
tmp = ((b * c) + (-4.0 * (t * a))) - t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = (4.0d0 * (x * i)) + (27.0d0 * (j * k))
if (y <= (-8.8d+256)) then
tmp = (18.0d0 * (t * (x * (y * z)))) - t_1
else
tmp = ((b * c) + ((-4.0d0) * (t * a))) - t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (4.0 * (x * i)) + (27.0 * (j * k));
double tmp;
if (y <= -8.8e+256) {
tmp = (18.0 * (t * (x * (y * z)))) - t_1;
} else {
tmp = ((b * c) + (-4.0 * (t * a))) - t_1;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = (4.0 * (x * i)) + (27.0 * (j * k)) tmp = 0 if y <= -8.8e+256: tmp = (18.0 * (t * (x * (y * z)))) - t_1 else: tmp = ((b * c) + (-4.0 * (t * a))) - t_1 return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(4.0 * Float64(x * i)) + Float64(27.0 * Float64(j * k))) tmp = 0.0 if (y <= -8.8e+256) tmp = Float64(Float64(18.0 * Float64(t * Float64(x * Float64(y * z)))) - t_1); else tmp = Float64(Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))) - t_1); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = (4.0 * (x * i)) + (27.0 * (j * k));
tmp = 0.0;
if (y <= -8.8e+256)
tmp = (18.0 * (t * (x * (y * z)))) - t_1;
else
tmp = ((b * c) + (-4.0 * (t * a))) - t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision] + N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -8.8e+256], N[(N[(18.0 * N[(t * N[(x * N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision], N[(N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := 4 \cdot \left(x \cdot i\right) + 27 \cdot \left(j \cdot k\right)\\
\mathbf{if}\;y \leq -8.8 \cdot 10^{+256}:\\
\;\;\;\;18 \cdot \left(t \cdot \left(x \cdot \left(y \cdot z\right)\right)\right) - t\_1\\
\mathbf{else}:\\
\;\;\;\;\left(b \cdot c + -4 \cdot \left(t \cdot a\right)\right) - t\_1\\
\end{array}
\end{array}
if y < -8.7999999999999997e256Initial program 71.4%
Simplified78.6%
Taylor expanded in a around 0 78.6%
Taylor expanded in b around 0 78.6%
if -8.7999999999999997e256 < y Initial program 84.8%
Simplified86.6%
Taylor expanded in y around 0 81.6%
Final simplification81.4%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. (FPCore (x y z t a b c i j k) :precision binary64 (if (or (<= x -1.08e+83) (not (<= x 5e+60))) (* x (+ (* (* 18.0 t) (* y z)) (* i -4.0))) (- (+ (* b c) (* -4.0 (* t a))) (* 27.0 (* j k)))))
assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((x <= -1.08e+83) || !(x <= 5e+60)) {
tmp = x * (((18.0 * t) * (y * z)) + (i * -4.0));
} else {
tmp = ((b * c) + (-4.0 * (t * a))) - (27.0 * (j * k));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if ((x <= (-1.08d+83)) .or. (.not. (x <= 5d+60))) then
tmp = x * (((18.0d0 * t) * (y * z)) + (i * (-4.0d0)))
else
tmp = ((b * c) + ((-4.0d0) * (t * a))) - (27.0d0 * (j * k))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((x <= -1.08e+83) || !(x <= 5e+60)) {
tmp = x * (((18.0 * t) * (y * z)) + (i * -4.0));
} else {
tmp = ((b * c) + (-4.0 * (t * a))) - (27.0 * (j * k));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if (x <= -1.08e+83) or not (x <= 5e+60): tmp = x * (((18.0 * t) * (y * z)) + (i * -4.0)) else: tmp = ((b * c) + (-4.0 * (t * a))) - (27.0 * (j * k)) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if ((x <= -1.08e+83) || !(x <= 5e+60)) tmp = Float64(x * Float64(Float64(Float64(18.0 * t) * Float64(y * z)) + Float64(i * -4.0))); else tmp = Float64(Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))) - Float64(27.0 * Float64(j * k))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if ((x <= -1.08e+83) || ~((x <= 5e+60)))
tmp = x * (((18.0 * t) * (y * z)) + (i * -4.0));
else
tmp = ((b * c) + (-4.0 * (t * a))) - (27.0 * (j * k));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[Or[LessEqual[x, -1.08e+83], N[Not[LessEqual[x, 5e+60]], $MachinePrecision]], N[(x * N[(N[(N[(18.0 * t), $MachinePrecision] * N[(y * z), $MachinePrecision]), $MachinePrecision] + N[(i * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.08 \cdot 10^{+83} \lor \neg \left(x \leq 5 \cdot 10^{+60}\right):\\
\;\;\;\;x \cdot \left(\left(18 \cdot t\right) \cdot \left(y \cdot z\right) + i \cdot -4\right)\\
\mathbf{else}:\\
\;\;\;\;\left(b \cdot c + -4 \cdot \left(t \cdot a\right)\right) - 27 \cdot \left(j \cdot k\right)\\
\end{array}
\end{array}
if x < -1.08e83 or 4.99999999999999975e60 < x Initial program 73.4%
Simplified82.0%
Taylor expanded in a around 0 73.8%
Taylor expanded in b around 0 64.3%
Taylor expanded in x around inf 73.6%
cancel-sign-sub-inv73.6%
associate-*r*73.6%
*-commutative73.6%
metadata-eval73.6%
*-commutative73.6%
Simplified73.6%
if -1.08e83 < x < 4.99999999999999975e60Initial program 91.4%
Simplified89.0%
Taylor expanded in x around 0 76.8%
Final simplification75.5%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* 4.0 (- a))))
(if (<= t -1.4e+128)
(* t (- t_1 (* (* x (* y z)) -18.0)))
(if (<= t 2.1e+54)
(- (* b c) (+ (* 4.0 (* x i)) (* 27.0 (* j k))))
(* t (- t_1 (* (* y x) (* z -18.0))))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = 4.0 * -a;
double tmp;
if (t <= -1.4e+128) {
tmp = t * (t_1 - ((x * (y * z)) * -18.0));
} else if (t <= 2.1e+54) {
tmp = (b * c) - ((4.0 * (x * i)) + (27.0 * (j * k)));
} else {
tmp = t * (t_1 - ((y * x) * (z * -18.0)));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = 4.0d0 * -a
if (t <= (-1.4d+128)) then
tmp = t * (t_1 - ((x * (y * z)) * (-18.0d0)))
else if (t <= 2.1d+54) then
tmp = (b * c) - ((4.0d0 * (x * i)) + (27.0d0 * (j * k)))
else
tmp = t * (t_1 - ((y * x) * (z * (-18.0d0))))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = 4.0 * -a;
double tmp;
if (t <= -1.4e+128) {
tmp = t * (t_1 - ((x * (y * z)) * -18.0));
} else if (t <= 2.1e+54) {
tmp = (b * c) - ((4.0 * (x * i)) + (27.0 * (j * k)));
} else {
tmp = t * (t_1 - ((y * x) * (z * -18.0)));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = 4.0 * -a tmp = 0 if t <= -1.4e+128: tmp = t * (t_1 - ((x * (y * z)) * -18.0)) elif t <= 2.1e+54: tmp = (b * c) - ((4.0 * (x * i)) + (27.0 * (j * k))) else: tmp = t * (t_1 - ((y * x) * (z * -18.0))) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(4.0 * Float64(-a)) tmp = 0.0 if (t <= -1.4e+128) tmp = Float64(t * Float64(t_1 - Float64(Float64(x * Float64(y * z)) * -18.0))); elseif (t <= 2.1e+54) tmp = Float64(Float64(b * c) - Float64(Float64(4.0 * Float64(x * i)) + Float64(27.0 * Float64(j * k)))); else tmp = Float64(t * Float64(t_1 - Float64(Float64(y * x) * Float64(z * -18.0)))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = 4.0 * -a;
tmp = 0.0;
if (t <= -1.4e+128)
tmp = t * (t_1 - ((x * (y * z)) * -18.0));
elseif (t <= 2.1e+54)
tmp = (b * c) - ((4.0 * (x * i)) + (27.0 * (j * k)));
else
tmp = t * (t_1 - ((y * x) * (z * -18.0)));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(4.0 * (-a)), $MachinePrecision]}, If[LessEqual[t, -1.4e+128], N[(t * N[(t$95$1 - N[(N[(x * N[(y * z), $MachinePrecision]), $MachinePrecision] * -18.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.1e+54], N[(N[(b * c), $MachinePrecision] - N[(N[(4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision] + N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t * N[(t$95$1 - N[(N[(y * x), $MachinePrecision] * N[(z * -18.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := 4 \cdot \left(-a\right)\\
\mathbf{if}\;t \leq -1.4 \cdot 10^{+128}:\\
\;\;\;\;t \cdot \left(t\_1 - \left(x \cdot \left(y \cdot z\right)\right) \cdot -18\right)\\
\mathbf{elif}\;t \leq 2.1 \cdot 10^{+54}:\\
\;\;\;\;b \cdot c - \left(4 \cdot \left(x \cdot i\right) + 27 \cdot \left(j \cdot k\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t \cdot \left(t\_1 - \left(y \cdot x\right) \cdot \left(z \cdot -18\right)\right)\\
\end{array}
\end{array}
if t < -1.39999999999999991e128Initial program 85.6%
Taylor expanded in t around -inf 80.7%
associate-*r*80.7%
neg-mul-180.7%
cancel-sign-sub-inv80.7%
metadata-eval80.7%
*-commutative80.7%
*-commutative80.7%
Simplified80.7%
if -1.39999999999999991e128 < t < 2.09999999999999986e54Initial program 84.9%
Simplified83.9%
Taylor expanded in t around 0 72.2%
if 2.09999999999999986e54 < t Initial program 79.5%
Taylor expanded in t around -inf 65.9%
associate-*r*65.9%
neg-mul-165.9%
cancel-sign-sub-inv65.9%
metadata-eval65.9%
*-commutative65.9%
*-commutative65.9%
Simplified65.9%
Taylor expanded in x around 0 65.9%
*-commutative65.9%
associate-*r*65.9%
associate-*l*65.9%
associate-*r*66.0%
Simplified66.0%
Final simplification72.3%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. (FPCore (x y z t a b c i j k) :precision binary64 (if (or (<= (* b c) -5.2e+87) (not (<= (* b c) 7.5e+151))) (+ (* b c) (* -4.0 (* t a))) (* -4.0 (+ (* x i) (* t a)))))
assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -5.2e+87) || !((b * c) <= 7.5e+151)) {
tmp = (b * c) + (-4.0 * (t * a));
} else {
tmp = -4.0 * ((x * i) + (t * a));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if (((b * c) <= (-5.2d+87)) .or. (.not. ((b * c) <= 7.5d+151))) then
tmp = (b * c) + ((-4.0d0) * (t * a))
else
tmp = (-4.0d0) * ((x * i) + (t * a))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -5.2e+87) || !((b * c) <= 7.5e+151)) {
tmp = (b * c) + (-4.0 * (t * a));
} else {
tmp = -4.0 * ((x * i) + (t * a));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if ((b * c) <= -5.2e+87) or not ((b * c) <= 7.5e+151): tmp = (b * c) + (-4.0 * (t * a)) else: tmp = -4.0 * ((x * i) + (t * a)) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if ((Float64(b * c) <= -5.2e+87) || !(Float64(b * c) <= 7.5e+151)) tmp = Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))); else tmp = Float64(-4.0 * Float64(Float64(x * i) + Float64(t * a))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if (((b * c) <= -5.2e+87) || ~(((b * c) <= 7.5e+151)))
tmp = (b * c) + (-4.0 * (t * a));
else
tmp = -4.0 * ((x * i) + (t * a));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[Or[LessEqual[N[(b * c), $MachinePrecision], -5.2e+87], N[Not[LessEqual[N[(b * c), $MachinePrecision], 7.5e+151]], $MachinePrecision]], N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-4.0 * N[(N[(x * i), $MachinePrecision] + N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -5.2 \cdot 10^{+87} \lor \neg \left(b \cdot c \leq 7.5 \cdot 10^{+151}\right):\\
\;\;\;\;b \cdot c + -4 \cdot \left(t \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;-4 \cdot \left(x \cdot i + t \cdot a\right)\\
\end{array}
\end{array}
if (*.f64 b c) < -5.19999999999999997e87 or 7.49999999999999977e151 < (*.f64 b c) Initial program 77.9%
distribute-rgt-out--81.0%
associate-*r*79.9%
*-commutative79.9%
associate-*l*79.9%
associate-*r*79.9%
Applied egg-rr79.9%
Taylor expanded in y around 0 79.8%
Taylor expanded in j around 0 78.5%
Taylor expanded in i around 0 70.8%
if -5.19999999999999997e87 < (*.f64 b c) < 7.49999999999999977e151Initial program 87.9%
distribute-rgt-out--89.8%
associate-*r*89.4%
*-commutative89.4%
associate-*l*89.4%
associate-*r*89.4%
Applied egg-rr89.4%
Taylor expanded in y around 0 78.2%
Taylor expanded in j around 0 56.7%
Taylor expanded in b around 0 53.6%
cancel-sign-sub-inv53.6%
metadata-eval53.6%
distribute-lft-out53.6%
*-commutative53.6%
*-commutative53.6%
Simplified53.6%
Final simplification60.2%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. (FPCore (x y z t a b c i j k) :precision binary64 (if (or (<= (* b c) -1.16e+244) (not (<= (* b c) 4.7e+232))) (* b c) (* -4.0 (+ (* x i) (* t a)))))
assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -1.16e+244) || !((b * c) <= 4.7e+232)) {
tmp = b * c;
} else {
tmp = -4.0 * ((x * i) + (t * a));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if (((b * c) <= (-1.16d+244)) .or. (.not. ((b * c) <= 4.7d+232))) then
tmp = b * c
else
tmp = (-4.0d0) * ((x * i) + (t * a))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -1.16e+244) || !((b * c) <= 4.7e+232)) {
tmp = b * c;
} else {
tmp = -4.0 * ((x * i) + (t * a));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if ((b * c) <= -1.16e+244) or not ((b * c) <= 4.7e+232): tmp = b * c else: tmp = -4.0 * ((x * i) + (t * a)) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if ((Float64(b * c) <= -1.16e+244) || !(Float64(b * c) <= 4.7e+232)) tmp = Float64(b * c); else tmp = Float64(-4.0 * Float64(Float64(x * i) + Float64(t * a))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if (((b * c) <= -1.16e+244) || ~(((b * c) <= 4.7e+232)))
tmp = b * c;
else
tmp = -4.0 * ((x * i) + (t * a));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[Or[LessEqual[N[(b * c), $MachinePrecision], -1.16e+244], N[Not[LessEqual[N[(b * c), $MachinePrecision], 4.7e+232]], $MachinePrecision]], N[(b * c), $MachinePrecision], N[(-4.0 * N[(N[(x * i), $MachinePrecision] + N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -1.16 \cdot 10^{+244} \lor \neg \left(b \cdot c \leq 4.7 \cdot 10^{+232}\right):\\
\;\;\;\;b \cdot c\\
\mathbf{else}:\\
\;\;\;\;-4 \cdot \left(x \cdot i + t \cdot a\right)\\
\end{array}
\end{array}
if (*.f64 b c) < -1.16000000000000003e244 or 4.69999999999999992e232 < (*.f64 b c) Initial program 75.0%
distribute-rgt-out--76.9%
associate-*r*75.0%
*-commutative75.0%
associate-*l*75.0%
associate-*r*75.0%
Applied egg-rr75.0%
Taylor expanded in b around inf 79.2%
if -1.16000000000000003e244 < (*.f64 b c) < 4.69999999999999992e232Initial program 86.4%
distribute-rgt-out--88.9%
associate-*r*88.5%
*-commutative88.5%
associate-*l*88.5%
associate-*r*88.5%
Applied egg-rr88.5%
Taylor expanded in y around 0 80.2%
Taylor expanded in j around 0 61.7%
Taylor expanded in b around 0 52.6%
cancel-sign-sub-inv52.6%
metadata-eval52.6%
distribute-lft-out52.6%
*-commutative52.6%
*-commutative52.6%
Simplified52.6%
Final simplification58.0%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= (* b c) -4.9e+88)
(+ (* b c) (* -4.0 (* t a)))
(if (<= (* b c) 4.8e+42)
(* -4.0 (+ (* x i) (* t a)))
(- (* b c) (* x (* 4.0 i))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((b * c) <= -4.9e+88) {
tmp = (b * c) + (-4.0 * (t * a));
} else if ((b * c) <= 4.8e+42) {
tmp = -4.0 * ((x * i) + (t * a));
} else {
tmp = (b * c) - (x * (4.0 * i));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if ((b * c) <= (-4.9d+88)) then
tmp = (b * c) + ((-4.0d0) * (t * a))
else if ((b * c) <= 4.8d+42) then
tmp = (-4.0d0) * ((x * i) + (t * a))
else
tmp = (b * c) - (x * (4.0d0 * i))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((b * c) <= -4.9e+88) {
tmp = (b * c) + (-4.0 * (t * a));
} else if ((b * c) <= 4.8e+42) {
tmp = -4.0 * ((x * i) + (t * a));
} else {
tmp = (b * c) - (x * (4.0 * i));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if (b * c) <= -4.9e+88: tmp = (b * c) + (-4.0 * (t * a)) elif (b * c) <= 4.8e+42: tmp = -4.0 * ((x * i) + (t * a)) else: tmp = (b * c) - (x * (4.0 * i)) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (Float64(b * c) <= -4.9e+88) tmp = Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))); elseif (Float64(b * c) <= 4.8e+42) tmp = Float64(-4.0 * Float64(Float64(x * i) + Float64(t * a))); else tmp = Float64(Float64(b * c) - Float64(x * Float64(4.0 * i))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if ((b * c) <= -4.9e+88)
tmp = (b * c) + (-4.0 * (t * a));
elseif ((b * c) <= 4.8e+42)
tmp = -4.0 * ((x * i) + (t * a));
else
tmp = (b * c) - (x * (4.0 * i));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[N[(b * c), $MachinePrecision], -4.9e+88], N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 4.8e+42], N[(-4.0 * N[(N[(x * i), $MachinePrecision] + N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(b * c), $MachinePrecision] - N[(x * N[(4.0 * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -4.9 \cdot 10^{+88}:\\
\;\;\;\;b \cdot c + -4 \cdot \left(t \cdot a\right)\\
\mathbf{elif}\;b \cdot c \leq 4.8 \cdot 10^{+42}:\\
\;\;\;\;-4 \cdot \left(x \cdot i + t \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot c - x \cdot \left(4 \cdot i\right)\\
\end{array}
\end{array}
if (*.f64 b c) < -4.9000000000000002e88Initial program 76.4%
distribute-rgt-out--81.6%
associate-*r*81.6%
*-commutative81.6%
associate-*l*81.6%
associate-*r*81.6%
Applied egg-rr81.6%
Taylor expanded in y around 0 79.6%
Taylor expanded in j around 0 81.7%
Taylor expanded in i around 0 71.9%
if -4.9000000000000002e88 < (*.f64 b c) < 4.7999999999999997e42Initial program 88.1%
distribute-rgt-out--89.5%
associate-*r*88.2%
*-commutative88.2%
associate-*l*88.3%
associate-*r*88.3%
Applied egg-rr88.3%
Taylor expanded in y around 0 79.3%
Taylor expanded in j around 0 57.0%
Taylor expanded in b around 0 55.5%
cancel-sign-sub-inv55.5%
metadata-eval55.5%
distribute-lft-out55.5%
*-commutative55.5%
*-commutative55.5%
Simplified55.5%
if 4.7999999999999997e42 < (*.f64 b c) Initial program 81.8%
Simplified85.4%
Taylor expanded in t around 0 70.7%
Taylor expanded in i around inf 60.8%
associate-*r*60.8%
*-commutative60.8%
Simplified60.8%
Final simplification60.3%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= (* b c) -6.6e+88)
(+ (* b c) (* -4.0 (* t a)))
(if (<= (* b c) 4.7e+147)
(* -4.0 (+ (* x i) (* t a)))
(- (* b c) (* 27.0 (* j k))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((b * c) <= -6.6e+88) {
tmp = (b * c) + (-4.0 * (t * a));
} else if ((b * c) <= 4.7e+147) {
tmp = -4.0 * ((x * i) + (t * a));
} else {
tmp = (b * c) - (27.0 * (j * k));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if ((b * c) <= (-6.6d+88)) then
tmp = (b * c) + ((-4.0d0) * (t * a))
else if ((b * c) <= 4.7d+147) then
tmp = (-4.0d0) * ((x * i) + (t * a))
else
tmp = (b * c) - (27.0d0 * (j * k))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((b * c) <= -6.6e+88) {
tmp = (b * c) + (-4.0 * (t * a));
} else if ((b * c) <= 4.7e+147) {
tmp = -4.0 * ((x * i) + (t * a));
} else {
tmp = (b * c) - (27.0 * (j * k));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if (b * c) <= -6.6e+88: tmp = (b * c) + (-4.0 * (t * a)) elif (b * c) <= 4.7e+147: tmp = -4.0 * ((x * i) + (t * a)) else: tmp = (b * c) - (27.0 * (j * k)) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (Float64(b * c) <= -6.6e+88) tmp = Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))); elseif (Float64(b * c) <= 4.7e+147) tmp = Float64(-4.0 * Float64(Float64(x * i) + Float64(t * a))); else tmp = Float64(Float64(b * c) - Float64(27.0 * Float64(j * k))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if ((b * c) <= -6.6e+88)
tmp = (b * c) + (-4.0 * (t * a));
elseif ((b * c) <= 4.7e+147)
tmp = -4.0 * ((x * i) + (t * a));
else
tmp = (b * c) - (27.0 * (j * k));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[N[(b * c), $MachinePrecision], -6.6e+88], N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 4.7e+147], N[(-4.0 * N[(N[(x * i), $MachinePrecision] + N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(b * c), $MachinePrecision] - N[(27.0 * N[(j * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -6.6 \cdot 10^{+88}:\\
\;\;\;\;b \cdot c + -4 \cdot \left(t \cdot a\right)\\
\mathbf{elif}\;b \cdot c \leq 4.7 \cdot 10^{+147}:\\
\;\;\;\;-4 \cdot \left(x \cdot i + t \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot c - 27 \cdot \left(j \cdot k\right)\\
\end{array}
\end{array}
if (*.f64 b c) < -6.6000000000000006e88Initial program 76.4%
distribute-rgt-out--81.6%
associate-*r*81.6%
*-commutative81.6%
associate-*l*81.6%
associate-*r*81.6%
Applied egg-rr81.6%
Taylor expanded in y around 0 79.6%
Taylor expanded in j around 0 81.7%
Taylor expanded in i around 0 71.9%
if -6.6000000000000006e88 < (*.f64 b c) < 4.7000000000000003e147Initial program 87.8%
distribute-rgt-out--89.8%
associate-*r*89.3%
*-commutative89.3%
associate-*l*89.3%
associate-*r*89.3%
Applied egg-rr89.3%
Taylor expanded in y around 0 78.1%
Taylor expanded in j around 0 57.1%
Taylor expanded in b around 0 53.9%
cancel-sign-sub-inv53.9%
metadata-eval53.9%
distribute-lft-out53.9%
*-commutative53.9%
*-commutative53.9%
Simplified53.9%
if 4.7000000000000003e147 < (*.f64 b c) Initial program 80.5%
Simplified80.5%
Taylor expanded in t around 0 72.5%
Taylor expanded in i around 0 67.9%
Final simplification60.2%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= (* b c) -5.5e+87)
(+ (* b c) (* -4.0 (* t a)))
(if (<= (* b c) 4.3e+144)
(* -4.0 (+ (* x i) (* t a)))
(+ (* b c) (* j (* k -27.0))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((b * c) <= -5.5e+87) {
tmp = (b * c) + (-4.0 * (t * a));
} else if ((b * c) <= 4.3e+144) {
tmp = -4.0 * ((x * i) + (t * a));
} else {
tmp = (b * c) + (j * (k * -27.0));
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if ((b * c) <= (-5.5d+87)) then
tmp = (b * c) + ((-4.0d0) * (t * a))
else if ((b * c) <= 4.3d+144) then
tmp = (-4.0d0) * ((x * i) + (t * a))
else
tmp = (b * c) + (j * (k * (-27.0d0)))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((b * c) <= -5.5e+87) {
tmp = (b * c) + (-4.0 * (t * a));
} else if ((b * c) <= 4.3e+144) {
tmp = -4.0 * ((x * i) + (t * a));
} else {
tmp = (b * c) + (j * (k * -27.0));
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if (b * c) <= -5.5e+87: tmp = (b * c) + (-4.0 * (t * a)) elif (b * c) <= 4.3e+144: tmp = -4.0 * ((x * i) + (t * a)) else: tmp = (b * c) + (j * (k * -27.0)) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (Float64(b * c) <= -5.5e+87) tmp = Float64(Float64(b * c) + Float64(-4.0 * Float64(t * a))); elseif (Float64(b * c) <= 4.3e+144) tmp = Float64(-4.0 * Float64(Float64(x * i) + Float64(t * a))); else tmp = Float64(Float64(b * c) + Float64(j * Float64(k * -27.0))); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if ((b * c) <= -5.5e+87)
tmp = (b * c) + (-4.0 * (t * a));
elseif ((b * c) <= 4.3e+144)
tmp = -4.0 * ((x * i) + (t * a));
else
tmp = (b * c) + (j * (k * -27.0));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[N[(b * c), $MachinePrecision], -5.5e+87], N[(N[(b * c), $MachinePrecision] + N[(-4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 4.3e+144], N[(-4.0 * N[(N[(x * i), $MachinePrecision] + N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(b * c), $MachinePrecision] + N[(j * N[(k * -27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -5.5 \cdot 10^{+87}:\\
\;\;\;\;b \cdot c + -4 \cdot \left(t \cdot a\right)\\
\mathbf{elif}\;b \cdot c \leq 4.3 \cdot 10^{+144}:\\
\;\;\;\;-4 \cdot \left(x \cdot i + t \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot c + j \cdot \left(k \cdot -27\right)\\
\end{array}
\end{array}
if (*.f64 b c) < -5.50000000000000022e87Initial program 76.4%
distribute-rgt-out--81.6%
associate-*r*81.6%
*-commutative81.6%
associate-*l*81.6%
associate-*r*81.6%
Applied egg-rr81.6%
Taylor expanded in y around 0 79.6%
Taylor expanded in j around 0 81.7%
Taylor expanded in i around 0 71.9%
if -5.50000000000000022e87 < (*.f64 b c) < 4.29999999999999984e144Initial program 87.8%
distribute-rgt-out--89.8%
associate-*r*89.3%
*-commutative89.3%
associate-*l*89.3%
associate-*r*89.3%
Applied egg-rr89.3%
Taylor expanded in y around 0 78.1%
Taylor expanded in j around 0 57.1%
Taylor expanded in b around 0 53.9%
cancel-sign-sub-inv53.9%
metadata-eval53.9%
distribute-lft-out53.9%
*-commutative53.9%
*-commutative53.9%
Simplified53.9%
if 4.29999999999999984e144 < (*.f64 b c) Initial program 80.5%
Simplified82.9%
Taylor expanded in b around inf 67.9%
Final simplification60.2%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. (FPCore (x y z t a b c i j k) :precision binary64 (if (or (<= (* b c) -2.3e+70) (not (<= (* b c) 2.5e+152))) (* b c) (* -4.0 (* x i))))
assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -2.3e+70) || !((b * c) <= 2.5e+152)) {
tmp = b * c;
} else {
tmp = -4.0 * (x * i);
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if (((b * c) <= (-2.3d+70)) .or. (.not. ((b * c) <= 2.5d+152))) then
tmp = b * c
else
tmp = (-4.0d0) * (x * i)
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -2.3e+70) || !((b * c) <= 2.5e+152)) {
tmp = b * c;
} else {
tmp = -4.0 * (x * i);
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if ((b * c) <= -2.3e+70) or not ((b * c) <= 2.5e+152): tmp = b * c else: tmp = -4.0 * (x * i) return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if ((Float64(b * c) <= -2.3e+70) || !(Float64(b * c) <= 2.5e+152)) tmp = Float64(b * c); else tmp = Float64(-4.0 * Float64(x * i)); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if (((b * c) <= -2.3e+70) || ~(((b * c) <= 2.5e+152)))
tmp = b * c;
else
tmp = -4.0 * (x * i);
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[Or[LessEqual[N[(b * c), $MachinePrecision], -2.3e+70], N[Not[LessEqual[N[(b * c), $MachinePrecision], 2.5e+152]], $MachinePrecision]], N[(b * c), $MachinePrecision], N[(-4.0 * N[(x * i), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -2.3 \cdot 10^{+70} \lor \neg \left(b \cdot c \leq 2.5 \cdot 10^{+152}\right):\\
\;\;\;\;b \cdot c\\
\mathbf{else}:\\
\;\;\;\;-4 \cdot \left(x \cdot i\right)\\
\end{array}
\end{array}
if (*.f64 b c) < -2.29999999999999994e70 or 2.5e152 < (*.f64 b c) Initial program 78.3%
distribute-rgt-out--81.3%
associate-*r*80.3%
*-commutative80.3%
associate-*l*80.3%
associate-*r*80.3%
Applied egg-rr80.3%
Taylor expanded in b around inf 56.4%
if -2.29999999999999994e70 < (*.f64 b c) < 2.5e152Initial program 87.8%
distribute-rgt-out--89.7%
associate-*r*89.2%
*-commutative89.2%
associate-*l*89.2%
associate-*r*89.2%
Applied egg-rr89.2%
Taylor expanded in i around inf 33.2%
*-commutative33.2%
Simplified33.2%
Final simplification42.2%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. (FPCore (x y z t a b c i j k) :precision binary64 (if (or (<= (* b c) -3.2e+102) (not (<= (* b c) 1.12e+21))) (* b c) (* (* j k) -27.0)))
assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -3.2e+102) || !((b * c) <= 1.12e+21)) {
tmp = b * c;
} else {
tmp = (j * k) * -27.0;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if (((b * c) <= (-3.2d+102)) .or. (.not. ((b * c) <= 1.12d+21))) then
tmp = b * c
else
tmp = (j * k) * (-27.0d0)
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -3.2e+102) || !((b * c) <= 1.12e+21)) {
tmp = b * c;
} else {
tmp = (j * k) * -27.0;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if ((b * c) <= -3.2e+102) or not ((b * c) <= 1.12e+21): tmp = b * c else: tmp = (j * k) * -27.0 return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if ((Float64(b * c) <= -3.2e+102) || !(Float64(b * c) <= 1.12e+21)) tmp = Float64(b * c); else tmp = Float64(Float64(j * k) * -27.0); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if (((b * c) <= -3.2e+102) || ~(((b * c) <= 1.12e+21)))
tmp = b * c;
else
tmp = (j * k) * -27.0;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[Or[LessEqual[N[(b * c), $MachinePrecision], -3.2e+102], N[Not[LessEqual[N[(b * c), $MachinePrecision], 1.12e+21]], $MachinePrecision]], N[(b * c), $MachinePrecision], N[(N[(j * k), $MachinePrecision] * -27.0), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -3.2 \cdot 10^{+102} \lor \neg \left(b \cdot c \leq 1.12 \cdot 10^{+21}\right):\\
\;\;\;\;b \cdot c\\
\mathbf{else}:\\
\;\;\;\;\left(j \cdot k\right) \cdot -27\\
\end{array}
\end{array}
if (*.f64 b c) < -3.1999999999999999e102 or 1.12e21 < (*.f64 b c) Initial program 80.1%
distribute-rgt-out--83.6%
associate-*r*83.6%
*-commutative83.6%
associate-*l*83.6%
associate-*r*83.6%
Applied egg-rr83.6%
Taylor expanded in b around inf 52.1%
if -3.1999999999999999e102 < (*.f64 b c) < 1.12e21Initial program 87.3%
Simplified88.2%
Taylor expanded in j around inf 29.2%
Final simplification39.4%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. (FPCore (x y z t a b c i j k) :precision binary64 (* b c))
assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return b * c;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = b * c
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return b * c;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): return b * c
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) return Float64(b * c) end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp = code(x, y, z, t, a, b, c, i, j, k)
tmp = b * c;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(b * c), $MachinePrecision]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
b \cdot c
\end{array}
Initial program 84.1%
distribute-rgt-out--86.4%
associate-*r*85.8%
*-commutative85.8%
associate-*l*85.8%
associate-*r*85.8%
Applied egg-rr85.8%
Taylor expanded in b around inf 25.2%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (+ (* a t) (* i x)) 4.0))
(t_2
(-
(- (* (* 18.0 t) (* (* x y) z)) t_1)
(- (* (* k j) 27.0) (* c b)))))
(if (< t -1.6210815397541398e-69)
t_2
(if (< t 165.68027943805222)
(+ (- (* (* 18.0 y) (* x (* z t))) t_1) (- (* c b) (* 27.0 (* k j))))
t_2))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = ((a * t) + (i * x)) * 4.0;
double t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b));
double tmp;
if (t < -1.6210815397541398e-69) {
tmp = t_2;
} else if (t < 165.68027943805222) {
tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j)));
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = ((a * t) + (i * x)) * 4.0d0
t_2 = (((18.0d0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0d0) - (c * b))
if (t < (-1.6210815397541398d-69)) then
tmp = t_2
else if (t < 165.68027943805222d0) then
tmp = (((18.0d0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0d0 * (k * j)))
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 c, double i, double j, double k) {
double t_1 = ((a * t) + (i * x)) * 4.0;
double t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b));
double tmp;
if (t < -1.6210815397541398e-69) {
tmp = t_2;
} else if (t < 165.68027943805222) {
tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j)));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j, k): t_1 = ((a * t) + (i * x)) * 4.0 t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b)) tmp = 0 if t < -1.6210815397541398e-69: tmp = t_2 elif t < 165.68027943805222: tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j))) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(Float64(a * t) + Float64(i * x)) * 4.0) t_2 = Float64(Float64(Float64(Float64(18.0 * t) * Float64(Float64(x * y) * z)) - t_1) - Float64(Float64(Float64(k * j) * 27.0) - Float64(c * b))) tmp = 0.0 if (t < -1.6210815397541398e-69) tmp = t_2; elseif (t < 165.68027943805222) tmp = Float64(Float64(Float64(Float64(18.0 * y) * Float64(x * Float64(z * t))) - t_1) + Float64(Float64(c * b) - Float64(27.0 * Float64(k * j)))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k) t_1 = ((a * t) + (i * x)) * 4.0; t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b)); tmp = 0.0; if (t < -1.6210815397541398e-69) tmp = t_2; elseif (t < 165.68027943805222) tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j))); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(N[(a * t), $MachinePrecision] + N[(i * x), $MachinePrecision]), $MachinePrecision] * 4.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(N[(18.0 * t), $MachinePrecision] * N[(N[(x * y), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision] - N[(N[(N[(k * j), $MachinePrecision] * 27.0), $MachinePrecision] - N[(c * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[t, -1.6210815397541398e-69], t$95$2, If[Less[t, 165.68027943805222], N[(N[(N[(N[(18.0 * y), $MachinePrecision] * N[(x * N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision] + N[(N[(c * b), $MachinePrecision] - N[(27.0 * N[(k * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(a \cdot t + i \cdot x\right) \cdot 4\\
t_2 := \left(\left(18 \cdot t\right) \cdot \left(\left(x \cdot y\right) \cdot z\right) - t\_1\right) - \left(\left(k \cdot j\right) \cdot 27 - c \cdot b\right)\\
\mathbf{if}\;t < -1.6210815397541398 \cdot 10^{-69}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t < 165.68027943805222:\\
\;\;\;\;\left(\left(18 \cdot y\right) \cdot \left(x \cdot \left(z \cdot t\right)\right) - t\_1\right) + \left(c \cdot b - 27 \cdot \left(k \cdot j\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
herbie shell --seed 2024144
(FPCore (x y z t a b c i j k)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, E"
:precision binary64
:alt
(! :herbie-platform default (if (< t -8105407698770699/5000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (- (- (* (* 18 t) (* (* x y) z)) (* (+ (* a t) (* i x)) 4)) (- (* (* k j) 27) (* c b))) (if (< t 8284013971902611/50000000000000) (+ (- (* (* 18 y) (* x (* z t))) (* (+ (* a t) (* i x)) 4)) (- (* c b) (* 27 (* k j)))) (- (- (* (* 18 t) (* (* x y) z)) (* (+ (* a t) (* i x)) 4)) (- (* (* k j) 27) (* c b))))))
(- (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) (* (* j 27.0) k)))