Updated EGVEctorField

cardmech/src/local_mpp/EGVectorFieldElement.cpp

@@ -43,12 +43,15 @@ EGVectorFieldElementT<TT, sDim, tDim>::VectorValue(int q, int j) const {
 template<typename TT, int sDim, int tDim>
 VectorFieldT<TT, sDim>
 EGVectorFieldElementT<TT, sDim, tDim>::VectorValue(int q, const Vector &u) const {
-  VectorFieldT<TT, sDim> V;
-  V += u(this->r(this->shape.size()), 0);
+  VectorFieldT<TT, sDim> V{};
+  for (int k = 0; k < dim; ++k) {
+    V[k] += u(this->r(this->shape.size()), k);
+  }
+
   for (int i = 0; i < this->shape.size(); ++i) {
     TT s = this->Value(q, i);
     for (int k = 0; k < dim; ++k) {
-      V[k] += s * u(this->r(i), 0);
+      V[k] += s * u(this->r(i), k);
     }
   }
   return V;
@@ -60,7 +63,7 @@ EGVectorFieldElementT<TT, sDim, tDim>::VectorValue(const PointT<TT, sDim, tDim>
                                                    const Vector &u) const {
   VectorFieldT<TT, sDim> V;
   for (int k = 0; k < dim; ++k) {
-    V[k] += this->Value(z, u);
+    V[k] += this->Value(z, u, k);
   }
   return V;
 }
@@ -93,10 +96,11 @@ TensorT<TT, sDim> EGVectorFieldElementT<TT, sDim, tDim>::VectorGradient(int q, i
 template<typename TT, int sDim, int tDim>
 TensorT<TT, sDim>
 EGVectorFieldElementT<TT, sDim, tDim>::VectorGradient(int q, const Vector &u) const {
-  TensorT<TT, sDim> T;
+  TensorT<TT, sDim> T{};
+
   for (int i = 0; i < this->shape.size(); ++i) {
     for (int k = 0; k < dim; ++k) {
-      VectorFieldT<TT, sDim> G = u(this->r(i), 0) * this->Derivative(q, i);
+      VectorFieldT<TT, sDim> G = u(this->r(i), k) * this->Derivative(q, i);
       for (int l = 0; l < dim; ++l)
         T[k][l] += G[l];
     }