ERF_Advance.cpp File Reference

#include <ERF.H>
#include <ERF_Utils.H>

Include dependency graph for ERF_Advance.cpp:

Functions
void	check_for_negative_theta (amrex::MultiFab &S_old)

Function Documentation

check_for_negative_theta()

void check_for_negative_theta

(

amrex::MultiFab &

S_old

)

Function that advances the solution at one level for a single time step – this does some preliminaries then calls erf_advance

Parameters

[in]	lev	level of refinement (coarsest level is 0)
[in]	time	start time for time advance
[in]	dt_lev	time step for this time advance

{
    // *****************************************************************************
    // Test for negative (rho theta)
    // *****************************************************************************
    for (MFIter mfi(S_old); mfi.isValid(); ++mfi)
    {
        Box bx = mfi.tilebox();
        const Array4<Real> &cell_data  = S_old.array(mfi);
        ParallelFor(bx, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
        {
#ifdef AMREX_USE_GPU
            if (cell_data(i,j,k,RhoTheta_comp) <= 0.) AMREX_DEVICE_PRINTF("BAD THETA AT %d %d %d %e %e \n",
                i,j,k,cell_data(i,j,k,RhoTheta_comp),cell_data(i,j,k+1,RhoTheta_comp));
#else
            if (cell_data(i,j,k,RhoTheta_comp) <= 0.) {
                printf("BAD THETA AT %d %d %d %e %e \n",
                i,j,k,cell_data(i,j,k,RhoTheta_comp),cell_data(i,j,k+1,RhoTheta_comp));
                amrex::Abort("Bad theta in check_for_negative_theta");
            }
#endif
            });
    } // mfi
}

Referenced by ERF::Advance().

Here is the caller graph for this function: