[Libmesh-users] Adaptative Mesh Refinement

[Vinicius C. R. <vin...@co...>]

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