00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00050 #ifndef __vtkQuadraticTetra_h
00051 #define __vtkQuadraticTetra_h
00052
00053 #include "vtkNonLinearCell.h"
00054
00055 class vtkPolyData;
00056 class vtkQuadraticEdge;
00057 class vtkQuadraticTriangle;
00058 class vtkTetra;
00059
00060 class VTK_COMMON_EXPORT vtkQuadraticTetra : public vtkNonLinearCell
00061 {
00062 public:
00063 static vtkQuadraticTetra *New();
00064 vtkTypeRevisionMacro(vtkQuadraticTetra,vtkNonLinearCell);
00065
00067
00069 int GetCellType() {return VTK_QUADRATIC_TETRA;}
00070 int GetCellDimension() {return 3;}
00071 int GetNumberOfEdges() {return 6;}
00072 int GetNumberOfFaces() {return 4;}
00073 vtkCell *GetEdge(int);
00074 vtkCell *GetFace(int);
00076
00077 int CellBoundary(int subId, float pcoords[3], vtkIdList *pts);
00078 void Contour(float value, vtkDataArray *cellScalars,
00079 vtkPointLocator *locator, vtkCellArray *verts,
00080 vtkCellArray *lines, vtkCellArray *polys,
00081 vtkPointData *inPd, vtkPointData *outPd,
00082 vtkCellData *inCd, vtkIdType cellId, vtkCellData *outCd);
00083 int EvaluatePosition(float x[3], float* closestPoint,
00084 int& subId, float pcoords[3],
00085 float& dist2, float *weights);
00086 void EvaluateLocation(int& subId, float pcoords[3], float x[3],
00087 float *weights);
00088 int Triangulate(int index, vtkIdList *ptIds, vtkPoints *pts);
00089 void Derivatives(int subId, float pcoords[3], float *values,
00090 int dim, float *derivs);
00091
00093
00095 void Clip(float value, vtkDataArray *cellScalars,
00096 vtkPointLocator *locator, vtkCellArray *tetras,
00097 vtkPointData *inPd, vtkPointData *outPd,
00098 vtkCellData *inCd, vtkIdType cellId, vtkCellData *outCd,
00099 int insideOut);
00101
00103
00105 int IntersectWithLine(float p1[3], float p2[3], float tol, float& t,
00106 float x[3], float pcoords[3], int& subId);
00108
00109
00111 int GetParametricCenter(float pcoords[3]);
00112
00114
00115 static void InterpolationFunctions(float pcoords[3], float weights[3]);
00116 static void InterpolationDerivs(float pcoords[3], float derivs[3]);
00118
00122 void JacobianInverse(float pcoords[3], double **inverse, float derivs[30]);
00123
00124 protected:
00125 vtkQuadraticTetra();
00126 ~vtkQuadraticTetra();
00127
00128 vtkQuadraticEdge *Edge;
00129 vtkQuadraticTriangle *Face;
00130 vtkTetra *Tetra;
00131 vtkFloatArray *Scalars;
00132
00133 private:
00134 vtkQuadraticTetra(const vtkQuadraticTetra&);
00135 void operator=(const vtkQuadraticTetra&);
00136 };
00137
00138 #endif
00139
00140