VTK  9.3.1
vtkHigherOrderTriangle.h
Go to the documentation of this file.
1 // SPDX-FileCopyrightText: Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
2 // SPDX-License-Identifier: BSD-3-Clause
19 #ifndef vtkHigherOrderTriangle_h
20 #define vtkHigherOrderTriangle_h
21 
22 #include <functional> //For std::function
23 
24 #include "vtkCommonDataModelModule.h" // For export macro
25 #include "vtkNew.h" // For member variable.
26 #include "vtkNonLinearCell.h"
27 #include "vtkSmartPointer.h" // For member variable.
28 
29 #include <vector> // For caching
30 
31 VTK_ABI_NAMESPACE_BEGIN
32 class vtkDoubleArray;
34 class vtkTriangle;
35 
36 class VTKCOMMONDATAMODEL_EXPORT vtkHigherOrderTriangle : public vtkNonLinearCell
37 {
38 public:
40  void PrintSelf(ostream& os, vtkIndent indent) override;
41 
42  int GetCellType() override = 0;
43  int GetCellDimension() override { return 2; }
44  int RequiresInitialization() override { return 1; }
45  int GetNumberOfEdges() override { return 3; }
46  int GetNumberOfFaces() override { return 0; }
47  vtkCell* GetEdge(int edgeId) override = 0;
48  void SetEdgeIdsAndPoints(int edgeId,
49  const std::function<void(const vtkIdType&)>& set_number_of_ids_and_points,
50  const std::function<void(const vtkIdType&, const vtkIdType&)>& set_ids_and_points);
51  vtkCell* GetFace(int) override { return nullptr; }
52 
53  void Initialize() override;
54 
55  int CellBoundary(int subId, const double pcoords[3], vtkIdList* pts) override;
56  int EvaluatePosition(const double x[3], double closestPoint[3], int& subId, double pcoords[3],
57  double& dist2, double weights[]) override;
58  void EvaluateLocation(int& subId, const double pcoords[3], double x[3], double* weights) override;
59  void Contour(double value, vtkDataArray* cellScalars, vtkIncrementalPointLocator* locator,
60  vtkCellArray* verts, vtkCellArray* lines, vtkCellArray* polys, vtkPointData* inPd,
61  vtkPointData* outPd, vtkCellData* inCd, vtkIdType cellId, vtkCellData* outCd) override;
62  void Clip(double value, vtkDataArray* cellScalars, vtkIncrementalPointLocator* locator,
63  vtkCellArray* polys, vtkPointData* inPd, vtkPointData* outPd, vtkCellData* inCd,
64  vtkIdType cellId, vtkCellData* outCd, int insideOut) override;
65  int IntersectWithLine(const double p1[3], const double p2[3], double tol, double& t, double x[3],
66  double pcoords[3], int& subId) override;
67  int Triangulate(int index, vtkIdList* ptIds, vtkPoints* pts) override;
68  void JacobianInverse(const double pcoords[3], double** inverse, double* derivs);
69  void Derivatives(
70  int subId, const double pcoords[3], const double* values, int dim, double* derivs) override;
71  void SetParametricCoords();
72  double* GetParametricCoords() override;
73 
74  int GetParametricCenter(double pcoords[3]) override;
75  double GetParametricDistance(const double pcoords[3]) override;
76 
77  void InterpolateFunctions(const double pcoords[3], double* weights) override = 0;
78  void InterpolateDerivs(const double pcoords[3], double* derivs) override = 0;
79 
80  vtkIdType GetOrder() const { return this->Order; }
81  vtkIdType ComputeOrder();
84  static bool PointCountSupportsUniformOrder(vtkIdType pointsPerTri);
85 
86  void ToBarycentricIndex(vtkIdType index, vtkIdType* bindex);
87  vtkIdType ToIndex(const vtkIdType* bindex);
88 
89  static void BarycentricIndex(vtkIdType index, vtkIdType* bindex, vtkIdType order);
90  static vtkIdType Index(const vtkIdType* bindex, vtkIdType order);
91 
92  static double Eta(vtkIdType n, vtkIdType chi, double sigma);
93  static double Deta(vtkIdType n, vtkIdType chi, double sigma);
94  virtual vtkHigherOrderCurve* GetEdgeCell() = 0;
95 
96 protected:
98  ~vtkHigherOrderTriangle() override;
99 
100  vtkIdType GetNumberOfSubtriangles() const { return this->NumberOfSubtriangles; }
101  vtkIdType ComputeNumberOfSubtriangles();
102 
103  // Description:
104  // Given the index of the subtriangle, compute the barycentric indices of
105  // the subtriangle's vertices.
106  void SubtriangleBarycentricPointIndices(vtkIdType cellIndex, vtkIdType (&pointBIndices)[3][3]);
107 
109  vtkDoubleArray* Scalars; // used to avoid New/Delete in contouring/clipping
113 
114  std::vector<vtkIdType> BarycentricIndexMap;
115  std::vector<vtkIdType> IndexMap;
116  std::vector<vtkIdType> SubtriangleIndexMap;
117 
118 private:
120  void operator=(const vtkHigherOrderTriangle&) = delete;
121 };
122 
123 VTK_ABI_NAMESPACE_END
124 #endif
vtkIdType GetNumberOfSubtriangles() const
represent and manipulate point attribute data
Definition: vtkPointData.h:29
int RequiresInitialization() override
Some cells require initialization prior to access.
std::vector< vtkIdType > IndexMap
virtual int GetCellType()=0
Return the type of cell.
vtkSmartPointer< vtkPoints > PointParametricCoordinates
vtkCell * GetFace(int) override
Return the face cell from the faceId of the cell.
represent and manipulate cell attribute data
Definition: vtkCellData.h:30
Abstract class in support of both point location and point insertion.
virtual int Triangulate(int index, vtkIdList *ptIds, vtkPoints *pts)=0
Generate simplices of proper dimension.
virtual void InterpolateFunctions(const double vtkNotUsed(pcoords)[3], double *vtkNotUsed(weight))
Compute the interpolation functions/derivatives (aka shape functions/derivatives) No-ops at this leve...
Definition: vtkCell.h:380
abstract superclass for non-linear cells
int vtkIdType
Definition: vtkType.h:315
std::vector< vtkIdType > BarycentricIndexMap
void PrintSelf(ostream &os, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
virtual int CellBoundary(int subId, const double pcoords[3], vtkIdList *pts)=0
Given parametric coordinates of a point, return the closest cell boundary, and whether the point is i...
virtual void InterpolateDerivs(const double vtkNotUsed(pcoords)[3], double *vtkNotUsed(derivs))
Definition: vtkCell.h:383
dynamic, self-adjusting array of double
virtual double GetParametricDistance(const double pcoords[3])
Return the distance of the parametric coordinate provided to the cell.
abstract class to specify cell behavior
Definition: vtkCell.h:49
virtual void EvaluateLocation(int &subId, const double pcoords[3], double x[3], double *weights)=0
Determine global coordinate (x[3]) from subId and parametric coordinates.
a simple class to control print indentation
Definition: vtkIndent.h:28
list of point or cell ids
Definition: vtkIdList.h:22
abstract superclass for arrays of numeric data
Definition: vtkDataArray.h:44
int GetNumberOfFaces() override
Return the number of faces in the cell.
int GetNumberOfEdges() override
Return the number of edges in the cell.
virtual void Clip(double value, vtkDataArray *cellScalars, vtkIncrementalPointLocator *locator, vtkCellArray *connectivity, vtkPointData *inPd, vtkPointData *outPd, vtkCellData *inCd, vtkIdType cellId, vtkCellData *outCd, int insideOut)=0
Cut (or clip) the cell based on the input cellScalars and the specified value.
int GetCellDimension() override
Return the topological dimensional of the cell (0,1,2, or 3).
virtual int EvaluatePosition(const double x[3], double closestPoint[3], int &subId, double pcoords[3], double &dist2, double weights[])=0
Given a point x[3] return inside(=1), outside(=0) cell, or (-1) computational problem encountered; ev...
object to represent cell connectivity
Definition: vtkCellArray.h:175
virtual vtkCell * GetEdge(int edgeId)=0
Return the edge cell from the edgeId of the cell.
a cell that represents a triangle
Definition: vtkTriangle.h:27
virtual void Contour(double value, vtkDataArray *cellScalars, vtkIncrementalPointLocator *locator, vtkCellArray *verts, vtkCellArray *lines, vtkCellArray *polys, vtkPointData *inPd, vtkPointData *outPd, vtkCellData *inCd, vtkIdType cellId, vtkCellData *outCd)=0
Generate contouring primitives.
A 2D cell that represents an arbitrary order HigherOrder triangle.
virtual void Derivatives(int subId, const double pcoords[3], const double *values, int dim, double *derivs)=0
Compute derivatives given cell subId and parametric coordinates.
std::vector< vtkIdType > SubtriangleIndexMap
virtual void Initialize()
Definition: vtkCell.h:104
virtual double * GetParametricCoords())
Return a contiguous array of parametric coordinates of the points defining this cell.
virtual int GetParametricCenter(double pcoords[3])
Return center of the cell in parametric coordinates.
virtual int IntersectWithLine(const double p1[3], const double p2[3], double tol, double &t, double x[3], double pcoords[3], int &subId)=0
Intersect with a ray.
represent and manipulate 3D points
Definition: vtkPoints.h:28