Function TestAMRBoxIntersection

Common/DataModel/Testing/Cxx/TestAMRBox.cxx:222–277

------------------------------------------------------------------------------

//------------------------------------------------------------------------------
int TestAMRBoxIntersection()
{
  int rc = 0;

  double h[3];
  int lo[3];
  int hi[3];

  vtkAMRBox A0, A, B, I;

  // Here is the initial AMR box
  h[0] = h[1] = h[2] = 1.0;
  lo[0] = lo[1] = lo[2] = 8;
  hi[0] = hi[1] = hi[2] = 16;
  Construct3DAMRBox(A, lo, hi);

  // Save the initial
  A0 = A;

  B = A;
  B.Shrink(2);
  bool doesIntersect = A.Intersect(B);
  if (!doesIntersect || (A != B))
  {
    std::cerr << "ERROR: Intersecting a fully encompassing box failed!\n";
    rc++;
  }

  A = A0;
  B = A;
  B.Shift(2, 2, 2);

  // Here is the expected box after intersecting
  h[0] = h[1] = h[2] = 1.0;
  lo[0] = lo[1] = lo[2] = 10;
  hi[0] = hi[1] = hi[2] = 16;
  Construct3DAMRBox(I, lo, hi);

  doesIntersect = A.Intersect(B);
  if (!doesIntersect || (A != I))
  {
    std::cerr << "ERROR: Intersecting a partially overlapping box failed!\n";
    rc++;
  }

  A = A0;
  B = A;
  B.Shift(10, 10, 10);
  doesIntersect = A.Intersect(B);
  if (doesIntersect)
  {
    std::cerr << "ERROR: Intersecting a non-overlapping box failed!\n";
    rc++;
  }
  return (rc);
}

//------------------------------------------------------------------------------