gsGaussRule.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsAssembler/gsQuadRule.h>
 
namespace gismo
{
 
template<class T>
class gsGaussRule GISMO_FINAL : public gsQuadRule<T>
{
public:
 
    typedef memory::unique_ptr<gsGaussRule> uPtr;
 
    gsGaussRule() { }
 
    gsGaussRule(gsVector<index_t> const & numNodes,
                const unsigned digits = 0 )
    {
        gsGaussRule::setNodes(numNodes, digits);
    }
 
    static uPtr make(gsVector<index_t> const & numNodes,
                        const unsigned digits = 0 )
    { return uPtr( new gsGaussRule(numNodes,digits) ); }
 
    gsGaussRule(index_t numNodes, const unsigned digits = 0 )
    {
        this->setNodes(numNodes, digits);
    }
 
    gsGaussRule(const gsBasis<T> & basis, const T quA, const index_t quB, short_t fixDir = -1);
    //const unsigned digits = std::numeric_limits<T>::digits10 );
 
    gsGaussRule(const gsBasis<T> & basis, const gsOptionList & options, short_t fixDir = -1);
    //const unsigned digits = std::numeric_limits<T>::digits10 );
 
    ~gsGaussRule() { }
    
public:
    // see gsQuadRule.h for documentation
    void setNodes( gsVector<index_t> const & numNodes, 
                   unsigned digits = 0 );
 
    using gsQuadRule<T>::setNodes; // unhide base
 
private:
 
    void init(const gsBasis<T> & basis, const T quA, const index_t quB, short_t fixDir);
    
    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 gsGaussRule
 
 
} // namespace gismo
 
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsGaussRule.hpp)
#endif