1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
/*************************************************************************
* 2021-2024 jerome.duriez@inrae.fr *
* This program is free software, see file LICENSE for details. *
*************************************************************************/
#ifdef YADE_LS_DEM
#pragma once
#include <core/IGeom.hpp>
#include <pkg/dem/ScGeom.hpp>
namespace yade {
class MultiScGeom : public IGeom {
public:
// hasNode and iteratorToNode below could certainly be const-declared but this may require a better handling of int / const int argument, maybe with respect to find in iteratorToNode which would expect a const int as last argument
bool hasNode(int); // to check whether some node (idx) is in nodesIds. Could there be sthg as simple as Python "in" ?
std::vector<int>::iterator iteratorToNode(int); // for locating some node (idx) in nodesIds
// clang-format off
YADE_CLASS_BASE_DOC_ATTRS_CTOR(MultiScGeom,IGeom,"A set of :yref:`ScGeom` for describing the kinematics of an interaction with multiple contact points between two :yref:`LevelSet` bodies, as a set of :yref:`ScGeom` items in :yref:`contacts<MultiScGeom.contacts>`. To combine with :yref:`MultiFrictPhys` and associated classes.",
((vector< shared_ptr<ScGeom> >,contacts,,,"The actual list of :yref:`ScGeom` items corresponding to the different contact points."))
((vector< int >,nodesIds,,,"List of :yref:`surface nodes<LevelSet.surfaceNodes>` (on id1 if that body is smaller -- or equal -- in volume, or id2 otherwise) making :yref:`contacts<MultiScGeom.contacts>`. Contact point for a node of index nodesIds[i] has kinematic properties stored in contacts[i]. Should be equal to :yref:`MultiFrictPhys.nodesIds` by design")) // do we need both ?"
, // existing Law2::go() will have to apply on a reference but contacts itself can not be a vector of references: https://stackoverflow.com/questions/922360/why-cant-i-make-a-vector-of-references
createIndex(); // this class will enter InteractionLoop dispatch, we need a create_index() here, and a REGISTER_*_INDEX below (https://yade-dem.org/doc/prog.html#indexing-dispatch-types)
);
// clang-format on
DECLARE_LOGGER;
REGISTER_CLASS_INDEX(MultiScGeom, IGeom); // see createIndex() remark
};
REGISTER_SERIALIZABLE(MultiScGeom);
class LSnodeGeom : public ScGeom {
public:
// clang-format off
YADE_CLASS_BASE_DOC_ATTRS_CTOR(LSnodeGeom,ScGeom,"Extends :yref:`ScGeom` for faster :yref:`LevelSet`-based simulations of convex bodies while storing previous closest node in :yref:`surfNodeIdx<LSnodeGeom.surfNodeIdx>`, for an optimal contact detection.",
((int,surfNodeIdx,-1,Attr::readonly,"Index (within :yref:`surfNodes<LevelSet.surfNodes>` of id2 if Volume(id1) > Volume(id2); else id1) of the surface node that shows the smallest interparticle distance (overlap or not). A negative value that would be encountered during an actual simulation would mean an inconsistency in implementation."))
((int,usedWorkflow,-1,Attr::readonly,"Tells for debugging purpose whether this interaction has just (in last Ig2 call) been handled by looping over full nodes because geom was just created (0) or because surfNodeIdx was directly detected as negative (1); or by looping just the neighbor subset (2, desired case), which could itself deteriorate into an additional full loop (3)."))
,
createIndex(); // same remark as MultiScGeom
);
// clang-format on
REGISTER_CLASS_INDEX(LSnodeGeom, ScGeom); // see createIndex() remark
};
REGISTER_SERIALIZABLE(LSnodeGeom);
} // namespace yade
#endif // YADE_LS_DEM
|
