00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018 #ifndef __itkFEMElement3DC0LinearTetrahedron_h
00019 #define __itkFEMElement3DC0LinearTetrahedron_h
00020
00021
00022 #include "itkFEMElementStd.h"
00023
00024 #include <vnl/vnl_matrix.h>
00025 #include <vnl/algo/vnl_matrix_inverse.h>
00026
00027 namespace itk {
00028 namespace fem {
00029
00030
00031
00032
00037 class Element3DC0LinearTetrahedron : public ElementStd<4,3>
00038 {
00039 typedef ElementStd<4,3> TemplatedParentClass;
00040 FEM_ABSTRACT_CLASS( Element3DC0LinearTetrahedron, TemplatedParentClass )
00041 public:
00042
00043
00045
00046
00047
00048
00049 virtual void GetIntegrationPointAndWeight(unsigned int i, VectorType& pt, Float& w, unsigned int order) const;
00050
00051 virtual unsigned int GetNumberOfIntegrationPoints(unsigned int order) const;
00052
00053
00054
00056
00057
00058
00059
00060 virtual VectorType ShapeFunctions( const VectorType& pt ) const;
00061
00062 virtual void ShapeFunctionDerivatives( const VectorType& pt, MatrixType& shapeD ) const;
00063
00064 virtual bool GetLocalFromGlobalCoordinates( const VectorType& globalPt, VectorType& localPt ) const;
00065
00069 #ifdef FEM_BUILD_VISUALIZATION
00070 void Draw(CDC* pDC, Solution::ConstPointer sol) const;
00071 #endif
00072
00073
00074 };
00075
00076
00077
00078
00079 }}
00080
00081 #endif