ITK  4.2.0
Insight Segmentation and Registration Toolkit
itkFEMElement2DC0LinearTriangular.h
Go to the documentation of this file.
1 /*=========================================================================
2  *
3  * Copyright Insight Software Consortium
4  *
5  * Licensed under the Apache License, Version 2.0 (the "License");
6  * you may not use this file except in compliance with the License.
7  * You may obtain a copy of the License at
8  *
9  * http://www.apache.org/licenses/LICENSE-2.0.txt
10  *
11  * Unless required by applicable law or agreed to in writing, software
12  * distributed under the License is distributed on an "AS IS" BASIS,
13  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14  * See the License for the specific language governing permissions and
15  * limitations under the License.
16  *
17  *=========================================================================*/
18 
19 #ifndef __itkFEMElement2DC0LinearTriangular_h
20 #define __itkFEMElement2DC0LinearTriangular_h
21 
22 #include "itkFEMElementStd.h"
23 
24 namespace itk
25 {
26 namespace fem
27 {
58 {
59 public:
60 
67 
70 
71  // ////////////////////////////////////////////////////////////////////////
77 
79  virtual void GetIntegrationPointAndWeight(unsigned int i, VectorType & pt, Float & w, unsigned int order) const;
80 
82  virtual unsigned int GetNumberOfIntegrationPoints(unsigned int order) const;
83 
84  // ////////////////////////////////////////////////////////////////////////
90  virtual VectorType ShapeFunctions(const VectorType & pt) const;
91 
93  virtual void ShapeFunctionDerivatives(const VectorType & pt, MatrixType & shapeD) const;
94 
96  virtual bool GetLocalFromGlobalCoordinates(const VectorType & globalPt, VectorType & localPt) const;
97 
98  // Since the Jacobian is not quadratic, we need to provide our
99  // own implementation of calculating the determinant and inverse.
100  virtual Float JacobianDeterminant(const VectorType & pt, const MatrixType *pJ = 0) const;
101 
103  virtual void JacobianInverse(const VectorType & pt, MatrixType & invJ, const MatrixType *pJ = 0) const;
104 
108  static const Float trigGaussRuleInfo[6][7][4];
109 
114  static const unsigned int Nip[6];
115 protected:
116  virtual void PopulateEdgeIds(void);
117 
118  virtual void PrintSelf(std::ostream& os, Indent indent) const;
119 
120 };
121 }
122 } // end namespace itk::fem
123 
124 #endif // #ifndef __itkFEMElement2DC0LinearTriangular_h
125