I was using flag_elements_by_nelem_target for an adaptive computation, and the assert assert(parent->has_children()); near the bottom of MeshRefinement::flag_elements_by_nelem_target failed. It looks to me like the problem is that active elements don't get a flag of -1 from the function MeshRefinement::create_parent_error_vector?
I added a line error_per_parent[elem->id()] = -1.0; in the first element loop in MeshRefinement::create_parent_error_vector (around line 186), and it worked OK. Does this sound like a correct alteration? - Dave ------------------------------------------------------------------------- This SF.net email is sponsored by: Microsoft Defy all challenges. Microsoft(R) Visual Studio 2008. http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/ _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
