VTK
|
00001 /*========================================================================= 00002 00003 Program: Visualization Toolkit 00004 Module: vtkHexagonalPrism.h 00005 00006 Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen 00007 All rights reserved. 00008 See Copyright.txt or http://www.kitware.com/Copyright.htm for details. 00009 00010 This software is distributed WITHOUT ANY WARRANTY; without even 00011 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 00012 PURPOSE. See the above copyright notice for more information. 00013 00014 =========================================================================*/ 00037 #ifndef __vtkHexagonalPrism_h 00038 #define __vtkHexagonalPrism_h 00039 00040 #include "vtkCommonDataModelModule.h" // For export macro 00041 #include "vtkCell3D.h" 00042 00043 class vtkLine; 00044 class vtkPolygon; 00045 class vtkQuad; 00046 00047 class VTKCOMMONDATAMODEL_EXPORT vtkHexagonalPrism : public vtkCell3D 00048 { 00049 public: 00050 static vtkHexagonalPrism *New(); 00051 vtkTypeMacro(vtkHexagonalPrism,vtkCell3D); 00052 void PrintSelf(ostream& os, vtkIndent indent); 00053 00055 00056 virtual void GetEdgePoints(int edgeId, int* &pts); 00057 virtual void GetFacePoints(int faceId, int* &pts); 00059 00061 00062 int GetCellType() {return VTK_HEXAGONAL_PRISM;}; 00063 int GetCellDimension() {return 3;}; 00064 int GetNumberOfEdges() {return 18;}; 00065 int GetNumberOfFaces() {return 8;}; 00066 vtkCell *GetEdge(int edgeId); 00067 vtkCell *GetFace(int faceId); 00068 int CellBoundary(int subId, double pcoords[3], vtkIdList *pts); 00070 00071 int EvaluatePosition(double x[3], double* closestPoint, 00072 int& subId, double pcoords[3], 00073 double& dist2, double *weights); 00074 void EvaluateLocation(int& subId, double pcoords[3], double x[3], 00075 double *weights); 00076 int IntersectWithLine(double p1[3], double p2[3], double tol, double& t, 00077 double x[3], double pcoords[3], int& subId); 00078 int Triangulate(int index, vtkIdList *ptIds, vtkPoints *pts); 00079 void Derivatives(int subId, double pcoords[3], double *values, 00080 int dim, double *derivs); 00081 double *GetParametricCoords(); 00082 00084 int GetParametricCenter(double pcoords[3]); 00085 00088 static void InterpolationFunctions(double pcoords[3], double weights[12]); 00091 static void InterpolationDerivs(double pcoords[3], double derivs[36]); 00093 00095 virtual void InterpolateFunctions(double pcoords[3], double weights[12]) 00096 { 00097 vtkHexagonalPrism::InterpolationFunctions(pcoords,weights); 00098 } 00099 virtual void InterpolateDerivs(double pcoords[3], double derivs[36]) 00100 { 00101 vtkHexagonalPrism::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[36]); 00116 00117 protected: 00118 vtkHexagonalPrism(); 00119 ~vtkHexagonalPrism(); 00120 00121 vtkLine *Line; 00122 vtkQuad *Quad; 00123 vtkPolygon *Polygon; 00124 00125 private: 00126 vtkHexagonalPrism(const vtkHexagonalPrism&); // Not implemented. 00127 void operator=(const vtkHexagonalPrism&); // Not implemented. 00128 }; 00129 00130 //---------------------------------------------------------------------------- 00131 inline int vtkHexagonalPrism::GetParametricCenter(double pcoords[3]) 00132 { 00133 pcoords[0] = pcoords[1] = 0.5; 00134 pcoords[2] = 0.5; 00135 return 0; 00136 } 00137 #endif 00138 00139