VoronoiBond.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_Prop_VoronoiBond_HPP
#define _LibDEM_Prop_VoronoiBond_HPP
 
#include "GranOO3/Core/Prop.hpp"
#include "GranOO3/DEM/Bond.hpp"
#include "GranOO3/DEM/ElementT.hpp"
#include "GranOO3/Shape/Polyhedron.hpp"
 
namespace GranOO3
{
  namespace DEM
  {
 
    template<class T> class ElementT;
 
    class VoronoiBond : public Core::Prop<Bond>,
            public Core::Register<VoronoiBond>
    {
 
      GRANOO_CLASS(DEM::VoronoiBond, Core::Prop<Bond>);
 
    public:
      // Constructors and destructor
      VoronoiBond();
      virtual ~VoronoiBond();
 
      void init();
      void update_element();
 
      GRANOO_ACCESS_GET(facet_center, Geom::Point, _facet_center);
                
      GRANOO_ACCESS(facet_vertex1, std::vector<double*>, _vertex1);
      GRANOO_ACCESS(facet_vertex2, std::vector<double*>, _vertex2);
      GRANOO_ACCESS(normal       , Geom::Vector       , _normal );
 
      GRANOO_ACCESS_REF      (el1, ElementT<Shape::Polyhedron>, *_el1);
      GRANOO_ACCESS_REF      (el2, ElementT<Shape::Polyhedron>, *_el2);
      GRANOO_ACCESS_REF_CONST(el1, ElementT<Shape::Polyhedron>, *_el1);
      GRANOO_ACCESS_REF_CONST(el2, ElementT<Shape::Polyhedron>, *_el2);
 
      GRANOO_ACCESS_GET(basefacet_index1, unsigned int, _basefacet_index1);
      GRANOO_ACCESS_GET(basefacet_index2, unsigned int, _basefacet_index2);
      
      Geom::Point compute_facet_center()  const;
      Geom::Point compute_facet_center1() const;
      Geom::Point compute_facet_center2() const;
 
      unsigned int get_basefacet_index_of(const Element&) const;
      bool is_inconsistent() const;
 
    private:
      friend class boost::serialization::access;
      template<class Archive> void serialize(Archive& ar, const unsigned int );
 
    private:
      ElementT<Shape::Polyhedron>* _el1;
      ElementT<Shape::Polyhedron>* _el2;
      unsigned int _basefacet_index1;
      unsigned int _basefacet_index2;
      Geom::Point _center;
      Geom::Point _facet_center;
      Geom::Vector _normal;
      std::vector<double*> _vertex1;
      std::vector<double*> _vertex2;
      bool _inconsistent;
    };
 
    template<class Archive> void
    VoronoiBond::serialize(Archive& ar, const unsigned int) {
      ar & boost::serialization::base_object<Core::Prop<Bond> >(*this);
      ar & _el1;
      ar & _el2;
      ar & _basefacet_index1;
      ar & _basefacet_index2;
      ar & _normal;
      ar & _center;
      ar & _facet_center;
    }
 
  }
}
 
namespace GranOO3
{
  GRANOO_CLASS_DECLARE_TPL(DEM::VoronoiBond);
}
 
 
 
#endif