AmrBoundary.h

Go to the documentation of this file.

#ifndef AMR_BOUNDARY_H
#define AMR_BOUNDARY_H

#include <vector>

#include "Amr/AmrDefs.h"

template <class Level>
class AmrBoundary {

public:
    typedef typename Level::umap_t      umap_t;
    typedef typename Level::lo_t        lo_t;
    typedef typename Level::go_t        go_t;
    typedef typename Level::scalar_t    scalar_t;
    typedef typename Level::basefab_t   basefab_t;
    typedef amr::AmrIntVect_t           AmrIntVect_t;
    
public:
    
    AmrBoundary(go_t nPoints) : nPoints_m(nPoints) { };
    
    
    bool isBoundary(const AmrIntVect_t& iv, const go_t* nr) const {
        return AMREX_D_TERM(   isBoundary(iv, 0, nr),
                            || isBoundary(iv, 1, nr),
                            || isBoundary(iv, 2, nr));
    }

    bool isBoundary(const AmrIntVect_t& iv,
                    const lo_t& dir,
                    const go_t* nr) const {
        return ( iv[dir] < 0 || iv[dir] >= nr[0] );
    }
    
    virtual void apply(const AmrIntVect_t& iv,
                       const lo_t& dir,
                       umap_t& map,
                       const scalar_t& value,
                       Level* mglevel,
                       const go_t* nr) = 0;
    
    const go_t& getNumberOfPoints() const {
        return nPoints_m;
    }
    
private:
    const go_t nPoints_m;    
};

#endif