exec/CellLocatorUniformBins.h

Go to the documentation of this file.

//============================================================================
//  The contents of this file are covered by the Viskores license. See
//  LICENSE.txt for details.
//
//  By contributing to this file, all contributors agree to the Developer
//  Certificate of Origin Version 1.1 (DCO 1.1) as stated in DCO.txt.
//============================================================================
 
//============================================================================
//  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.
//============================================================================
#ifndef viskores_exec_CellLocatorUniformBins_h
#define viskores_exec_CellLocatorUniformBins_h
 
#include <viskores/exec/CellInside.h>
#include <viskores/exec/ParametricCoordinates.h>
 
#include <viskores/cont/ArrayHandle.h>
#include <viskores/cont/ArrayHandleGroupVecVariable.h>
#include <viskores/cont/CoordinateSystem.h>
 
#include <viskores/TopologyElementTag.h>
#include <viskores/VecFromPortalPermute.h>
#include <viskores/VecTraits.h>
 
namespace viskores
{
namespace exec
{
 
//--------------------------------------------------------------------
 
template <typename CellStructureType>
class VISKORES_ALWAYS_EXPORT CellLocatorUniformBins
{
  template <typename T>
  using ReadPortal = typename viskores::cont::ArrayHandle<T>::ReadPortalType;
 
  using CoordsPortalType =
    typename viskores::cont::CoordinateSystem::MultiplexerArrayType::ReadPortalType;
 
  using CellIdArrayType = viskores::cont::ArrayHandle<viskores::Id>;
  using CellIdOffsetArrayType = viskores::cont::ArrayHandle<viskores::Id>;
  using CellIdReadPortal =
    typename viskores::cont::ArrayHandleGroupVecVariable<CellIdArrayType,
                                                         CellIdOffsetArrayType>::ReadPortalType;
 
public:
  template <typename CellSetType>
  VISKORES_CONT CellLocatorUniformBins(
    const viskores::Id3& cellDims,
    const viskores::Vec3f& origin,
    const viskores::Vec3f& maxPoint,
    const viskores::Vec3f& invSpacing,
    const viskores::Id3& maxCellIds,
    const viskores::cont::ArrayHandleGroupVecVariable<CellIdArrayType, CellIdOffsetArrayType>&
      cellIds,
    const CellSetType& cellSet,
    const viskores::cont::CoordinateSystem& coords,
    viskores::cont::DeviceAdapterId device,
    viskores::cont::Token& token)
    : CellDims(cellDims)
    , Origin(origin)
    , MaxPoint(maxPoint)
    , InvSpacing(invSpacing)
    , MaxCellIds(maxCellIds)
    , CellIds(cellIds.PrepareForInput(device, token))
    , CellSet(cellSet.PrepareForInput(device,
                                      viskores::TopologyElementTagCell{},
                                      viskores::TopologyElementTagPoint{},
                                      token))
    , Coords(coords.GetDataAsMultiplexer().PrepareForInput(device, token))
  {
  }
 
  struct LastCell
  {
    viskores::Id CellId = -1;
    viskores::Id BinIdx = -1;
  };
 
  VISKORES_EXEC viskores::ErrorCode FindCell(const viskores::Vec3f& point,
                                             viskores::Id& cellId,
                                             viskores::Vec3f& parametric) const
  {
    LastCell lastCell;
    return this->FindCellImpl(point, cellId, parametric, lastCell);
  }
 
  VISKORES_EXEC viskores::ErrorCode FindCell(const viskores::Vec3f& point,
                                             viskores::Id& cellId,
                                             viskores::Vec3f& parametric,
                                             LastCell& lastCell) const
  {
    viskores::Id binIdx = this->FindBinIdx(point);
 
    if (binIdx == -1)
    {
      lastCell.CellId = -1;
      lastCell.BinIdx = -1;
      cellId = -1;
      return viskores::ErrorCode::CellNotFound;
    }
    //See if the point is still in the same bin.
    else if (binIdx == lastCell.BinIdx && this->LastCellValid(lastCell))
    {
      viskores::Vec3f pc;
      //Check the last cell first.
      if (this->PointInCell(point, lastCell.CellId, pc))
      {
        parametric = pc;
        cellId = lastCell.CellId;
        return viskores::ErrorCode::Success;
      }
      //Otherwise, check cells in the bin, but skip lastCell.CellId
      else if (this->PointInBin(point, lastCell.BinIdx, cellId, pc, lastCell.CellId))
      {
        parametric = pc;
        return viskores::ErrorCode::Success;
      }
    }
    //if cell still not found, drop to the general find cell below.
 
    //LastCell not initialized, or not in the same bin: do a full test.
    //Since already computed the binIdx, re-use it.
    return this->FindCellImpl(point, cellId, parametric, lastCell, binIdx);
  }
 
  VISKORES_DEPRECATED(1.6, "Locators are no longer pointers. Use . operator.")
  VISKORES_EXEC CellLocatorUniformBins* operator->() { return this; }
  VISKORES_DEPRECATED(1.6, "Locators are no longer pointers. Use . operator.")
  VISKORES_EXEC const CellLocatorUniformBins* operator->() const { return this; }
 
private:
  VISKORES_EXEC viskores::Id FindBinIdx(const viskores::Vec3f& point) const
  {
    if (!this->IsInside(point))
      return -1;
 
    viskores::Vec3f temp;
    temp = point - this->Origin;
    temp = temp * this->InvSpacing;
 
    //make sure that if we border the upper edge, we sample the correct cell
    viskores::Id3 logicalCell = viskores::Min(viskores::Id3(temp), this->MaxCellIds);
 
    viskores::Id binIdx =
      (logicalCell[2] * this->CellDims[1] + logicalCell[1]) * this->CellDims[0] + logicalCell[0];
 
    return binIdx;
  }
 
  VISKORES_EXEC bool LastCellValid(const LastCell& lastCell) const
  {
    return lastCell.BinIdx >= 0 && lastCell.BinIdx < this->CellIds.GetNumberOfValues() &&
      lastCell.CellId >= 0 && lastCell.CellId < this->CellSet.GetNumberOfElements();
  }
 
  VISKORES_EXEC bool IsInside(const viskores::Vec3f& point) const
  {
    if (point[0] < this->Origin[0] || point[0] > this->MaxPoint[0])
      return false;
    if (point[1] < this->Origin[1] || point[1] > this->MaxPoint[1])
      return false;
    if (point[2] < this->Origin[2] || point[2] > this->MaxPoint[2])
      return false;
 
    return true;
  }
 
  VISKORES_EXEC
  viskores::ErrorCode FindCellImpl(const viskores::Vec3f& point,
                                   viskores::Id& cellId,
                                   viskores::Vec3f& parametric,
                                   LastCell& lastCell,
                                   viskores::Id ptBinIdx = -1) const
  {
    lastCell.CellId = -1;
    lastCell.BinIdx = -1;
 
    //if ptBinIdx is set, use it. Otherwise, compute the bin idx.
    viskores::Id binIdx = -1;
    if (ptBinIdx == -1)
      binIdx = this->FindBinIdx(point);
    else
      binIdx = ptBinIdx;
 
    //point not in a bin. return not found.
    if (binIdx == -1)
    {
      cellId = -1;
      return viskores::ErrorCode::CellNotFound;
    }
 
    //search cells in the bin.
    viskores::Vec3f pc;
    if (this->PointInBin(point, binIdx, cellId, pc))
    {
      parametric = pc;
      lastCell.CellId = cellId;
      lastCell.BinIdx = binIdx;
      return viskores::ErrorCode::Success;
    }
 
    return viskores::ErrorCode::CellNotFound;
  }
 
  template <typename PointsVecType>
  VISKORES_EXEC viskores::Bounds ComputeCellBounds(const PointsVecType& points) const
  {
    auto numPoints = viskores::VecTraits<PointsVecType>::GetNumberOfComponents(points);
 
    viskores::Bounds bounds;
    for (viskores::IdComponent i = 0; i < numPoints; ++i)
      bounds.Include(points[i]);
 
    return bounds;
  }
 
  // TODO: This function may return false positives for non 3D cells as the
  // tests are done on the projection of the point on the cell. Extra checks
  // should be added to test if the point actually falls on the cell.
  template <typename CellShapeTag, typename CoordsType>
  VISKORES_EXEC viskores::ErrorCode PointInsideCell(viskores::Vec3f point,
                                                    CellShapeTag cellShape,
                                                    CoordsType cellPoints,
                                                    viskores::Vec3f& parametricCoordinates,
                                                    bool& inside) const
  {
    auto bounds = this->ComputeCellBounds(cellPoints);
    if (bounds.Contains(point))
    {
      VISKORES_RETURN_ON_ERROR(viskores::exec::WorldCoordinatesToParametricCoordinates(
        cellPoints, point, cellShape, parametricCoordinates));
      inside = viskores::exec::CellInside(parametricCoordinates, cellShape);
    }
    else
    {
      inside = false;
    }
    // Return success error code even point is not inside this cell
    return viskores::ErrorCode::Success;
  }
 
  VISKORES_EXEC
  bool PointInBin(const viskores::Vec3f& point,
                  const viskores::Id& binIdx,
                  viskores::Id& cellId,
                  viskores::Vec3f& parametric,
                  const viskores::Id& skipCellId = -1) const
  {
    auto binIds = this->CellIds.Get(binIdx);
 
    viskores::Vec3f pc;
    for (viskores::IdComponent i = 0; i < binIds.GetNumberOfComponents(); i++)
    {
      viskores::Id cid = binIds[i];
      if (cid != skipCellId && this->PointInCell(point, cid, pc))
      {
        cellId = cid;
        parametric = pc;
        return true;
      }
    }
 
    return false;
  }
 
  VISKORES_EXEC
  bool PointInCell(const viskores::Vec3f& point,
                   const viskores::Id& cid,
                   viskores::Vec3f& parametric) const
  {
    auto indices = this->CellSet.GetIndices(cid);
    auto pts = viskores::make_VecFromPortalPermute(&indices, this->Coords);
    viskores::Vec3f pc;
    bool inside;
    auto status = this->PointInsideCell(point, this->CellSet.GetCellShape(cid), pts, pc, inside);
    if (status == viskores::ErrorCode::Success && inside)
    {
      parametric = pc;
      return true;
    }
 
    return false;
  }
 
  viskores::Id3 CellDims;
  viskores::Vec3f Origin;
  viskores::Vec3f MaxPoint;
  viskores::Vec3f InvSpacing;
  viskores::Id3 MaxCellIds;
 
  CellIdReadPortal CellIds;
 
  CellStructureType CellSet;
  CoordsPortalType Coords;
};
 
}
} // viskores::exec
 
#endif //viskores_exec_CellLocatorUniformBins_h