Method construct_new_gravity

Source/gravity/Castro_gravity.cpp:142–298

}

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

    const Real strt_time = ParallelDescriptor::second();

    MultiFab& grav_new = get_new_data(Gravity_Type);
    MultiFab& phi_new = get_new_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_new.setVal(0.0);
    }

    if (!do_grav) {

        grav_new.setVal(0.0);

        return;

    }

    // If we're doing Poisson gravity, do the new-time level or composite solve here.

    if (gravity->get_gravity_type() == "PoissonGrav")
    {
        if (level == 0 && parent->subcyclingMode() == "None") {
            if (castro::verbose > 0) {
                amrex::Print() << "\n... new-time composite Poisson gravity solve from level " << level << " to level " << parent->finestLevel() << std::endl << std::endl;
            }

            // Use the "old" phi from the current time step as a guess for this solve.

            for (int lev = level; lev <= parent->finestLevel(); ++lev) {
                MultiFab& lev_phi_old = getLevel(lev).get_old_data(PhiGrav_Type);
                MultiFab& lev_phi_new = getLevel(lev).get_new_data(PhiGrav_Type);

                MultiFab::Copy(lev_phi_new, lev_phi_old, 0, 0, 1, lev_phi_new.nGrow());
            }

            gravity->multilevel_solve_for_new_phi(level, parent->finestLevel());
        }
        else if (parent->subcyclingMode() != "None") {
            // Use the "old" phi from the current time step as a guess for this solve.

            MultiFab& phi_old = get_old_data(PhiGrav_Type);

            MultiFab::Copy(phi_new, phi_old, 0, 0, 1, phi_new.nGrow());

            // Subtract off the (composite - level) contribution for the purposes
            // of the level solve. We'll add it back later.

            if (gravity->DoCompositeCorrection() && level < parent->finestLevel() && level <= gravity->get_max_solve_level()) {
                phi_new.minus(comp_minus_level_phi, 0, 1, 0);
            }