Hello, I am trying to compute the L2/H1 error norms for an elliptic problem defined on two different grids (not necessarily embedded, different h-p refinements) for convergence studies. In the presence of an exact solution, I can use the VectorTools::integrate_difference routine to accomplish what I need but am having a hard time doing it cleanly without resorting to a separate projection of the solution on to the target mesh and a consequent wrapping of the resultant solution using a Function to use the VectorTools routine like above in the case when I have a reference solution on a very finely refined grid.
This idea perhaps seems roundabout and was wondering if there was a better API that I have overlooked to accomplish this task. Any pointers to examples or routines that might help out will be much appreciated. Thank you, Vijay _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
