initialize *vel_old to zero to avoid math on unitialized data

seahorce-scidac · Sep 6, 2024 · 3accc1c · 3accc1c
1 parent 54fef6d
commit 3accc1c
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/Source/REMORA.cpp b/Source/REMORA.cpp
@@ -600,6 +600,10 @@ REMORA::init_only (int lev, Real time)
     yvel_new[lev]->setVal(0.0_rt);
     zvel_new[lev]->setVal(0.0_rt);
 
+    xvel_old[lev]->setVal(0.0_rt);
+    yvel_old[lev]->setVal(0.0_rt);
+    zvel_old[lev]->setVal(0.0_rt);
+
     vec_ru[lev]->setVal(0.0_rt);
     vec_rv[lev]->setVal(0.0_rt);