Hi, I'm trying to perform a h-refinement and a subdomain evaluation/reevaluation for each refinement iteration, prior to initializing an EquationSystems object. Currently when call init() I get a out_of_range exception. I suspect it is thrown from within its node or the element loop. Everything works fine if only I leave the subdomain evaluation prior to the refinements (i.e., I comment the subdomain reevaluation inside the refinement iterations).
The refinement is made by setting the refine flag directly in each selected element and calling refine_elements() in the MeshRefinement object. I'm not using any ErrorEstimator because there's no error to be estimated a the point the refinement is performed. Can someone point me what might be the wrong assumption(s) I'm making? Is it possible for a parent element to have childrens in different subdomains? If I iterate from *MeshBase->active_local_subdomain_elements_begin ([MySubDomain])* through *MeshBase->active_local_subdomain_elements_begin ([MySubDomain])* right after a refinement step will the iterators be correctly pointing to the new subdomain's range? I'm using libmesh compiled with Petsc and MPI support, nevertheless I'm running it in a desktop and with only one node (-np1). If it helps to contextualize, I'm working on an embedded domain problem. I'd appreciate any help! Thanks, VinÃcius da Costa Reis Master's Candidate Department of Civil Engineering COPPE/Federal University of Rio de Janeiro ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
