The problem the I see with using a system with LAGRANGE variables is that
the code for assembling the element stiffness matrix becomes a bit messy
and less flexible to switch from 2 to 3 dimensions and vice versa.
For instance, in the example
https://github.com/libMesh/libmesh/blob/master/example
On Wed, 3 Aug 2016, Salazar De Troya, Miguel wrote:
where coef is the vector with the element components of the
solution. The problem that I see here is that for the LAGRANGE
element, dphi[l][qp] is a RealGradient with mostly zeros and it’s
the same information for shape functions on the same n
Hello
I have noticed that dphi from FEVectorBase has a structure like this:
std::vector > dphi;
dphi[shape_funct][qp](i,j)
shape_funct goes from zero to the number of degrees of freedom per element, qp
is the gaussian point. This makes the interpolation operation for the gradient
of a field “u
