Mihai, > for (typename Triangulation<dim>::active_cell_iterator > cell = triangulation.begin_active(triangulation.n_levels()-2); > cell != triangulation.end(); ++cell) > cell->clear_coarsen_flag ();
Here is the problem: you start with the first active cell on level 2 (the third one) and continue to the last active cell (on the final level), thus you clear all coarsening flags, not only those on the desired level. In your particular case you should end the loop with cell != triangulation.begin_active(triangulation.n_levels()-1) Best regards, Tobias _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
