Leave density-weighted mass matrix update in updateDensity fcn

pecos · Dec 3, 2024 · c5c636f · c5c636f
1 parent 461cea5
commit c5c636f
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/src/reactingFlow.cpp b/src/reactingFlow.cpp
@@ -1522,10 +1522,6 @@ void ReactingFlow::step() {
 
   updateDensity(1.0);
   updateDiffusivity();
-  Array<int> empty;
-  MsRho_form_->Update();
-  MsRho_form_->Assemble();
-  MsRho_form_->FormSystemMatrix(empty, MsRho_);
 
   // PART I: Form and solve implicit systems for species and temperature.
   //
@@ -2545,6 +2541,10 @@ void ReactingFlow::updateDensity(double tStep) {
   }
   rn_gf_.SetFromTrueDofs(rn_);
 
+  MsRho_form_->Update();
+  MsRho_form_->Assemble();
+  MsRho_form_->FormSystemMatrix(empty, MsRho_);
+
   // project to p-space in case not same as vel-temp
   R0PM0_gf_.SetFromTrueDofs(rn_);
 }