> 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.
--
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
