Hi all,

After I solve my PDE for the solution u, I want to compute the integral

\int_\Gamma ( \partial{u} / \partial{n} )^2 ds

where \Gamma is the domain boundary. I have the discrete solution u^h as a 
vector of values; is there a function in the library that helps with this? I 
have looked in VectorTools but could not find something that calculates 
integrals of directional derivatives (it does not fit the scope of the class I 
guess). How would a quadrature calculation look like? I would have to iterate 
over each DG element that has a face on the boundary, get the quadrature 
points/weights, and then sum-up the boundary integral by hand? 


Could I find something helpful in one of the tutorial examples?

Thanks!
  -- M


_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to