Update probe filter

In the code path for probing on to uniform points, the code to compute the points
inside a cell has been made more accurate.
This commit is contained in:
Sujin Philip 2018-02-06 14:21:01 -05:00
parent e70e6ebb82
commit f7f08d3db2

@ -105,8 +105,10 @@ public:
// Compute points inside cell bounds
auto portal = points.GetPortal();
auto minp = static_cast<vtkm::Id3>((cbmin - portal.GetOrigin()) / portal.GetSpacing());
auto maxp = static_cast<vtkm::Id3>((cbmax - portal.GetOrigin()) / portal.GetSpacing());
auto minp =
static_cast<vtkm::Id3>(vtkm::Ceil((cbmin - portal.GetOrigin()) / portal.GetSpacing()));
auto maxp =
static_cast<vtkm::Id3>(vtkm::Floor((cbmax - portal.GetOrigin()) / portal.GetSpacing()));
// clamp
minp = vtkm::Max(minp, vtkm::Id3(0));