99 int row,
col, pp, length, ok ;
105 printf (
"ccolamd %d-by-%d input matrix:\n",
A_NROW,
A_NCOL) ;
106 for (col = 0 ; col <
A_NCOL ; col++)
108 length = p [col+1] - p [
col] ;
109 printf (
"Column %d, with %d entries:\n", col, length) ;
110 for (pp = p [col] ; pp < p [col+1] ; pp++)
113 printf (
" row %d\n", row) ;
126 printf (
"ccolamd error!\n") ;
134 printf (
"ccolamd column ordering:\n") ;
135 printf (
"1st column: %d\n", p [0]) ;
136 printf (
"2nd column: %d\n", p [1]) ;
137 printf (
"3rd column: %d\n", p [2]) ;
138 printf (
"4th column: %d\n", p [3]) ;
144 printf (
"\n\ncsymamd %d-by-%d input matrix:\n",
B_N,
B_N) ;
145 printf (
"Entries in strictly lower triangular part:\n") ;
146 for (col = 0 ; col <
B_N ; col++)
148 length = q [col+1] - q [
col] ;
149 printf (
"Column %d, with %d entries:\n", col, length) ;
150 for (pp = q [col] ; pp < q [col+1] ; pp++)
153 printf (
" row %d\n", row) ;
161 ok =
csymamd (B_N, B, q, perm, (
double *) NULL, stats, &calloc, &free,
167 printf (
"csymamd error!\n") ;
175 printf (
"csymamd column ordering:\n") ;
176 printf (
"1st row/column: %d\n", perm [0]) ;
177 printf (
"2nd row/column: %d\n", perm [1]) ;
178 printf (
"3rd row/column: %d\n", perm [2]) ;
179 printf (
"4th row/column: %d\n", perm [3]) ;
180 printf (
"5th row/column: %d\n", perm [4]) ;
void csymamd_report(int stats [CCOLAMD_STATS])
void ccolamd_report(int stats [CCOLAMD_STATS])
idx_t idx_t idx_t idx_t idx_t * perm
EIGEN_DEVICE_FUNC const Scalar & q
int csymamd(int n, int A [], int p [], int perm [], double knobs [CCOLAMD_KNOBS], int stats [CCOLAMD_STATS], void *(*allocate)(size_t, size_t), void(*release)(void *), int cmember [], int stype)
int ccolamd(int n_row, int n_col, int Alen, int A [], int p [], double knobs [CCOLAMD_KNOBS], int stats [CCOLAMD_STATS], int cmember [])