XQFT6DACFOMNNXDYUZRGRCMEPWF34KW5PRZ7RVMVJHYCS76JIZSAC
DFELEQUVC3SUCWUKMBHDVCIDMCHWCUAJBC47TA5NNRUYEIYA47GAC
XU5HOJREK4XY4NBCJINLZPKQNKSYOLUDTWR47REFSNQKDOSNDXLQC
2DKSL6DKZAIYQUJGDULORCKU5K4Z5Z3W4RIKQYDSLKMCNQNDZFBAC
722HZ7UFINNE3YKSYKP2NHZ5XEG5QQLQHSKC7PREJZR3EX6RDYUAC
M5R6KQLXLGYSVKHVAX5AJKD6NYE6IM5Z6WVTR3BTKPJDNNKF3ARAC
BJDGFYBMECTJG7BHLNHLSCUCBVYHAY6OGY37FIJP6JDGNDXQNQVAC
33IC3UHCEPZLGS5ACS2JXHGT6CRU5LXU6PM6RDHCERPOIELVRVXQC
RTNZAS3UPI6GG3KY4Z5WVXJ4R2YF5427BB6WAV3GHRS5W7XPOSUQC
KG47IF4CPCUT3BHS34WDRHTH5HYMBTY4OSTB3X7APR2E5ZJ32IYQC
MSBBCXVGD3GRLE5KAI6BKAFRV7SQUWI2SNN43AJAUD3ISRCEXY6QC
TVBD244E7Q7WV44CRBTFST535NUP3JAZH6OLL4IKDR3OWEXSU7HAC
JHHVFKSN2L3YMG4JKXUDS6EBBYCDUDO56LNAZMG4QYS2FFPZF3SQC
3HPNY5XHNPWR6RVX4LYEVDFJQWSB2CMS7J6XPAKATQ2UPFTXD33QC
BVR7DVINVPQG7PA6Z7QYVYNQ43YZL7XCC6AOMSMWMGAAB2Q43STAC
BPRNUTY7MHK7LK4EY5MY5OFFG3ABOL7LWXD574L35M74YSQPULFAC
NE2O3IMQBR3OAPLJC22SPJ7AX2ABULA7XNSIZYFHCQLWC45T6GJAC
Z2NTVAM6CXOOV34ZXN6A5WNAUOPWPFTSH5723LZZEEGYEWH36WWAC
JD6PQOJ6YYNQYEEWEXO2NM7NVYNBUI6V7ZU6Q3FNHGAT2VYOF5WAC
UUGQGVC4WEKN64WAP7F5QPS2UHGQB5ZLMFRIYNWKMIEBDO3QRX4AC
FEMASUBNU32NSG4DNXZX54CGCA57PVRGYO46L3A6F2EJ4BCSJ3SAC
5XGIB7XMEZNBLA5ZLQTXRTC3AZ5CTRGMXWBPVMWXG4DPHKWDF4ZAC
int maxnumlevels = ghext->amrcore->maxLevel() + 1;
for (int level = 0; level < maxnumlevels; ++level) {
CCTK_VINFO("Geometry level %d:", level);
cout << ghext->amrcore->Geom(level) << "\n";
#pragma omp critical
{
const int maxnumlevels = ghext->amrcore->maxLevel() + 1;
for (int level = 0; level < maxnumlevels; ++level) {
CCTK_VINFO("Geometry level %d:", level);
cout << ghext->amrcore->Geom(level) << "\n";
}
for (int level = 0; level < int(ghext->leveldata.size()); ++level)
CCTK_VINFO("norm_inf rhs[%d]: %g", level,
double(rhss.at(level)->norminf(0, 0, false, true)));
for (int level = 0; level < int(ghext->leveldata.size()); ++level)
CCTK_VINFO("norm_inf sol[%d]: %g", level,
double(sols.at(level)->norminf(0, 0, false, true)));
for (int level = 0; level < int(ghext->leveldata.size()); ++level)
CCTK_VINFO("norm_inf res[%d]: %g", level,
double(ress.at(level)->norminf(0, 0, false, true)));
#pragma omp critical
{
for (int level = 0; level < int(ghext->leveldata.size()); ++level)
CCTK_VINFO("norm_inf rhs[%d]: %g", level,
double(rhss.at(level)->norminf(0, 0, false, true)));
for (int level = 0; level < int(ghext->leveldata.size()); ++level)
CCTK_VINFO("norm_inf sol[%d]: %g", level,
double(sols.at(level)->norminf(0, 0, false, true)));
for (int level = 0; level < int(ghext->leveldata.size()); ++level)
CCTK_VINFO("norm_inf res[%d]: %g", level,
double(ress.at(level)->norminf(0, 0, false, true)));
}
nan_imin[2], nan_imax[0], nan_imax[1], nan_imax[2], double(nan_xmin[0]),
double(nan_xmin[1]), double(nan_xmin[2]), double(nan_xmax[0]),
double(nan_xmax[1]), double(nan_xmax[2]),
nan_imin[2], nan_imax[0], nan_imax[1], nan_imax[2],
double(nan_xmin[0]), double(nan_xmin[1]), double(nan_xmin[2]),
double(nan_xmax[0]), double(nan_xmax[1]), double(nan_xmax[2]),
CCTK_VINFO(" xmin=[%g,%g,%g]", x0[0], x0[1], x0[2]);
CCTK_VINFO(" xmax=[%g,%g,%g]", x1[0], x1[1], x1[2]);
CCTK_VINFO(" base dx=[%g,%g,%g]", dx[0], dx[1], dx[2]);
CCTK_VINFO(" xmin=[%.17g,%.17g,%.17g]", double(x0[0]), double(x0[1]),
double(x0[2]));
CCTK_VINFO(" xmax=[%.17g,%.17g,%.17g]", double(x1[0]), double(x1[1]),
double(x1[2]));
CCTK_VINFO(" base dx=[%.17g,%.17g,%.17g]", double(dx[0]),
double(dx[1]), double(dx[2]));