Hi Guido, > in the latest version, how can I get the number of subfaces on a given > face? I would have expected that TriaObjectAccessor or CellAccessor has > a function telling me this, but I cannot find it.
For a general face, what exactly is the number of subfaces on a given face? Before we can implement a function which gives us a number we have to answer that question. I cannot. However, the most relevant cases are explained in the step-30 tutorial, but I can try to summarize them here: a) The maximum number of direct children a face can have is given by GeometryInfo<dim>::max_children_per_face. This is the number of children of an isotropically refined face. b) The actual number of direct offsprings of a face can be obtained via face->n_children(). c) If you are on an active cell and and consider one of the faces of that cell, then face->number_of_children() gives the total number of active offsprings of that face, which (for a refined face) is always two in 2d, but could be two, three or four in 3d. This is the number of subfaces to integrate over in DG methods, for example. Best, Tobias _______________________________________________
