If you plan to use any domain that is not a square (or an affine transformation), you have to make sure you integrate exactly the product of two polynomials of order degree and of the determinant of the Jacobian. This last term is constant only for simple meshes, but it is the square root of a polynomial of order (degree-1) in more complicated cases.
2*fe_degree is ok for most cases, but I would not use this for serious calculations. I prefer to be on the safe side… :) L. > On 11 Apr 2019, at 19:34, Robert Spartus <bobspar...@gmail.com> wrote: > > Dear Luca, > > Thanks for your suggestion. Unfortunately, it did not solve the problem. I am > sending the modified version, as well as the output of the program. > > Out of curiosity, what is the reason to use (2*fe_degree + 1)? Checking > step-8, I notice that there a quadrature degree one larger than the > polynomial degree is also used. > > Bests, > Bob > -- 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
