Method construct_old_gravity

Source/gravity/Castro_gravity.cpp:12–140

using namespace amrex;

void
Castro::construct_old_gravity (Real time)
{
    BL_PROFILE("Castro::construct_old_gravity()");

    const Real strt_time = ParallelDescriptor::second();

    MultiFab& grav_old = get_old_data(Gravity_Type);
    MultiFab& phi_old = get_old_data(PhiGrav_Type);

    // Always set phi to zero initially since some gravity modes
    // don't use it and we want to have valid data.

    if (gravity->get_gravity_type() != "PoissonGrav") {
        phi_old.setVal(0.0);
    }

    if (!do_grav) {

        grav_old.setVal(0.0);

        return;

    }

    // Do level solve at beginning of time step in order to compute the
    // difference between the multilevel and the single level solutions.
    // Note that we don't need to do this solve for single-level runs,
    // since the solution at the end of the last timestep won't have changed.
    // Similarly, we can skip this if we aren't subcycling on this level
    // or all levels above this level, since the new-time composite solve
    // at the new time in the last step will still be valid.

    bool do_old_solve = true;

    if (parent->subcyclingMode() == "None" || parent->finestLevel() == 0) {
        do_old_solve = false;
    }

    if (gravity->get_gravity_type() == "PoissonGrav" && do_old_solve)
    {
        // Create a copy of the current (composite) data on this level.

        MultiFab comp_phi;
        Vector<std::unique_ptr<MultiFab> > comp_gphi(AMREX_SPACEDIM);

        if (gravity->DoCompositeCorrection() && level < parent->finestLevel() && level <= gravity->get_max_solve_level()) {

            comp_phi.define(phi_old.boxArray(), phi_old.DistributionMap(), phi_old.nComp(), phi_old.nGrow());
            MultiFab::Copy(comp_phi, phi_old, 0, 0, phi_old.nComp(), phi_old.nGrow());

            for (int n = 0; n < AMREX_SPACEDIM; ++n) {
                comp_gphi[n] = std::make_unique<MultiFab>(getEdgeBoxArray(n), dmap, 1, 0);
                MultiFab::Copy(*comp_gphi[n], *gravity->get_grad_phi_prev(level)[n], 0, 0, 1, 0);
            }

        }