Method flame_width_properties

Exec/science/flame/Prob.cpp:7–57

using namespace amrex;

void
Castro::flame_width_properties (amrex::Real time,
                                amrex::Real& T_max, amrex::Real& T_min,
                                amrex::Real& grad_T_max)
{
    BL_PROFILE("Castro::flame_width_properties()");

    const auto dx = geom.CellSizeArray();

    auto mf = derive("Temp",time,1);

    BL_ASSERT(mf != nullptr);

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

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

        const auto temp = (*mf)[mfi].array();

        reduce_op.eval(box, reduce_data,
        [=] AMREX_GPU_HOST_DEVICE (int i, int j, int k) -> ReduceTuple
        {
            // Assumes 1D simulation, right now. Also assumes that
            // we have at least one ghost cell in the x dimension.

            amrex::Real T = temp(i,j,k);
            amrex::Real grad_T = std::abs(temp(i+1,j,k) - temp(i-1,j,k)) / (2.0_rt * dx[0]);

            // Ignore problem zones where we have a negative temperature

            if (T < 0.0_rt) {
                T = 0.0_rt;
                grad_T = 0.0_rt;
            }

            return {T, T, grad_T};
        });
    }

    ReduceTuple hv = reduce_data.value();
    T_max = amrex::max(T_max, amrex::get<0>(hv));
    T_min = amrex::min(T_min, amrex::get<1>(hv));
    grad_T_max = amrex::max(grad_T_max, amrex::get<2>(hv));
}