void F4GB::make_matrix()
{
/* loop through all spairs, process,
then while there are any columns to process, do so,
then process rows.
Is this the best order to do it in? Maybe not...
*/
spair *p;
while (p = S->get_next_pair())
{
process_s_pair(p);
}
while (next_col_to_process < mat->columns.size())
process_column(next_col_to_process++);
// DEBUGGING:
if (M2_gbTrace >= 2) {
fprintf(stderr, "--matrix--%ld by %ld\n",
(long)mat->rows.size(), (long)mat->columns.size());
if (using_syz)
fprintf(stderr, "-syzygies-%ld by %ld\n",
(long)syz->rows.size(), (long)syz->columns.size());
}
// show_row_info();
// show_column_info();
// show_matrix();
// Now we reorder the columns, possibly rows?
reorder_columns();
//reorder_rows(); // This is only here so we can see what we are doing...?
}
void
GaussJordan::uptriangle ()
{
for (int k = 0; k < size; k++)
{
find_pivot (k);
process_column (k);
}
}
