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
|
/*=========================================================================
Program: VMTK
Module: $RCSfile: vtkvmtkFEShapeFunctions.h,v $
Language: C++
Date: $Date: 2006/04/06 16:46:44 $
Version: $Revision: 1.4 $
Copyright (c) Luca Antiga, David Steinman. All rights reserved.
See LICENCE file for details.
Portions of this code are covered under the VTK copyright.
See VTKCopyright.txt or http://www.kitware.com/VTKCopyright.htm
for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
// .NAME vtkvmtkFEShapeFunctions - ..
// .SECTION Description
// ..
#ifndef __vtkvmtkFEShapeFunctions_h
#define __vtkvmtkFEShapeFunctions_h
#include "vtkObject.h"
#include "vtkvmtkWin32Header.h"
#include "vtkCell.h"
#include "vtkDoubleArray.h"
class VTK_VMTK_DIFFERENTIAL_GEOMETRY_EXPORT vtkvmtkFEShapeFunctions : public vtkObject
{
public:
vtkTypeRevisionMacro(vtkvmtkFEShapeFunctions,vtkObject);
static vtkvmtkFEShapeFunctions* New();
void Initialize(vtkCell* cell, vtkDoubleArray* pcoords);
double GetPhi(vtkIdType id, vtkIdType i)
{ return this->Phi->GetValue(id*this->NumberOfCellPoints+i); }
double* GetDPhi(vtkIdType id, vtkIdType i)
{ return this->DPhi->GetTuple(id*this->NumberOfCellPoints+i); }
void GetDPhi(vtkIdType id, vtkIdType i, double* dphi)
{ this->DPhi->GetTuple(id*this->NumberOfCellPoints+i,dphi); }
double GetDPhi(vtkIdType id, vtkIdType i, int c)
{ return this->DPhi->GetComponent(id*this->NumberOfCellPoints+i,c); }
double GetJacobian(vtkIdType i)
{ return this->Jacobians->GetValue(i); }
static void GetInterpolationFunctions(vtkCell* cell, double* pcoords, double* sf);
static void GetInterpolationDerivs(vtkCell* cell, double* pcoords, double* derivs);
static double ComputeJacobian(vtkCell* cell, double* pcoords);
protected:
vtkvmtkFEShapeFunctions();
~vtkvmtkFEShapeFunctions();
static void ComputeInverseJacobianMatrix2D(vtkCell* cell, double* pcoords, double inverseJacobianMatrix[2][3]);
static void ComputeInverseJacobianMatrix3D(vtkCell* cell, double* pcoords, double inverseJacobianMatrix[3][3]);
vtkDoubleArray* Phi;
vtkDoubleArray* DPhi;
vtkDoubleArray* Jacobians;
vtkIdType NumberOfCellPoints;
private:
vtkvmtkFEShapeFunctions(const vtkvmtkFEShapeFunctions&); // Not implemented.
void operator=(const vtkvmtkFEShapeFunctions&); // Not implemented.
};
#endif
|
