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