improved structure in LagrangeElasticity

src/elasticity/assemble/LagrangeElasticity.cpp

@@ -278,12 +278,10 @@ void LagrangeElasticity::Jacobi(const cell &c, const Vector &U, Matrix &A) const
         double localPressure = -eProblem.Pressure(Time(), FE.QPoint(q), E.Bnd(face));
         for (int i = 0; i < FE.NodalPoints(); ++i) {
           for (int k = 0; k < c.dim(); ++k) {
-            auto u_ik = FE.VectorValue(q, i, k);
-            Tensor F_ik = FE.VectorGradient(q, i, k);
+            auto u_ik = FE.VectorComponentValue(q, i, k);
             for (int j = 0; j < FE.NodalPoints(); ++j) {
               for (int l = 0; l < c.dim(); ++l) {
-                auto u_jl = FE.VectorValue(q, j, l);
-                auto F_jl = FE.VectorGradient(q, j, l);
+                auto F_jl = FE.VectorRowGradient(q, j, l);
                 A_c(i, j, k, l) -= w * (localPressure) * (Cross(F, F_jl) * N) * u_ik;
               }
             }