2017-07-11 12:33:56 +00:00
|
|
|
//============================================================================
|
|
|
|
// Copyright (c) Kitware, Inc.
|
|
|
|
// All rights reserved.
|
|
|
|
// See LICENSE.txt for details.
|
2019-04-15 23:24:21 +00:00
|
|
|
//
|
2017-07-11 12:33:56 +00:00
|
|
|
// 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.
|
|
|
|
//============================================================================
|
|
|
|
|
|
|
|
#include <vtkm/worklet/DispatcherPointNeighborhood.h>
|
|
|
|
#include <vtkm/worklet/WorkletPointNeighborhood.h>
|
|
|
|
|
|
|
|
#include <vtkm/worklet/ScatterIdentity.h>
|
|
|
|
#include <vtkm/worklet/ScatterUniform.h>
|
|
|
|
|
|
|
|
#include <vtkm/Math.h>
|
|
|
|
#include <vtkm/VecAxisAlignedPointCoordinates.h>
|
|
|
|
|
|
|
|
#include <vtkm/cont/DataSet.h>
|
2019-04-03 20:28:31 +00:00
|
|
|
#include <vtkm/cont/DeviceAdapterTag.h>
|
2017-07-11 12:33:56 +00:00
|
|
|
|
|
|
|
#include <vtkm/cont/testing/MakeTestDataSet.h>
|
|
|
|
#include <vtkm/cont/testing/Testing.h>
|
|
|
|
|
|
|
|
namespace test_pointneighborhood
|
|
|
|
{
|
|
|
|
|
2018-11-18 17:41:13 +00:00
|
|
|
struct MaxNeighborValue : public vtkm::worklet::WorkletPointNeighborhood
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
|
|
|
|
2019-01-10 18:59:25 +00:00
|
|
|
using ControlSignature = void(FieldInNeighborhood neighbors, CellSetIn, FieldOut maxV);
|
2017-07-11 12:33:56 +00:00
|
|
|
|
2018-12-03 19:44:17 +00:00
|
|
|
using ExecutionSignature = void(Boundary, _1, _3);
|
2017-07-11 12:33:56 +00:00
|
|
|
//verify input domain can be something other than first parameter
|
2018-02-22 13:29:13 +00:00
|
|
|
using InputDomain = _2;
|
2017-07-11 12:33:56 +00:00
|
|
|
|
|
|
|
template <typename FieldIn, typename FieldOut>
|
2018-12-03 19:44:17 +00:00
|
|
|
VTKM_EXEC void operator()(const vtkm::exec::BoundaryState& boundary,
|
|
|
|
const vtkm::exec::FieldNeighborhood<FieldIn>& inputField,
|
2017-07-11 12:33:56 +00:00
|
|
|
FieldOut& output) const
|
|
|
|
{
|
|
|
|
using ValueType = typename FieldIn::ValueType;
|
|
|
|
|
|
|
|
auto* nboundary = inputField.Boundary;
|
|
|
|
|
2019-07-19 14:42:29 +00:00
|
|
|
if (!(nboundary->IsRadiusInXBoundary(1) == boundary.IsRadiusInXBoundary(1)))
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
|
|
|
this->RaiseError("Got invalid XPos boundary state");
|
|
|
|
}
|
|
|
|
|
2019-07-19 14:42:29 +00:00
|
|
|
if (!(nboundary->IsRadiusInYBoundary(1) == boundary.IsRadiusInYBoundary(1)))
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
|
|
|
this->RaiseError("Got invalid YPos boundary state");
|
|
|
|
}
|
|
|
|
|
2019-07-19 14:42:29 +00:00
|
|
|
if (!(nboundary->IsRadiusInZBoundary(1) == boundary.IsRadiusInZBoundary(1)))
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
|
|
|
this->RaiseError("Got invalid ZPos boundary state");
|
|
|
|
}
|
|
|
|
|
2019-07-19 14:42:29 +00:00
|
|
|
if (!(nboundary->IsRadiusInBoundary(1) == boundary.IsRadiusInBoundary(1)))
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
2018-11-18 17:41:13 +00:00
|
|
|
this->RaiseError("Got invalid boundary state");
|
2017-07-11 12:33:56 +00:00
|
|
|
}
|
|
|
|
|
2019-07-19 14:42:29 +00:00
|
|
|
if (nboundary->IsRadiusInXBoundary(1) !=
|
|
|
|
(boundary.IsNeighborInXBoundary(-1) && boundary.IsNeighborInXBoundary(1)))
|
|
|
|
{
|
|
|
|
this->RaiseError("Neighbor/Radius boundary mismatch in X dimension.");
|
|
|
|
}
|
|
|
|
|
|
|
|
if (nboundary->IsRadiusInYBoundary(1) !=
|
|
|
|
(boundary.IsNeighborInYBoundary(-1) && boundary.IsNeighborInYBoundary(1)))
|
|
|
|
{
|
|
|
|
this->RaiseError("Neighbor/Radius boundary mismatch in Y dimension.");
|
|
|
|
}
|
|
|
|
|
|
|
|
if (nboundary->IsRadiusInZBoundary(1) !=
|
|
|
|
(boundary.IsNeighborInZBoundary(-1) && boundary.IsNeighborInZBoundary(1)))
|
|
|
|
{
|
|
|
|
this->RaiseError("Neighbor/Radius boundary mismatch in Z dimension.");
|
|
|
|
}
|
|
|
|
|
|
|
|
if (nboundary->IsRadiusInBoundary(1) !=
|
|
|
|
(boundary.IsNeighborInBoundary({ -1 }) && boundary.IsNeighborInBoundary({ 1 })))
|
|
|
|
{
|
|
|
|
this->RaiseError("Neighbor/Radius boundary mismatch.");
|
|
|
|
}
|
|
|
|
|
2017-07-11 12:33:56 +00:00
|
|
|
|
2018-11-18 17:41:13 +00:00
|
|
|
auto minNeighbors = boundary.MinNeighborIndices(1);
|
|
|
|
auto maxNeighbors = boundary.MaxNeighborIndices(1);
|
2017-07-11 12:33:56 +00:00
|
|
|
|
|
|
|
ValueType maxV = inputField.Get(0, 0, 0); //our value
|
2018-11-18 17:41:13 +00:00
|
|
|
for (vtkm::IdComponent k = minNeighbors[2]; k <= maxNeighbors[2]; ++k)
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
2018-11-18 17:41:13 +00:00
|
|
|
for (vtkm::IdComponent j = minNeighbors[1]; j <= maxNeighbors[1]; ++j)
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
2018-11-18 17:41:13 +00:00
|
|
|
for (vtkm::IdComponent i = minNeighbors[0]; i <= maxNeighbors[0]; ++i)
|
|
|
|
{
|
|
|
|
maxV = vtkm::Max(maxV, inputField.Get(i, j, k));
|
|
|
|
}
|
2017-07-11 12:33:56 +00:00
|
|
|
}
|
|
|
|
}
|
|
|
|
output = static_cast<FieldOut>(maxV);
|
|
|
|
}
|
|
|
|
};
|
|
|
|
|
2018-11-18 17:41:13 +00:00
|
|
|
struct ScatterIdentityNeighbor : public vtkm::worklet::WorkletPointNeighborhood
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
2019-01-10 18:59:25 +00:00
|
|
|
using ControlSignature = void(CellSetIn topology, FieldIn pointCoords);
|
2018-05-25 21:18:41 +00:00
|
|
|
using ExecutionSignature =
|
|
|
|
void(_2, WorkIndex, InputIndex, OutputIndex, ThreadIndices, VisitIndex);
|
2017-07-11 12:33:56 +00:00
|
|
|
|
|
|
|
VTKM_CONT
|
|
|
|
ScatterIdentityNeighbor() {}
|
|
|
|
|
|
|
|
template <typename T>
|
|
|
|
VTKM_EXEC void operator()(
|
|
|
|
const vtkm::Vec<T, 3>& vtkmNotUsed(coords),
|
|
|
|
const vtkm::Id& workIndex,
|
|
|
|
const vtkm::Id& inputIndex,
|
|
|
|
const vtkm::Id& outputIndex,
|
2018-11-18 17:41:13 +00:00
|
|
|
const vtkm::exec::arg::ThreadIndicesPointNeighborhood& vtkmNotUsed(threadIndices),
|
2017-07-11 12:33:56 +00:00
|
|
|
const vtkm::Id& visitIndex) const
|
|
|
|
{
|
|
|
|
if (workIndex != inputIndex)
|
|
|
|
{
|
|
|
|
this->RaiseError("Got wrong input value.");
|
|
|
|
}
|
|
|
|
if (outputIndex != workIndex)
|
|
|
|
{
|
|
|
|
this->RaiseError("Got work and output index don't match.");
|
|
|
|
}
|
|
|
|
if (visitIndex != 0)
|
|
|
|
{
|
|
|
|
this->RaiseError("Got wrong visit value1.");
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
using ScatterType = vtkm::worklet::ScatterIdentity;
|
|
|
|
};
|
|
|
|
|
2018-11-18 17:41:13 +00:00
|
|
|
struct ScatterUniformNeighbor : public vtkm::worklet::WorkletPointNeighborhood
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
2019-01-10 18:59:25 +00:00
|
|
|
using ControlSignature = void(CellSetIn topology, FieldIn pointCoords);
|
2018-05-25 21:18:41 +00:00
|
|
|
using ExecutionSignature =
|
|
|
|
void(_2, WorkIndex, InputIndex, OutputIndex, ThreadIndices, VisitIndex);
|
2017-07-11 12:33:56 +00:00
|
|
|
|
|
|
|
VTKM_CONT
|
|
|
|
ScatterUniformNeighbor() {}
|
|
|
|
|
|
|
|
template <typename T>
|
|
|
|
VTKM_EXEC void operator()(
|
|
|
|
const vtkm::Vec<T, 3>& vtkmNotUsed(coords),
|
|
|
|
const vtkm::Id& workIndex,
|
|
|
|
const vtkm::Id& inputIndex,
|
|
|
|
const vtkm::Id& outputIndex,
|
2018-11-18 17:41:13 +00:00
|
|
|
const vtkm::exec::arg::ThreadIndicesPointNeighborhood& vtkmNotUsed(threadIndices),
|
2017-07-11 12:33:56 +00:00
|
|
|
const vtkm::Id& visitIndex) const
|
|
|
|
{
|
|
|
|
if ((workIndex / 3) != inputIndex)
|
|
|
|
{
|
|
|
|
this->RaiseError("Got wrong input value.");
|
|
|
|
}
|
|
|
|
if (outputIndex != workIndex)
|
|
|
|
{
|
|
|
|
this->RaiseError("Got work and output index don't match.");
|
|
|
|
}
|
|
|
|
if ((workIndex % 3) != visitIndex)
|
|
|
|
{
|
|
|
|
this->RaiseError("Got wrong visit value2.");
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
|
Move Scatter from Worklet to Dispatcher
Previously, when a Worklet needed a scatter, the scatter object was
stored in the Worklet object. That was problematic because that means
the Scatter, which is a control object, was shoved into the execution
environment.
To prevent that, move the Scatter into the Dispatcher object. The
worklet still declares a ScatterType alias, but no longer has a
GetScatter method. Instead, the Dispatcher now takes a Scatter object in
its constructor. If using the default scatter (ScatterIdentity), the
default constructor is used. If using another type of Scatter that
requires data to set up its state, then the caller of the worklet needs
to provide that to the dispatcher. For convenience, worklets are
encouraged to have a MakeScatter method to help construct a proper
scatter object.
2018-04-27 04:43:51 +00:00
|
|
|
using ScatterType = vtkm::worklet::ScatterUniform<3>;
|
2017-07-11 12:33:56 +00:00
|
|
|
};
|
|
|
|
}
|
|
|
|
|
|
|
|
namespace
|
|
|
|
{
|
|
|
|
|
|
|
|
static void TestMaxNeighborValue();
|
|
|
|
static void TestScatterIdentityNeighbor();
|
|
|
|
static void TestScatterUnfiormNeighbor();
|
|
|
|
|
2018-10-17 19:17:29 +00:00
|
|
|
void TestWorkletPointNeighborhood(vtkm::cont::DeviceAdapterId id)
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
2018-10-17 19:17:29 +00:00
|
|
|
std::cout << "Testing Point Neighborhood Worklet on device adapter: " << id.GetName()
|
|
|
|
<< std::endl;
|
2017-07-11 12:33:56 +00:00
|
|
|
|
|
|
|
TestMaxNeighborValue();
|
|
|
|
TestScatterIdentityNeighbor();
|
|
|
|
TestScatterUnfiormNeighbor();
|
|
|
|
}
|
|
|
|
|
|
|
|
static void TestMaxNeighborValue()
|
|
|
|
{
|
|
|
|
std::cout << "Testing MaxPointOfCell worklet" << std::endl;
|
|
|
|
|
|
|
|
|
|
|
|
vtkm::cont::testing::MakeTestDataSet testDataSet;
|
|
|
|
|
|
|
|
vtkm::worklet::DispatcherPointNeighborhood<::test_pointneighborhood::MaxNeighborValue> dispatcher;
|
|
|
|
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::Float32> output;
|
|
|
|
|
|
|
|
vtkm::cont::DataSet dataSet3D = testDataSet.Make3DUniformDataSet0();
|
2021-01-25 17:36:56 +00:00
|
|
|
dispatcher.Invoke(dataSet3D.GetField("pointvar")
|
|
|
|
.GetData()
|
|
|
|
.ResetTypes<vtkm::TypeListFieldScalar, VTKM_DEFAULT_STORAGE_LIST>(),
|
|
|
|
dataSet3D.GetCellSet(),
|
|
|
|
output);
|
2017-07-11 12:33:56 +00:00
|
|
|
|
|
|
|
vtkm::Float32 expected3D[18] = { 110.3f, 120.3f, 120.3f, 110.3f, 120.3f, 120.3f,
|
|
|
|
170.5f, 180.5f, 180.5f, 170.5f, 180.5f, 180.5f,
|
|
|
|
170.5f, 180.5f, 180.5f, 170.5f, 180.5f, 180.5f };
|
|
|
|
for (int i = 0; i < 18; ++i)
|
|
|
|
{
|
2020-01-28 19:14:32 +00:00
|
|
|
VTKM_TEST_ASSERT(test_equal(output.ReadPortal().Get(i), expected3D[i]),
|
2017-07-11 12:33:56 +00:00
|
|
|
"Wrong result for MaxNeighborValue worklet");
|
|
|
|
}
|
|
|
|
|
|
|
|
vtkm::cont::DataSet dataSet2D = testDataSet.Make2DUniformDataSet1();
|
2021-01-25 17:36:56 +00:00
|
|
|
dispatcher.Invoke(dataSet2D.GetField("pointvar")
|
|
|
|
.GetData()
|
|
|
|
.ResetTypes<vtkm::TypeListFieldScalar, VTKM_DEFAULT_STORAGE_LIST>(),
|
|
|
|
dataSet2D.GetCellSet(),
|
|
|
|
output);
|
2017-07-11 12:33:56 +00:00
|
|
|
|
|
|
|
vtkm::Float32 expected2D[25] = { 100.0f, 100.0f, 78.0f, 49.0f, 33.0f, 100.0f, 100.0f,
|
|
|
|
78.0f, 50.0f, 48.0f, 94.0f, 94.0f, 91.0f, 91.0f,
|
|
|
|
91.0f, 52.0f, 52.0f, 91.0f, 91.0f, 91.0f, 12.0f,
|
|
|
|
51.0f, 91.0f, 91.0f, 91.0f };
|
|
|
|
|
|
|
|
for (int i = 0; i < 25; ++i)
|
|
|
|
{
|
2020-01-28 19:14:32 +00:00
|
|
|
VTKM_TEST_ASSERT(test_equal(output.ReadPortal().Get(i), expected2D[i]),
|
2017-07-11 12:33:56 +00:00
|
|
|
"Wrong result for MaxNeighborValue worklet");
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
static void TestScatterIdentityNeighbor()
|
|
|
|
{
|
|
|
|
std::cout << "Testing identity scatter with PointNeighborhood" << std::endl;
|
|
|
|
|
|
|
|
vtkm::cont::testing::MakeTestDataSet testDataSet;
|
|
|
|
|
|
|
|
vtkm::worklet::DispatcherPointNeighborhood<::test_pointneighborhood::ScatterIdentityNeighbor>
|
|
|
|
dispatcher;
|
|
|
|
|
|
|
|
vtkm::cont::DataSet dataSet3D = testDataSet.Make3DUniformDataSet0();
|
|
|
|
dispatcher.Invoke(dataSet3D.GetCellSet(), dataSet3D.GetCoordinateSystem());
|
|
|
|
|
|
|
|
vtkm::cont::DataSet dataSet2D = testDataSet.Make2DUniformDataSet0();
|
|
|
|
dispatcher.Invoke(dataSet2D.GetCellSet(), dataSet2D.GetCoordinateSystem());
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
static void TestScatterUnfiormNeighbor()
|
|
|
|
{
|
|
|
|
std::cout << "Testing uniform scatter with PointNeighborhood" << std::endl;
|
|
|
|
|
|
|
|
vtkm::cont::testing::MakeTestDataSet testDataSet;
|
|
|
|
|
|
|
|
vtkm::worklet::DispatcherPointNeighborhood<::test_pointneighborhood::ScatterUniformNeighbor>
|
|
|
|
dispatcher;
|
|
|
|
|
|
|
|
vtkm::cont::DataSet dataSet3D = testDataSet.Make3DUniformDataSet0();
|
|
|
|
dispatcher.Invoke(dataSet3D.GetCellSet(), dataSet3D.GetCoordinateSystem());
|
|
|
|
|
|
|
|
vtkm::cont::DataSet dataSet2D = testDataSet.Make2DUniformDataSet0();
|
|
|
|
dispatcher.Invoke(dataSet2D.GetCellSet(), dataSet2D.GetCoordinateSystem());
|
|
|
|
}
|
|
|
|
|
|
|
|
} // anonymous namespace
|
|
|
|
|
2018-10-17 19:17:29 +00:00
|
|
|
int UnitTestWorkletMapPointNeighborhood(int argc, char* argv[])
|
2017-07-11 12:33:56 +00:00
|
|
|
{
|
2018-10-17 19:17:29 +00:00
|
|
|
return vtkm::cont::testing::Testing::RunOnDevice(TestWorkletPointNeighborhood, argc, argv);
|
2017-07-11 12:33:56 +00:00
|
|
|
}
|