gsPointLoads.h

Go to the documentation of this file.

#pragma once

namespace gismo
{



template<class T>
struct point_load
{
    point_load(const gsVector<T> & _point, 
               const T             _value, 
               int _patch = 0, 
               bool _parametric = true)
    :
    patch(_patch), value(_value), point(1), parametric(_parametric)
    { 
        point[0] = _value;
    }

    point_load(const gsVector<T> & _point, 
               const gsVector<T> & _value, 
               int _patch = 0, 
               bool _parametric = true)
    :
    patch(_patch), value(_value), point(_point), parametric(_parametric)
    { }

    int patch;

    gsVector<T> value;

    gsVector<T> point;

    bool parametric;
};



template<class T>
class gsPointLoads
{
public:
    typedef point_load<T> pLoad;
    typedef typename std::vector<pLoad> plContainer;

    typedef typename std::vector<pLoad>::iterator iterator;

    typedef typename std::vector<pLoad>::const_iterator const_iterator;

public:

    std::ostream & print(std::ostream &os) const
    { 
        os << "gsPointLoads: "<<m_pointLoads.size()<<"\n";
        return os; 
    }

public:

    gsPointLoads() 
    { }

    ~gsPointLoads() // Destructor
    { }
    
public:

    void clear()
    {
        m_pointLoads.clear();
    }

    inline pLoad   operator [] (size_t i) const { return m_pointLoads[i]; }
    inline pLoad & operator [] (size_t i) { return m_pointLoads[i]; }

    void addLoad(const gsVector<T> & _point, 
                 const gsVector<T> & _value, 
                 int _patch = 0, 
                 bool _parametric = true)
    {
        m_pointLoads.push_back( pLoad(_point,_value,_patch,_parametric) );
    }

    void addLoad(const gsVector<T> & _point, 
                 const T             _value, 
                 int _patch = 0, 
                 bool _parametric = true)
    {
        m_pointLoads.push_back( pLoad(_point,_value,_patch,_parametric) );
    }

    size_t numLoads() const { return  m_pointLoads.size(); }

private:

    plContainer  m_pointLoads; 

}; // class gsPointLoads

template<class T>
std::ostream &operator<<(std::ostream &os, const gsPointLoads<T>& pls)
{return pls.print(os); }

#ifdef GISMO_WITH_PYBIND11

  void pybind11_init_gsPointLoads(pybind11::module &m);

#endif // GISMO_WITH_PYBIND11

} // namespace gismo