VTK
dox/Common/DataModel/vtkHexagonalPrism.h
Go to the documentation of this file.
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