On Thu, 10 Nov 2016, Michael Povolotskyi wrote: > Dear Libmesh developers and users, > I need to solve the following problem > 1) Create a mesh > 2) refine part of it for computing one 3D integral and deactivate > elements beyond that part > 3) compute integral on the active elements > 4) revert the mesh to the non-refined state > 5) refine a different part of the mesh, deactivate elements beyond that part > 6) compute another 3D integral. > > My question is how to do the steps 2 and 4. > The code know which elements needs to be refined and deactivated based > on the element center coordinates.
After step 1, create your system with per-subdomain variables, only active on one subdomain id. For step 2 and 4, loop over all elements, setting refinement flags for the elements you want to refine and setting subdomain ids for elements you want active vs inactive. You can use coarsen_uniformly() in step 4 rather than manually setting coarsen flags. --- Roy ------------------------------------------------------------------------------ Developer Access Program for Intel Xeon Phi Processors Access to Intel Xeon Phi processor-based developer platforms. With one year of Intel Parallel Studio XE. Training and support from Colfax. Order your platform today. http://sdm.link/xeonphi _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
