> given a cell in my DG finite element triangulation, how can I find the > entries in the solution vector that correspond to the DoFs defined on > the cell? Also, is there a way to construct the inverse mapping (from > solution entries to DoFs)?
You can ask the cell for the degrees of freedom defined on it like so: std::vector<unsigned int> local_dof_indices (fe.dofs_per_cell); cell->get_dof_indices (local_dof_indices); You should be able to construct the reverse mapping from this forward mapping. W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
