ITK  4.9.0
Insight Segmentation and Registration Toolkit
itkFEMElement2DC0QuadraticTriangular.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 itkFEMElement2DC0QuadraticTriangular_h
20 #define itkFEMElement2DC0QuadraticTriangular_h
21 
22 #include "itkFEMElementStd.h"
23 #include "ITKFEMExport.h"
24 
25 namespace itk
26 {
27 namespace fem
28 {
58 class ITKFEM_EXPORT Element2DC0QuadraticTriangular : public ElementStd<6, 2>
59 {
60 public:
61 
68 
71 
72  // ////////////////////////////////////////////////////////////////////////
77  enum { DefaultIntegrationOrder = 2 };
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 
106 protected:
107  virtual void PrintSelf(std::ostream& os, Indent indent) const ITK_OVERRIDE;
108 
109  virtual void PopulateEdgeIds(void) ITK_OVERRIDE;
110 
111 };
112 }
113 } // end namespace itk::fem
114 
115 #endif // #ifndef itkFEMElement2DC0QuadraticTriangular_h
Implements standard node management in the element classes.
6-noded, quadratic, C0 continuous finite element in 2D space that defines a triangle.
Control indentation during Print() invocation.
Definition: itkIndent.h:49