Hi there, I actually figured out a problem with the function n_global_active_cells from parallel::distributed::Triangulation.
If I call this function after a simple global refinement, the numbercache will not be updated (maybe not). I figured out, that there is a function update_number_cache(), but it is private, so I can't do it directly. So the output of this function will be the number of locally owned cells only. A workaround is maybe to use the MPI::sum overall locally_owned cells. Maybe someone should fix this. (I think calling update_number_cache in the function n_global_active_cell should be the fix). Note: after using the function execute_coarsening_and_refinement() the number cache will be updated, and I think the n_global_active_cells shall give the expected result. Best Uwe -- Dipl.-Ing. Uwe Köcher Hamburg, Helmut Schmidt University University of the German Federal Armed Forces Email: [email protected] Web: http://www.hsu-hh.de/mb-mathe _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
