chombo-discharge
Loading...
Searching...
No Matches
CD_MultiIndex.H
Go to the documentation of this file.
1/*
2 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
3 *
4 * SPDX-License-Identifier: GPL-3.0-or-later
5 */
6
13#ifndef CD_MULTIINDEX_H
14#define CD_MULTIINDEX_H
15
16// Std includes
17#include <map>
18
19// Chombo includes
20#include <IntVect.H>
21#include <RealVect.H>
22
23// Our includes
24#include <CD_NamespaceHeader.H>
25
31{
32public:
36 MultiIndex() = delete;
37
42 MultiIndex(const MultiIndex& a_index) = default;
43
48 MultiIndex(const int a_order);
49
53 ~MultiIndex() = default;
54
59 void
60 define(const int a_order);
61
66 int
67 factorial() const;
68
74 getCurrentIndex() const;
75
80 int
81 getOrder() const;
82
88 int
90
97 getMultiIndex(const int a_linearIndex) const;
98
105 operator=(const MultiIndex& a_other) = default;
106
112 int
113 operator[](const int a_dir) const;
114
118 void
119 operator++();
120
124 void
125 reset();
126
131 bool
132 ok() const;
133
139 int
140 norm() const;
141
146 int
147 getNumIndices() const;
148
154 Real
155 pow(const RealVect& a_vec);
156
157protected:
162
167
171 std::vector<IntVect>::iterator m_iter;
172
177
182
187
193 int
194 factorial(const int a_n) const;
195
199 void
200 makeIndices();
201
205 void
206 makeMaps();
207
208private:
212 static int
213 norm(const IntVect a_iv);
214};
215
216#include <CD_NamespaceFooter.H>
217
218#endif
SpaceDim multi-index type for use with higher order Taylor series.
Definition CD_MultiIndex.H:31
void operator++()
Increment multi-index lexicographically.
Definition CD_MultiIndex.cpp:152
void makeMaps()
Fill maps.
Definition CD_MultiIndex.cpp:195
MultiIndex()=delete
Weak construction is bad.
int factorial() const
Factorial function.
Definition CD_MultiIndex.cpp:111
int norm() const
Get norm of multi-index.
Definition CD_MultiIndex.cpp:122
std::vector< IntVect >::iterator m_iter
Pointer to index.
Definition CD_MultiIndex.H:171
void makeIndices()
Fill m_indices.
Definition CD_MultiIndex.cpp:158
IntVect m_curIndex
Current idnex.
Definition CD_MultiIndex.H:166
void reset()
Reset multi-index to point to the first index.
Definition CD_MultiIndex.cpp:34
IntVect getCurrentIndex() const
Get index.
Definition CD_MultiIndex.cpp:40
~MultiIndex()=default
Destructor.
int m_order
Maximum order representable by this multi-index.
Definition CD_MultiIndex.H:161
int operator[](const int a_dir) const
Get one index.
Definition CD_MultiIndex.cpp:88
IntVect getMultiIndex(const int a_linearIndex) const
Get a multi-index corresponding to the linear index (i.e. column number)
Definition CD_MultiIndex.cpp:73
MultiIndex(const MultiIndex &a_index)=default
Copy constructor.
MultiIndex & operator=(const MultiIndex &a_other)=default
Assignment operator.
std::map< int, IntVect > m_mapToMultiIndex
Map going from linear index to multi-index.
Definition CD_MultiIndex.H:186
bool ok() const
Check that multi-index is ok.
Definition CD_MultiIndex.cpp:94
std::map< IntVect, int > m_mapToLinearIndex
Map for going from lexicographically ordered multi-index to linear index.
Definition CD_MultiIndex.H:181
int getLinearIndex(const IntVect a_MultiIndex) const
Get the linear index (i.e. column number) for a given multi-index.
Definition CD_MultiIndex.cpp:58
int getOrder() const
Get highest order for this index.
Definition CD_MultiIndex.cpp:46
void define(const int a_order)
Define function.
Definition CD_MultiIndex.cpp:23
std::vector< IntVect > m_indices
All indices.
Definition CD_MultiIndex.H:176
Real pow(const RealVect &a_vec)
Power function for multi-indices.
Definition CD_MultiIndex.cpp:139
int getNumIndices() const
Return number of indices in MultiIndex range.
Definition CD_MultiIndex.cpp:52
Base class for a tracer particle solver. This solver can advance particles in a pre-defined velocity ...
Definition CD_TracerParticleSolver.H:38
TracerParticleSolver()
Default constructor.
Definition CD_TracerParticleSolverImplem.H:26