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=========================================================================*/
30
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
44class vtkDoubleArray;
46class vtkTriangle;
47
48class VTKCOMMONDATAMODEL_EXPORT vtkHigherOrderTriangle : public vtkNonLinearCell
49{
50public:
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);
82 int subId, const double pcoords[3], const double* values, int dim, double* derivs) override;
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; }
94
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")
110
111protected:
114
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
133private:
135 void operator=(const vtkHigherOrderTriangle&) = delete;
136};
137
138#endif
object to represent cell connectivity
represent and manipulate cell attribute data
Definition vtkCellData.h:42
dynamic, self-adjusting array of double
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)
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.
double * GetParametricCoords() override
Return a contiguous array of parametric coordinates of the points defining this 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.
static vtkIdType Index(const vtkIdType *bindex, vtkIdType order)
void InterpolateFunctions(const double pcoords[3], double *weights) override=0
std::vector< vtkIdType > SubtriangleIndexMap
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
vtkCell * GetFace(int) override
Return the face cell from the faceId of the cell.
virtual vtkHigherOrderCurve * getEdgeCell()
std::vector< vtkIdType > IndexMap
double GetParametricDistance(const double pcoords[3]) override
Return the distance of the parametric coordinate provided to the cell.
void ToBarycentricIndex(vtkIdType index, vtkIdType *bindex)
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.
static double d_eta(vtkIdType n, vtkIdType chi, double sigma)
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.
virtual vtkHigherOrderCurve * GetEdgeCell()=0
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)
vtkCell * GetEdge(int edgeId) override=0
Return the edge cell from the edgeId of the cell.
static double Eta(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:34
Abstract class in support of both point location and point insertion.
a simple class to control print indentation
Definition vtkIndent.h:40
represent and manipulate point attribute data
represent and manipulate 3D points
Definition vtkPoints.h:40
Hold a reference to a vtkObjectBase instance.
a cell that represents a triangle
Definition vtkTriangle.h:39
#define vtkDataArray
#define VTK_DEPRECATED_IN_9_1_0(reason)
int vtkIdType
Definition vtkType.h:332