Failure in cell locator with explicit grids
I'm having trouble with the cell locator for explicit grids. I've created this branch that illustrates the bug: https://gitlab.kitware.com/dpugmire/vtk-m/tree/cellLocatorBug
I get the following output when running UnitTestCellLocator
Any thoughts on this?
./bin/UnitTests_vtkm_cont_testing UnitTestCellLocator Explicit Dataset DataSet: CoordSystems[1] Coordinate System coordinates assoc= Points valueType=N4vtkm3VecIfLi3EEE storageType=N4vtkm4cont17StorageTagVirtualE numValues=64 bytes=768 [(0,0,0) (1,0,0) (2,0,0) ... (1,3,3) (2,3,3) (3,3,3)] CellSets[1] ExplicitCellSet: cells PointToCell: Shapes: valueType=h storageType=N4vtkm4cont15StorageTagBasicE numValues=27 bytes=27 [12 12 12 ... 12 12 12] NumIndices: valueType=i storageType=N4vtkm4cont15StorageTagBasicE numValues=27 bytes=108 [3 3 3 ... 3 3 3] Connectivity: valueType=x storageType=N4vtkm4cont15StorageTagBasicE numValues=216 bytes=1728 [0 1 4 ... 0 0 0] IndexOffsets: Not Allocated CellToPoint: Not Allocated Fields[0] UnitTests_vtkm_cont_testing: /mnt/disk2TB/proj/vtkm/cellLocatorBug/vtk-m/vtkm/exec/CellInterpolate.h:376: typename FieldVecType::ComponentType vtkm::exec::CellInterpolate(const FieldVecType&, const vtkm::Vec<ParametricCoordType, 3>&, vtkm::CellShapeTagHexahedron, const vtkm::exec::FunctorBase&) [with FieldVecType = vtkm::VecVariable<vtkm::Vec<float, 3>, 8>; ParametricCoordType = float; typename FieldVecType::ComponentType = vtkm::Vec<float, 3>]: Assertion `field.GetNumberOfComponents() == 8' failed. zsh: abort (core dumped) ./bin/UnitTests_vtkm_cont_testing UnitTestCellLocator