gsLobattoRule.h

Go to the documentation of this file.

#pragma once

namespace gismo
{

template<class T>
class gsLobattoRule GISMO_FINAL : public gsQuadRule<T>
{
public:

    typedef memory::unique_ptr<gsLobattoRule> uPtr;

    gsLobattoRule() { }

    gsLobattoRule(gsVector<index_t> const & numNodes,
                const unsigned digits = 0 )
    {
        gsLobattoRule::setNodes(numNodes, digits);
    }

    static uPtr make(gsVector<index_t> const & numNodes,
                        const unsigned digits = 0 )
    { return uPtr( new gsLobattoRule(numNodes,digits) ); }

    gsLobattoRule(index_t numNodes, const unsigned digits = 0 )
    {
        this->setNodes(numNodes, digits);
    }

    ~gsLobattoRule() { }

public:
    // see gsQuadRule.h for documentation
    void setNodes( gsVector<index_t> const & numNodes,
                   unsigned digits = 0 );

    using gsQuadRule<T>::setNodes;// unhide base

private:

    static void computeReference(index_t n, gsVector<T> & x, gsVector<T> & w,
                                 unsigned digits = 0 );

    static bool lookupReference (index_t n, gsVector<T> & x, gsVector<T> & w);

}; // class gsLobattoRule


} // namespace gismo


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