gsGridHierarchy.h

Go to the documentation of this file.

#pragma once

#include <gsCore/gsDofMapper.h>
#include <gsNurbs/gsKnotVector.h>
#include <gsIO/gsOptionList.h>

namespace gismo
{

template< typename T >
class gsGridHierarchy
{

public:

    static gsGridHierarchy buildByRefinement(
        gsMultiBasis<T> mBasis,
        const gsBoundaryConditions<T>& boundaryConditions,
        const gsOptionList& assemblerOptions,
        index_t levels,
        index_t numberOfKnotsToBeInserted = 1,
        index_t multiplicityOfKnotsToBeInserted = 1,
        index_t unk = 0
        );

    static gsGridHierarchy buildByRefinement(
        gsMultiBasis<T> mBasis,
        const gsBoundaryConditions<T>& boundaryConditions,
        const gsOptionList& options
        )
    {
        return gsGridHierarchy::buildByRefinement(
            give(mBasis),
            boundaryConditions,
            options,
            options.askInt( "Levels", 3 ),
            options.askInt( "NumberOfKnotsToBeInserted", 1 ),
            options.askInt( "MultiplicityOfKnotsToBeInserted", 1 ),
            0
        );
    }


    static gsGridHierarchy buildByCoarsening(
        gsMultiBasis<T> mBasis,
        const gsBoundaryConditions<T>& boundaryConditions,
        const gsOptionList& assemblerOptions,
        index_t levels,
        index_t degreesOfFreedom = 0,
        index_t unk = 0
        );

    static gsGridHierarchy buildByCoarsening(
        gsMultiBasis<T> mBasis,
        const gsBoundaryConditions<T>& boundaryConditions,
        const gsOptionList& options
    )
    {
        return gsGridHierarchy::buildByCoarsening(
            give(mBasis),
            boundaryConditions,
            options,
            options.askInt( "Levels", 3 ),
            options.askInt( "DegreesOfFreedom", 0 ),
            0
        );
    }

    static gsOptionList defaultOptions()
    {
        gsOptionList opt;
        opt.addInt( "DirichletStrategy", "Method for enforcement of Dirichlet BCs [11..14]", 11 );
        opt.addInt( "InterfaceStrategy", "Method of treatment of patch interfaces [0..3]", 1  );
        opt.addInt( "Levels", "Number of levels to be constructed in the grid hierarchy", 3 );
        opt.addInt( "DegreesOfFreedom",   "Number of dofs in the coarsest grid in the grid hierarchy (only buildByCoarsening)", 0 );
        opt.addInt( "NumberOfKnotsToBeInserted", "The number of knots to be inserted (only buildByRefinement)", 1 );
        opt.addInt( "MultiplicityOfKnotsToBeInserted",   "The multiplicity of the knots to be inserted (only buildByRefinement)", 1 );
        return opt;
    }

    void clear()
    {
        m_mBases.clear();
        m_transferMatrices.clear();
    }

    const std::vector< gsMultiBasis<T> >& getMultiBases() const
    { return m_mBases; }

    gsGridHierarchy& moveMultiBasesTo( std::vector< gsMultiBasis<T> >& o )
    { o = give(m_mBases); return *this; }

    const std::vector< gsSparseMatrix<T, RowMajor> >& getTransferMatrices() const
    { return m_transferMatrices; }

    gsGridHierarchy& moveTransferMatricesTo( std::vector< gsSparseMatrix<T, RowMajor> >& o )
    { o = give(m_transferMatrices); return *this; }

private:
    std::vector< gsMultiBasis<T> > m_mBases;
    std::vector< gsSparseMatrix<T, RowMajor> > m_transferMatrices;
};

} // namespace gismo

#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsGridHierarchy.hpp)
#endif