begin, end,
dof_info, info_box,
integrator, assembler);
+
+ matrix.compress(VectorOperation::add);
}
MatrixIntegrator<dim> integrator;
MeshWorker::integration_loop<dim, dim> (
- begin, end,
+ begin, end,
dof_info, info_box,
integrator, assembler);
+
+ for (unsigned int level=mg_matrix.min_level(); level <= mg_matrix.max_level();++level)
+ {
+ mg_matrix[level].compress(VectorOperation::add);
+ if (level > mg_matrix.min_level())
+ {
+ mg_matrix_dg_up[level].compress(VectorOperation::add);
+ mg_matrix_dg_down[level].compress(VectorOperation::add);
+ }
+ }
}
dof_info, info_box,
integrator, assembler);
+ right_hand_side.compress(VectorOperation::add);
right_hand_side *= -1.;
}
In the beginning the Universe was created. This has made a lot of
people very angry and has been widely regarded as a bad move.
Douglas Adams