[deal.II] Re: Why fe.shape_value( const unsigned int i, const Point< dim > & p ) return a negative number?

2018-05-05 Thread Zhao Yidong
Thanks a lot, it's really clear to understand by the help of the website you share! Indeed, my problem is this: 1. I have fe which is a FiniteElement class, and I have fe_values which is a FEValues class. I also have QGauss quadrature to create fe_values. 2. I iterate on every quadrature point

[deal.II] Re: Why fe.shape_value( const unsigned int i, const Point< dim > & p ) return a negative number?

2018-05-05 Thread Denis Davydov
Hi, it's perfectly fine for a general FE shape functions to be negative at some points, that's clearly the case for quadratic http://hplgit.github.io/INF5620/doc/pub/sphinx-fem/._main_fem003.html#fem-approx-fe-fig-p2 Denis. On Saturday, May 5, 2018 at 3:11:24 AM UTC+2, Zhao Yidong wrote: > >