make_matrix()

void F4GB::make_matrix ( )

private

Definition at line 469 of file f4.cpp.

{
  /* 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;
  std::pair<bool,spair> p;
  p = mSPairSet.get_next_pair();
  while (p.first)
    {
      process_s_pair(&p.second);
      p = mSPairSet.get_next_pair();
    }
 
  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());
    }
  //  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...?
 
  // For debugging: let's compute the number of non-zero elements above the diagonal in A matrix.
  // Also, let's compute the amount of space used for:
  //  column info
  //  row info
  //  sum of values in each row
  // Maybe so A B C D matrices separately.
  macaulayMatrixStats();
}

References M2_gbTrace, macaulayMatrixStats(), mat, mSPairSet, next_col_to_process, p, process_column(), process_s_pair(), and reorder_columns().

Referenced by do_spairs().