//============================================================================
//  Copyright (c) Kitware, Inc.
//  All rights reserved.
//  See LICENSE.txt 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.
//
//  Copyright 2014 Sandia Corporation.
//  Copyright 2014 UT-Battelle, LLC.
//  Copyright 2014 Los Alamos National Security.
//
//  Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
//  the U.S. Government retains certain rights in this software.
//
//  Under the terms of Contract DE-AC52-06NA25396 with Los Alamos National
//  Laboratory (LANL), the U.S. Government retains certain rights in
//  this software.
//============================================================================

#include <vtkm/cont/testing/MakeTestDataSet.h>
#include <vtkm/cont/testing/Testing.h>

#include <vtkm/filter/ExternalFaces.h>

using vtkm::cont::testing::MakeTestDataSet;

namespace {
void TestExternalFacesExplicitGrid()
{
  vtkm::cont::DataSet ds;

  const int nVerts = 8; //A cube that is tetrahedralized
  typedef vtkm::Vec<vtkm::Float32,3> CoordType;
  CoordType coordinates[nVerts] = {
    CoordType(0, 0, 0),
    CoordType(1, 0, 0),
    CoordType(1, 1, 0),
    CoordType(0, 1, 0),
    CoordType(0, 0, 1),
    CoordType(1, 0, 1),
    CoordType(1, 1, 1),
    CoordType(0, 1, 1)
    };

  ds.AddCoordinateSystem(
    vtkm::cont::CoordinateSystem("coordinates", coordinates, nVerts));

  //Construct the VTK-m shapes and numIndices connectivity arrays
  const int nCells = 6;  //The tetrahedrons of the cube
  vtkm::IdComponent cellVerts[nCells][4] = {
                                            {4,7,6,3}, {4,6,3,2}, {4,0,3,2},
                                            {4,6,5,2}, {4,5,0,2}, {1,0,5,2}
                                           };
  vtkm::cont::CellSetExplicit<> cellSet("cells");

  vtkm::cont::ArrayHandle<vtkm::UInt8>       shapes;
  vtkm::cont::ArrayHandle<vtkm::IdComponent> numIndices;
  vtkm::cont::ArrayHandle<vtkm::Id>          conn;
  shapes.Allocate(static_cast<vtkm::Id>(nCells));
  numIndices.Allocate(static_cast<vtkm::Id>(nCells));
  conn.Allocate(static_cast<vtkm::Id>(4 * nCells));

  int index = 0;
  for(int j = 0; j < nCells; j++)
  {
    shapes.GetPortalControl().Set(j, static_cast<vtkm::UInt8>(vtkm::CELL_SHAPE_TETRA));
    numIndices.GetPortalControl().Set(j, 4);
    for(int k = 0; k < 4; k++)
      conn.GetPortalControl().Set(index++, cellVerts[j][k]);
  }

  cellSet.Fill(nVerts, shapes, numIndices, conn);

  //Add the VTK-m cell set
  ds.AddCellSet(cellSet);

  //Run the External Faces filter
  vtkm::filter::ExternalFaces externalFaces;
  vtkm::filter::ResultDataSet result = externalFaces.Execute(ds);

  //Validate the number of external faces (output) returned by the worklet
  VTKM_TEST_ASSERT(result.IsValid(), "Results should be valid");

  vtkm::cont::CellSetExplicit<> &new_cellSet =
      result.GetDataSet().GetCellSet(0).Cast<vtkm::cont::CellSetExplicit<> >();

  const vtkm::Id numExtFaces_out = new_cellSet.GetNumberOfCells();
  const vtkm::Id numExtFaces_actual = 12;
  VTKM_TEST_ASSERT(numExtFaces_out == numExtFaces_actual, "Number of External Faces mismatch");
}

void TestExternalFacesFilter()
{
  TestExternalFacesExplicitGrid();
}

} // anonymous namespace


int UnitTestExternalFacesFilter(int, char *[])
{
  return vtkm::cont::testing::Testing::Run(TestExternalFacesFilter);
}
