typedef LaplaceOperator<dim,fe_degree,n_q_points_1d,number> LevelMatrixType;
MGLevelObject<LevelMatrixType> mg_matrices;
- mg_matrices.resize(0, dof.get_triangulation().n_levels()-1);
- for (unsigned int level = 0; level<dof.get_triangulation().n_levels(); ++level)
+ mg_matrices.resize(0, dof.get_triangulation().n_global_levels()-1);
+ for (unsigned int level = 0; level<dof.get_triangulation().n_global_levels(); ++level)
{
mg_matrices[level].initialize(mapping, dof, dirichlet_boundaries, level);
}
mg_smoother;
MGLevelObject<typename SMOOTHER::AdditionalData> smoother_data;
- smoother_data.resize(0, dof.get_triangulation().n_levels()-1);
- for (unsigned int level = 0; level<dof.get_triangulation().n_levels(); ++level)
+ smoother_data.resize(0, dof.get_triangulation().n_global_levels()-1);
+ for (unsigned int level = 0; level<dof.get_triangulation().n_global_levels(); ++level)
{
smoother_data[level].smoothing_range = 15.;
smoother_data[level].degree = 5;