/**
* This function calls the compress()
* command of all matrices after
- * the sparsity pattern has been
- * generated or the assembly is
+ * the assembly is
* completed. Note that all MPI
* processes need to call this
- * command (in contrast to sparsity
- * pattern generation, which is
- * probably only called on one
- * single process) before any
+ * command (whereas the individual
+ * assembly routines will most probably
+ * only be called on each processor
+ * individually) before any
* can complete it.
*/
void compress ();
* subobjects. You *must* call
* this function each time after
* you have changed the size of
- * the sub-objects.
+ * the sub-objects. Note that
+ * this is a collective
+ * operation, i.e., it needs to
+ * be called on all MPI processes.
*/
void collect_sizes ();
matrix = std::auto_ptr<Epetra_FECrsMatrix>
(new Epetra_FECrsMatrix(Copy, row_map, col_map,
&n_entries_per_row[0], true));
-
-// Assert (matrix->NumGlobalCols() == col_map.NumGlobalElements(),
-// ExcInternalError());
-// Assert (matrix->NumGlobalRows() == row_map.NumGlobalElements(),
-// ExcInternalError());
-// Assert (matrix->NumGlobalCols() == sparsity_pattern.n_cols(),
-// ExcInternalError());
-// Assert (matrix->NumGlobalRows() == sparsity_pattern.n_rows(),
-// ExcInternalError());
const unsigned int n_max_entries_per_row = *std::max_element (
&n_entries_per_row[0], &n_entries_per_row[n_rows-1]);