ITK  4.9.0
Insight Segmentation and Registration Toolkit
itkFEMElement3DC0LinearTriangular.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 itkFEMElement3DC0LinearTriangular_h
20 #define itkFEMElement3DC0LinearTriangular_h
21 
22 #include "itkFEMElementStd.h"
23 #include "ITKFEMExport.h"
24 
25 namespace itk
26 {
27 namespace fem
28 {
52 class ITKFEM_EXPORT Element3DC0LinearTriangular : public ElementStd<3, 3>
53 {
54 public:
61 
64 
65 // ////////////////////////////////////////////////////////////////////////
66 /*
67  * Methods related to numeric integration
68  */
69 
70  enum { DefaultIntegrationOrder = 1 };
71 
73  virtual void GetIntegrationPointAndWeight(unsigned int i, VectorType & pt, Float & w, unsigned int order) const ITK_OVERRIDE;
74 
76  virtual unsigned int GetNumberOfIntegrationPoints(unsigned int order) const ITK_OVERRIDE;
77 
78  // ////////////////////////////////////////////////////////////////////////
79  /*
80  * Methods related to the geometry of an element
81  */
82 
84  virtual VectorType ShapeFunctions(const VectorType & pt) const ITK_OVERRIDE;
85 
87  virtual void ShapeFunctionDerivatives(const VectorType & pt, MatrixType & shapeD) const ITK_OVERRIDE;
88 
90  virtual bool GetLocalFromGlobalCoordinates(const VectorType & globalPt, VectorType & localPt) const ITK_OVERRIDE;
91 
93  virtual Float JacobianDeterminant(const VectorType & pt, const MatrixType *pJ = ITK_NULLPTR) const ITK_OVERRIDE;
94 
96  virtual void JacobianInverse(const VectorType & pt, MatrixType & invJ, const MatrixType *pJ = ITK_NULLPTR) const ITK_OVERRIDE;
97 
99  virtual void PopulateEdgeIds() ITK_OVERRIDE;
100 
104  void ComputeNormalDirection(const VectorType & v1, const VectorType & v2, const VectorType & v3,
105  VectorType & n) const;
106 
110  void GeneralizedProjectPoint(const VectorType & x, const VectorType & origin, const VectorType & normal,
111  VectorType & xproj) const;
112 
114  itk::fem::Element::Float Determinant2x2(const VectorType & c1, const VectorType & c2) const;
115 
119  static const Float trigGaussRuleInfo[6][7][4];
120 
125  static const unsigned int Nip[6];
126 
127 protected:
128  virtual void PrintSelf(std::ostream& os, Indent indent) const ITK_OVERRIDE;
129 
130 };
131 }
132 } // end namespace itk::fem
133 
134 #endif // #ifndef itkFEMElement3DC0LinearTriangular_h
Implements standard node management in the element classes.
Abstract base element class.
3-noded, linear, C0 continuous finite element in 2D space.
Control indentation during Print() invocation.
Definition: itkIndent.h:49