Merge pull request #264 from hklion/init_old_vels

initialize *vel_old to zero to avoid math on unitialized data
seahorce-scidac · Sep 6, 2024 · e9aba33 · e9aba33
2 parents 54fef6d + 3accc1c
commit e9aba33
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);