contact_solver_factory.hpp

Go to the documentation of this file.

#pragma once
 
#include "bond_solver_pp.hpp"
#include "bond_solver_pw.hpp"
#include "collision_solver_pp.hpp"
#include "collision_solver_pw.hpp"
#include "utils_math.hpp"
#include <memory>
#include <unordered_map>
 
namespace netdem {
 
class ContactSolverSettings {
public:
  enum SolverType { gjk, sdf, automatic };
 
  SolverType solver_type{SolverType::automatic};
 
  bool gjk_use_erosion{false};
 
  double gjk_erosion_ratio_initial{0.01};
 
  double gjk_erosion_ratio_increment{0.01};
 
  bool gjk_enable_logging{false};
 
  int sdf_potential_type{0};
 
  bool sdf_solve_two_sides{false};
 
  bool sdf_use_equivalent_stiffness{true};
};
 
class ContactSolverFactory {
public:
  ContactSolverSettings settings;
 
private:
  BondSolverPP bond_solver_pp;
 
  BondSolverPW bond_solver_pw;
 
  VecNT<VecNT<int, Shape::Type::num_shapes>, Shape::Type::num_shapes>
      solver_id_pp_list, solver_id_pw_list;
 
  std::unordered_map<std::pair<int, int>, int, PairHash> solver_id_customized;
 
  VecXT<CollisionSolverPP *> solver_pp_pool;
 
  VecXT<CollisionSolverPW *> solver_pw_pool;
 
public:
  ContactSolverFactory();
 
  ContactSolverFactory(ContactSolverFactory &tmp_factory);
 
  ContactSolverFactory(ContactSolverFactory &&tmp_factory);
 
  ContactSolverFactory &operator=(ContactSolverFactory &tmp_factory);
 
  ContactSolverFactory &operator=(ContactSolverFactory &&tmp_factory);
 
  BondSolverPP *GetBondSolver(Particle *const p1, Particle *const p2);
 
  BondSolverPW *GetBondSolver(Particle *const p, Wall *const w);
 
  CollisionSolverPP *GetCollisionSolver(Particle *const p1, Particle *const p2);
 
  CollisionSolverPW *GetCollisionSolver(Particle *const p, Wall *const w);
 
  int InsertSolver(CollisionSolverPP *const cnt_solver);
 
  int InsertSolver(CollisionSolverPW *const cnt_solver);
 
  void CustomizeSolverPP(int shape_id1, int shape_id2, int solver_id);
 
  void CustomizeSolverPW(int shape_id1, int shape_id2, int solver_id);
 
  ~ContactSolverFactory();
 
private:
  CollisionSolverPP *NewCollisionSolver(Particle *const p1, Particle *const p2);
 
  CollisionSolverPW *NewCollisionSolver(Particle *const p, Wall *const w);
};
 
} // namespace netdem