VTK  9.0.2
vtkQuadraticPyramid.h
Go to the documentation of this file.
1 /*=========================================================================
2 
3  Program: Visualization Toolkit
4  Module: vtkQuadraticPyramid.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 =========================================================================*/
39 #ifndef vtkQuadraticPyramid_h
40 #define vtkQuadraticPyramid_h
41 
42 #include "vtkCommonDataModelModule.h" // For export macro
43 #include "vtkNonLinearCell.h"
44 
45 class vtkQuadraticEdge;
46 class vtkQuadraticQuad;
48 class vtkTetra;
49 class vtkPyramid;
50 class vtkDoubleArray;
51 
52 class VTKCOMMONDATAMODEL_EXPORT vtkQuadraticPyramid : public vtkNonLinearCell
53 {
54 public:
57  void PrintSelf(ostream& os, vtkIndent indent) override;
58 
60 
64  int GetCellType() override { return VTK_QUADRATIC_PYRAMID; }
65  int GetCellDimension() override { return 3; }
66  int GetNumberOfEdges() override { return 8; }
67  int GetNumberOfFaces() override { return 5; }
68  vtkCell* GetEdge(int edgeId) override;
69  vtkCell* GetFace(int faceId) override;
71 
72  int CellBoundary(int subId, const double pcoords[3], vtkIdList* pts) override;
73  void Contour(double value, vtkDataArray* cellScalars, vtkIncrementalPointLocator* locator,
74  vtkCellArray* verts, vtkCellArray* lines, vtkCellArray* polys, vtkPointData* inPd,
75  vtkPointData* outPd, vtkCellData* inCd, vtkIdType cellId, vtkCellData* outCd) override;
76  int EvaluatePosition(const double x[3], double closestPoint[3], int& subId, double pcoords[3],
77  double& dist2, double weights[]) override;
78  void EvaluateLocation(int& subId, const double pcoords[3], double x[3], double* weights) override;
79  int Triangulate(int index, vtkIdList* ptIds, vtkPoints* pts) override;
81  int subId, const double pcoords[3], const double* values, int dim, double* derivs) override;
82  double* GetParametricCoords() override;
83 
89  void Clip(double value, vtkDataArray* cellScalars, vtkIncrementalPointLocator* locator,
90  vtkCellArray* tets, vtkPointData* inPd, vtkPointData* outPd, vtkCellData* inCd,
91  vtkIdType cellId, vtkCellData* outCd, int insideOut) override;
92 
97  int IntersectWithLine(const double p1[3], const double p2[3], double tol, double& t, double x[3],
98  double pcoords[3], int& subId) override;
99 
103  int GetParametricCenter(double pcoords[3]) override;
104 
105  static void InterpolationFunctions(const double pcoords[3], double weights[13]);
106  static void InterpolationDerivs(const double pcoords[3], double derivs[39]);
108 
112  void InterpolateFunctions(const double pcoords[3], double weights[13]) override
113  {
115  }
116  void InterpolateDerivs(const double pcoords[3], double derivs[39]) override
117  {
119  }
121 
122 
129  static const vtkIdType* GetEdgeArray(vtkIdType edgeId);
130  static const vtkIdType* GetFaceArray(vtkIdType faceId);
132 
138  void JacobianInverse(const double pcoords[3], double** inverse, double derivs[39]);
139 
140 protected:
143 
152  vtkDoubleArray* Scalars; // used to avoid New/Delete in contouring/clipping
153 
155 
161  void Subdivide(
162  vtkPointData* inPd, vtkCellData* inCd, vtkIdType cellId, vtkDataArray* cellScalars);
164 
165 
171  void ResizeArrays(vtkIdType newSize);
173 
174 private:
175  vtkQuadraticPyramid(const vtkQuadraticPyramid&) = delete;
176  void operator=(const vtkQuadraticPyramid&) = delete;
177 };
178 //----------------------------------------------------------------------------
179 // Return the center of the quadratic pyramid in parametric coordinates.
180 //
181 inline int vtkQuadraticPyramid::GetParametricCenter(double pcoords[3])
182 {
183  pcoords[0] = pcoords[1] = 6.0 / 13.0;
184  pcoords[2] = 3.0 / 13.0;
185  return 0;
186 }
187 
188 #endif
object to represent cell connectivity
Definition: vtkCellArray.h:180
represent and manipulate cell attribute data
Definition: vtkCellData.h:33
abstract class to specify cell behavior
Definition: vtkCell.h:57
virtual int GetParametricCenter(double pcoords[3])
Return center of the cell in parametric coordinates.
abstract superclass for arrays of numeric data
Definition: vtkDataArray.h:50
dynamic, self-adjusting array of double
list of point or cell ids
Definition: vtkIdList.h:31
Abstract class in support of both point location and point insertion.
a simple class to control print indentation
Definition: vtkIndent.h:34
abstract superclass for non-linear cells
represent and manipulate point attribute data
Definition: vtkPointData.h:32
represent and manipulate 3D points
Definition: vtkPoints.h:34
a 3D cell that represents a linear pyramid
Definition: vtkPyramid.h:44
cell represents a parabolic, isoparametric edge
cell represents a parabolic, 13-node isoparametric pyramid
int GetCellType() override
Implement the vtkCell API.
vtkDoubleArray * Scalars
vtkQuadraticQuad * Face
double * GetParametricCoords() override
Return a contiguous array of parametric coordinates of the points defining this cell.
vtkCell * GetFace(int faceId) override
Return the face cell from the faceId of the cell.
void ResizeArrays(vtkIdType newSize)
Resize the superclasses' member arrays to newSize where newSize should either be 13 or 14.
int GetParametricCenter(double pcoords[3]) override
Return the center of the quadratic pyramid in parametric coordinates.
void Clip(double value, vtkDataArray *cellScalars, vtkIncrementalPointLocator *locator, vtkCellArray *tets, vtkPointData *inPd, vtkPointData *outPd, vtkCellData *inCd, vtkIdType cellId, vtkCellData *outCd, int insideOut) override
Clip this quadratic triangle using scalar value provided.
int GetCellDimension() override
Return the topological dimensional of the cell (0,1,2, or 3).
static const vtkIdType * GetEdgeArray(vtkIdType edgeId)
Return the ids of the vertices defining edge/face (edgeId/‘faceId’).
int Triangulate(int index, vtkIdList *ptIds, vtkPoints *pts) override
Generate simplices of proper dimension.
vtkQuadraticTriangle * TriangleFace
vtkCell * GetEdge(int edgeId) override
Return the edge cell from the edgeId of the cell.
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.
~vtkQuadraticPyramid() override
void Derivatives(int subId, const double pcoords[3], const double *values, int dim, double *derivs) override
Compute derivatives given cell subId and parametric coordinates.
static void InterpolationDerivs(const double pcoords[3], double derivs[39])
void InterpolateFunctions(const double pcoords[3], double weights[13]) override
Compute the interpolation functions/derivatives (aka shape functions/derivatives)
void JacobianInverse(const double pcoords[3], double **inverse, double derivs[39])
Given parametric coordinates compute inverse Jacobian transformation matrix.
int GetNumberOfFaces() override
Return the number of faces in the cell.
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...
static vtkQuadraticPyramid * New()
static void InterpolationFunctions(const double pcoords[3], double weights[13])
void EvaluateLocation(int &subId, const double pcoords[3], double x[3], double *weights) override
Determine global coordinate (x[3]) from subId and parametric coordinates.
static const vtkIdType * GetFaceArray(vtkIdType faceId)
vtkQuadraticEdge * Edge
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 PrintSelf(ostream &os, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
void InterpolateDerivs(const double pcoords[3], double derivs[39]) override
void Subdivide(vtkPointData *inPd, vtkCellData *inCd, vtkIdType cellId, vtkDataArray *cellScalars)
This method adds in a point at the center of the quadrilateral face and then interpolates values to t...
int IntersectWithLine(const double p1[3], const double p2[3], double tol, double &t, double x[3], double pcoords[3], int &subId) override
Line-edge intersection.
int GetNumberOfEdges() override
Return the number of edges in the cell.
vtkDoubleArray * CellScalars
cell represents a parabolic, 8-node isoparametric quad
cell represents a parabolic, isoparametric triangle
a 3D cell that represents a tetrahedron
Definition: vtkTetra.h:42
@ value
Definition: vtkX3D.h:226
@ index
Definition: vtkX3D.h:252
@ VTK_QUADRATIC_PYRAMID
Definition: vtkCellType.h:72
int vtkIdType
Definition: vtkType.h:338