Method problem_post_init

Exec/gravity_tests/uniform_cube/Prob.cpp:15–159

using namespace amrex;

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

    gravity->multilevel_solve_for_new_phi(0, parent->finestLevel());

    const int norm_power = 2;

    ReduceOps<ReduceOpSum, ReduceOpSum> reduce_op;
    ReduceData<Real, Real> reduce_data(reduce_op);
    using ReduceTuple = typename decltype(reduce_data)::Type;

    for (int lev = 0; lev <= parent->finestLevel(); lev++)
    {
        const auto dx = getLevel(lev).geom.CellSizeArray();
        const auto problo = getLevel(lev).geom.ProbLoArray();

        const Real time = getLevel(lev).state[State_Type].curTime();

        auto phiGrav = getLevel(lev).derive("phiGrav", time, 0);

        if (lev < parent->finestLevel())
        {
            const MultiFab& mask = getLevel(lev+1).build_fine_mask();
            MultiFab::Multiply(*phiGrav, mask, 0, 0, 1, 0);
        }

#ifdef _OPENMP
#pragma omp parallel
#endif
        for (MFIter mfi(*phiGrav, TilingIfNotGPU()); mfi.isValid(); ++mfi)
        {
            const Box& box = mfi.tilebox();

            auto phi = (*phiGrav)[mfi].array();
            auto vol = getLevel(lev).Volume()[mfi].array();

            // Compute the norm of the difference between the calculated potential
            // and the analytical solution.

            reduce_op.eval(box, reduce_data,
            [=] AMREX_GPU_HOST_DEVICE (int i, int j, int k) -> ReduceTuple
            {
                Real radius = 0.5_rt * problem::diameter;
                Real mass = (4.0_rt / 3.0_rt) * M_PI * radius * radius * radius * problem::density;

                Real xx = problo[0] + dx[0] * (static_cast<Real>(i) + 0.5_rt) - problem::center[0];

#if AMREX_SPACEDIM >= 2
                Real yy = problo[1] + dx[1] * (static_cast<Real>(j) + 0.5_rt) - problem::center[1];
#else
                Real yy = 0.0_rt;
#endif

#if AMREX_SPACEDIM == 3
                Real zz = problo[2] + dx[2] * (static_cast<Real>(k) + 0.5_rt) - problem::center[2];
#else
                Real zz = 0.0_rt;