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); } }