gsDomainIterator.h

Go to the documentation of this file.

#pragma once

//#include <gsCore/gsBasis.h> // todo: remove
#include <gsCore/gsDomain.h>
#include <gsCore/gsBoundary.h>

namespace gismo
{

template <class T>
class gsDomainIterator
{
public:
    typedef memory::shared_ptr< gsDomainIterator > Ptr;
    typedef memory::unique_ptr< gsDomainIterator > uPtr;

public:

    gsDomainIterator( ) : m_basis(NULL), m_isGood( true ) { }

    gsDomainIterator( const gsBasis<T>& basisParam, const boxSide & s = boundary::none)
        : center( gsVector<T>::Zero(basisParam.dim()) ), m_basis( &basisParam ),
          m_isGood( true ), m_side(s)
    { }

    virtual ~gsDomainIterator() { }

public:

    virtual bool next() = 0;

    virtual bool next(index_t increment) = 0;

    virtual void reset()
    {
        GISMO_NO_IMPLEMENTATION
    }

public:

    size_t id() const   { return m_id; }

    bool good() const   { return m_isGood; }

    short_t dim() const   { return center.size(); }

    virtual void adjacent( const gsVector<bool> & ,
                           gsDomainIterator & )
    {
        GISMO_NO_IMPLEMENTATION
    }

    const gsVector<T>& centerPoint () const
    { return center; }

    virtual const gsVector<T>& lowerCorner() const
    {
        GISMO_NO_IMPLEMENTATION
    }

    virtual const gsVector<T>& upperCorner() const
    {
        GISMO_NO_IMPLEMENTATION
    }

    virtual const T getPerpendicularCellSize() const
    {
        GISMO_NO_IMPLEMENTATION
    }

    virtual bool isBoundaryElement() const
    {
        GISMO_NO_IMPLEMENTATION
    }

    T getCellSize() const
    {
        return (upperCorner() - lowerCorner()).norm();
    }

    T getMinCellLength() const
    {
        return (upperCorner() - lowerCorner()).minCoeff();
    }

    T getMaxCellLength() const
    {
        return (upperCorner() - lowerCorner()).maxCoeff();
    }

    T volume() const
    { return (upperCorner() - lowerCorner()).prod(); }

    virtual size_t numElements() const
    {
        //\todo Remove this implementation. Probably using a shallow
        //copy, "reset" and "next" would do this better.

        // Buggy, and probably a terrible implementation,
        // but needed and therefore can be useful
        // sometimes.
        typename gsBasis<T>::domainIter domIter = m_basis->makeDomainIterator(m_side);

        size_t numEl = 0;
        for (; domIter->good(); domIter->next(), numEl++){}

        return numEl;
    }


    inline boxSide side() const {return m_side;}

public:

    gsVector<T> center;

protected:
    const gsBasis<T> * m_basis;

    bool m_isGood;

    boxSide m_side;

    size_t m_id;

private:
    // disable copying
    gsDomainIterator( const gsDomainIterator& );
    gsDomainIterator& operator= ( const gsDomainIterator& );
}; // class gsDomainIterator


//template<class T>
//std::ostream &operator<<(std::ostream &os, const gsDomainIterator<T>& b)
//{return b.print(os); }


} // namespace gismo