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;
    };

}