domain.hh

Go to the documentation of this file.

#ifndef DUNE_SPGRID_DOMAIN_HH
#define DUNE_SPGRID_DOMAIN_HH
 
#include <dune/grid/spgrid/cube.hh>
#include <dune/grid/spgrid/topology.hh>
#include <dune/grid/spgrid/refinement.hh>
 
namespace Dune
{
 
  // SPDomain
  // --------
 
  template< class ct, int dim >
  class SPDomain
  {
    typedef SPDomain< ct, dim > This;
 
  public:
    typedef SPCube< ct, dim > Cube;
    typedef SPTopology< dim > Topology;
 
    typedef typename Cube::ctype ctype;
 
    static const int dimension = Cube::dimension;
 
    typedef typename Cube::GlobalVector GlobalVector;
 
    SPDomain ( const GlobalVector &a, const GlobalVector &b );
 
    SPDomain ( const std::vector< Cube > &cubes, const Topology &topology );
 
    const Cube &cube () const { return cube_; }
 
    const Topology &topology () const { return topology_; }
 
    bool contains ( const GlobalVector &x ) const { return cube().contains( x ); }
 
    static This unitCube ();
 
  private:
    Cube cube_;
    Topology topology_;
  };
 
 
 
  // Implementation of SPDomain
  // --------------------------
 
  template< class ct, int dim >
  inline SPDomain< ct, dim >
    ::SPDomain ( const GlobalVector &a, const GlobalVector &b )
  : cube_( a, b ),
    topology_()
  {}
 
 
  template< class ct, int dim >
  inline SPDomain< ct, dim >
    ::SPDomain ( const std::vector< Cube > &cubes, const Topology &topology )
  : cube_( cubes[ 0 ] ),
    topology_( topology )
  {}
 
 
  template< class ct, int dim >
  inline typename SPDomain< ct, dim >::This
  SPDomain< ct, dim >::unitCube ()
  {
    GlobalVector a, b;
    for( int i = 0; i < dimension; ++i )
    {
      a = ctype( 0 );
      b = ctype( 1 );
    }
    return This( a, b );
  }
 
} // namespace Dune
 
#endif // #ifndef DUNE_SPGRID_DOMAIN_HH