Hello all, I am still trying to find a good solution to this problem.
I now do the following: const dealii::FE_DGQ<dim> fe_cell; dealii::DoFHandler<dim> dh_cell; fe_cell is of degree zero, dh_cell is the corresponding dof handler; I have one dof per cell. I want to compute a shock indicator on each cell: dealii::Vector<double> shock_indicator; Using MeshWorker, when I am in the face integration function, I calculate a jump term for the face. I now want to add the jump term to the shock_indicator variable on the two cells sharing the face. My problem is how to get the dof indices for fe_cell on the current cell, inside the integrate face/boundary functions. Note that the MeshWorker is based on a different fe, not on fe_cell. Thanks praveen On Sun, Mar 27, 2011 at 5:54 PM, Wolfgang Bangerth <[email protected]>wrote: > > > yes, every cell has a index. You can access it by cell->index (). > > But this only gives the index of the cell *within the current mesh level* > and > the set of indices may also have holes. > > Best > W. > >
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
