Check is wrong inside prolongate_3d_o1_rf2.cc

Create issue
Issue #470 closed
Steven R. Brandt created an issue

The problem materialized when creating wave equation tutorials for the Teragrid conference. I was unable to run an FMR application in parallel. The check comprises lines 76-80 of the file:

if (not regbbox.expand(offsetlo, offsethi).is_contained_in(srcbbox) or
    not regbbox                           .is_contained_in(dstbbox))
{
  CCTK_WARN (0, "Internal error: region extent is not contained in array extent");
}

Keyword:

Comments (5)

  1. Erik Schnetter
    • removed comment

    It seems that offsethi is set to 1 even if the interpolation falls onto a grid point, where no offset is required. Other parts of Carpet handle this correctly and pass in an interpolation region that is large enough, but smaller than expected here.

  2. Roland Haas
    • removed comment

    This is fixed in the mercurial version of Carpet which has become the official verion in the ET Maxwell release. Should we close the ticket since it only affects an "old" version of Carpet?

  3. Log in to comment