Is it stil right if value of nodes on two different boundary sufaces/lines are should same value as a constraint for periodic boundary condition in the domain?
I need to have a node number 1 10 grid x--------------------------------x subject to u(1) = u(10) or possibly u(1) = u(10) - g (for inhomogenousity) On Mon, August 3, 2009 5:39 pm, Guido Kanschat wrote: > 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
