Hi Tobias, Thanks a lot - your fix worked!
-- Mihai ________________________________ Von: "[email protected]" <[email protected]> An: [email protected]; [email protected] Gesendet: Donnerstag, den 18. November 2010, 13:45:32 Uhr Betreff: AW: [deal.II] controlled coarsening and refinement 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
