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 ________________________________ From: dealii@googlegroups.com <dealii@googlegroups.com> on behalf of Praveen C <cprav...@gmail.com> Sent: Tuesday, June 1, 2021 5:34 AM To: Deal. II Googlegroup <dealii@googlegroups.com> Subject: [deal.II] order of dofs within a cell in a FESystem Is it possible to change the order of dofs within a cell in a FESystem ? I have situation like this: FESystem<dim> fe; degree = 1 nvar = 4 // number of unknowns at each support point fe(FE_DGQArbitrarynodes(QGauss<1>(degree+1), nvar) I print for i = 0; i<fe.dofs_per_cell; ++i comp_i = fe.system_to_component_index(i).first indx_i = fe.system_to_component_index(i).second i comp_i indx_i 0 0 0 1 0 1 2 0 2 3 0 3 4 1 0 5 1 1 6 1 2 7 1 3 8 2 0 9 2 1 10 2 2 11 2 3 12 3 0 13 3 1 14 3 2 15 3 3 All 0'th component are consecutive 1'th component are consecutive etc. Is it possible to have the other ordering i comp_i indx_i 0 0 0 1 1 0 2 2 0 3 3 0 4 0 1 5 1 1 6 2 1 7 3 1 8 0 2 9 1 2 10 2 2 11 3 2 12 0 3 13 1 3 14 2 3 15 3 3 Thanks praveen -- 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CC6DB04A-6A8C-4CEE-B895-1B36D151EAFB%40gmail.com. -- 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/BN7PR03MB435685C9E893AF40F59FC257ED3E9%40BN7PR03MB4356.namprd03.prod.outlook.com.
