diff options
Diffstat (limited to 'Source/WarpX.cpp')
| -rw-r--r-- | Source/WarpX.cpp | 13 |
1 files changed, 13 insertions, 0 deletions
diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp index 951f566ba..85dfe9333 100644 --- a/Source/WarpX.cpp +++ b/Source/WarpX.cpp @@ -69,6 +69,8 @@ long WarpX::particle_pusher_algo; int WarpX::maxwell_fdtd_solver_id; long WarpX::load_balance_costs_update_algo; int WarpX::do_dive_cleaning = 0; +int WarpX::em_solver_medium; +int WarpX::macroscopic_solver_algo; long WarpX::n_rz_azimuthal_modes = 1; long WarpX::ncomps = 1; @@ -212,6 +214,13 @@ WarpX::WarpX () pml.resize(nlevs_max); costs.resize(nlevs_max); + + if (em_solver_medium == MediumForEM::Macroscopic) { + // create object for macroscopic solver + m_macroscopic_properties = std::unique_ptr<MacroscopicProperties> (new MacroscopicProperties()); + } + + // Set default values for particle and cell weights for costs update; // Default values listed here for the case AMREX_USE_GPU are determined // from single-GPU tests on Summit. @@ -618,6 +627,10 @@ WarpX::ReadParameters () l_lower_order_in_v = false; } load_balance_costs_update_algo = GetAlgorithmInteger(pp, "load_balance_costs_update"); + em_solver_medium = GetAlgorithmInteger(pp, "em_solver_medium"); + if (em_solver_medium == MediumForEM::Macroscopic ) { + macroscopic_solver_algo = GetAlgorithmInteger(pp,"macroscopic_sigma_method"); + } pp.query("costs_heuristic_cells_wt", costs_heuristic_cells_wt); pp.query("costs_heuristic_particles_wt", costs_heuristic_particles_wt); } |