domain.h

#ifndef DOMAIN_H
#define DOMAIN_H
 
#include <vector>
 
#include "../svector/svector.h"
#include "fields.h"
 
namespace domain {
 
enum Face
{
    iMinus,
    jMinus,
    kMinus,
    iPlus,
    jPlus,
    kPlus
};
 
constexpr Face getOppositeFace(const Face& f)
{
    return static_cast<Face>((f + 3) % 6);
};
 
struct Domain {
    Domain(const int& ni, const int& nj, const int& nk, const int& nn);
 
    const int n[3];
 
    const int& ni = n[0];
 
    const int& nj = n[1];
 
    const int& nk = n[2];
 
    const int nn;
 
    std::vector<fields::Owner<svec::SVector>> s;
 
    std::vector<fields::Owner<svec::NormalizedSVector>> c;
 
    fields::Helper<svec::SVector> getGhost(const Face& f, const int& n);
 
    fields::Helper<svec::SVector> getEdge(const Face& f, const int& n);
 
    constexpr bool hasNeighbor(Face f) const;
 
    template <class T>
    T getMax(const T& in) const;
};
 
constexpr bool domain::Domain::hasNeighbor(Face f) const
{
    return false;
}
 
template <class T>
inline T Domain::getMax(const T& in) const
{
    return in;
}
 
} // namespace domain
 
#endif