Method problem_post_init

Exec/gravity_tests/uniform_sphere/Prob.cpp:15–203

using namespace amrex;

void Castro::problem_post_init()
{
    BL_ASSERT(level == 0);

    // Add up the mass on the domain and then update the density
    // in the sphere so that it has the 'correct' amount of total mass.

    Real actual_mass = 0.0;

    bool local_flag = true;
    Real time = state[State_Type].curTime();

    for (int lev = 0; lev <= parent->finestLevel(); ++lev) {
        actual_mass += getLevel(lev).volWgtSum("density", time, local_flag);
    }

    ParallelDescriptor::ReduceRealSum(actual_mass);

    // The correct amount of mass is the mass of a sphere
    // with the given diameter and density.

    Real target_mass = problem::density * (1.0e0 / 6.0e0) * M_PI * std::pow(problem::diameter, 3);

    Real update_factor = target_mass / actual_mass;

    // Now update the density given this factor.

    amrex::Print() << "\n";
    amrex::Print() << "  Updating density by the factor " << update_factor << " to ensure total mass matches target mass.\n";
    amrex::Print() << "\n";

    problem::density = problem::density * update_factor;

    for (int lev = 0; lev <= parent->finestLevel(); lev++)
    {
        MultiFab& state = getLevel(lev).get_new_data(State_Type);

        const auto dx = getLevel(lev).geom.CellSizeArray();
        const auto problo = getLevel(lev).geom.ProbLoArray();

#ifdef _OPENMP
#pragma omp parallel
#endif
        for (MFIter mfi(state, TilingIfNotGPU()); mfi.isValid(); ++mfi) {

            const Box& box = mfi.tilebox();

            auto u = state[mfi].array();

            // Update the density field. This ensures that the sum of the
            // mass on the domain is what we intend it to be.

            amrex::ParallelFor(box,
            [=] AMREX_GPU_HOST_DEVICE (int i, int j, int k)
            {
                Real xx = problo[0] + dx[0] * (static_cast<Real>(i) + 0.5_rt) - problem::center[0];

#if AMREX_SPACEDIM >= 2