ITK  4.2.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 
24 namespace itk
25 {
26 namespace fem
27 {
52 {
53 public:
60 
63 
64 // ////////////////////////////////////////////////////////////////////////
65 /*
66  * Methods related to numeric integration
67  */
68 
70 
72  virtual void GetIntegrationPointAndWeight(unsigned int i, VectorType & pt, Float & w, unsigned int order) const;
73 
75  virtual unsigned int GetNumberOfIntegrationPoints(unsigned int order) const;
76 
77  // ////////////////////////////////////////////////////////////////////////
78  /*
79  * Methods related to the geometry of an element
80  */
81 
83  virtual VectorType ShapeFunctions(const VectorType & pt) const;
84 
86  virtual void ShapeFunctionDerivatives(const VectorType & pt, MatrixType & shapeD) const;
87 
89  virtual bool GetLocalFromGlobalCoordinates(const VectorType & globalPt, VectorType & localPt) const;
90 
92  virtual Float JacobianDeterminant(const VectorType & pt, const MatrixType *pJ = 0) const;
93 
95  virtual void JacobianInverse(const VectorType & pt, MatrixType & invJ, const MatrixType *pJ = 0) const;
96 
98  virtual void PopulateEdgeIds();
99 
103  void ComputeNormalDirection(const VectorType & v1, const VectorType & v2, const VectorType & v3,
104  VectorType & n) const;
105 
109  void GeneralizedProjectPoint(const VectorType & x, const VectorType & origin, const VectorType & normal,
110  VectorType & xproj) const;
111 
113  itk::fem::Element::Float Determinant2x2(const VectorType & c1, const VectorType & c2) const;
114 
118  static const Float trigGaussRuleInfo[6][7][4];
119 
124  static const unsigned int Nip[6];
125 protected:
126  virtual void PrintSelf(std::ostream& os, Indent indent) const;
127 
128 };
129 }
130 } // end namespace itk::fem
131 
132 #endif // #ifndef __itkFEMElement3DC0LinearTriangular_h
133