/usr/include/ITK-4.5/itkFEMElement3DC0LinearTriangular.h is in libinsighttoolkit4-dev 4.5.0-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 | /*=========================================================================
*
* Copyright Insight Software Consortium
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
#ifndef __itkFEMElement3DC0LinearTriangular_h
#define __itkFEMElement3DC0LinearTriangular_h
#include "itkFEMElementStd.h"
namespace itk
{
namespace fem
{
/**
* \class Element3DC0LinearTriangular
* \brief 3-noded, linear, C0 continuous finite element in 2D space.
*
* The ordering of the nodes is counter clockwise. That is the nodes
* should be defined in the following order:
*
* (0,1)
* 2
* *
* |\
* | \
* | \
* | \
* | \
* | \
* *------*
* 0 1
* (0,0) (0,1)
*
* \ingroup ITKFEM
*/
class Element3DC0LinearTriangular : public ElementStd<3, 3>
{
public:
/** Standard class typedefs. */
typedef Element3DC0LinearTriangular Self;
typedef ElementStd<3, 3> TemplatedParentClass;
typedef TemplatedParentClass Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Run-time type information (and related methods). */
itkTypeMacro(Element3DC0LinearTriangular, TemplatedParentClass);
// ////////////////////////////////////////////////////////////////////////
/*
* Methods related to numeric integration
*/
enum { DefaultIntegrationOrder = 1 };
/** Get the Integration point and weight */
virtual void GetIntegrationPointAndWeight(unsigned int i, VectorType & pt, Float & w, unsigned int order) const;
/** Get the number of integration points */
virtual unsigned int GetNumberOfIntegrationPoints(unsigned int order) const;
// ////////////////////////////////////////////////////////////////////////
/*
* Methods related to the geometry of an element
*/
/** Return the shape functions used to interpolate across the element */
virtual VectorType ShapeFunctions(const VectorType & pt) const;
/** Return the shape functions derivatives in the shapeD matrix */
virtual void ShapeFunctionDerivatives(const VectorType & pt, MatrixType & shapeD) const;
/** Convert from global to local coordinates */
virtual bool GetLocalFromGlobalCoordinates(const VectorType & globalPt, VectorType & localPt) const;
/** Compute the determinate of the Jacobian matrix */
virtual Float JacobianDeterminant(const VectorType & pt, const MatrixType *pJ = 0) const;
/** Compute the inverse of the Jacobian matrix */
virtual void JacobianInverse(const VectorType & pt, MatrixType & invJ, const MatrixType *pJ = 0) const;
/** Define the edges and nodes that correspond to the edges */
virtual void PopulateEdgeIds();
/**
* Normal of the triangle element
*/
void ComputeNormalDirection(const VectorType & v1, const VectorType & v2, const VectorType & v3,
VectorType & n) const;
/**
* Project the point x onto the plane containing the triangle element
*/
void GeneralizedProjectPoint(const VectorType & x, const VectorType & origin, const VectorType & normal,
VectorType & xproj) const;
/** Return the determinate of a 2x2 matrix */
itk::fem::Element::Float Determinant2x2(const VectorType & c1, const VectorType & c2) const;
/**
* Constants for integration rules.
*/
static const Float trigGaussRuleInfo[6][7][4];
/**
* Array that holds number of integration point for each order
* of numerical integration.
*/
static const unsigned int Nip[6];
protected:
virtual void PrintSelf(std::ostream& os, Indent indent) const;
};
}
} // end namespace itk::fem
#endif // #ifndef __itkFEMElement3DC0LinearTriangular_h
|