gsDomain.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsCore/gsLinearAlgebra.h>
 
namespace gismo
{
 
template<class T>
class gsDomain
{
public:
 
    virtual ~gsDomain() { }
 
#if EIGEN_HAS_RVALUE_REFERENCES && EIGEN_GNUC_AT_MOST(4,7) && !EIGEN_COMP_PGI
    // defaulted declaration required at least in Gcc 4.7.2
    gsDomain() = default;
    gsDomain(const gsDomain&) = default;
    gsDomain(gsDomain&&) = default;
    gsDomain & operator=(const gsDomain&) = default;
    gsDomain & operator=(gsDomain&&) = default;
#endif
 
public:
 
    virtual short_t dim() const
    { gsWarn << "gsDomain: dimension() not defined at "<< *this << "\n"; return 0; }
 
    virtual gsMatrix<T> boundingBox()
    {GISMO_NO_IMPLEMENTATION}
 
    virtual gsMatrix<T> elements()
    {GISMO_NO_IMPLEMENTATION}
 
    virtual gsMatrix<T> mesh()
    {GISMO_NO_IMPLEMENTATION }
 
    virtual T minMeshSize()
    {GISMO_NO_IMPLEMENTATION}
 
    virtual std::vector<T> breaks() const
    {GISMO_NO_IMPLEMENTATION}
    
    virtual gsDomain* clone() const = 0;
 
    virtual std::ostream &print(std::ostream &os) const = 0;
 
    /*
      Member functions that may be implemented or not in the derived class
    */
 
    virtual std::vector<T> unique() const 
    { gsWarn<<"gsDomain: unique() was not defined at"<< *this <<"\n"; return std::vector<T>() ;}
 
    void merge(gsDomain<T> * other ) 
    { gsWarn<<"gsDomain: merge(..) was not defined in "<< *this <<"\n"; return;}
 
}; // class gsDomain
 
template<class T>
std::ostream &operator<<(std::ostream &os, const gsDomain<T>& b)
{ return b.print(os); }
 
 
} // namespace gismo