Program Listing for File broadphase-callbacks.hpp

Return to documentation for file (include/pinocchio/collision/broadphase-callbacks.hpp)

//
// Copyright (c) 2022 INRIA
//

#ifndef __pinocchio_collision_broadphase_callback_hpp__
#define __pinocchio_collision_broadphase_callback_hpp__

#include <hpp/fcl/broadphase/broadphase_callbacks.h>

#include "pinocchio/multibody/fcl.hpp"
#include "pinocchio/multibody/geometry.hpp"

#include "pinocchio/collision/collision.hpp"

namespace pinocchio
{

  struct CollisionCallBackBase : hpp::fcl::CollisionCallBackBase
  {
    CollisionCallBackBase(const GeometryModel & geometry_model, GeometryData & geometry_data)
    : geometry_model_ptr(&geometry_model)
    , geometry_data_ptr(&geometry_data)
    , collision(false)
    , accumulate(false)
    {
    }

    const GeometryModel & getGeometryModel() const
    {
      return *geometry_model_ptr;
    }
    const GeometryData & getGeometryData() const
    {
      return *geometry_data_ptr;
    }
    GeometryData & getGeometryData()
    {
      return *geometry_data_ptr;
    }

    virtual bool stop() const = 0;

    virtual void done() {};

  protected:
    const GeometryModel * geometry_model_ptr;

    GeometryData * geometry_data_ptr;

  public:
    bool collision;

    bool accumulate;
  };

  struct CollisionCallBackDefault : CollisionCallBackBase
  {
    CollisionCallBackDefault(
      const GeometryModel & geometry_model,
      GeometryData & geometry_data,
      bool stopAtFirstCollision = false)
    : CollisionCallBackBase(geometry_model, geometry_data)
    , stopAtFirstCollision(stopAtFirstCollision)
    , count(0)
    //  , visited(Eigen::MatrixXd::Zero(geometry_model.ngeoms,geometry_model.ngeoms))
    {
    }

    void init()
    {
      if (accumulate) // skip reseting of the parameters
        return;

      count = 0;
      collision = false;
      collisionPairIndex = std::numeric_limits<PairIndex>::max();
      //    visited.setZero();
    }

    bool collide(hpp::fcl::CollisionObject * o1, hpp::fcl::CollisionObject * o2)
    {

      assert(!stop() && "must never happened");

      CollisionObject & co1 = reinterpret_cast<CollisionObject &>(*o1);
      CollisionObject & co2 = reinterpret_cast<CollisionObject &>(*o2);

      const Eigen::DenseIndex go1_index = (Eigen::DenseIndex)co1.geometryObjectIndex;
      const Eigen::DenseIndex go2_index = (Eigen::DenseIndex)co2.geometryObjectIndex;

      const GeometryModel & geometry_model = *geometry_model_ptr;

      PINOCCHIO_CHECK_INPUT_ARGUMENT(
        go1_index < (Eigen::DenseIndex)geometry_model.ngeoms && go1_index >= 0);
      PINOCCHIO_CHECK_INPUT_ARGUMENT(
        go2_index < (Eigen::DenseIndex)geometry_model.ngeoms && go2_index >= 0);

      const int pair_index = geometry_model.collisionPairMapping(go1_index, go2_index);
      if (pair_index == -1)
        return false;

      const GeometryData & geometry_data = *geometry_data_ptr;
      const CollisionPair & cp = geometry_model.collisionPairs[(PairIndex)pair_index];
      const bool do_collision_check =
        geometry_data.activeCollisionPairs[(PairIndex)pair_index]
        && !(
          geometry_model.geometryObjects[cp.first].disableCollision
          || geometry_model.geometryObjects[cp.second].disableCollision);
      if (!do_collision_check)
        return false;

      count++;

      fcl::CollisionRequest collision_request(
        geometry_data_ptr->collisionRequests[size_t(pair_index)]);
      collision_request.gjk_variant = fcl::GJKVariant::NesterovAcceleration;
      //    collision_request.gjk_initial_guess = fcl::GJKInitialGuess::BoundingVolumeGuess;

      if (
        co1.collisionGeometry().get()
          != geometry_model.geometryObjects[size_t(go1_index)].geometry.get()
        || co2.collisionGeometry().get()
             != geometry_model.geometryObjects[size_t(go2_index)].geometry.get())
        PINOCCHIO_THROW_PRETTY(
          std::logic_error, "go1: " << go1_index << " or go2: " << go2_index
                                    << " have not been updated and have missmatching pointers.");
      //    if(!(co1.collisionGeometry()->aabb_local.volume() < 0 ||
      //    co2.collisionGeometry()->aabb_local.volume() <0)) { // TODO(jcarpent): check potential
      //    bug
      //      collision_request.gjk_initial_guess = fcl::GJKInitialGuess::BoundingVolumeGuess;
      //    }

      bool res;
      try
      {
        res = computeCollision(
          *geometry_model_ptr, *geometry_data_ptr, (PairIndex)pair_index, collision_request);
      }
      catch (std::logic_error & e)
      {
        PINOCCHIO_THROW_PRETTY(
          std::logic_error, "Geometries with index go1: "
                              << go1_index << " or go2: " << go2_index
                              << " have produced an internal error within HPP-FCL.\n what:\n"
                              << e.what());
      }

      if (res && !collision)
      {
        collision = true;
        collisionPairIndex = (PairIndex)pair_index;
      }

      if (!stopAtFirstCollision)
        return false;
      else
        return res;
    }

    bool stop() const final
    {
      if (stopAtFirstCollision && collision)
        return true;

      return false;
    }

    void done() final
    {
      if (collision)
        geometry_data_ptr->collisionPairIndex = collisionPairIndex;
    }

    bool stopAtFirstCollision;

    PairIndex collisionPairIndex;

    size_t count;

    //  Eigen::MatrixXd visited;
  };

} // namespace pinocchio

/* --- Details -------------------------------------------------------------------- */
#include "pinocchio/collision/broadphase-callbacks.hxx"

#endif // ifndef __pinocchio_collision_broadphase_callback_hpp__