ITK
4.1.0
Insight Segmentation and Registration Toolkit
|
00001 /*========================================================================= 00002 * 00003 * Copyright Insight Software Consortium 00004 * 00005 * Licensed under the Apache License, Version 2.0 (the "License"); 00006 * you may not use this file except in compliance with the License. 00007 * You may obtain a copy of the License at 00008 * 00009 * http://www.apache.org/licenses/LICENSE-2.0.txt 00010 * 00011 * Unless required by applicable law or agreed to in writing, software 00012 * distributed under the License is distributed on an "AS IS" BASIS, 00013 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 00014 * See the License for the specific language governing permissions and 00015 * limitations under the License. 00016 * 00017 *=========================================================================*/ 00018 00019 #ifndef __itkFEMElement3DC0LinearTriangular_h 00020 #define __itkFEMElement3DC0LinearTriangular_h 00021 00022 #include "itkFEMElementStd.h" 00023 00024 namespace itk 00025 { 00026 namespace fem 00027 { 00051 class Element3DC0LinearTriangular : public ElementStd<3, 3> 00052 { 00053 public: 00055 typedef Element3DC0LinearTriangular Self; 00056 typedef ElementStd<3, 3> TemplatedParentClass; 00057 typedef TemplatedParentClass Superclass; 00058 typedef SmartPointer<Self> Pointer; 00059 typedef SmartPointer<const Self> ConstPointer; 00060 00062 itkTypeMacro(Element3DC0LinearTriangular, TemplatedParentClass); 00063 00064 // //////////////////////////////////////////////////////////////////////// 00065 /* 00066 * Methods related to numeric integration 00067 */ 00068 00069 enum { DefaultIntegrationOrder = 1 }; 00070 00072 virtual void GetIntegrationPointAndWeight(unsigned int i, VectorType & pt, Float & w, unsigned int order) const; 00073 00075 virtual unsigned int GetNumberOfIntegrationPoints(unsigned int order) const; 00076 00077 // //////////////////////////////////////////////////////////////////////// 00078 /* 00079 * Methods related to the geometry of an element 00080 */ 00081 00083 virtual VectorType ShapeFunctions(const VectorType & pt) const; 00084 00086 virtual void ShapeFunctionDerivatives(const VectorType & pt, MatrixType & shapeD) const; 00087 00089 virtual bool GetLocalFromGlobalCoordinates(const VectorType & globalPt, VectorType & localPt) const; 00090 00092 virtual Float JacobianDeterminant(const VectorType & pt, const MatrixType *pJ = 0) const; 00093 00095 virtual void JacobianInverse(const VectorType & pt, MatrixType & invJ, const MatrixType *pJ = 0) const; 00096 00098 virtual void PopulateEdgeIds(); 00099 00103 void ComputeNormalDirection(const VectorType & v1, const VectorType & v2, const VectorType & v3, 00104 VectorType & n) const; 00105 00109 void GeneralizedProjectPoint(const VectorType & x, const VectorType & origin, const VectorType & normal, 00110 VectorType & xproj) const; 00111 00113 itk::fem::Element::Float Determinant2x2(const VectorType & c1, const VectorType & c2) const; 00114 00118 static const Float trigGaussRuleInfo[6][7][4]; 00119 00124 static const unsigned int Nip[6]; 00125 protected: 00126 virtual void PrintSelf(std::ostream& os, Indent indent) const; 00127 00128 }; 00129 } 00130 } // end namespace itk::fem 00131 00132 #endif // #ifndef __itkFEMElement3DC0LinearTriangular_h 00133