VTK  9.3.0
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);
70  int subId, const double pcoords[3], const double* values, int dim, double* derivs) override;
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; }
84  static bool PointCountSupportsUniformOrder(vtkIdType pointsPerTri);
85 
87  vtkIdType ToIndex(const vtkIdType* bindex);
88 
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);
95 
96 protected:
99 
100  vtkIdType GetNumberOfSubtriangles() const { return this->NumberOfSubtriangles; }
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
object to represent cell connectivity
Definition: vtkCellArray.h:185
represent and manipulate cell attribute data
Definition: vtkCellData.h:40
abstract class to specify cell behavior
Definition: vtkCell.h:59
abstract superclass for arrays of numeric data
Definition: vtkDataArray.h:54
dynamic, self-adjusting array of double
A 2D cell that represents an arbitrary order HigherOrder triangle.
void SetEdgeIdsAndPoints(int edgeId, const std::function< void(const vtkIdType &)> &set_number_of_ids_and_points, const std::function< void(const vtkIdType &, const vtkIdType &)> &set_ids_and_points)
~vtkHigherOrderTriangle() override
int GetNumberOfFaces() override
Return the number of faces in the cell.
void SubtriangleBarycentricPointIndices(vtkIdType cellIndex, vtkIdType(&pointBIndices)[3][3])
vtkIdType GetNumberOfSubtriangles() const
int GetNumberOfEdges() override
Return the number of edges in the cell.
vtkCell * GetFace(int) override
Return the face cell from the faceId of the cell.
int GetParametricCenter(double pcoords[3]) override
Return center of the cell in parametric coordinates.
int Triangulate(int index, vtkIdList *ptIds, vtkPoints *pts) override
Generate simplices of proper dimension.
virtual vtkHigherOrderCurve * GetEdgeCell()=0
static vtkIdType Index(const vtkIdType *bindex, vtkIdType order)
void InterpolateFunctions(const double pcoords[3], double *weights) override=0
std::vector< vtkIdType > SubtriangleIndexMap
double * GetParametricCoords() override
Return a contiguous array of parametric coordinates of the points defining this cell.
void PrintSelf(ostream &os, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
void EvaluateLocation(int &subId, const double pcoords[3], double x[3], double *weights) override
Determine global coordinate (x[3]) from subId and parametric coordinates.
int IntersectWithLine(const double p1[3], const double p2[3], double tol, double &t, double x[3], double pcoords[3], int &subId) override
Intersect with a ray.
void Derivatives(int subId, const double pcoords[3], const double *values, int dim, double *derivs) override
Compute derivatives given cell subId and parametric coordinates.
void InterpolateDerivs(const double pcoords[3], double *derivs) override=0
std::vector< vtkIdType > IndexMap
double GetParametricDistance(const double pcoords[3]) override
Return the distance of the parametric coordinate provided to the cell.
vtkCell * GetEdge(int edgeId) override=0
Return the edge cell from the edgeId of the cell.
void ToBarycentricIndex(vtkIdType index, vtkIdType *bindex)
static bool PointCountSupportsUniformOrder(vtkIdType pointsPerTri)
Return true if the number of points supports a cell of uniform degree along each axis.
int GetCellDimension() override
Return the topological dimensional of the cell (0,1,2, or 3).
int RequiresInitialization() override
Some cells require initialization prior to access.
vtkSmartPointer< vtkPoints > PointParametricCoordinates
int CellBoundary(int subId, const double pcoords[3], vtkIdList *pts) override
Given parametric coordinates of a point, return the closest cell boundary, and whether the point is i...
void Contour(double value, vtkDataArray *cellScalars, vtkIncrementalPointLocator *locator, vtkCellArray *verts, vtkCellArray *lines, vtkCellArray *polys, vtkPointData *inPd, vtkPointData *outPd, vtkCellData *inCd, vtkIdType cellId, vtkCellData *outCd) override
Generate contouring primitives.
vtkIdType ComputeNumberOfSubtriangles()
vtkIdType ToIndex(const vtkIdType *bindex)
static void BarycentricIndex(vtkIdType index, vtkIdType *bindex, vtkIdType order)
void Initialize() override
int GetCellType() override=0
Return the type of cell.
void Clip(double value, vtkDataArray *cellScalars, vtkIncrementalPointLocator *locator, vtkCellArray *polys, vtkPointData *inPd, vtkPointData *outPd, vtkCellData *inCd, vtkIdType cellId, vtkCellData *outCd, int insideOut) override
Cut (or clip) the cell based on the input cellScalars and the specified value.
void JacobianInverse(const double pcoords[3], double **inverse, double *derivs)
std::vector< vtkIdType > BarycentricIndexMap
static double Deta(vtkIdType n, vtkIdType chi, double sigma)
static double Eta(vtkIdType n, vtkIdType chi, double sigma)
vtkIdType ComputeOrder()
int EvaluatePosition(const double x[3], double closestPoint[3], int &subId, double pcoords[3], double &dist2, double weights[]) override
Given a point x[3] return inside(=1), outside(=0) cell, or (-1) computational problem encountered; ev...
list of point or cell ids
Definition: vtkIdList.h:32
Abstract class in support of both point location and point insertion.
a simple class to control print indentation
Definition: vtkIndent.h:38
abstract superclass for non-linear cells
represent and manipulate point attribute data
Definition: vtkPointData.h:39
represent and manipulate 3D points
Definition: vtkPoints.h:38
a cell that represents a triangle
Definition: vtkTriangle.h:37
@ order
Definition: vtkX3D.h:440
@ function
Definition: vtkX3D.h:249
@ value
Definition: vtkX3D.h:220
@ index
Definition: vtkX3D.h:246
int vtkIdType
Definition: vtkType.h:315