bndprojection.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_ALU_BNDPROJECTION_HH
#define DUNE_ALU_BNDPROJECTION_HH
 
namespace Dune {
 
  template <class GridImp, class ctype = double >
  class ALUGridBoundaryProjection
#ifdef ALUGRID_VERTEX_PROJECTION
    : public GridImp :: ALUGridVertexProjectionType
#endif
  {
    typedef GridImp GridType;
    // type of double coordinate vector
    typedef ctype coord_t[ GridType :: dimensionworld ];
  protected:
 
    const GridType& grid_;
  public:
    typedef typename GridType :: DuneBoundaryProjectionType DuneBoundaryProjectionType;
 
    typedef typename DuneBoundaryProjectionType :: CoordinateType CoordinateType;
 
    ALUGridBoundaryProjection(const GridType& grid)
      : grid_( grid )
    {}
 
    int operator () (const coord_t &orig,
                     coord_t &prj) const
    {
      return this->operator()( orig, 0, prj);
    }
 
    int operator () (const coord_t &orig,
                     const int segmentIndex,
                     coord_t &prj) const
    {
#ifdef ALUGRID_VERTEX_PROJECTION
      // get boundary projection
      const DuneBoundaryProjectionType* bndPrj =
        grid_.boundaryProjection( segmentIndex );
 
      // if pointer is zero we do nothing, i.e. identity mapping
      if( bndPrj )
      {
        // call projection operator
        reinterpret_cast<CoordinateType &> (* (&prj[0])) =
          (*bndPrj)( reinterpret_cast<const CoordinateType &> (* (&orig[0])) );
      }
#endif
 
      // return 1 for success
      return 1;
    }
  };
 
} // end namespace Dune
#endif