19 #include "monomial_ordering.hpp" 25 unsigned long **A =
new unsigned long * [ROWDIM];
26 for (
unsigned i = 0; i < ROWDIM; ++i)
27 A[i] =
new unsigned long [COLDIM];
32 A[0][0] = A[0][1] = 0; A[0][2] = A[0][3] = 1;
33 A[1][0] = A[1][1] = A[1][2] = 1; A[1][3] = 2;
34 A[2][0] = 1; A[2][1] = 2; A[2][2] = A[2][3] = 3;
35 A[3][0] = A[3][1] = A[3][2] = 3; A[3][3] = 6;
36 std::cout <<
"Nonsingular? " <<
nonsingular(ROWDIM, COLDIM, const_cast<const unsigned long **>(A)) << std::endl;
bool nonsingular(NVAR_TYPE m, NVAR_TYPE n, const WT_TYPE **A)
verifies that a matrix supplied for an ordering is not nonsingular