index.h

Go to the documentation of this file.

#ifndef INDEX_H
#define INDEX_H

#include <vector>

namespace Parallel {
    class Index
    {
    private:
        static std::vector<unsigned int> m_N;
        static std::vector<long long> m_NTot;
    public:
        Index();
        ~Index();

        // Index getter
        static inline unsigned long cubeIndex(unsigned long int i, unsigned long int j, unsigned long int k, unsigned long int Ni, unsigned long int Nj) {
            return i + Ni*(j + Nj*k);
        }

        static inline unsigned long getIndex(unsigned int i, unsigned int j, unsigned int k, unsigned int l) {
            return i + m_N[0]*(j + m_N[1]*(k + m_N[2]*l));
        }

        static inline long long getGlobalIndex(long long i, long long j, long long k, long long l) {
            return i + m_NTot[0]*(j + m_NTot[1]*(k + m_NTot[2]*l)); // column-major
        }

        // Setters
        static void setN(std::vector<unsigned int> N);
        static void setNTot(unsigned int NSpatial, unsigned int NTemporal);
    };
}
#endif // INDEX_H