Bug fixes for ustruct Guccione model and follower pressure load (SimV…

…ascular#201)

* Bug fixes for ustruct Gucci and follower pressure

* Reverting ustruct/LV_Guccione_active max_iterations

* Updating LV_Guccione_active reference solution

* Replacing Sb = Sb + with Sb +=

Code/Source/svFSI/mat_models.cpp

@@ -690,10 +690,10 @@ void get_pk2cc_dev(const ComMod& com_mod, const CepMod& cep_mod, const dmnType&
            2.0*g3*(Es(0,1)*RmRm.slice(3) + Es(0,2)*RmRm.slice(5));
 
       auto CCb = 2.0*ten_dyad_prod(Sb, Sb, nsd);
-      Sb += Sb * r2;
+      Sb = Sb * r2;
 
       // Fiber reinforcement/active stress
-      Sb += Sb + Tfa*mat_dyad_prod(fl.col(0), fl.col(0), nsd);
+      Sb += Tfa*mat_dyad_prod(fl.col(0), fl.col(0), nsd);
 
       double r1 = J2d*mat_ddot(C, Sb, nsd) / nd;
       S = J2d*Sb - r1*Ci;

Code/Source/svFSI/ustruct.cpp

@@ -142,7 +142,7 @@ void b_ustruct_3d(const ComMod& com_mod, const int eNoN, const double w, const V
   for (int a = 0; a  < eNoN; a++) {
     NxFi(0,a) = Nx(0,a)*Fi(0,0) + Nx(1,a)*Fi(1,0) + Nx(2,a)*Fi(2,0);
     NxFi(1,a) = Nx(0,a)*Fi(0,1) + Nx(1,a)*Fi(1,1) + Nx(2,a)*Fi(2,1);
-    NxFi(2,a) = Nx(1,a)*Fi(1,2) + Nx(1,a)*Fi(1,2) + Nx(2,a)*Fi(2,2);
+    NxFi(2,a) = Nx(0,a)*Fi(0,2) + Nx(1,a)*Fi(1,2) + Nx(2,a)*Fi(2,2);
   }
 
   nFi(0) = nV(0)*Fi(0,0) + nV(1)*Fi(1,0) + nV(2)*Fi(2,0);