00001 /*========================================================================= 00002 00003 Program: Insight Segmentation & Registration Toolkit 00004 Module: $RCSfile: itkFEMElement3DC0LinearTetrahedron.h,v $ 00005 Language: C++ 00006 Date: $Date: 2004-07-22 15:52:10 $ 00007 Version: $Revision: 1.10 $ 00008 00009 Copyright (c) 2002 Insight Consortium. All rights reserved. 00010 See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details. 00011 00012 This software is distributed WITHOUT ANY WARRANTY; without even 00013 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 00014 PURPOSE. See the above copyright notices for more information. 00015 00016 =========================================================================*/ 00017 00018 #ifndef __itkFEMElement3DC0LinearTetrahedron_h 00019 #define __itkFEMElement3DC0LinearTetrahedron_h 00020 00021 00022 #include "itkFEMElementStd.h" 00023 //to make some checks in GetLocalFromGlobalCoordinates 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 * Methods related to numeric integration 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 * Methods related to the geometry of an element 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 }} // end namespace itk::fem 00080 00081 #endif // #ifndef __itkFEMElement3DC0LinearTetrahedron_h 00082