GhostPointGenerator.cpp

#include "GhostPointGenerator.hpp"

namespace
{
        size_t IsBoundaryEdge(Edge const& edge, int npoints)
        {
                if (edge.neighbors.first >= npoints)
                        return 1;
                if (edge.neighbors.second >= npoints)
                        return 2;
                else
                        return 0;
        }
}

GhostPointGenerator::~GhostPointGenerator(void) {}

vector<std::pair<size_t, size_t> > GhostPointGenerator::GetOuterEdgesIndeces(Tessellation const& tess)const
{
        vector<Edge> const& edges = tess.getAllEdges();
        vector<std::pair<size_t, size_t> > res;
        int npoints = tess.GetPointNo();
        for (size_t i = 0; i < edges.size(); ++i)
        {
                const size_t ghostindex = IsBoundaryEdge(edges[i], npoints);
                if (ghostindex == 1)
                {
                        if(tess.GetOriginalIndex(edges[i].neighbors.first)<npoints)
                                res.push_back(std::pair<size_t, size_t>(i, 1));
                }
                if (ghostindex == 2)
                {
                        if (tess.GetOriginalIndex(edges[i].neighbors.second)<npoints)
                                res.push_back(std::pair<size_t, size_t>(i, 2));
                }
        }
        return res;
}