Use self_fields_verbosity value for magnetostatic solver (#5560)

roelof-groenewald · web-flow · commit 961d9a6af1f6 · 2025-01-21T17:25:33.000-08:00
This is another temporary fix (similar to #5517) to close #5530.
diff --git a/Source/FieldSolver/MagnetostaticSolver/MagnetostaticSolver.cpp b/Source/FieldSolver/MagnetostaticSolver/MagnetostaticSolver.cpp
@@ -136,13 +136,12 @@ WarpX::AddMagnetostaticFieldLabFrame()
     else {
         required_precision = 1e-11;
     }
-    const int verbosity = 2;
 
     computeVectorPotential(
         m_fields.get_mr_levels_alldirs(FieldType::current_fp, finest_level),
         m_fields.get_mr_levels_alldirs(FieldType::vector_potential_fp_nodal, finest_level),
         required_precision, absolute_tolerance, magnetostatic_solver_max_iters,
-        verbosity
+        magnetostatic_solver_verbosity
     );
 }
 
diff --git a/Source/WarpX.H b/Source/WarpX.H
@@ -903,6 +903,7 @@ public:
     // Magnetostatic Solver Interface
     MagnetostaticSolver::VectorPoissonBoundaryHandler m_vector_poisson_boundary_handler;
     int magnetostatic_solver_max_iters = 200;
+    int magnetostatic_solver_verbosity = 2;
     void ComputeMagnetostaticField ();
     void AddMagnetostaticFieldLabFrame ();
     void computeVectorPotential (ablastr::fields::MultiLevelVectorField const& curr,
diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp
@@ -692,6 +692,7 @@ WarpX::ReadParameters ()
         "To use the FFT Poisson solver, compile with WARPX_USE_FFT=ON.");
 #endif
         utils::parser::queryWithParser(pp_warpx, "self_fields_max_iters", magnetostatic_solver_max_iters);
+        utils::parser::queryWithParser(pp_warpx, "self_fields_verbosity", magnetostatic_solver_verbosity);
 
         WARPX_ALWAYS_ASSERT_WITH_MESSAGE(
         (

Original file line number	Diff line number	Diff line change
`@@ -136,13 +136,12 @@ WarpX::AddMagnetostaticFieldLabFrame()`
`136`	`136`	`else {`
`137`	`137`	`required_precision = 1e-11;`
`138`	`138`	`}`
`139`		`- const int verbosity = 2;`
`140`	`139`
`141`	`140`	`computeVectorPotential(`
`142`	`141`	`m_fields.get_mr_levels_alldirs(FieldType::current_fp, finest_level),`
`143`	`142`	`m_fields.get_mr_levels_alldirs(FieldType::vector_potential_fp_nodal, finest_level),`
`144`	`143`	`required_precision, absolute_tolerance, magnetostatic_solver_max_iters,`
`145`		`- verbosity`
	`144`	`+ magnetostatic_solver_verbosity`
`146`	`145`	`);`
`147`	`146`	`}`
`148`	`147`
Original file line number	Diff line number	Diff line change
`@@ -692,6 +692,7 @@ WarpX::ReadParameters ()`
`692`	`692`	`"To use the FFT Poisson solver, compile with WARPX_USE_FFT=ON.");`
`693`	`693`	`#endif`
`694`	`694`	`utils::parser::queryWithParser(pp_warpx, "self_fields_max_iters", magnetostatic_solver_max_iters);`
	`695`	`+ utils::parser::queryWithParser(pp_warpx, "self_fields_verbosity", magnetostatic_solver_verbosity);`
`695`	`696`
`696`	`697`	`WARPX_ALWAYS_ASSERT_WITH_MESSAGE(`
`697`	`698`	`(`