VTK
vtkLagrangeHexahedron.h
Go to the documentation of this file.
1 /*=========================================================================
2 
3  Program: Visualization Toolkit
4  Module: vtkLagrangeHexahedron.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 vtkLagrangeHexahedron
16 // .SECTION Description
17 // .SECTION See Also
18 
19 #ifndef vtkLagrangeHexahedron_h
20 #define vtkLagrangeHexahedron_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 vtkHexahedron;
31 class vtkIdList;
32 class vtkLagrangeCurve;
35 class vtkPointData;
36 class vtkPoints;
37 class vtkVector3d;
38 class vtkVector3i;
39 
40 class VTKCOMMONDATAMODEL_EXPORT vtkLagrangeHexahedron : public vtkNonLinearCell
41 {
42 public:
43  static vtkLagrangeHexahedron* New();
45  void PrintSelf(ostream& os, vtkIndent indent) override;
46 
47  int GetCellType() override { return VTK_LAGRANGE_HEXAHEDRON; }
48  int GetCellDimension() override { return 3; }
49  int RequiresInitialization() override { return 1; }
50  int GetNumberOfEdges() override { return 12; }
51  int GetNumberOfFaces() override { return 6; }
52  vtkCell* GetEdge(int edgeId) override;
53  vtkCell* GetFace(int faceId) override;
54 
55  void Initialize() override;
56 
57  int CellBoundary(int subId, double pcoords[3], vtkIdList* pts) override;
58  int EvaluatePosition(
59  double x[3], double* closestPoint,
60  int& subId, double pcoords[3],
61  double& dist2, double* weights) override;
62  void EvaluateLocation(
63  int& subId, double pcoords[3], double x[3],
64  double* weights) override;
65  void Contour(
66  double value, vtkDataArray* cellScalars,
68  vtkCellArray* lines, vtkCellArray* polys,
69  vtkPointData* inPd, vtkPointData* outPd,
70  vtkCellData* inCd, vtkIdType cellId, vtkCellData* outCd) override;
71  void Clip(
72  double value, vtkDataArray* cellScalars,
74  vtkPointData* inPd, vtkPointData* outPd,
75  vtkCellData* inCd, vtkIdType cellId, vtkCellData* outCd,
76  int insideOut) override;
78  double p1[3], double p2[3], double tol, double& t,
79  double x[3], double pcoords[3], int& subId) override;
80  int Triangulate(int index, vtkIdList* ptIds, vtkPoints* pts) override;
81  void Derivatives(
82  int subId, double pcoords[3], double* values,
83  int dim, double* derivs) override;
84  double* GetParametricCoords() override;
85  int GetParametricCenter(double center[3]) override;
86 
87  double GetParametricDistance(double pcoords[3]) override;
88 
89  const int* GetOrder();
90  int GetOrder(int i) { return this->GetOrder()[i]; }
91 
92  void InterpolateFunctions(double pcoords[3], double* weights) override;
93  void InterpolateDerivs(double pcoords[3], double* derivs) override;
94 
95  bool SubCellCoordinatesFromId(vtkVector3i& ijk, int subId);
96  bool SubCellCoordinatesFromId(int& i, int& j, int& k, int subId);
97  static int PointIndexFromIJK(int i, int j, int k, const int* order);
98  int PointIndexFromIJK(int i, int j, int k);
99  bool TransformApproxToCellParams(int subCell, double* pcoords);
100  bool TransformFaceToCellParams(int bdyFace, double* pcoords);
101 
102 protected:
104  ~vtkLagrangeHexahedron() override;
105 
106  vtkHexahedron* GetApprox();
107  void PrepareApproxData(vtkPointData* pd, vtkCellData* cd, vtkIdType cellId, vtkDataArray* cellScalars);
108  vtkHexahedron* GetApproximateHex(
109  int subId, vtkDataArray* scalarsIn = nullptr, vtkDataArray* scalarsOut = nullptr);
110 
111  int Order[4];
123 
124 private:
126  void operator=(const vtkLagrangeHexahedron&) = delete;
127 };
128 
130 {
131  center[0] = center[1] = center[2] = 0.5;
132  return 0;
133 }
134 
135 #endif // vtkLagrangeHexahedron_h
represent and manipulate point attribute data
Definition: vtkPointData.h:31
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
vtkNew< vtkLagrangeInterpolation > Interp
represent and manipulate cell attribute data
Definition: vtkCellData.h:32
vtkSmartPointer< vtkPoints > PointParametricCoordinates
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...
abstract superclass for non-linear cells
vtkNew< vtkPoints > TmpPts
vtkNew< vtkLagrangeQuadrilateral > FaceCell
int vtkIdType
Definition: vtkType.h:345
vtkSmartPointer< vtkPointData > ApproxPD
void PrintSelf(ostream &os, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
int GetCellType() override
Return the type of cell.
vtkNew< vtkDoubleArray > CellScalars
dynamic, self-adjusting array of double
vtkSmartPointer< vtkCellData > ApproxCD
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.
int RequiresInitialization() override
Some cells require initialization prior to access.
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< vtkIdList > TmpIds
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.
a cell that represents a linear 3D hexahedron
Definition: vtkHexahedron.h:41
int GetNumberOfEdges() override
Return the number of edges in the cell.
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.
vtkSmartPointer< vtkHexahedron > Approx
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
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 GetCellDimension() override
Return the topological dimensional of the cell (0,1,2, or 3).
vtkNew< vtkDoubleArray > Scalars
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...
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.
int GetNumberOfFaces() override
Return the number of faces in the cell.
int GetParametricCenter(double center[3]) override
Return center of the cell in parametric coordinates.
vtkNew< vtkLagrangeCurve > EdgeCell
represent and manipulate 3D points
Definition: vtkPoints.h:33