ITK  4.9.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 #include "ITKFEMExport.h"
24 
25 namespace itk
26 {
27 namespace fem
28 {
58 class ITKFEM_EXPORT Element2DC0LinearTriangular : public ElementStd<3, 2>
59 {
60 public:
61 
68 
71 
72  // ////////////////////////////////////////////////////////////////////////
77  enum { DefaultIntegrationOrder = 1 };
78 
80  virtual void GetIntegrationPointAndWeight(unsigned int i, VectorType & pt, Float & w, unsigned int order) const ITK_OVERRIDE;
81 
83  virtual unsigned int GetNumberOfIntegrationPoints(unsigned int order) const ITK_OVERRIDE;
84 
85  // ////////////////////////////////////////////////////////////////////////
91  virtual VectorType ShapeFunctions(const VectorType & pt) const ITK_OVERRIDE;
92 
94  virtual void ShapeFunctionDerivatives(const VectorType & pt, MatrixType & shapeD) const ITK_OVERRIDE;
95 
97  virtual bool GetLocalFromGlobalCoordinates(const VectorType & globalPt, VectorType & localPt) const ITK_OVERRIDE;
98 
99  // Since the Jacobian is not quadratic, we need to provide our
100  // own implementation of calculating the determinant and inverse.
101  virtual Float JacobianDeterminant(const VectorType & pt, const MatrixType *pJ = ITK_NULLPTR) const ITK_OVERRIDE;
102 
104  virtual void JacobianInverse(const VectorType & pt, MatrixType & invJ, const MatrixType *pJ = ITK_NULLPTR) const ITK_OVERRIDE;
105 
109  static const Float trigGaussRuleInfo[6][7][4];
110 
115  static const unsigned int Nip[6];
116 
117 protected:
118  virtual void PopulateEdgeIds(void) ITK_OVERRIDE;
119 
120  virtual void PrintSelf(std::ostream& os, Indent indent) const ITK_OVERRIDE;
121 
122 };
123 }
124 } // end namespace itk::fem
125 
126 #endif // #ifndef itkFEMElement2DC0LinearTriangular_h
3-noded, linear, C0 continuous finite element in 2D space.The ordering of the nodes is counter clockw...
Implements standard node management in the element classes.
Control indentation during Print() invocation.
Definition: itkIndent.h:49