00001 
00002 
00003 
00004 
00005 
00006 
00007 
00008 
00009 
00010 
00011 
00012 
00013 
00014 
00035 #ifndef __vtkQuadraticTriangle_h
00036 #define __vtkQuadraticTriangle_h
00037 
00038 #include "vtkNonLinearCell.h"
00039 
00040 class vtkQuadraticEdge;
00041 class vtkTriangle;
00042 class vtkDoubleArray;
00043 
00044 class VTK_FILTERING_EXPORT vtkQuadraticTriangle : public vtkNonLinearCell
00045 {
00046 public:
00047   static vtkQuadraticTriangle *New();
00048   vtkTypeMacro(vtkQuadraticTriangle,vtkNonLinearCell);
00049   void PrintSelf(ostream& os, vtkIndent indent);
00050 
00052 
00054   int GetCellType() {return VTK_QUADRATIC_TRIANGLE;};
00055   int GetCellDimension() {return 2;}
00056   int GetNumberOfEdges() {return 3;}
00057   int GetNumberOfFaces() {return 0;}
00058   vtkCell *GetEdge(int edgeId);
00059   vtkCell *GetFace(int) {return 0;}
00061 
00062   int CellBoundary(int subId, double pcoords[3], vtkIdList *pts);
00063   void Contour(double value, vtkDataArray *cellScalars,
00064                vtkIncrementalPointLocator *locator, vtkCellArray *verts,
00065                vtkCellArray *lines, vtkCellArray *polys,
00066                vtkPointData *inPd, vtkPointData *outPd,
00067                vtkCellData *inCd, vtkIdType cellId, vtkCellData *outCd);
00068   int EvaluatePosition(double x[3], double* closestPoint,
00069                        int& subId, double pcoords[3],
00070                        double& dist2, double *weights);
00071   void EvaluateLocation(int& subId, double pcoords[3], double x[3],
00072                         double *weights);
00073   int Triangulate(int index, vtkIdList *ptIds, vtkPoints *pts);
00074   void Derivatives(int subId, double pcoords[3], double *values,
00075                    int dim, double *derivs);
00076   virtual double *GetParametricCoords();
00077 
00079 
00082   void Clip(double value, vtkDataArray *cellScalars,
00083             vtkIncrementalPointLocator *locator, vtkCellArray *polys,
00084             vtkPointData *inPd, vtkPointData *outPd,
00085             vtkCellData *inCd, vtkIdType cellId, vtkCellData *outCd,
00086             int insideOut);
00088 
00090 
00092   int IntersectWithLine(double p1[3], double p2[3], double tol, double& t,
00093                         double x[3], double pcoords[3], int& subId);
00095 
00096 
00099   int GetParametricCenter(double pcoords[3]);
00100 
00103   double GetParametricDistance(double pcoords[3]);
00104 
00106 
00108   static void InterpolationFunctions(double pcoords[3], double weights[6]);
00109   
00110   
00111   static void InterpolationDerivs(double pcoords[3], double derivs[12]);
00112   
00113   
00114   
00115   virtual void InterpolateFunctions(double pcoords[3], double weights[6])
00116     {
00117     vtkQuadraticTriangle::InterpolationFunctions(pcoords,weights);
00118     }
00119   virtual void InterpolateDerivs(double pcoords[3], double derivs[12])
00120     {
00121     vtkQuadraticTriangle::InterpolationDerivs(pcoords,derivs);
00122     }
00124 
00125 protected:
00126   vtkQuadraticTriangle();
00127   ~vtkQuadraticTriangle();
00128 
00129   vtkQuadraticEdge *Edge;
00130   vtkTriangle      *Face;
00131   vtkDoubleArray    *Scalars; 
00132 
00133 private:
00134   vtkQuadraticTriangle(const vtkQuadraticTriangle&);  
00135   void operator=(const vtkQuadraticTriangle&);  
00136 };
00137 
00138 inline int vtkQuadraticTriangle::GetParametricCenter(double pcoords[3])
00139 {
00140   pcoords[0] = pcoords[1] = 1./3;
00141   pcoords[2] = 0.0;
00142   return 0;
00143 }
00144 
00145 
00146 #endif
00147 
00148