gsBasisRefs.h

Go to the documentation of this file.

#pragma once

namespace gismo
{

template <class T>
class gsBasisRefs
{
private:
    typedef typename std::vector<gsMultiBasis<T> >::const_iterator mbIterator;

public:

    inline gsBasisRefs(std::vector<gsMultiBasis<T> > & bases, const size_t k)
    {
        GISMO_ASSERT(bases.size()>0, "Cannot construct empty list of gsBasis.");
        m_refs.reserve(bases.size());
        for(mbIterator  it = bases.begin(); it != bases.end(); ++it)
            m_refs.push_back( &(*it)[k] );
    }

    inline const gsBasis<T> & operator[](size_t i) const
    { return *m_refs[i]; }

    inline operator const gsBasis<T> &() const 
    { return *m_refs.front(); }

    inline const gsBasis<T> & front() const 
    { return *m_refs.front(); }

    inline const gsBasis<T> & back () const 
    { return *m_refs.back(); }

    inline size_t size () const
    { return m_refs.size(); }

    inline short_t dim() const { return m_refs.front()->dim();}

private:
    std::vector<const gsBasis<T>*> m_refs;

private:
    gsBasisRefs();
    gsBasisRefs(const gsBasisRefs &);
};

} // namespace gismo