Nektar++
NodalTetEvenlySpaced.h
Go to the documentation of this file.
1///////////////////////////////////////////////////////////////////////////////
2//
3// File: NodalTetEvenlySpaced.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// Permission is hereby granted, free of charge, to any person obtaining a
14// copy of this software and associated documentation files (the "Software"),
15// to deal in the Software without restriction, including without limitation
16// the rights to use, copy, modify, merge, publish, distribute, sublicense,
17// and/or sell copies of the Software, and to permit persons to whom the
18// Software is furnished to do so, subject to the following conditions:
19//
20// The above copyright notice and this permission notice shall be included
21// in all copies or substantial portions of the Software.
22//
23// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
24// OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
25// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
26// THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
27// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
28// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
29// DEALINGS IN THE SOFTWARE.
30//
31// Description: Header file of 2D Nodal Tetrahedron Evenly Spaced Points
32//
33///////////////////////////////////////////////////////////////////////////////
34
35#ifndef NODALTETEVENLYSPACED_H
36#define NODALTETEVENLYSPACED_H
37
38#include <memory>
39
42
44{
45class NodalTetEvenlySpaced : public Points<NekDouble>
46{
47public:
49 {
50 }
51
53 {
54 }
55
56 LIB_UTILITIES_EXPORT static std::shared_ptr<PointsBaseType> Create(
57 const PointsKey &key);
58
59protected:
60 const MatrixSharedPtrType v_GetI(const PointsKey &pkey) override
61 {
62 ASSERTL0(pkey.GetPointsDim() == 3,
63 "NodalTetEvenlySpaced Points can only interp to other "
64 "3d point distributions");
66 PointsManager()[pkey]->GetPoints(x, y, z);
67 return GetI(x, y, z);
68 }
69
73 const Array<OneD, const NekDouble> &z) override
74 {
75 size_t numpoints = x.size();
76 size_t np = GetTotNumPoints();
77
78 Array<OneD, NekDouble> interp(GetTotNumPoints() * numpoints);
79 CalculateInterpMatrix(x, y, z, interp);
80
81 NekDouble *d = interp.data();
82 return MemoryManager<NekMatrix<NekDouble>>::AllocateSharedPtr(numpoints,
83 np, d);
84 }
85
86private:
87 static bool initPointsManager[];
88
89 std::shared_ptr<NodalUtilTetrahedron> m_util;
90
93
95
96 void v_CalculatePoints() final;
97 void v_CalculateWeights() final;
98 void v_CalculateDerivMatrix() final;
99
100 void CalculateInterpMatrix(const Array<OneD, const NekDouble> &xi,
101 const Array<OneD, const NekDouble> &yi,
102 const Array<OneD, const NekDouble> &zi,
103 Array<OneD, NekDouble> &interp);
104}; // end of NodalTetEvenlySpaced
105} // namespace Nektar::LibUtilities
106
107#endif // NODALTETEVENLYSPACED_H
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:208
#define LIB_UTILITIES_EXPORT
std::shared_ptr< NodalUtilTetrahedron > m_util
const MatrixSharedPtrType v_GetI(const PointsKey &pkey) override
void CalculateInterpMatrix(const Array< OneD, const NekDouble > &xi, const Array< OneD, const NekDouble > &yi, const Array< OneD, const NekDouble > &zi, Array< OneD, NekDouble > &interp)
const MatrixSharedPtrType v_GetI(const Array< OneD, const NekDouble > &x, const Array< OneD, const NekDouble > &y, const Array< OneD, const NekDouble > &z) override
NodalTetEvenlySpaced(const NodalTetEvenlySpaced &points)=delete
static std::shared_ptr< PointsBaseType > Create(const PointsKey &key)
Stores a set of points of datatype DataT, defined by a PointKey.
Definition: Points.h:235
std::shared_ptr< NekMatrix< DataType > > MatrixSharedPtrType
Definition: Points.h:238
const MatrixSharedPtrType GetI(const PointsKey &key)
Definition: Points.h:317
Defines a specification for a set of points.
Definition: Points.h:50
size_t GetPointsDim() const
Definition: Points.h:128
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
PointsManagerT & PointsManager(void)
std::vector< double > z(NPUPPER)
std::vector< double > d(NPUPPER *NPUPPER)
double NekDouble