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