Skip to content

Fix an issue with structured (non-uniform) gradients

Kenneth Moreland requested to merge kmorel/vtk-m:structured-gradient-fix into master

This fixes a bug introduced in commit 93d87e06

Part of the commit moved fromusing a BoundaryState object to keeping boundary flags in a vtkm::Vec<bool,3>. However, the meaning of the flag got inverted. The BoundaryState::IsRadiusIn?Boundary methods returned true if the current location was inside the bounds (i.e. not on the boundary), but the new flags are true if the current location is on the boundary (the inverse).

The updating of the code for the StructuredPointGradient::Jacobian method got the flag somewhat wrong and scaled the wrong case.

Merge request reports