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.
For DG, it probably wouldn't be very difficult to implement something like
this, but it doesn't exist at the moment.
Can I ask why you want this particular property? Why is it not enough to have
all DoFs associated with a cell contiguous in memory? And why do you care
about the cell-local order of DoFs, and not the global DoF indices?
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
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/53eb39b8-c2c4-ca4d-74f3-3481b55f19c8%40colostate.edu.
