VTK
vtkLagrangeWedge.h
Go to the documentation of this file.
1 /*=========================================================================
2 
3  Program: Visualization Toolkit
4  Module: vtkLagrangeWedge.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 =========================================================================*/
15 // .NAME vtkLagrangeWedge
16 // .SECTION Description
17 // .SECTION See Also
18 
19 #ifndef vtkLagrangeWedge_h
20 #define vtkLagrangeWedge_h
21 
22 #include "vtkCommonDataModelModule.h" // For export macro
23 #include "vtkNonLinearCell.h"
24 #include "vtkSmartPointer.h" // For member variable.
25 #include "vtkCellType.h" // For GetCellType.
26 #include "vtkNew.h" // For member variable.
27 
28 class vtkCellData;
29 class vtkDoubleArray;
30 class vtkWedge;
31 class vtkIdList;
32 class vtkPointData;
33 class vtkPoints;
34 class vtkVector3d;
35 class vtkVector3i;
36 class vtkLagrangeCurve;
40 
41 class VTKCOMMONDATAMODEL_EXPORT vtkLagrangeWedge : public vtkNonLinearCell
42 {
43 public:
44  static vtkLagrangeWedge* New();
46  void PrintSelf(ostream& os, vtkIndent indent) override;
47 
48  int GetCellType() override { return VTK_LAGRANGE_WEDGE; }
49  int GetCellDimension() override { return 3; }
50  int RequiresInitialization() override { return 1; }
51  int GetNumberOfEdges() override { return 9; }
52  int GetNumberOfFaces() override { return 5; }
53  vtkCell* GetEdge(int edgeId) override;
54  vtkCell* GetFace(int faceId) override;
55 
56  void Initialize() override;
57 
58  int CellBoundary(int subId, double pcoords[3], vtkIdList* pts) override;
59  int EvaluatePosition(
60  double x[3], double* closestPoint,
61  int& subId, double pcoords[3],
62  double& dist2, double* weights) override;
63  void EvaluateLocation(
64  int& subId, double pcoords[3], double x[3],
65  double* weights) override;
66  void Contour(
67  double value, vtkDataArray* cellScalars,
69  vtkCellArray* lines, vtkCellArray* polys,
70  vtkPointData* inPd, vtkPointData* outPd,
71  vtkCellData* inCd, vtkIdType cellId, vtkCellData* outCd) override;
72  void Clip(
73  double value, vtkDataArray* cellScalars,
75  vtkPointData* inPd, vtkPointData* outPd,
76  vtkCellData* inCd, vtkIdType cellId, vtkCellData* outCd,
77  int insideOut) override;
79  double p1[3], double p2[3], double tol, double& t,
80  double x[3], double pcoords[3], int& subId) override;
81  int Triangulate(int index, vtkIdList* ptIds, vtkPoints* pts) override;
82  void Derivatives(
83  int subId, double pcoords[3], double* values,
84  int dim, double* derivs) override;
85  double* GetParametricCoords() override;
86  int GetParametricCenter(double center[3]) override;
87 
88  double GetParametricDistance(double pcoords[3]) override;
89 
90  const int* GetOrder();
91  int GetOrder(int i) { return this->GetOrder()[i]; }
92 
93  void InterpolateFunctions(double pcoords[3], double* weights) override;
94  void InterpolateDerivs(double pcoords[3], double* derivs) override;
95 
96  bool SubCellCoordinatesFromId(vtkVector3i& ijk, int subId);
97  bool SubCellCoordinatesFromId(int& i, int& j, int& k, int subId);
98  static int PointIndexFromIJK(int i, int j, int k, const int* order);
99  int PointIndexFromIJK(int i, int j, int k);
100  bool TransformApproxToCellParams(int subCell, double* pcoords);
101  bool TransformFaceToCellParams(int bdyFace, double* pcoords);
102 
103  static int GetNumberOfApproximatingWedges(const int* order);
105  { return vtkLagrangeWedge::GetNumberOfApproximatingWedges(this->GetOrder()); }
106 
107 protected:
109  ~vtkLagrangeWedge() override;
110 
111  vtkWedge* GetApprox();
112  void PrepareApproxData(vtkPointData* pd, vtkCellData* cd, vtkIdType cellId, vtkDataArray* cellScalars);
113  vtkWedge* GetApproximateWedge(
114  int subId, vtkDataArray* scalarsIn = nullptr, vtkDataArray* scalarsOut = nullptr);
115 
116  vtkLagrangeTriangle* GetTriangularFace(int iAxis, int k);
117  vtkLagrangeQuadrilateral* GetQuadrilateralFace(int di, int dj);
118 
119  int Order[4];
132 
133 private:
134  vtkLagrangeWedge(const vtkLagrangeWedge&) = delete;
135  void operator=(const vtkLagrangeWedge&) = delete;
136 };
137 
139 {
140  center[0] = center[1] = 1./3.;
141  center[2] = 0.5;
142  return 0;
143 }
144 
145 #endif // vtkLagrangeWedge_h
int GetParametricCenter(double center[3]) override
Return center of the cell in parametric coordinates.
represent and manipulate point attribute data
Definition: vtkPointData.h:31
vtkNew< vtkLagrangeCurve > BdyEdge
virtual void InterpolateFunctions(double vtkNotUsed(pcoords)[3], double *vtkNotUsed(weight))
Compute the interpolation functions/derivatives (aka shape functions/derivatives) No-ops at this leve...
Definition: vtkCell.h:354
represent and manipulate cell attribute data
Definition: vtkCellData.h:32
vtkSmartPointer< vtkWedge > Approx
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 EvaluateLocation(int &subId, double pcoords[3], double x[3], double *weights)=0
Determine global coordinate (x[3]) from subId and parametric coordinates.
virtual int EvaluatePosition(double x[3], double *closestPoint, 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...
int GetCellType() override
Return the type of cell.
abstract superclass for non-linear cells
int vtkIdType
Definition: vtkType.h:345
void PrintSelf(ostream &os, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
int GetNumberOfApproximatingWedges()
int GetCellDimension() override
Return the topological dimensional of the cell (0,1,2, or 3).
dynamic, self-adjusting array of double
vtkNew< vtkDoubleArray > Scalars
abstract class to specify cell behavior
Definition: vtkCell.h:56
virtual double GetParametricDistance(double pcoords[3])
Return the distance of the parametric coordinate provided to the cell.
vtkNew< vtkDoubleArray > CellScalars
a simple class to control print indentation
Definition: vtkIndent.h:33
list of point or cell ids
Definition: vtkIdList.h:30
virtual void Derivatives(int subId, double pcoords[3], double *values, int dim, double *derivs)=0
Compute derivatives given cell subId and parametric coordinates.
vtkNew< vtkLagrangeTriangle > BdyTri
abstract superclass for arrays of numeric data
Definition: vtkDataArray.h:48
virtual int IntersectWithLine(double p1[3], double p2[3], double tol, double &t, double x[3], double pcoords[3], int &subId)=0
Intersect with a ray.
virtual void InterpolateDerivs(double vtkNotUsed(pcoords)[3], double *vtkNotUsed(derivs))
Definition: vtkCell.h:357
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.
vtkNew< vtkLagrangeInterpolation > Interp
vtkSmartPointer< vtkCellData > ApproxCD
virtual vtkCell * GetFace(int faceId)=0
Return the face cell from the faceId of the cell.
object to represent cell connectivity
Definition: vtkCellArray.h:44
int RequiresInitialization() override
Some cells require initialization prior to access.
virtual vtkCell * GetEdge(int edgeId)=0
Return the edge cell from the edgeId of the cell.
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.
int GetNumberOfFaces() override
Return the number of faces in the cell.
vtkNew< vtkIdList > TmpIds
vtkSmartPointer< vtkPointData > ApproxPD
virtual int CellBoundary(int subId, 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 Initialize()
Definition: vtkCell.h:111
static vtkObject * New()
Create an object with Debug turned off, modified time initialized to zero, and reference counting on...
vtkSmartPointer< vtkPoints > PointParametricCoordinates
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.
vtkNew< vtkPoints > TmpPts
int GetNumberOfEdges() override
Return the number of edges in the cell.
vtkNew< vtkLagrangeQuadrilateral > BdyQuad
a 3D cell that represents a linear wedge
Definition: vtkWedge.h:43
represent and manipulate 3D points
Definition: vtkPoints.h:33