QNV4LD7UGYSSNYDXGJC6SRP6Y3USUVDQLEERBMBWRC7LILWREB7AC
PG2P3JQPLWZRLCAJ7JY2B7G7AM5DHOE2CJUKIPWREI3NAUKZF3TAC
RF3SDQ3ELEBCGJSFQL2VOSHNN34RAH6KM3LEITEYOHKJMMHHKMFAC
2DKSL6DKZAIYQUJGDULORCKU5K4Z5Z3W4RIKQYDSLKMCNQNDZFBAC
722HZ7UFINNE3YKSYKP2NHZ5XEG5QQLQHSKC7PREJZR3EX6RDYUAC
YIQN7NJTGEVKW7JZHL6CTH6EPCIXCNBYNURIGXPYZAOUX3VAJQMAC
FEMASUBNU32NSG4DNXZX54CGCA57PVRGYO46L3A6F2EJ4BCSJ3SAC
33IC3UHCEPZLGS5ACS2JXHGT6CRU5LXU6PM6RDHCERPOIELVRVXQC
UZAKARMGORRQG733ZUPJOEGL5FG243I32NCC2SRSFDCZKUQ5A52QC
BVR7DVINVPQG7PA6Z7QYVYNQ43YZL7XCC6AOMSMWMGAAB2Q43STAC
// Set grid functions to nan
auto mfitinfo = MFItInfo().SetDynamic(true).EnableTiling(
{max_tile_size_x, max_tile_size_y, max_tile_size_z});
if (poison_undefined_points) {
// Set grid functions to nan
auto mfitinfo = MFItInfo().SetDynamic(true).EnableTiling(
{max_tile_size_x, max_tile_size_y, max_tile_size_z});
for (MFIter mfi(*leveldata.mfab0, mfitinfo); mfi.isValid(); ++mfi) {
GridPtrDesc grid(leveldata, mfi);
const Array4<CCTK_REAL> &vars = groupdata.mfab.at(tl)->array(mfi);
for (int vi = 0; vi < groupdata.numvars; ++vi) {
CCTK_REAL *restrict const ptr = grid.ptr(vars, vi);
grid.loop_all(groupdata.indextype, [&](const Loop::PointDesc &p) {
ptr[p.idx] = 0.0 / 0.0;
});
for (MFIter mfi(*leveldata.mfab0, mfitinfo); mfi.isValid(); ++mfi) {
GridPtrDesc grid(leveldata, mfi);
const Array4<CCTK_REAL> &vars = groupdata.mfab.at(tl)->array(mfi);
for (int vi = 0; vi < groupdata.numvars; ++vi) {
CCTK_REAL *restrict const ptr = grid.ptr(vars, vi);
grid.loop_all(groupdata.indextype, [&](const Loop::PointDesc &p) {
ptr[p.idx] = 0.0 / 0.0;
});
}
// Set grid functions to nan
auto mfitinfo = MFItInfo().SetDynamic(true).EnableTiling(
{max_tile_size_x, max_tile_size_y, max_tile_size_z});
if (poison_undefined_points) {
// Set grid functions to nan
auto mfitinfo = MFItInfo().SetDynamic(true).EnableTiling(
{max_tile_size_x, max_tile_size_y, max_tile_size_z});
for (MFIter mfi(*leveldata.mfab0, mfitinfo); mfi.isValid(); ++mfi) {
GridPtrDesc grid(leveldata, mfi);
const Array4<CCTK_REAL> &vars = mfab->array(mfi);
for (int vi = 0; vi < groupdata.numvars; ++vi) {
CCTK_REAL *restrict const ptr = grid.ptr(vars, vi);
grid.loop_all(groupdata.indextype, [&](const Loop::PointDesc &p) {
ptr[p.idx] = 0.0 / 0.0;
});
for (MFIter mfi(*leveldata.mfab0, mfitinfo); mfi.isValid(); ++mfi) {
GridPtrDesc grid(leveldata, mfi);
const Array4<CCTK_REAL> &vars = mfab->array(mfi);
for (int vi = 0; vi < groupdata.numvars; ++vi) {
CCTK_REAL *restrict const ptr = grid.ptr(vars, vi);
grid.loop_all(groupdata.indextype, [&](const Loop::PointDesc &p) {
ptr[p.idx] = 0.0 / 0.0;
});
}
// Set grid functions to nan
auto mfitinfo = MFItInfo().SetDynamic(true).EnableTiling(
{max_tile_size_x, max_tile_size_y, max_tile_size_z});
if (poison_undefined_points) {
// Set grid functions to nan
const int tl = 0;
auto mfitinfo = MFItInfo().SetDynamic(true).EnableTiling(
{max_tile_size_x, max_tile_size_y, max_tile_size_z});
for (MFIter mfi(*leveldata.mfab0, mfitinfo); mfi.isValid(); ++mfi) {
GridPtrDesc grid(leveldata, mfi);
const int tl = 0;
const Array4<CCTK_REAL> &vars = groupdata.mfab.at(tl)->array(mfi);
for (int vi = 0; vi < groupdata.numvars; ++vi) {
CCTK_REAL *restrict const ptr = grid.ptr(vars, vi);
grid.loop_all(groupdata.indextype, [&](const Loop::PointDesc &p) {
ptr[p.idx] = 0.0 / 0.0;
});
for (MFIter mfi(*leveldata.mfab0, mfitinfo); mfi.isValid(); ++mfi) {
GridPtrDesc grid(leveldata, mfi);
const Array4<CCTK_REAL> &vars = groupdata.mfab.at(tl)->array(mfi);
for (int vi = 0; vi < groupdata.numvars; ++vi) {
CCTK_REAL *restrict const ptr = grid.ptr(vars, vi);
grid.loop_all(groupdata.indextype, [&](const Loop::PointDesc &p) {
ptr[p.idx] = 0.0 / 0.0;
});
}