Method construct_old_react_source

Source/sdc/Castro_sdc.cpp:296–384

#ifdef REACTIONS
void
Castro::construct_old_react_source(MultiFab& U_state,
                                   MultiFab& R_source,
                                   const bool input_is_average)
{

    // note: there is a side effect here.  If we are fourth-order and
    // coming in with input_is_average = true, then we make a backup
    // of R_source when it is on cell-centers and store it in Sburn.

    BL_PROFILE("Castro::construct_old_react_source()");

    auto domain_lo = geom.Domain().loVect3d();
    auto domain_hi = geom.Domain().hiVect3d();

    if (sdc_order == 4 && input_is_average)
    {
        // we have cell-averages
        // Note: we cannot tile these operations

        FArrayBox U_center(The_Async_Arena());
        FArrayBox R_center(The_Async_Arena());
        FArrayBox tmp(The_Async_Arena());

        for (MFIter mfi(U_state); mfi.isValid(); ++mfi)
        {

            const Box& bx = mfi.tilebox();
            const Box& obx = mfi.growntilebox(1);

            // Convert to centers
            U_center.resize(obx, NUM_STATE);
            auto const U_center_arr = U_center.array();

            make_cell_center(obx, U_state.array(mfi), U_center_arr, domain_lo, domain_hi);

            // sometimes the Laplacian can make the species go
            // negative near discontinuities
            amrex::ParallelFor(obx,
            [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
            {
                normalize_species_sdc(i, j, k, U_center_arr);
            });

            // burn, including one ghost cell
            R_center.resize(obx, NUM_STATE);
            auto const R_center_arr = R_center.array();

            amrex::ParallelFor(obx,
            [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
            {
                instantaneous_react(i, j, k, U_center_arr, R_center_arr);
            });

            // at this point, we have the reaction term on centers,
            // including a ghost cell.  Save this into Sburn so we can
            // use it later for the plotfile filling
            Sburn[mfi].copy(R_center, obx, 0, obx, 0, NUM_STATE);