periodic.h

Go to the documentation of this file.

 
#ifndef PERIODIC_H
#define PERIODIC_H
 
#include <array>
#include <cstdint>
#include <map>
#include <string>
#include <vector>
 
#include "mesh.h"
 
struct PeriodicNodes {
    std::vector<std::string> zones;  // the zones joined by the pairs
    std::vector<std::array<rtype, N_DIM>> translations;
    std::vector<uint32_t> key;
    std::vector<std::array<int8_t, 3>> lattice;
};
 
std::vector<size_t> periodic_directions(const std::vector<Mesh::PeriodicPair> & pairs, PeriodicNodes & result);
 
rtype periodic_tolerance(rtype h);
 
class PeriodicGrid {
    public:
        using Point = std::array<rtype, N_DIM>;
        using Cell = std::array<int64_t, N_DIM>;
        static constexpr int N_NEIGHBORS = N_DIM == 2 ? 9 : 27;
 
        explicit PeriodicGrid(rtype spacing) : h(spacing), tol(periodic_tolerance(spacing)) {}
 
        Cell cell_of(const Point & x) const;
        static Cell neighbor(const Cell & c, int m);
        void insert(uint32_t id, const Point & x) { grid[cell_of(x)].push_back({id, x}); }
        uint32_t find(const Point & y, uint32_t & match) const;
 
    private:
        rtype h, tol;
        std::map<Cell, std::vector<std::pair<uint32_t, Point>>> grid;
};
 
std::string periodic_point_string(const std::array<rtype, N_DIM> & x);
 
PeriodicNodes match_periodic_nodes(const std::vector<std::array<rtype, N_DIM>> & nodes,
                                   const std::vector<Mesh::BoundaryFace> & boundary_faces,
                                   const std::vector<Mesh::PeriodicPair> & pairs);
 
#endif // PERIODIC_H