VTK  9.2.6
vtkHigherOrderTriangle.h
Go to the documentation of this file.
1 /*=========================================================================
2 
3  Program: Visualization Toolkit
4  Module: vtkHigherOrderTriangle.h
5 
6  Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
7  All rights reserved.
8  See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
9 
10  This software is distributed WITHOUT ANY WARRANTY; without even
11  the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
12  PURPOSE. See the above copyright notice for more information.
13 
14 =========================================================================*/
31 #ifndef vtkHigherOrderTriangle_h
32 #define vtkHigherOrderTriangle_h
33 
34 #include <functional> //For std::function
35 
36 #include "vtkCommonDataModelModule.h" // For export macro
37 #include "vtkDeprecation.h" // For deprecation macros
38 #include "vtkNew.h" // For member variable.
39 #include "vtkNonLinearCell.h"
40 #include "vtkSmartPointer.h" // For member variable.
41 
42 #include <vector> // For caching
43 
44 class vtkDoubleArray;
46 class vtkTriangle;
47 
48 class VTKCOMMONDATAMODEL_EXPORT vtkHigherOrderTriangle : public vtkNonLinearCell
49 {
50 public:
52  void PrintSelf(ostream& os, vtkIndent indent) override;
53 
54  int GetCellType() override = 0;
55  int GetCellDimension() override { return 2; }
56  int RequiresInitialization() override { return 1; }
57  int GetNumberOfEdges() override { return 3; }
58  int GetNumberOfFaces() override { return 0; }
59  vtkCell* GetEdge(int edgeId) override = 0;
60  void SetEdgeIdsAndPoints(int edgeId,
61  const std::function<void(const vtkIdType&)>& set_number_of_ids_and_points,
62  const std::function<void(const vtkIdType&, const vtkIdType&)>& set_ids_and_points);
63  vtkCell* GetFace(int) override { return nullptr; }
64 
65  void Initialize() override;
66 
67  int CellBoundary(int subId, const double pcoords[3], vtkIdList* pts) override;
68  int EvaluatePosition(const double x[3], double closestPoint[3], int& subId, double pcoords[3],
69  double& dist2, double weights[]) override;
70  void EvaluateLocation(int& subId, const double pcoords[3], double x[3], double* weights) override;
71  void Contour(double value, vtkDataArray* cellScalars, vtkIncrementalPointLocator* locator,
72  vtkCellArray* verts, vtkCellArray* lines, vtkCellArray* polys, vtkPointData* inPd,
73  vtkPointData* outPd, vtkCellData* inCd, vtkIdType cellId, vtkCellData* outCd) override;
74  void Clip(double value, vtkDataArray* cellScalars, vtkIncrementalPointLocator* locator,
75  vtkCellArray* polys, vtkPointData* inPd, vtkPointData* outPd, vtkCellData* inCd,
76  vtkIdType cellId, vtkCellData* outCd, int insideOut) override;
77  int IntersectWithLine(const double p1[3], const double p2[3], double tol, double& t, double x[3],
78  double pcoords[3], int& subId) override;
79  int Triangulate(int index, vtkIdList* ptIds, vtkPoints* pts) override;
80  void JacobianInverse(const double pcoords[3], double** inverse, double* derivs);
81  void Derivatives(
82  int subId, const double pcoords[3], const double* values, int dim, double* derivs) override;
83  void SetParametricCoords();
84  double* GetParametricCoords() override;
85 
86  int GetParametricCenter(double pcoords[3]) override;
87  double GetParametricDistance(const double pcoords[3]) override;
88 
89  void InterpolateFunctions(const double pcoords[3], double* weights) override = 0;
90  void InterpolateDerivs(const double pcoords[3], double* derivs) override = 0;
91 
92  vtkIdType GetOrder() const { return this->Order; }
93  vtkIdType ComputeOrder();
94 
95  void ToBarycentricIndex(vtkIdType index, vtkIdType* bindex);
96  vtkIdType ToIndex(const vtkIdType* bindex);
97 
98  static void BarycentricIndex(vtkIdType index, vtkIdType* bindex, vtkIdType order);
99  static vtkIdType Index(const vtkIdType* bindex, vtkIdType order);
100 
101  VTK_DEPRECATED_IN_9_1_0("renamed to Eta")
102  static double eta(vtkIdType n, vtkIdType chi, double sigma);
103  static double Eta(vtkIdType n, vtkIdType chi, double sigma);
104  VTK_DEPRECATED_IN_9_1_0("renamed to Deta")
105  static double d_eta(vtkIdType n, vtkIdType chi, double sigma);
106  static double Deta(vtkIdType n, vtkIdType chi, double sigma);
107  VTK_DEPRECATED_IN_9_1_0("renamed to GetEdgeCell")
108  virtual vtkHigherOrderCurve* getEdgeCell();
109  virtual vtkHigherOrderCurve* GetEdgeCell() = 0;
110 
111 protected:
113  ~vtkHigherOrderTriangle() override;
114 
115  vtkIdType GetNumberOfSubtriangles() const { return this->NumberOfSubtriangles; }
116  vtkIdType ComputeNumberOfSubtriangles();
117 
118  // Description:
119  // Given the index of the subtriangle, compute the barycentric indices of
120  // the subtriangle's vertices.
121  void SubtriangleBarycentricPointIndices(vtkIdType cellIndex, vtkIdType (&pointBIndices)[3][3]);
122 
124  vtkDoubleArray* Scalars; // used to avoid New/Delete in contouring/clipping
128 
129  std::vector<vtkIdType> BarycentricIndexMap;
130  std::vector<vtkIdType> IndexMap;
131  std::vector<vtkIdType> SubtriangleIndexMap;
132 
133 private:
135  void operator=(const vtkHigherOrderTriangle&) = delete;
136 };
137 
138 #endif
#define VTK_DEPRECATED_IN_9_1_0(reason)
represent and manipulate point attribute data
Definition: vtkPointData.h:41
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:41
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:391
abstract superclass for non-linear cells
int vtkIdType
Definition: vtkType.h:332
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:394
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:60
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:39
list of point or cell ids
Definition: vtkIdList.h:33
abstract superclass for arrays of numeric data
Definition: vtkDataArray.h:55
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:186
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:38
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:115
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:39