> I am using discontinuous elements :) In that case you'd update your deal.II copy from svn and use the function DoFTools::extract_dofs_with_support_on_boundary that I just checked in ;-)
Alternatively, you can also simply copy it from here: http://www.dealii.org/svn/dealii/trunk/deal.II/deal.II/source/dofs/dof_tools.cc If you search for the function, you'll find a version for the 1d case and one for the 2/3d case. That said, there was also a bug in the FE_DGQ class that you'll need to get the new function to work correctly. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
