SupportShape.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_SupportShape_tpp_
#define _libDEM_SupportShape_tpp_
 
#include <string>
 
#include "GranOO3/Core/SetOf.hpp"
#include "GranOO3/Physic/Ground.hpp"
 
#include "GranOO3/DEM/DiscreteElement.hpp"
#include "GranOO3/DEM/SupportShapeBoundary.hpp"
#include "GranOO3/Shape/Volume.hpp"
#include "GranOO3/Core/XmlParser.hpp"
 
 
namespace GranOO3
{
 
  extern template class Core::SetOf<DEM::DiscreteElement>;
 
  namespace DEM
  {
 
    class SupportShape : public Physic::Ground,
                         public Core::Register< SupportShape >
    {
 
      GRANOO_CLASS(DEM::SupportShape, Physic::Ground);
      GRANOO_CLASS_DEFAULT_COLOR(.2, .8, .4, .2);
 
    public:
      static SupportShape& get(const std::string& ID);
      static bool exist(const std::string& ID);
      static SupportShape& glob(const std::string& ID);
      static SupportShape* new_object(const TiXmlElement* el);
 
    public:
      //CONSTRUCTORS & DESTRUCTORS
      SupportShape(const std::string& ID, const Geom::Frame &f);
      virtual ~SupportShape();
 
      //INTERNAL USE (Not documented)
      virtual void read_xml_element(const TiXmlElement* el);
      virtual void add_glob(const std::string&) = 0;
 
      //USEFULL
      virtual void draw();
      void update_average_discrete_element_radius();
      double get_continuous_volume() const; 
      double get_volumic_fraction() const;
      double set_continuous_density(double density);
      virtual std::string info() const;
 
      virtual const Shape::Volume& to_volume_shape() const = 0;
      virtual Shape::Volume& to_volume_shape() = 0;
      virtual std::string shape_class_ID() const = 0;
 
      //DISCRETE SHAPE MANAGEMENT
      bool is_attach() const;
      virtual void attach_to(Core::SetOf<DiscreteElement>& set, bool update = true, double dist_factor = 1.);
      virtual void recoverer_boundary();
      virtual void update() = 0;
      double get_average_discrete_element_radius() const;
      Core::SetOf< DiscreteElement >& get_all_discrete_element();
      const Core::SetOf< DiscreteElement >& get_all_discrete_element() const;
      double get_discrete_volume() const;
 
      SupportShapeBoundary& get_boundary();
      const SupportShapeBoundary& get_boundary() const;
 
      // IO
      virtual std::ostream& write_ascii (std::ostream& out) const;
      virtual std::istream& read_ascii (std::istream& in);
      virtual std::ostream& write_ascii_shape (std::ostream& out) const = 0; // for .agdd compatibility
      virtual std::istream& read_ascii_shape (std::istream& in) = 0;         // for .agdd compatibility
 
    private:
      SupportShape(const SupportShape  &) = delete;              //Not allowed
      SupportShape  & operator=(const SupportShape  &) = delete; //Not allowed
 
      //SERIALIZATION
      friend class boost::serialization::access;
      template<class Archive> void save(Archive&, const unsigned int ) const;
      template<class Archive> void load(Archive&, const unsigned int);
      BOOST_SERIALIZATION_SPLIT_MEMBER();
 
    protected:
      SupportShape();
 
    private:
      bool _attached;
      double _average_discrete_element_radius;
      Core::SetOf< DiscreteElement >* _all;
      SupportShapeBoundary _boundary;
    };
 
 
    inline SupportShape& SupportShape::get(const std::string& id) {
      Core::SetOf< SupportShape >& set = Core::SetOf< SupportShape >::get_global_set();
      for (unsigned int i = 0; i<set.size(); ++i) {
        if (set[i]->get_ID()==id)
          return *set[i];
      }
      UserAssert(0, "can't find SupportShape named " + id);
      return *set[0];
    }
 
    inline bool SupportShape::exist(const std::string& id) {
      Core::SetOf< SupportShape >& set = Core::SetOf< SupportShape >::get_global_set();
      for (unsigned int i = 0; i<set.size(); ++i) {
        if (set[i]->get_ID()==id)
          return true;
      }
      return false;
    }
 
    inline std::string
    SupportShape::info() const {
      std::ostringstream os;
      os << Core::Base::info() << "User Id : " << _ID << '\n';
      return os.str();
    }
 
    inline SupportShapeBoundary&
    SupportShape::get_boundary() {
      AssertMsg(_attached == true, "The support shape must be attached first");
      return _boundary;
    }
 
    inline const SupportShapeBoundary&
    SupportShape::get_boundary() const {
      AssertMsg(_attached == true, "The support shape must be attached first");
      return _boundary;
    }
 
    template<class Archive>
    void
    SupportShape::save(Archive& ar, const unsigned int) const {
      ar << boost::serialization::base_object<Physic::Ground>(*this);
      ar << _attached;
      ar << _average_discrete_element_radius;
      ar << _all;
      ar << _boundary;
    }
 
    template<class Archive>
    void
    SupportShape::load(Archive& ar, const unsigned int) {
      ar >> boost::serialization::base_object<Physic::Ground>(*this);
      ar >> _attached;
      ar >> _average_discrete_element_radius;
      ar >> _all;
      ar >> _boundary;
      local_frame().update_glMatrix();
    }
 
  }
}
 
namespace GranOO3
{
  GRANOO_CLASS_DECLARE_TPL(DEM::SupportShape);
}
 
 
 
#endif