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 vtkTypeRevisionMacro(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 vtkPointLocator *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 vtkPointLocator *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