Proto_MBBoundaryData.H

Go to the documentation of this file.

#pragma once
#ifndef _PROTO_MB_BOUNDARY_DATA_
#define _PROTO_MB_BOUNDARY_DATA__

#include "Proto.H"
#include "Proto_MBDisjointBoxLayout.H"

namespace Proto
{
    template<typename T, unsigned int C, MemType MEM>
    struct MBBoundaryData
    {
        MBBoundaryData();

        MBBoundaryData(
                        MBDisjointBoxLayout a_layout,
            MBIndex a_localIndex,
            MBIndex a_adjIndex,
            Box a_localBox,
            Box a_adjBox,
            CoordPermutation a_adjToLocal,
            Point a_ghost = Point::Zeros());
        
        inline void define(
                        MBDisjointBoxLayout a_layout,
            MBIndex a_localIndex,
            MBIndex a_adjIndex,
            Box a_localBox,
            Box a_adjBox,
            CoordPermutation a_adjToLocal,
            Point a_ghost = Point::Zeros());
                inline Face localFace() const;
                inline Face adjFace() const;
        inline void fill(const BoxData<T, C, MEM>& a_data);

        inline Box localBox() const  {return localData->box().grow(-localGhost);}
        inline Box adjBox() const {return adjData->box().grow(-adjGhost);}

        inline void print();

                MBDisjointBoxLayout layout;
        MBIndex localIndex;
        MBIndex adjIndex;
                Point localDir;
                Point adjDir;
        CoordPermutation adjToLocal;
        Point localGhost;
        Point adjGhost;
        std::shared_ptr<BoxData<T, C, MEM>> localData;
        std::shared_ptr<BoxData<T, C, MEM>> adjData;
    };
    
    template<typename T, unsigned int C, MemType MEM>
    MBBoundaryData<T, C, MEM>::MBBoundaryData() {}

    template<typename T, unsigned int C, MemType MEM>
    MBBoundaryData<T, C, MEM>::MBBoundaryData(
                        MBDisjointBoxLayout a_layout,
            MBIndex a_localIndex,
            MBIndex a_adjIndex,
            Box a_localBox,
            Box a_adjBox,
            CoordPermutation a_adjToLocal,
            Point a_ghost)
    {
        define(a_layout, a_localIndex, a_adjIndex, a_localBox, a_adjBox, a_adjToLocal, a_ghost);
    }
    
    template<typename T, unsigned int C, MemType MEM>
    void MBBoundaryData<T, C, MEM>::define(
            MBDisjointBoxLayout a_layout,
                        MBIndex a_localIndex,
            MBIndex a_adjIndex,
            Box a_localBox,
            Box a_adjBox,
            CoordPermutation a_adjToLocal,
            Point a_ghost)
    {
                layout = a_layout;
        localGhost = a_ghost;
        localIndex = a_localIndex;
        adjIndex = a_adjIndex;
        Box localBox = a_localBox.grow(a_ghost);
        localData = std::make_shared<BoxData<T, C, MEM>>(localBox, 7);
        adjGhost = a_adjToLocal.inverse()(a_ghost).abs();
        Box adjBox = a_adjBox.grow(adjGhost);
        adjData = std::make_shared<BoxData<T, C, MEM>>(adjBox, 7);
        adjToLocal = a_adjToLocal;
                localDir = layout.connectivity(localIndex, adjIndex);
                adjDir = layout.connectivity(adjIndex, localIndex);
    }
    
    template<typename T, unsigned int C, MemType MEM>
    void MBBoundaryData<T, C, MEM>::fill(
            const BoxData<T, C, MEM>& a_data)
    {
        a_data.copyTo(*adjData);
        adjData->copyTo(*localData, adjToLocal);
    }
    
        template<typename T, unsigned int C, MemType MEM>
    Face MBBoundaryData<T, C, MEM>::localFace() const
        {
                PROTO_ASSERT(localDir.codim() == 1,
                        "MBBoundaryData::localFace | Error: This boundary is not a codimension 1 face");
                for (int ii = 0; ii < DIM; ii++)
                {
                        if (localDir[ii] != 0)
                        {
                                return Face(ii, localDir[ii] < 0 ? Side::Lo : Side::Hi);
                        }
                }
        return Face(0,Side::Lo); // dummy return to suppress a compilation warning
        }
         
        template<typename T, unsigned int C, MemType MEM>
    Face MBBoundaryData<T, C, MEM>::adjFace() const
        {
                PROTO_ASSERT(adjDir.codim() == 1,
                        "MBBoundaryData::localFace | Error: This boundary is not a codimension 1 face");
                for (int ii = 0; ii < DIM; ii++)
                {
                        if (adjDir[ii] != 0)
                        {
                                return Face(ii, adjDir[ii] < 0 ? Side::Lo : Side::Hi);
                        }
                }
        return Face(0,Side::Lo); // dummy return to suppress a compilation warning
        }
        
    template<typename T, unsigned int C, MemType MEM>
    void MBBoundaryData<T, C, MEM>::print()
    {
        pr_out() << " | localBoundary: " << localData->box();
        pr_out() << " | adjBoundary: " << adjData->box() << std::endl;
    }
} // end namespace Proto
#endif //end include guard