Main Page   Groups   Namespace List   Class Hierarchy   Alphabetical List   Compound List   File List   Namespace Members   Compound Members   File Members   Concepts

itkFEMElement2DC0LinearTriangular.h

Go to the documentation of this file.
00001 /*========================================================================= 00002 00003 Program: Insight Segmentation & Registration Toolkit 00004 Module: $RCSfile: itkFEMElement2DC0LinearTriangular.h,v $ 00005 Language: C++ 00006 Date: $Date: 2003/09/10 14:29:39 $ 00007 Version: $Revision: 1.6 $ 00008 00009 Copyright (c) Insight Software 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 __itkFEMElement2DC0LinearTriangular_h 00019 #define __itkFEMElement2DC0LinearTriangular_h 00020 00021 #include "itkFEMElementStd.h" 00022 00023 namespace itk { 00024 namespace fem { 00025 00026 00027 00028 00033 class Element2DC0LinearTriangular : public ElementStd<3,2> 00034 { 00035 typedef ElementStd<3,2> TemplatedParentClass; 00036 FEM_ABSTRACT_CLASS( Element2DC0LinearTriangular, TemplatedParentClass ) 00037 public: 00038 00039 00041 /* 00042 * Methods related to numeric integration 00043 */ 00044 00045 enum { DefaultIntegrationOrder = 1 }; 00046 00047 virtual void GetIntegrationPointAndWeight(unsigned int i, VectorType& pt, Float& w, unsigned int order) const; 00048 00049 virtual unsigned int GetNumberOfIntegrationPoints(unsigned int order) const; 00050 00051 00052 00054 /* 00055 * Methods related to the geometry of an element 00056 */ 00057 00058 virtual VectorType ShapeFunctions( const VectorType& pt ) const; 00059 00060 virtual void ShapeFunctionDerivatives( const VectorType& pt, MatrixType& shapeD ) const; 00061 00062 // FIXME: Write a proper implementation 00063 virtual bool GetLocalFromGlobalCoordinates( const VectorType& globalPt , VectorType& localPt) const; 00064 00065 // Since the Jacobian is not quadratic, we need to provide our 00066 // own implementation of calculating the determinant and inverse. 00067 virtual Float JacobianDeterminant( const VectorType& pt, const MatrixType* pJ = 0 ) const; 00068 virtual void JacobianInverse( const VectorType& pt, MatrixType& invJ, const MatrixType* pJ = 0 ) const; 00069 00073 #ifdef FEM_BUILD_VISUALIZATION 00074 void Draw(CDC* pDC, Solution::ConstPointer sol) const; 00075 #endif 00076 00080 static const Float trigGaussRuleInfo[6][7][4]; 00081 00086 static const unsigned int Nip[6]; 00087 00088 }; 00089 00090 00091 00092 00093 }} // end namespace itk::fem 00094 00095 #endif // #ifndef __itkFEMElement2DC0LinearTriangular_h

Generated at Sat Mar 31 02:18:09 2007 for ITK by doxygen 1.3.8 written by Dimitri van Heesch, © 1997-2000