ENH: Added test for AMR Blanking

Added a simple test to ensure that the visibility is computed correctly
for 2D (XY,XZ,YZ) and 3D AMR datasets.

Filtering/Testing/Cxx/CMakeLists.txt

@@ -22,6 +22,7 @@ CREATE_TEST_SOURCELIST(Tests ${KIT}CxxTests.cxx
   quadCellConsistency.cxx
   quadraticEvaluation.cxx
   TestAMRBox.cxx
+  TestAMRBlanking.cxx
   TestInterpolationFunctions.cxx
   TestInterpolationDerivs.cxx
   TestImageDataFindCell.cxx
@@ -44,7 +45,7 @@ CREATE_TEST_SOURCELIST(Tests ${KIT}CxxTests.cxx
 )
 
 ADD_EXECUTABLE(${KIT}CxxTests ${Tests})
-TARGET_LINK_LIBRARIES(${KIT}CxxTests vtkIO vtkFiltering vtkGraphics)
+TARGET_LINK_LIBRARIES(${KIT}CxxTests vtkIO vtkFiltering vtkGraphics vtkParallel)
 
 SET (TestsToRun ${Tests})
 REMOVE (TestsToRun ${KIT}CxxTests.cxx)

Filtering/Testing/Cxx/TestAMRBlanking.cxx

+/*=========================================================================
+
+  Program:   Visualization Toolkit
+  Module:    TestAMRBlanking.cxx
+
+  Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
+  All rights reserved.
+  See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
+
+     This software is distributed WITHOUT ANY WARRANTY; without even
+     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+     PURPOSE.  See the above copyright notice for more information.
+
+=========================================================================*/
+// .NAME TestAMRBlanking.cxx -- Simple test for AMR blanking (visibility)
+//
+// .SECTION Description
+//  Test blanking for Berger-Collela AMR datasets.
+
+#include "vtkOverlappingAMR.h"
+#include "vtkUniformGrid.h"
+#include "vtkAMRUtilities.h"
+#include "vtkStructuredData.h"
+#include "vtkXMLHierarchicalBoxDataWriter.h"
+
+#include <cassert>
+#include <sstream>
+
+namespace AMRVisibilityTests {
+
+//------------------------------------------------------------------------------
+vtkUniformGrid* GetGrid(double origin[3], double h[3], int ndims[3] )
+{
+  vtkUniformGrid* grid = vtkUniformGrid::New();
+  grid->Initialize();
+  grid->SetOrigin( origin );
+  grid->SetSpacing( h );
+  grid->SetDimensions( ndims );
+  return( grid );
+}
+
+//------------------------------------------------------------------------------
+vtkOverlappingAMR* GetAMRDataSet(const int description)
+{
+  double origin[3];
+  double spacing[3];
+  int ndims[3];
+
+  vtkOverlappingAMR *amrDataSet = vtkOverlappingAMR::New();
+  vtkUniformGrid *gridPtr = NULL;
+
+  switch( description )
+    {
+    case VTK_XY_PLANE:
+     // Root block
+     ndims[2]   = 1;
+     origin[0]  = origin[1]  = origin[2]  = 0.0;
+     spacing[0] = spacing[1] = spacing[2] = 1.0;
+     ndims[0]   = ndims[1]   = 4;
+     gridPtr    = GetGrid(origin,spacing,ndims);
+     amrDataSet->SetDataSet(0,0,gridPtr);
+     gridPtr->Delete();
+
+     // Refined patch that covers entire root domain
+     origin[0]  = origin[1]  = origin[2]  = 0.0;
+     spacing[0] = spacing[1] = spacing[2] = 0.5;
+     ndims[0]   = ndims[1]   = 6;
+     gridPtr    = GetGrid(origin,spacing,ndims);
+     amrDataSet->SetDataSet(1,0,gridPtr);
+     gridPtr->Delete();
+     break;
+    case VTK_XZ_PLANE:
+      // Root block
+      ndims[1]   = 1;
+      origin[0]  = origin[1]  = origin[2]  = 0.0;
+      spacing[0] = spacing[1] = spacing[2] = 1.0;
+      ndims[0]   = ndims[2]   = 4;
+      gridPtr    = GetGrid(origin,spacing,ndims);
+      amrDataSet->SetDataSet(0,0,gridPtr);
+      gridPtr->Delete();
+
+      // Refined patch that covers entire root domain
+      origin[0]  = origin[1]  = origin[2]  = 0.0;
+      spacing[0] = spacing[1] = spacing[2] = 0.5;
+      ndims[0]   = ndims[2]   = 6;
+      gridPtr    = GetGrid(origin,spacing,ndims);
+      amrDataSet->SetDataSet(1,0,gridPtr);
+      gridPtr->Delete();
+      break;
+    case VTK_YZ_PLANE:
+      // Root block
+      ndims[0]   = 1;
+      origin[0]  = origin[1]  = origin[2]  = 0.0;
+      spacing[0] = spacing[1] = spacing[2] = 1.0;
+      ndims[1]   = ndims[2]   = 4;
+      gridPtr    = GetGrid(origin,spacing,ndims);
+      amrDataSet->SetDataSet(0,0,gridPtr);
+      gridPtr->Delete();
+
+      // Refined patch that covers entire root domain
+      origin[0]  = origin[1]  = origin[2]  = 0.0;
+      spacing[0] = spacing[1] = spacing[2] = 0.5;
+      ndims[1]   = ndims[2]   = 6;
+      gridPtr    = GetGrid(origin,spacing,ndims);
+      amrDataSet->SetDataSet(1,0,gridPtr);
+      gridPtr->Delete();
+      break;
+    case VTK_XYZ_GRID:
+      // Root block
+      origin[0]  = origin[1]  = origin[2]  = 0.0;
+      spacing[0] = spacing[1] = spacing[2] = 1.0;
+      ndims[0]   = ndims[1]   = ndims[2]   = 4;
+      gridPtr    = GetGrid(origin,spacing,ndims);
+      amrDataSet->SetDataSet(0,0,gridPtr);
+      gridPtr->Delete();
+
+      // Refined patch that covers entire root domain
+      origin[0]  = origin[1]  = origin[2]  = 0.0;
+      spacing[0] = spacing[1] = spacing[2] = 0.5;
+      ndims[0]   = ndims[1]   = ndims[2]   = 6;
+      gridPtr    = GetGrid(origin,spacing,ndims);
+      amrDataSet->SetDataSet(1,0,gridPtr);
+      gridPtr->Delete();
+      break;
+    default:
+     assert(
+     "ERROR: Unhandled data description! Code should not reach here!" && false);
+    }
+  vtkAMRUtilities::GenerateMetaData(amrDataSet,NULL);
+  amrDataSet->GenerateVisibilityArrays();
+  return( amrDataSet );
+}
+
+//------------------------------------------------------------------------------
+int TestAMRVisibility( const int dataDescription )
+{
+  int rc = 0;
+  vtkOverlappingAMR *myAMR = GetAMRDataSet( dataDescription );
+  assert("pre: AMR dataset is NULL" && (myAMR != NULL) );
+
+  vtkUniformGrid *root = myAMR->GetDataSet(0,0);
+  assert("ERROR: root AMR block is NULL!" && (root != NULL) );
+
+  vtkIdType cellIdx = 0;
+  for(; cellIdx < root->GetNumberOfCells(); ++cellIdx )
+    {
+    if( root->IsCellVisible(cellIdx) )
+      {
+      rc++;
+      }
+    } // END for all cells
+
+  if( rc != 0 )
+    {
+    std::cerr << rc << "/" << root->GetNumberOfCells() << " are visible!\n";
+    }
+
+  myAMR->Delete();
+  return( rc );
+}
+
+//------------------------------------------------------------------------------
+void CheckTestStatus( int rc, std::string TestName )
+{
+  std::cout << "Test " << TestName << "...";
+  std::cout.flush();
+  if( rc == 0 )
+    {
+    std::cout << "PASSED!\n";
+    std::cout.flush();
+    }
+  else
+    {
+    std::cout << "FAILED!\n";
+    std::cout.flush();
+    }
+}
+
+} // END namespace AMRVisibilityTests
+
+//------------------------------------------------------------------------------
+int TestAMRBlanking(int, char *[])
+{
+  int rc     = 0;
+  int status = 0;
+
+  // XYZ test
+  status     = AMRVisibilityTests::TestAMRVisibility(VTK_XYZ_GRID);
+  AMRVisibilityTests::CheckTestStatus(status,"TestAMRVisibility-VTK_XYZ_GRID");
+  rc += status;
+
+  // XY PLANE
+  status = AMRVisibilityTests::TestAMRVisibility(VTK_XY_PLANE);
+  AMRVisibilityTests::CheckTestStatus(status,"TestAMRVisibility-VTK_XY_PLANE");
+  rc += status;
+
+  // XZ PLANE
+  status = AMRVisibilityTests::TestAMRVisibility(VTK_XZ_PLANE);
+  AMRVisibilityTests::CheckTestStatus(status,"TestAMRVisibility-VTK_XZ_PLANE");
+  rc += status;
+
+  // YZ PLANE
+  status = AMRVisibilityTests::TestAMRVisibility(VTK_YZ_PLANE);
+  AMRVisibilityTests::CheckTestStatus(status,"TestAMRVisibility-VTK_YZ_PLANE");
+  rc += status;
+  return( rc );
+}
+