Shuenhoy/vtk-m
1
0
Fork 0
mirror of https://gitlab.kitware.com/vtk/vtk-m synced 2024-09-08 21:33:55 +00:00
Code Issues 2 Actions Packages Projects Releases Wiki Activity

Adding fixes to identify degenerate cells

Browse Source
This commit is contained in:
Abhishek Yenpure 2022-01-22 21:06:28 -08:00
parent 4556eed085
commit 16f15c94b9
2 changed files with 6 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
vtkm/exec/CellMeasure.h
View File

@ -101,6 +101,12 @@ VTKM_EXEC OutType CellMeasure(const vtkm::IdComponent& numPts,
pts[0] - pts[3],
};
if (vtkm::MagnitudeSquared(edges[0]) == OutType(0.0) ||
vtkm::MagnitudeSquared(edges[1]) == OutType(0.0) ||
vtkm::MagnitudeSquared(edges[2]) == OutType(0.0) ||
vtkm::MagnitudeSquared(edges[3]) == OutType(0.0))
return OutType(0.0);
typename PointCoordVecType::ComponentType cornerNormals[4] = {
Cross(edges[3], edges[0]),
Cross(edges[0], edges[1]),

10
vtkm/worklet/MeshQuality.h
View File

@ -80,16 +80,6 @@ public:
const PointCoordVecType& pts,
OutType& metricValue) const
{
if (index == 2120)
{
std::cout << "Start Debug " << std::endl;
std::cout << "Num Points : " << numPoints << std::endl;
for (int i = 0; i < numPoints; i++)
{
std::cout << pts[i] << std::endl;
}
std::cout << "End Debug " << std::endl;
}
vtkm::UInt8 thisId = shape.Id;
if (shape.Id == vtkm::CELL_SHAPE_POLYGON)
{