00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00037 #ifndef __vtkPentagonalPrism_h
00038 #define __vtkPentagonalPrism_h
00039
00040 #include "vtkCell3D.h"
00041
00042 class vtkLine;
00043 class vtkPolygon;
00044 class vtkQuad;
00045
00046 class VTK_FILTERING_EXPORT vtkPentagonalPrism : public vtkCell3D
00047 {
00048 public:
00049 static vtkPentagonalPrism *New();
00050 vtkTypeRevisionMacro(vtkPentagonalPrism,vtkCell3D);
00051 void PrintSelf(ostream& os, vtkIndent indent);
00052
00054
00055 virtual void GetEdgePoints(int edgeId, int* &pts);
00056 virtual void GetFacePoints(int faceId, int* &pts);
00058
00060
00061 int GetCellType() {return VTK_PENTAGONAL_PRISM;};
00062 int GetCellDimension() {return 3;};
00063 int GetNumberOfEdges() {return 15;};
00064 int GetNumberOfFaces() {return 7;};
00065 vtkCell *GetEdge(int edgeId);
00066 vtkCell *GetFace(int faceId);
00067 int CellBoundary(int subId, double pcoords[3], vtkIdList *pts);
00069
00070 int EvaluatePosition(double x[3], double* closestPoint,
00071 int& subId, double pcoords[3],
00072 double& dist2, double *weights);
00073 void EvaluateLocation(int& subId, double pcoords[3], double x[3],
00074 double *weights);
00075 int IntersectWithLine(double p1[3], double p2[3], double tol, double& t,
00076 double x[3], double pcoords[3], int& subId);
00077 int Triangulate(int index, vtkIdList *ptIds, vtkPoints *pts);
00078 void Derivatives(int subId, double pcoords[3], double *values,
00079 int dim, double *derivs);
00080 double *GetParametricCoords();
00081
00083 int GetParametricCenter(double pcoords[3]);
00084
00086
00088 static void InterpolationFunctions(double pcoords[3], double weights[10]);
00089
00090
00091 static void InterpolationDerivs(double pcoords[3], double derivs[30]);
00092
00093
00094
00095 virtual void InterpolateFunctions(double pcoords[3], double weights[10])
00096 {
00097 vtkPentagonalPrism::InterpolationFunctions(pcoords, weights);
00098 }
00099 virtual void InterpolateDerivs(double pcoords[3], double derivs[30])
00100 {
00101 vtkPentagonalPrism::InterpolationDerivs(pcoords, derivs);
00102 }
00104
00106
00108 static int *GetEdgeArray(int edgeId);
00109 static int *GetFaceArray(int faceId);
00111
00115 void JacobianInverse(double pcoords[3], double **inverse, double derivs[30]);
00116
00117 protected:
00118 vtkPentagonalPrism();
00119 ~vtkPentagonalPrism();
00120
00121 vtkLine *Line;
00122 vtkQuad *Quad;
00123 vtkPolygon *Polygon;
00124
00125 private:
00126 vtkPentagonalPrism(const vtkPentagonalPrism&);
00127 void operator=(const vtkPentagonalPrism&);
00128 };
00129
00130
00131 inline int vtkPentagonalPrism::GetParametricCenter(double pcoords[3])
00132 {
00133 pcoords[0] = pcoords[1] = 0.5;
00134 pcoords[2] = 0.5;
00135
00136 return 0;
00137 }
00138 #endif
00139
00140