gsContainerBasisBase.h

Go to the documentation of this file.

 
/*
    TO DO
*/
 
#pragma once
 
#include<gsIO/gsOptionList.h>
#include<gsUnstructuredSplines/src/gsContainerBasis.h>
 
namespace gismo
{
 
    template<short_t d,class T>
    class gsContainerBasisBase
    {
 
    private:
        typedef gsContainerBasis<d,T> basis;
        typedef typename std::vector<basis> basisContainer;
 
        typedef memory::shared_ptr< gsContainerBasisBase > Ptr;
 
        typedef memory::unique_ptr< gsContainerBasisBase > uPtr;
 
    public:
 
        virtual ~gsContainerBasisBase() {};
 
        gsContainerBasisBase(gsMultiPatch<T> & mp, gsMultiBasis<T> & mb)
                : m_patches(mp), m_multiBasis(mb)
        {
 
        }
 
    public:
 
        virtual gsOptionList & options() { return m_options; }
        virtual void defaultOptions() { };
 
        virtual void setOptions(gsOptionList opt) {m_options.update(opt, gsOptionList::addIfUnknown); };
 
        virtual void init() = 0;
        virtual void compute() = 0;
        
 
        void getMultiBasis(gsMultiBasis<T> & multiBasis_result)
        {
            multiBasis_result.clear();
 
            std::vector<gsBasis<T> *> basis_temp = std::vector<gsBasis<T> *>(m_patches.nPatches());
            for (size_t np = 0; np < m_patches.nPatches(); np++) {
                gsContainerBasis<2, real_t>::uPtr basis = gsContainerBasis<d, T>::make(m_bases[np]);
                basis_temp[np] = static_cast<gsBasis<T> *>(basis.release());
            }
 
            multiBasis_result = gsMultiBasis<T>(basis_temp, m_patches.topology());
        };
 
        gsSparseMatrix<T> & getSystem() { return m_matrix; };
        void setSystem(gsSparseMatrix<T> & system) { m_matrix = system; };
 
    private:
        // Helper functions
 
    protected:
        // Data members
 
        // Put here the members of the shared functions
        gsMultiPatch<T> & m_patches;
        gsMultiBasis<T> & m_multiBasis;
 
        gsOptionList m_options;
 
        basisContainer m_bases;
 
        gsSparseMatrix<T> m_matrix;
    };
 
}