Nektar++
Main Page
Related Pages
Namespaces
Classes
Files
File List
File Members
•
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
library
LibUtilities
Foundations
FourierPoints.h
Go to the documentation of this file.
1
///////////////////////////////////////////////////////////////////////////////
2
//
3
// File FourierPoints.h
4
//
5
// For more information, please see: http://www.nektar.info
6
//
7
// The MIT License
8
//
9
// Copyright (c) 2006 Division of Applied Mathematics, Brown University (USA),
10
// Department of Aeronautics, Imperial College London (UK), and Scientific
11
// Computing and Imaging Institute, University of Utah (USA).
12
//
13
// License for the specific language governing rights and limitations under
14
// Permission is hereby granted, free of charge, to any person obtaining a
15
// copy of this software and associated documentation files (the "Software"),
16
// to deal in the Software without restriction, including without limitation
17
// the rights to use, copy, modify, merge, publish, distribute, sublicense,
18
// and/or sell copies of the Software, and to permit persons to whom the
19
// Software is furnished to do so, subject to the following conditions:
20
//
21
// The above copyright notice and this permission notice shall be included
22
// in all copies or substantial portions of the Software.
23
//
24
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
25
// OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
26
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
27
// THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
28
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
29
// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
30
// DEALINGS IN THE SOFTWARE.
31
//
32
// Description: Header file of 1D Evenly-Spaced Point Definitions
33
//
34
///////////////////////////////////////////////////////////////////////////////
35
36
#ifndef FOURIERPOINTS_H
37
#define FOURIERPOINTS_H
38
39
#include <
LibUtilities/Foundations/FoundationsFwd.hpp
>
40
#include <
LibUtilities/Foundations/Points.h
>
41
#include <boost/shared_ptr.hpp>
42
#include <
LibUtilities/Foundations/Foundations.hpp
>
43
//#include <LibUtilities/BasicUtils/BasicUtilsFwd.hpp> // for NekManager
44
#include <
LibUtilities/BasicUtils/NekManager.hpp
>
// for NekManager
45
46
47
namespace
Nektar
48
{
49
namespace
LibUtilities
50
{
51
class
FourierPoints
:
public
Points
<NekDouble>
52
{
53
public
:
54
virtual
~FourierPoints
()
55
{
56
}
57
58
LIB_UTILITIES_EXPORT
static
boost::shared_ptr< PointsBaseType >
Create
(
const
PointsKey
&key);
59
LIB_UTILITIES_EXPORT
boost::shared_ptr< NekMatrix<NekDouble> >
CreateMatrix
(
const
PointsKey
&pkey);
60
61
LIB_UTILITIES_EXPORT
const
MatrixSharedPtrType
GetI
(
const
PointsKey
&pkey);
62
LIB_UTILITIES_EXPORT
const
MatrixSharedPtrType
GetI
(
const
Array<OneD, const NekDouble>& x);
63
LIB_UTILITIES_EXPORT
const
MatrixSharedPtrType
GetI
(
unsigned
int
numpoints,
const
Array<OneD, const NekDouble>& x);
64
65
FourierPoints
(
const
PointsKey
&key):
PointsBaseType
(key)
66
{
67
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussGaussLegendre
),
68
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
69
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussRadauMLegendre
),
70
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
71
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussRadauPLegendre
),
72
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
73
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussLobattoLegendre
),
74
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
75
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussGaussChebyshev
),
76
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
77
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussRadauMChebyshev
),
78
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
79
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussRadauPChebyshev
),
80
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
81
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussLobattoChebyshev
),
82
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
83
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussRadauMAlpha0Beta1
),
84
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
85
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussRadauMAlpha0Beta2
),
86
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
87
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussRadauMAlpha1Beta0
),
88
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
89
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eGaussRadauMAlpha2Beta0
),
90
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
91
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
ePolyEvenlySpaced
),
92
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
93
m_InterpManager
.
RegisterCreator
(
PointsKey
(0,
eFourierEvenlySpaced
),
94
boost::bind(&
FourierPoints::CreateMatrix
,
this
, _1));
95
}
96
97
private
:
98
/// Default constructor should not be called except by Create method.
99
FourierPoints
();
100
101
/// Copy constructor should not be called.
102
FourierPoints
(
const
FourierPoints
&points);
103
104
void
CalculatePoints
();
105
void
CalculateWeights
();
106
void
CalculateDerivMatrix
();
107
108
void
CalculateInterpMatrix
(
unsigned
int
npts
,
const
Array<OneD, const NekDouble>& xpoints,
109
Array<OneD, NekDouble>& interp);
110
NekDouble
PeriodicSincFunction
(
const
NekDouble
x,
const
NekDouble
h);
111
};
// class FourierPoints
112
}
// end of namespace
113
}
// end of namespace
114
115
116
#endif //FOURIERPOINTS_H
Generated on Wed Mar 4 2015 16:33:14 for Nektar++ by
1.8.1.2