> I have a vector fe like fe(FE_DGQ<dim>(degree), n_components).
>
> On each cell I want to separate the solution dofs into support points like
>
> W[support_point_index][component_index]
>
> I then want to use W to compute a flux at every support point. This flux
> depends on the n_components at the support point.
>
> Note that I dont want to compute the support point coordinates.
If you're not really interested in the actual location of the DoF but only its
number within a single component, then you could write this as follows:
W[fe.system_to_component_index(i).second]
[fe.system_to_component_index(i).first]
See FiniteElement::system_to_component_index.
Best
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
