Thanks David for the clarification. I am using nodal DG basis.
I would like all components located at a support point to be consecutive in memory. I could achieve this by implementing a renumber_node_wise function. However, when I get local dofs on a cell cell->get_dof_indices(local_dof_indices) the indices pointed by local_dof_indices are not consecutive in memory, since this order is fixed by FESystem. I can make my own looping without using local_dof_indices which reads consecutively from a Vector but it would be nice if this was already supported by some FE class. A space like FE_DGTensor( fe_dg_1d, n_components) where n_components are together in memory, seems very useful. Best praveen > On 01-Jun-2021, at 8:42 PM, Wells, David <[email protected]> wrote: > > Hi Praveen, > > I don't think its possible to change the internal numbering used by FESystem > itself since it must support the case where not all finite elements have the > same number of degrees of freedom. If you really need to accomplish this you > will have to write a new finite element class. Another possibility (if all > the base elements are the same) would be to only have one finite element and > do the block-based numbering manually since each component will have the same > number of dofs. > > Best, > David -- 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/072850D0-76DA-4F82-ACC4-B38333CC4F3B%40gmail.com.
