Hi Greg, I guess what you could do is to treat the vector as vector associated with a DoFHandler with FE_DGQ(0). During each refinement you "interpolate" the vector to the new mesh and decrease the value by one. This approach naturally extends to a parallel setting.
Alternately, you could create a single vector associated the coarse cells and while looping over active cells of the refined mesh you would recursively call the parent cell until you reach the coarsest cell whose id you can use to access the vector. In the parallel setting, the vector associated to the coarse cells would need to replicated between processes. Hope this help, PM On Thursday, 9 November 2023 at 03:32:42 UTC+1 [email protected] wrote: > Dear all, > > This might be a somewhat odd feature request but I was wondering whether > the functionality described below can be readily achieved in deal.II > already: > > Say we are given a coarse mesh with every cell never refined, and a vector > of non-negative integers associated with each cell. The integer is the > level of isotropic refinement to be applied to each cell. The number can be > 0, which means not doing anything here; it can be 1, which corresponds to > what set_refine_flag() does. But it can also be any number bigger than 1. > Say we have a quadrilateral cell and the number associated with it is 2, > then the cell is to be refined isotropically once into four children and > then every child is refined again, resulting in 16 grandchildren. > > If we were to hard code this, would it be sound to simply rely on > set_refine_flag() and execute the refinement one level after another for > M levels, with M being the maximum number of the given vector? And perhaps > to do this we would need iterators triangulation.begin(n), > triangulation.end(n) and isotropic_child() involved somehow? > > Any help and/or pro tips would be greatly appreciated! > > Best regards, > Yuan > > -- 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/78e80d9e-8b6f-4a99-beee-5f2df247383dn%40googlegroups.com.
