Jaekwang, On Tuesday, November 1, 2016 at 5:43:36 PM UTC-4, Jaekwang Kim wrote: > > I think that I didn't make any mistake on my code lines because at least > my numerical solution is converging for the case p=1. > > but I couldn't get such convergence for the case p=2, and 3. I have tried > to change tolerance number in iterative scheme, but it didn't work. > > > Are there other things that I failed to consider? > > Or did I make wrong way to test the code? > What you did looks good. You should check that you use the right quadrature everywhere. If somewhere in your code you use a low order quadrature, it will work for p=1 but you will have a larger error when you increase the quadrature. You should also look at the solution for p=2 and p=3. This might help you understand what is going on.
Best, Bruno -- 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]. For more options, visit https://groups.google.com/d/optout.
