In the viscous tensor solve, do better for the initial guess given

to the solver
cgilet · Apr 18, 2024 · d4a8d9f · d4a8d9f
1 parent 308603c
commit d4a8d9f
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/Source/Diffusion.cpp b/Source/Diffusion.cpp
@@ -827,6 +827,9 @@ Diffusion::diffuse_tensor_velocity (Real                   dt,
              if ( has_delta_rhs ) {
                 rhs(i,j,k,n) += deltarhs(i,j,k,n) * dt;
              }
+             // Put unew back since it's a reference to the MF that
+             // ultimately gets used as an initial guess for the solve
+             unew(i,j,k,n) /= rho(i,j,k);
           }
        });
     }