Good evening everybody, I could really use some generous advice. Distributed code, PETSc wrappers, adaptive refinements. We get to the point of refining:
Vector<float> error_indicators(this->triangulation->n_active_cells()); ... some strategy to compute error_indicators ... parallel::distributed::GridRefinement::refine_and_coarsen_fixed_fraction( *this->triangulation, error_indicators, 0.8, 0.02); cout<<"Rank "<<this->this_mpi_process<<" starts refining and coarsening" <<std::endl; until here all fine this->triangulation->execute_coarsening_and_refinement(); But here it breaks... sometimes. Let me explain my bewilderment: It may go well for 6 refinement cycles and break at the 7th. Does anybody have any hint at what I should investigate to figure out the problem? Thanks in advance! Best, Matteo Malvestiti -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion visit https://groups.google.com/d/msgid/dealii/c26cb4d1-c804-4741-9e0f-33a5859873b6n%40googlegroups.com.
