Method Gravity

Source/gravity/Gravity.cpp:64–102

AMREX_GPU_MANAGED Array2D<Real, 0, multipole::lnum_max, 0, multipole::lnum_max> multipole::parity_qC_qS;

Gravity::Gravity(Amr* Parent, int _finest_level, BCRec* _phys_bc, int _density)
  :
    parent(Parent),
    LevelData(MAX_LEV),
    grad_phi_curr(MAX_LEV),
    grad_phi_prev(MAX_LEV),
    grids(Parent->boxArray()),
    dmap(Parent->DistributionMap()),
    abs_tol(MAX_LEV),
    rel_tol(MAX_LEV),
    level_solver_resnorm(MAX_LEV),
    volume(MAX_LEV),
    area(MAX_LEV),
    phys_bc(_phys_bc)
{

     amrex::ignore_unused(_finest_level);

     AMREX_ALWAYS_ASSERT(parent->maxLevel() < MAX_LEV);

     Density = _density;
     read_params();
     finest_level_allocated = -1;

     radial_grav_old.resize(MAX_LEV);
     radial_grav_new.resize(MAX_LEV);
     radial_mass.resize(MAX_LEV);
     radial_vol.resize(MAX_LEV);
#ifdef GR_GRAV
     radial_pres.resize(MAX_LEV);
#endif

     if (gravity::gravity_type == "PoissonGrav") {
         make_mg_bc();
         init_multipole_grav();
     }
     max_rhs = 0.0;
     numpts_at_level = -1;
}

Gravity::~Gravity() = default;