Skip to content

Fixes in PolynomialBasis for derivatives in higher dimension

Maik Porrmann requested to merge maik.porrmann/dune-localfunctions:fix_PolynomialBasis into master

I found several bugs in the PolynomialBasis class when using (higher order) partial derivatives in three dimensions. The Errors can be reproduced by changing line 28 (or any of the corresponding lines below) of localfunctions/test/raviartthomaselementtest.cc to TEST_FE3(rt0simplex3dlfem, DisableNone, 2);

Fixes:

  • Missing template parameter in Base typedef of PolynomialBasisWithMatrix
  • Transforming multiindex for partial now works when differentiating twice in the same direction
  • Change order of derivatives when assembling hessian from vector with unique partial derivatives
  • q<hsize assertion moved before increment
  • resize vector at beginning of partial
Edited by Maik Porrmann

Merge request reports