Shuenhoy/vtk-m
1
0
Fork 0
mirror of https://gitlab.kitware.com/vtk/vtk-m synced 2024-07-30 18:50:59 +00:00
Code Issues 2 Actions Packages Projects Releases Wiki Activity

Fix divide-by-zero in UnitTestCellInterpolate

Browse Source
This commit is contained in:
Kenneth Moreland 2021-07-12 07:08:33 -06:00
parent e4ae3cee1b
commit 9ce97352d9
1 changed files with 11 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
vtkm/exec/testing/UnitTestCellInterpolate.cxx
View File

@ -59,6 +59,11 @@ struct TestInterpolateFunctor
void DoTestWithField(CellShapeTag shape, const FieldVecType& fieldValues) const
{
vtkm::IdComponent numPoints = fieldValues.GetNumberOfComponents();
if (numPoints < 1)
{
return;
}
FieldType averageValue = vtkm::TypeTraits<FieldType>::ZeroInitialization();
for (vtkm::IdComponent pointIndex = 0; pointIndex < numPoints; pointIndex++)
{
@ -77,16 +82,13 @@ struct TestInterpolateFunctor
"Interpolation at point not point value.");
}
if (numPoints > 0)
{
vtkm::Vec3f pcoord;
CHECK_CALL(vtkm::exec::ParametricCoordinatesCenter(numPoints, shape, pcoord));
FieldType interpolatedValue;
CHECK_CALL(vtkm::exec::CellInterpolate(fieldValues, pcoord, shape, interpolatedValue));
vtkm::Vec3f pcoord;
CHECK_CALL(vtkm::exec::ParametricCoordinatesCenter(numPoints, shape, pcoord));
FieldType interpolatedValue;
CHECK_CALL(vtkm::exec::CellInterpolate(fieldValues, pcoord, shape, interpolatedValue));
VTKM_TEST_ASSERT(test_equal(averageValue, interpolatedValue),
"Interpolation at center not average value.");
}
VTKM_TEST_ASSERT(test_equal(averageValue, interpolatedValue),
"Interpolation at center not average value.");
}
template <typename CellShapeTag>