More easily, you can use FETools::compute_block_renumbering(...)
Guido
actually I want to access to the global index of an individual dof of an individual type in an individual block in vector value problems; namely I want to know for example what is the global index of the 4th v-velocity in the 34th cell of the domain; in another note I am looking for the global index of cells[34]->block(1).member[4]; he I assume bilinear FE for pressure and biquadratic FE for both u and v; hence each cell would have a matrix of three block within each block we have finite number of dofs depnding upon the type of FE we've employed.You can do this as follows: cell = dof_handler.begin_active(); for (i=0; i<34; ++i) ++cell std::vector<uint> local_dof_indices (fe.dofs_per_cell); cell->get_dof_indices(local_dof_indices); and the global dof index you want is then local_dof_indices[fe.component_to_system_index(1,4)]; Best W.
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
