ElementT.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 _LibDEM_ElementT_TPP
#define _LibDEM_ElementT_TPP
 
 
 
#include <boost/serialization/version.hpp>
#include <boost/serialization/map.hpp>
 
 
#include "GranOO3/DEM/Element.hpp"
 
#include "GranOO3/Shape/Sphere.hpp"
#include "GranOO3/Shape/Box.hpp"
#include "GranOO3/Shape/Cylinder.hpp"
#include "GranOO3/Shape/Cone.hpp"
#include "GranOO3/Shape/ConeTruncated.hpp"
#include "GranOO3/Shape/Polyhedron.hpp"
 
 
namespace GranOO3
{
  namespace DEM
  {
 
    template<class shape>
    class ElementT : public Element,
             public shape,
             public Core::PropClass< ElementT<shape> >,
             public Core::Register< ElementT<shape> >
    {
 
      GRANOO_CLASS_T(DEM::ElementT, shape, Element);
 
    public:
      // Constructors and destructor
      ElementT(const Geom::Point& pos, double radius, double density); //Just for spherical DE
      ElementT(const Geom::Point& pos, double density);
      ElementT(const Geom::Point& pos);
      ElementT();
 
      virtual ~ElementT();
 
      GRANOO_ACCESS_REF_CONST(local_frame, Geom::Frame,      Node::local_frame());
      GRANOO_ACCESS_GET      (center     , Geom::Point     , Node::get_center() );
      GRANOO_ACCESS          (quaternion , Geom::Quaternion, Node::quaternion() );
 
      virtual void draw();
      virtual void draw_edge();
      virtual std::string info() const;
      std::string shape_class_ID() const;
 
 
      // For avoiding confusion with base classes
      double get_volume() const;
      const Geom::Vector& get_position() const;
 
      // cloning
      virtual Core::Base& clone();
      void make_equal_to(const ElementT& el);
 
      // Shape casting
      virtual bool is_spherical_element() const;
      virtual bool is_polyhedral_element() const;
 
      // IO
      virtual std::ostream& write_ascii (std::ostream& os) const;
      virtual std::istream& read_ascii (std::istream& is);
 
 
    private:
      ElementT(const ElementT&) = delete;
      ElementT & operator=(const ElementT &) = delete;
 
      //BOOST SERIALIZATION
      friend class boost::serialization::access;
      template<class Archive> void serialize(Archive&, const unsigned int);
 
    };
 
 
    template<class shape>
    ElementT<shape>::ElementT(const Geom::Point& pos, double density)
      : Element(*this, pos),
    shape(Element::local_frame()) {
      //neighbour_(*this)
      Shape::Volume::update_bounding_shape();
      Body::set_density(density);
    }
 
    template<class shape>
    ElementT<shape>::ElementT(const Geom::Point& pos)
      : Element(*this, pos),
      shape(Element::local_frame()) {
      //neighbour_(*this)
    }
 
    template<class shape>
    ElementT<shape>::ElementT()
      : ElementT<shape>(Geom::global::center) {
    }
 
    template<class shape>
    ElementT<shape>::~ElementT() {
    }
 
    template<class shape> std::string
    ElementT<shape>::shape_class_ID() const {
      return shape::class_ID();
    }
 
    template<class shape>
    double ElementT<shape>::get_volume() const {
      return shape::get_volume();
    }
 
    template<class shape> const Geom::Vector&
    ElementT<shape>::get_position() const {
      return Physic::Node::get_position();
    }
 
 
    template<class shape>
    Core::Base& ElementT<shape>::clone() {
      ElementT<shape>* el = new ElementT<shape>();
      el->make_equal_to(*this);
      return *el;
    }
    
    template<class shape>
    void ElementT<shape>::make_equal_to(const ElementT& el) {
      this->shape::make_equal_to(el);
      Element::make_equal_to(el);
    }
 
    template<class shape> void
    ElementT<shape>::draw() {
#ifndef SERVER
      Drawable::apply_color();
      glPushMatrix();
      glMultMatrixd(local_frame().get_glMatrix());
      shape::draw_gl();
      glPopMatrix();
#endif
    }
 
    template<class shape> void
    ElementT<shape>::draw_edge() {
#ifndef SERVER
      Drawable::apply_edge_color();
      glPushMatrix();
      glMultMatrixd(local_frame().get_glMatrix());
      shape::draw_gl_edge();
      glPopMatrix();
#endif
    }
 
    template<class shape> std::string
    ElementT<shape>::info() const {
      std::ostringstream os;
      os << Element::info() << '\n'
         << Body::info() << '\n'
         << shape::info();
      return os.str();
    }
 
    template<class shape> std::ostream&
    ElementT<shape>::write_ascii (std::ostream& out) const {
      Element::write_ascii(out);
      shape::write_ascii(out);
      return out;
    }
 
    template<class shape> std::istream&
    ElementT<shape>::read_ascii (std::istream& in) {
      Element::read_ascii(in);
      shape::read_ascii(in);
      set_density(1000.);
      return in;
    }
 
 
    template<class shape> bool
    ElementT<shape>::is_spherical_element() const {
      return false;
    }
 
    template<class shape> bool
    ElementT<shape>::is_polyhedral_element() const {
      return false;
    }
 
 
    template<class shape> template<class Archive> void
    ElementT<shape>::serialize(Archive & ar, const unsigned int version) {
      ar & boost::serialization::base_object<Element>(*this);
      ar & boost::serialization::base_object<shape>(*this);
    }
 
  }
}
 
namespace GranOO3
{
  namespace Core
  {
 
    template<class shape>
    struct Read<DEM::ElementT<shape> >
    {
      static Base* new_object(const std::string& classID, std::istream& in) {
        DEM::ElementT<shape>* pt = new DEM::ElementT<shape>();
        pt->read_ascii(in);
        return pt;
      }
 
    };
 
  }
}
 
 
 
namespace GranOO3
{
  GRANOO_CLASS_DECLARE_TPL(DEM::ElementT<Shape::Sphere>);
  GRANOO_CLASS_DECLARE_TPL(DEM::ElementT<Shape::Box>);
  GRANOO_CLASS_DECLARE_TPL(DEM::ElementT<Shape::Cylinder>);
  GRANOO_CLASS_DECLARE_TPL(DEM::ElementT<Shape::Cone>);
  GRANOO_CLASS_DECLARE_TPL(DEM::ElementT<Shape::Polyhedron>);
 
  extern template class Core::PropClass<DEM::ElementT<Shape::Sphere> >;
  extern template class Core::PropClass<DEM::ElementT<Shape::Box> >;
  extern template class Core::PropClass<DEM::ElementT<Shape::Cylinder> >;
  extern template class Core::PropClass<DEM::ElementT<Shape::Cone> >;
  extern template class Core::PropClass<DEM::ElementT<Shape::Polyhedron> >;
}
 
 
#endif