try to fix memory issue.

LLNL · Sep 13, 2024 · c96517d · c96517d
1 parent 34da99e
commit c96517d
Showing 1 changed file with 5 additions and 5 deletions.
diff --git a/src/serac/physics/solid_mechanics_contact.hpp b/src/serac/physics/solid_mechanics_contact.hpp
@@ -121,7 +121,7 @@ class SolidMechanicsContact<order, dim, Parameters<parameter_space...>,
       r_blk = res;
       mfem::Vector uPlusShapeDisp(u.Size());
       uPlusShapeDisp = u;
-      uPlusShapeDisp += shape_displacement_;
+      uPlusShapeDisp.Add(1.0, shape_displacement_);
 
       contact_.residualFunction(uPlusShapeDisp, r);
       r_blk.SetSubVector(bcs_.allEssentialTrueDofs(), 0.0);
@@ -143,7 +143,7 @@ class SolidMechanicsContact<order, dim, Parameters<parameter_space...>,
 
             mfem::Vector uPlusShapeDisp(u.Size());
             uPlusShapeDisp = u;
-            uPlusShapeDisp += shape_displacement_;
+            uPlusShapeDisp.Add(1.0, shape_displacement_);
 
             // create block operator holding jacobian contributions
             J_constraint_ = contact_.jacobianFunction(uPlusShapeDisp, J_.release());
@@ -183,7 +183,7 @@ class SolidMechanicsContact<order, dim, Parameters<parameter_space...>,
 
             mfem::Vector uPlusShapeDisp(u.Size());
             uPlusShapeDisp = u;
-            uPlusShapeDisp += shape_displacement_;
+            uPlusShapeDisp.Add(1.0, shape_displacement_);
 
             // get 11-block holding jacobian contributions
             auto block_J         = contact_.jacobianFunction(uPlusShapeDisp, J_.release());
@@ -281,7 +281,7 @@ class SolidMechanicsContact<order, dim, Parameters<parameter_space...>,
 
     mfem::Vector uPlusShapeDisp(displacement_.Size());
     uPlusShapeDisp = displacement_;
-    uPlusShapeDisp += shape_displacement_;
+    uPlusShapeDisp.Add(1.0, shape_displacement_);
 
     auto block_J         = contact_.jacobianFunction(uPlusShapeDisp, jacobian.release());
     block_J->owns_blocks = false;
@@ -308,7 +308,7 @@ class SolidMechanicsContact<order, dim, Parameters<parameter_space...>,
 
     mfem::Vector uPlusShapeDisp(displacement_.Size());
     uPlusShapeDisp = displacement_;
-    uPlusShapeDisp += shape_displacement_;
+    uPlusShapeDisp.Add(1.0, shape_displacement_);
 
     auto block_J         = contact_.jacobianFunction(uPlusShapeDisp, drdshape_mat.release());
     block_J->owns_blocks = false;