Shungho, if you do DG, you do not want to eliminate constraints. You rather want to compute fluxes between the corresponding boundary cells.
Best, Guido > I am doing with DG (based on tutorial 11) + periodic boundary via > Constrant Matrix as given in deal.II library. I think that it surely needs > to eliminate(or condense) the constraint lines after assembling global > matrix. However, I could not find the relevant function for this one in > deal.II docoment. The type of matrix is now > "TrilinosWrappers::BlockSparseMatrix". > _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
