Skip to content

Commit

Permalink
Use self_fields_verbosity value for magnetostatic solver
Browse files Browse the repository at this point in the history
  • Loading branch information
roelof-groenewald committed Jan 15, 2025
1 parent 57703f8 commit bd1dc53
Show file tree
Hide file tree
Showing 3 changed files with 3 additions and 2 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -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
);
}

Expand Down
1 change: 1 addition & 0 deletions Source/WarpX.H
Original file line number Diff line number Diff line change
Expand Up @@ -902,6 +902,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,
Expand Down
1 change: 1 addition & 0 deletions Source/WarpX.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -684,6 +684,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(
(
Expand Down

0 comments on commit bd1dc53

Please sign in to comment.