Polyhedron.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 _GranOO_LibShape__Polyhedron_Corehh_
#define _GranOO_LibShape__Polyhedron_Corehh_
 
#include <iostream>
#include <fstream>
#include <string>
 
#include <boost/serialization/string.hpp>
#include <boost/serialization/vector.hpp>
 
#include <vtkUnstructuredGrid.h>
 
 
#include "GranOO3/Algo/PolyhedronCore.hpp"
#include "GranOO3/Geom/Point.hpp"
#include "GranOO3/Shape/Sphere.hpp"
#include "GranOO3/Shape/Volume.hpp"
#include "GranOO3/Core/ObjectFactory.hpp"
#include "GranOO3/Core/String.hpp"
 
namespace GranOO3
{
  namespace Shape
  {
    class Polyhedron : public Volume, public Core::ObjectFactoryInterface<Base, Polyhedron>
    {
      GRANOO_OBJECT_FACTORY(Polyhedron);
 
    public:
      static std::string class_ID() {return std::string("Polyhedron");}
      static std::string desc();
 
    public:
      //CONSTRUCTORS & DESTRUCTORS
      Polyhedron();
      Polyhedron(std::string fileName, const Geom::Frame &frame);
      Polyhedron(const Polyhedron &polyhedron, const Geom::Frame &frame);
      Polyhedron(const Geom::Frame &frame);
      Polyhedron(const Polyhedron &s);
      virtual ~Polyhedron();
 
      //IO
      void read_xml_element(const TiXmlElement *element);
      virtual std::ostream &write_ascii(std::ostream &out) const;
      virtual std::istream &read_ascii(std::istream &in);
      virtual void to_vtk(vtkUnstructuredGrid *data);
      std::string info() const;
      void to_vtk(const unsigned int shift, unsigned int* conn, unsigned int* cell_type, double* point) const;
      unsigned int vtk_conn_number() const;
 
      //CLONING
      void make_equal_to(const Polyhedron&);
 
      //BOUNDING SHAPE
      void update_bounding_box();
      void update_bounding_sphere();
      double get_greatest_dimension() const;
      AABB aabb() const;
 
      //USEFULL
      void set_input_file(std::stringstream &stringStream, bool update = true, bool translate_to_centroid = true);
      void set_input_file(const std::string &offFile, bool update = true, bool translate_to_centroid = true);
      void set_output_file(const std::string &offFile);
      void set_input_offset_file(const std::string &offFile);
      double get_volume() const;
      void compute_inertia_tensor(double density, Geom::Tensor &) const;
      Geom::Point farthest_point_along(const Geom::Vector &) const;
      bool basefacet_from_normal(const Geom::Vector &normal, unsigned int &basefacetIndex) const;
      Geom::Point basefacet_center(const unsigned int basefacetIndex) const;
 
      double area_of_facet_crossed_by_vector(const Geom::Point &, const Geom::Vector &, unsigned int &, Geom::Vector &) const;
      Geom::Point initial_centroid() const;
      PolyhedronCore& get_polyhedron_core();
      const PolyhedronCore& get_polyhedron_core() const;
 
      const std::map<std::string, std::vector<int>> &get_facet_map() const { return _map_facet; };
      
      
      void draw_gl() const;
      void draw_gl_edge() const;
      void scale(double ratio);
 
      unsigned int get_number_of_vertex() const;
      unsigned int get_number_of_facet() const;
      
      //CONTACT DETECTION
      Interference locate(const Geom::Point &) const;
      Interference locate(const Geom::Point &, double radius) const;
      bool overlap(const Sphere &, double margin, std::vector<std::string> &facet) const;
      double compute_distance(const Geom::Point &);
 
      //CONTACT DETECTION v2
      GRANOO_DECLARE_SHAPE_COLLISION();
      bool collide_inside(const Geom::Point& p1, double radius, Collision::Data& info) const;
      bool collide_outside(const Geom::Point& p1, double radius, Collision::Data& info) const;
 
      std::ostream& to_povray(std::ostream& out, const Core::Color& ) const;
 
    private:
      Polyhedron & operator = (const Polyhedron &); // copy constructor forbidden
 
      //SERIALIZATION
      friend class boost::serialization::access;
      template<class Archive> void save(Archive& ar, const unsigned int version) const;
      template<class Archive> void load(Archive& ar, const unsigned int version);
      BOOST_SERIALIZATION_SPLIT_MEMBER();
 
    private:
      PolyhedronCore _polyhedron_core;
      std::map<std::string, std::vector<int> > _map_facet;
    };
 
 
    inline PolyhedronCore&
    Polyhedron::get_polyhedron_core() {
      return _polyhedron_core;
    }
 
    inline const
    PolyhedronCore&
    Polyhedron::get_polyhedron_core() const {
      return _polyhedron_core;
    }
 
    template<class Archive>
    inline void
    Polyhedron::save(Archive& ar, const unsigned int version) const {
      ar << boost::serialization::base_object<Volume>(*this);
 
      std::stringstream out;
      _polyhedron_core.to_off_format(out);
 
      std::string str = out.str();
      ar << str;
    }
 
    template<class Archive>
    inline void
    Polyhedron::load(Archive& ar, const unsigned int version) {
      ar >> boost::serialization::base_object<Volume>(*this);
 
      std::string str;
      ar >> str;
      std::stringstream stream(str);
      _polyhedron_core = PolyhedronCore(stream);
      update_bounding_shape();
    }
 
  }
}
 
//Don't forget to export for dynamic serialization of child class
#include <boost/serialization/export.hpp>
BOOST_CLASS_EXPORT_KEY(GranOO3::Shape::Polyhedron)
 
#include <boost/serialization/version.hpp>
BOOST_CLASS_VERSION(GranOO3::Shape::Polyhedron, 0)
GRANOO_SERIALIZE_SHAPE(Polyhedron)
 
 
namespace GranOO3
{
  extern template class Core::ObjectFactoryInterface<Shape::Base, Shape::Polyhedron>;
}
 
 
#endif