00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00057 #ifndef __vtkPyramid_h
00058 #define __vtkPyramid_h
00059
00060 #include "vtkCell3D.h"
00061 #include "vtkLine.h"
00062 #include "vtkTriangle.h"
00063 #include "vtkQuad.h"
00064
00065 class vtkUnstructuredGrid;
00066
00067 class VTK_COMMON_EXPORT vtkPyramid : public vtkCell3D
00068 {
00069 public:
00070 vtkTypeMacro(vtkPyramid,vtkCell);
00071
00073 static vtkPyramid *New();
00074
00076
00077 virtual void GetEdgePoints(int edgeId, int* &pts);
00078 virtual void GetFacePoints(int faceId, int* &pts);
00080
00082
00083 vtkCell *MakeObject();
00084 int GetCellType() {return VTK_PYRAMID;}
00085 int GetCellDimension() {return 3;}
00086 int GetNumberOfEdges() {return 8;}
00087 int GetNumberOfFaces() {return 5;}
00088 vtkCell *GetEdge(int edgeId);
00089 vtkCell *GetFace(int faceId);
00090 int CellBoundary(int subId, float pcoords[3], vtkIdList *pts);
00091 void Contour(float value, vtkDataArray *cellScalars,
00092 vtkPointLocator *locator, vtkCellArray *verts,
00093 vtkCellArray *lines, vtkCellArray *polys,
00094 vtkPointData *inPd, vtkPointData *outPd,
00095 vtkCellData *inCd, vtkIdType cellId, vtkCellData *outCd);
00096 int EvaluatePosition(float x[3], float* closestPoint,
00097 int& subId, float pcoords[3],
00098 float& dist2, float *weights);
00099 void EvaluateLocation(int& subId, float pcoords[3], float x[3],
00100 float *weights);
00101 int IntersectWithLine(float p1[3], float p2[3], float tol, float& t,
00102 float x[3], float pcoords[3], int& subId);
00103 int Triangulate(int index, vtkIdList *ptIds, vtkPoints *pts);
00104 void Derivatives(int subId, float pcoords[3], float *values,
00105 int dim, float *derivs);
00107
00109 int GetParametricCenter(float pcoords[3]);
00110
00112
00114 static void InterpolationFunctions(float pcoords[3], float weights[5]);
00115 static void InterpolationDerivs(float pcoords[3], float derivs[15]);
00116 int JacobianInverse(float pcoords[3], double **inverse, float derivs[15]);
00117 static int *GetEdgeArray(int edgeId);
00118 static int *GetFaceArray(int faceId);
00120
00121 protected:
00122 vtkPyramid();
00123 ~vtkPyramid();
00124
00125 vtkLine *Line;
00126 vtkTriangle *Triangle;
00127 vtkQuad *Quad;
00128
00129 private:
00130 vtkPyramid(const vtkPyramid&);
00131 void operator=(const vtkPyramid&);
00132 };
00133
00134 inline int vtkPyramid::GetParametricCenter(float pcoords[3])
00135 {
00136 pcoords[0] = pcoords[1] = 0.5;
00137 pcoords[2] = 0.333333;
00138 return 0;
00139 }
00140
00141 #endif
00142
00143
00144