ERF_Tagging.cpp File Reference

#include <ERF.H>
#include <ERF_Derive.H>

Include dependency graph for ERF_Tagging.cpp:

Functions
void	tag_on_distance_from_eye (const Geometry &cgeom, TagBoxArray *tags, const Real eye_x, const Real eye_y, const Real rad_tag)

Function Documentation

tag_on_distance_from_eye()

void tag_on_distance_from_eye	(	const Geometry &	cgeom,
		TagBoxArray *	tags,
		const Real	eye_x,
		const Real	eye_y,
		const Real	rad_tag
	)

Function to tag cells for refinement – this overrides the pure virtual function in AmrCore

Parameters

[in]	levc	level of refinement at which we tag cells (0 is coarsest level)
[out]	tags	array of tagged cells
[in]	time	current time

{
    const auto dx      = cgeom.CellSizeArray();
    const auto prob_lo = cgeom.ProbLoArray();
 
    for (MFIter mfi(*tags); mfi.isValid(); ++mfi) {
        TagBox& tag = (*tags)[mfi];
        auto tag_arr = tag.array();  // Get device-accessible array
 
        const Box& tile_box = mfi.tilebox(); // The box for this tile
 
        ParallelFor(tile_box, [=] AMREX_GPU_DEVICE(int i, int j, int k) {
            // Compute cell center coordinates
            Real x = prob_lo[0] + (i + 0.5) * dx[0];
            Real y = prob_lo[1] + (j + 0.5) * dx[1];
 
            Real dist = std::sqrt((x - eye_x)*(x - eye_x) + (y - eye_y)*(y - eye_y));
 
            if (dist < rad_tag) {
                tag_arr(i,j,k) = TagBox::SET;
            } else {
                tag_arr(i,j,k) = TagBox::CLEAR;
            }
        });
    }
}

Referenced by ERF::ErrorEst(), and ERF::HurricaneTracker().

Here is the caller graph for this function: