dem_solver.hpp

Go to the documentation of this file.

#pragma once
 
#include "contact_solver_factory.hpp"
#include "dem_profiler.hpp"
 
namespace netdem {
 
class Simulation;
class MPIManager;
class ModifierManager;
class Scene;
 
class DEMSolver {
public:
  enum CyclePoint {
    pre,   
    mid_0, 
    mid_1, 
    mid_2, 
    mid_3, 
    mid_4, 
    post,  
    num_cycle_points 
  };
 
  double timestep{1.0e-4};
 
  DEMProfiler dem_profiler;
 
  ContactSolverFactory contact_solver_factory;
 
  DEMSolver();
 
  void Init(Simulation *sim);
 
  void UpdatePreModifiers();
 
  void UpdateLinkedList();
 
  void UpdateContacts();
 
  void UpdateParticleContactForce();
 
  void UpdateWallContactForce();
 
  void UpdateParticleMotion();
 
  void UpdateWallMotion();
 
  void UpdatePostModifiers();
 
  void Cycle(int num_cycles);
 
  void Solve(double time);
 
private:
  Simulation *sim{nullptr};
 
  void UpdateMidModifiers(CyclePoint cyc_point);
 
  void Cycle();
 
  void DryCycle();
 
  void SolveContactPP(Particle *const p_ii, NeighPofP *const neigh_tuple,
                      ContactSolverFactory *const solver_factory);
 
  void SolveContactPW(Particle *const p_ii, NeighWofP *const neigh_tuple,
                      ContactSolverFactory *const solver_factory);
};
 
} // namespace netdem