MArray.hpp

Go to the documentation of this file.

// This file is part of GranOO, a workbench for DEM simulation.
//
// Author(s)    : - Damien Andre       IRCER/UNILIM, Limoges France
//                  <damien.andre@unilim.fr>
//                - Jean-luc Charles   Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <jean-luc.charles@ensam.eu>
//                - Jeremie Girardot   Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <jeremie.girardot@ensam.eu>
//                - Cedric Hubert      LAMIH/UPHF, Valenciennes France
//                  <cedric.hubert@uphf.fr>
//                - Ivan Iordanoff     Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <ivan.iordanoff@ensam.eu>
//
// Copyright (C) 2008-2019 D. Andre, JL. Charles, J. Girardot, C. Hubert, I. Iordanoff
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program.  If not, see <http://www.gnu.org/licenses/>.
 
 
#ifndef _libCore_MArray_hpp_
#define _libCore_MArray_hpp_
 
#include <boost/archive/text_oarchive.hpp>
#include <boost/archive/text_iarchive.hpp>
 
#include "GranOO3/Core/SetOf.hpp"
#include "GranOO3/3rdParty/Eigen/Dense"
 
namespace GranOO3
{
  namespace Core
  {
 
    // MPtr is a small utility class able to switch between a legacy object and a pointer
    template < typename T>
    class MPtr
    {
    public:
      template<typename... Args> explicit MPtr(Args... args);
      explicit MPtr(const T&);
 
      
      const T* operator->() const;
      T* operator->();
      
      T& operator* ();
      const T& operator* () const;
      
      void set_ptr(T* pt);
      T*&  ptr();
      const T*&  get_ptr() const;
      
      friend class boost::serialization::access;
      template<class Archive> void serialize(Archive&, const unsigned int);
 
      T&    operator=(const T& v   ) {*_p= v; return *_p;};
      MPtr& operator=(const MPtr& v) {*_p=*v; return *this;};
      template<typename T1> void operator+= (T1 v) {*_p+=v;}
      template<typename T1> void operator-= (T1 v) {*_p-=v;}
      template<typename T1> void operator*= (T1 v) {*_p*=v;}
      template<typename T1> void operator/= (T1 v) {*_p/=v;}
 
    private:
      MPtr(const MPtr&) = delete;
      
      
    private:
      T  _t;
      T* _p;
    };
 
 
 
    // MArray class is able to initialize an Eigen::Vector from a collection of C objects
    // It changes on-the-fly the MPtr objects for pointing to the good position in the Eigen::Vector
    // MArray should not be used as this, this an intermediate class
    class MArray
    {
    public:
      static constexpr int N();
 
      template <typename C, typename T, typename... Rest>
      static constexpr int N(Core::MPtr<T>& (C::*f)(),  Rest... rest);
 
 
    public:
      MArray();
      
      
      template <typename C, typename T, typename... Rest>
      void init(const std::vector<C*>& set, Core::MPtr<T>& (C::*f)(),  Rest... rest);
      
      template <typename C>
      void fill(size_t index, C* item);
      
      template <typename C, typename T, typename... Rest>
      void fill(size_t index, C* item, Core::MPtr<T>& (C::*f)(),  Rest... rest);
      
      Eigen::VectorXd& array(); 
      const Eigen::VectorXd& array() const;
      
    private:
      Eigen::VectorXd _array;
      bool            _init; 
    };
 
 
    
 
  }
}
 
 
#include "GranOO3/Core/MArray.tpp"
 
 
 
#endif