Felix, the usual approach for verifying your code is the method of manufactured solution. Choose functions that are elements of your ansatz space for all solution variables (possibly chaningyour grid to a cartesian one) and compute the corresponding right-hand side. Check that the error vanishes. If that works, choose smooth functions not part of your ansatz space and check the rate of convergence when refining your mesh. You may also want to have a look at https://www.dealii.org/8.5.0/doxygen/deal.II/step_7.html for this.
Best, Daniel -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/eca2f947-b72f-4385-bba4-1c7e7b7baede%40googlegroups.com.
