Between2SetOf.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_ContactDetection_Between2SetOf_hpp_
#define _libDEM_ContactDetection_Between2SetOf_hpp_
 
#include "GranOO3/Core/ObjectFactory.hpp"
 
#include "GranOO3/Core/NeedSetOf.hpp"
#include "GranOO3/DEM/ContactDetection.hpp"
#include "GranOO3/DEM/DiscreteElement.hpp"
 
//
// Base class for detection method
//
 
namespace GranOO3
{
  namespace DEM
  {
 
    class ContactDetection_Between2SetOf : public ContactDetection<DEM::DiscreteElement>,
                       public Core::NeedSetOf< DiscreteElement, 2>
    {
 
      GRANOO_OBJECT_FACTORY(ContactDetection_Between2SetOf)
 
    public:
      static std::string class_ID() {return "Between2SetOf";}
 
    public:
      //CONSTRUCTORS & DESTRUCTORS
      ContactDetection_Between2SetOf();
      virtual ~ContactDetection_Between2SetOf();
 
      //USEFULL
      virtual void detect_contact();
      virtual void parse_xml();
 
    private:
      ContactDetection_Between2SetOf(const ContactDetection_Between2SetOf&) = delete;
      ContactDetection_Between2SetOf& operator=(const ContactDetection_Between2SetOf&) = delete; //Not allowed
 
    private:
      unsigned int _contact_number;
    };
 
 
    inline
    ContactDetection_Between2SetOf::ContactDetection_Between2SetOf()
      : ContactDetection(),
      Core::NeedSetOf< DiscreteElement,2 >(true) {
    }
 
    inline
    ContactDetection_Between2SetOf::~ContactDetection_Between2SetOf() {
    }
 
    inline void
    ContactDetection_Between2SetOf::parse_xml() {
      ContactDetection::parse_xml();
      Core::NeedSetOf< DiscreteElement,2 >::parse_xml();
 
      Core::SetOf<DiscreteElement>& set1 = Core::NeedSetOf< DiscreteElement,2 >::get_setof(0);
      Core::SetOf<DiscreteElement>& set2 = Core::NeedSetOf< DiscreteElement,2 >::get_setof(1);
 
 
      AssertMsg(&set1 != &set2, "you can't assign the same SetOf for both");
 
    }
 
    inline void
    ContactDetection_Between2SetOf::detect_contact() {
      Core::SetOfBase<DiscreteElement>& set1 = Core::NeedSetOf<DiscreteElement, 2>::get_setof(0);
      Core::SetOfBase<DiscreteElement>& set2 = Core::NeedSetOf<DiscreteElement, 2>::get_setof(1);
      for (unsigned int i = 0; i<set1.size(); i++) {
        for (unsigned int j = 0; j<set2.size(); j++) {
          DEM::DiscreteElement& part1 = *set1[i];
          DEM::DiscreteElement& part2 = *set2[j];
          const Geom::Vector vDist = (part1.get_position() - part2.get_position());
          const double dist = vDist.norm();
          if (dist < part1.get_radius() + part2.get_radius()) {
            ContactDetection::Process(part1, part2, vDist/dist, part1.get_radius() + part2.get_radius()-dist);
          }
        }
      }
    }
 
  }
}
 
 
#endif