const unsigned int n_levels = triangulation.n_levels();
mg_interface_in.resize(0, n_levels-1);
- mg_interface_in.clear ();
+ mg_interface_in.clear_elements ();
mg_interface_out.resize(0, n_levels-1);
- mg_interface_out.clear ();
+ mg_interface_out.clear_elements ();
mg_matrices.resize(0, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
mg_sparsity_patterns.resize(0, n_levels-1);
// Now, we have to provide a matrix on each level. To this end, we first
setup_time = 0;
system_matrix.clear();
- mg_matrices.clear();
- mg_constraints.clear();
+ mg_matrices.clear_elements();
+ mg_constraints.clear_elements();
dof_handler.distribute_dofs (fe);
dof_handler.distribute_mg_dofs (fe);
// The global system is set up, now we attend to the level matrices. We
// resize all matrix objects to hold one matrix per level.
mg_matrix.resize(0, n_levels-1);
- mg_matrix.clear();
+ mg_matrix.clear_elements();
mg_matrix_dg_up.resize(0, n_levels-1);
- mg_matrix_dg_up.clear();
+ mg_matrix_dg_up.clear_elements();
mg_matrix_dg_down.resize(0, n_levels-1);
- mg_matrix_dg_down.clear();
+ mg_matrix_dg_down.clear_elements();
// It is important to update the sparsity patterns after <tt>clear()</tt>
// was called for the level matrices, since the matrices lock the sparsity
// pattern through the SmartPointer and Subscriptor mechanism.
const unsigned int n_levels = triangulation.n_global_levels();
mg_interface_matrices.resize(0, n_levels-1);
- mg_interface_matrices.clear ();
+ mg_interface_matrices.clear_elements ();
mg_matrices.resize(0, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
// Now, we have to provide a matrix on each
// level. To this end, we first use the
inline void
SmootherRelaxation<RelaxationType, VectorType>::clear ()
{
- MGLevelObject<RelaxationType>::clear();
+ MGLevelObject<RelaxationType>::clear_elements();
}
inline void
MGSmootherRelaxation<MatrixType, RelaxationType, VectorType>::clear ()
{
- smoothers.clear();
+ smoothers.clear_elements();
unsigned int i=matrices.min_level(),
max_level=matrices.max_level();
inline void
MGSmootherPrecondition<MatrixType, PreconditionerType, VectorType>::clear ()
{
- smoothers.clear();
+ smoothers.clear_elements();
unsigned int i=matrices.min_level(),
max_level=matrices.max_level();
SparseMatrixCollection<number>::resize(const unsigned int minlevel, const unsigned int maxlevel)
{
matrix.resize(minlevel, maxlevel);
- matrix.clear();
+ matrix.clear_elements();
matrix_up.resize(minlevel+1, maxlevel);
- matrix_up.clear();
+ matrix_up.clear_elements();
matrix_down.resize(minlevel+1, maxlevel);
- matrix_down.clear();
+ matrix_down.clear_elements();
matrix_in.resize(minlevel, maxlevel);
- matrix_in.clear();
+ matrix_in.clear_elements();
matrix_out.resize(minlevel, maxlevel);
- matrix_out.clear();
+ matrix_out.clear_elements();
sparsity.resize(minlevel, maxlevel);
sparsity_edge.resize(minlevel, maxlevel);
}
void LaplaceProblem<dim>::setup_system ()
{
system_matrix.clear();
- mg_matrices.clear();
- mg_constraints.clear();
+ mg_matrices.clear_elements();
+ mg_constraints.clear_elements();
dof_handler.distribute_dofs (fe);
dof_handler.distribute_mg_dofs (fe);
const unsigned int n_levels = triangulation.n_global_levels();
mg_matrix.resize(0, n_levels-1);
- mg_matrix.clear();
+ mg_matrix.clear_elements();
mg_matrix_dg_up.resize(0, n_levels-1);
- mg_matrix_dg_up.clear();
+ mg_matrix_dg_up.clear_elements();
mg_matrix_dg_down.resize(0, n_levels-1);
- mg_matrix_dg_down.clear();
+ mg_matrix_dg_down.clear_elements();
for (unsigned int level=mg_matrix.min_level();
level<=mg_matrix.max_level(); ++level)
const unsigned int n_levels = triangulation.n_global_levels();
mg_interface_matrices.resize(0, n_levels-1);
- mg_interface_matrices.clear ();
+ mg_interface_matrices.clear_elements ();
mg_matrices.resize(0, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
for (unsigned int level=0; level<n_levels; ++level)
{
const unsigned int n_levels = triangulation.n_global_levels();
mg_interface_matrices.resize(0, n_levels-1);
- mg_interface_matrices.clear ();
+ mg_interface_matrices.clear_elements ();
mg_matrices.resize(0, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
for (unsigned int level=0; level<n_levels; ++level)
{
const unsigned int n_levels = triangulation.n_global_levels();
mg_matrix.resize(0, n_levels-1);
- mg_matrix.clear();
+ mg_matrix.clear_elements();
mg_matrix_dg_up.resize(0, n_levels-1);
- mg_matrix_dg_up.clear();
+ mg_matrix_dg_up.clear_elements();
mg_matrix_dg_down.resize(0, n_levels-1);
- mg_matrix_dg_down.clear();
+ mg_matrix_dg_down.clear_elements();
for (unsigned int level=mg_matrix.min_level();
level<=mg_matrix.max_level(); ++level)
}
mg_matrices.resize(0, nlevels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
mg_matrices_renumbered.resize(0, nlevels-1);
- mg_matrices_renumbered.clear ();
+ mg_matrices_renumbered.clear_elements ();
mg_sparsity.resize(0, nlevels-1);
mg_sparsity_renumbered.resize(0, nlevels-1);
const unsigned int n_levels = triangulation.n_levels();
mg_interface_in.resize(0, n_levels-1);
- mg_interface_in.clear ();
+ mg_interface_in.clear_elements ();
mg_interface_out.resize(0, n_levels-1);
- mg_interface_out.clear ();
+ mg_interface_out.clear_elements ();
mg_matrices.resize(0, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
mg_sparsity_patterns.resize(0, n_levels-1);
for (unsigned int level=0; level<n_levels; ++level)
const unsigned int n_levels = triangulation.n_levels();
mg_interface_matrices.resize(min_level, n_levels-1);
- mg_interface_matrices.clear ();
+ mg_interface_matrices.clear_elements ();
mg_matrices.resize(min_level, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
mg_sparsity_patterns.resize(min_level, n_levels-1);
for (unsigned int level=min_level; level<n_levels; ++level)
const unsigned int n_levels = triangulation.n_levels();
mg_interface_matrices.resize(min_level, n_levels-1);
- mg_interface_matrices.clear ();
+ mg_interface_matrices.clear_elements ();
mg_matrices.resize(min_level, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
mg_sparsity_patterns.resize(min_level, n_levels-1);
for (unsigned int level=min_level; level<n_levels; ++level)
const unsigned int n_levels = triangulation.n_global_levels();
mg_interface_matrices.resize(0, n_levels-1);
- mg_interface_matrices.clear ();
+ mg_interface_matrices.clear_elements ();
mg_matrices.resize(0, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
for (unsigned int level=0; level<n_levels; ++level)
{
const unsigned int n_levels = triangulation.n_levels();
mg_interface_matrices.resize(0, n_levels-1);
- mg_interface_matrices.clear ();
+ mg_interface_matrices.clear_elements ();
mg_matrices.resize(0, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
mg_sparsity_patterns.resize(0, n_levels-1);
for (unsigned int level=0; level<n_levels; ++level)
const unsigned int n_levels = triangulation.n_levels();
mg_interface_in.resize(0, n_levels-1);
- mg_interface_in.clear ();
+ mg_interface_in.clear_elements ();
mg_interface_out.resize(0, n_levels-1);
- mg_interface_out.clear ();
+ mg_interface_out.clear_elements ();
mg_matrices.resize(0, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
mg_sparsity_patterns.resize(0, n_levels-1);
for (unsigned int level=0; level<n_levels; ++level)
const unsigned int n_levels = triangulation.n_levels();
mg_interface_matrices.resize(0, n_levels-1);
- mg_interface_matrices.clear ();
+ mg_interface_matrices.clear_elements ();
mg_matrices.resize(0, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
mg_sparsity_patterns.resize(0, n_levels-1);
for (unsigned int level=0; level<n_levels; ++level)
const unsigned int n_levels = triangulation.n_levels();
mg_matrix.resize(0, n_levels-1);
- mg_matrix.clear();
+ mg_matrix.clear_elements();
mg_matrix_dg_up.resize(0, n_levels-1);
- mg_matrix_dg_up.clear();
+ mg_matrix_dg_up.clear_elements();
mg_matrix_dg_down.resize(0, n_levels-1);
- mg_matrix_dg_down.clear();
+ mg_matrix_dg_down.clear_elements();
mg_matrix_in_out.resize(0, n_levels-1);
- mg_matrix_in_out.clear();
+ mg_matrix_in_out.clear_elements();
mg_sparsity.resize(0, n_levels-1);
mg_sparsity_dg_interface.resize(0, n_levels-1);
const unsigned int n_levels = triangulation.n_levels();
mg_matrix.resize(0, n_levels-1);
- mg_matrix.clear();
+ mg_matrix.clear_elements();
mg_matrix_dg_up.resize(0, n_levels-1);
- mg_matrix_dg_up.clear();
+ mg_matrix_dg_up.clear_elements();
mg_matrix_dg_down.resize(0, n_levels-1);
- mg_matrix_dg_down.clear();
+ mg_matrix_dg_down.clear_elements();
mg_matrix_in_out.resize(0, n_levels-1);
- mg_matrix_in_out.clear();
+ mg_matrix_in_out.clear_elements();
mg_sparsity.resize(0, n_levels-1);
mg_sparsity_dg_interface.resize(0, n_levels-1);
const unsigned int n_levels = triangulation.n_levels();
mg_matrix.resize(0, n_levels-1);
- mg_matrix.clear();
+ mg_matrix.clear_elements();
mg_matrix_dg_up.resize(0, n_levels-1);
- mg_matrix_dg_up.clear();
+ mg_matrix_dg_up.clear_elements();
mg_matrix_dg_down.resize(0, n_levels-1);
- mg_matrix_dg_down.clear();
+ mg_matrix_dg_down.clear_elements();
mg_matrix_in_out.resize(0, n_levels-1);
- mg_matrix_in_out.clear();
+ mg_matrix_in_out.clear_elements();
mg_sparsity.resize(0, n_levels-1);
mg_sparsity_dg_interface.resize(0, n_levels-1);
const unsigned int n_levels = triangulation.n_levels();
mg_matrix.resize(0, n_levels-1);
- mg_matrix.clear();
+ mg_matrix.clear_elements();
mg_matrix_dg_up.resize(0, n_levels-1);
- mg_matrix_dg_up.clear();
+ mg_matrix_dg_up.clear_elements();
mg_matrix_dg_down.resize(0, n_levels-1);
- mg_matrix_dg_down.clear();
+ mg_matrix_dg_down.clear_elements();
mg_sparsity.resize(0, n_levels-1);
mg_sparsity_dg_interface.resize(0, n_levels-1);
const unsigned int n_levels = triangulation.n_global_levels();
mg_interface_matrices.resize(0, n_levels-1);
- mg_interface_matrices.clear ();
+ mg_interface_matrices.clear_elements ();
mg_matrices.resize(0, n_levels-1);
- mg_matrices.clear ();
+ mg_matrices.clear_elements ();
for (unsigned int level=0; level<n_levels; ++level)
{