Add runtime parameter tensor_diffusion.use_gauss_seidel

The default is true. One can use Jacobi smoothers for tensor, mac and nodal solves with tensor.use_gauss_seidel = false mac_proj.use_gauss_seidel = false nodal_proj.use_gauss_seidel = false
AMReX-Fluids · Aug 18, 2024 · 167bfea · 167bfea
1 parent 94aa2ba
commit 167bfea
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 0 deletions.
diff --git a/src/diffusion/DiffusionTensorOp.H b/src/diffusion/DiffusionTensorOp.H
@@ -66,6 +66,7 @@ private:
     amrex::Real m_mg_rtol = amrex::Real(1.0e-11);
     amrex::Real m_mg_atol = amrex::Real(1.0e-14);
 #endif
+    bool m_mg_use_gauss_seidel = true;
     std::string m_bottom_solver = "bicgstab";
 };
 

diff --git a/src/diffusion/DiffusionTensorOp.cpp b/src/diffusion/DiffusionTensorOp.cpp
@@ -59,6 +59,7 @@ DiffusionTensorOp::DiffusionTensorOp (incflo* a_incflo)
                                                 m_incflo->boxArray(0,finest_level),
                                                 m_incflo->DistributionMap(0,finest_level),
                                                 info_solve);
+            m_reg_solve_op->setGaussSeidel(m_mg_use_gauss_seidel);
             m_reg_solve_op->setMaxOrder(m_mg_maxorder);
             m_reg_solve_op->setDomainBC(m_incflo->get_diffuse_tensor_bc(Orientation::low),
                                         m_incflo->get_diffuse_tensor_bc(Orientation::high));
@@ -96,6 +97,8 @@ DiffusionTensorOp::readParameters ()
 
     pp.query("num_pre_smooth", m_num_pre_smooth);
     pp.query("num_post_smooth", m_num_post_smooth);
+
+    pp.query("use_gauss_seidel", m_mg_use_gauss_seidel);
 }
 
 void