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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
/*************************************************************************
* Copyright (C) 2004 by Janek Kozicki *
* cosurgi@berlios.de *
* *
* This program is free software; it is licensed under the terms of the *
* GNU General Public License v2 or later. See file LICENSE for details. *
*************************************************************************/
#include <pkg/common/Collider.hpp>
namespace yade { // Cannot have #include directive inside.
YADE_PLUGIN((Collider));
bool Collider::mayCollide(
const Body* b1,
const Body* b2
#ifdef YADE_MPI
,
Body::id_t subdomain
#endif
)
{
return
#ifdef YADE_MPI //skip interactions outside subdomain, and between the subdomain and its own bodies
((subdomain == b1->subdomain or subdomain == b2->subdomain)
and not(b1->subdomain == b2->subdomain and (b1->getIsSubdomain() or b2->getIsSubdomain())))
&& handleFluidDomainCollision(b1, b2) &&
#endif
// might be called with deleted bodies, i.e. NULL pointers
(b1 != NULL && b2 != NULL) &&
// only collide if at least one particle is standalone or they belong to different clumps
(b1->isStandalone() || b2->isStandalone()
|| b1->clumpId != b2->clumpId
#ifdef YADE_SPH
// If SPH-mode enabled, we do not skip interactions between clump-members
// to get the correct calculation of density b->rho
|| true
#endif
)
&&
// do not collide clumps, since they are just containers, never interact
!b1->isClump() && !b2->isClump() &&
// masks must have at least 1 bit in common
b1->maskCompatible(b2->groupMask) &&
// avoid contact between particles having the same mask compatible with the avoidSelfInteractionMask.
!((b1->groupMask == b2->groupMask) && b1->maskCompatible(avoidSelfInteractionMask));
}
void Collider::pyHandleCustomCtorArgs(boost::python::tuple& t, boost::python::dict& /*d*/)
{
if (boost::python::len(t) == 0) return; // nothing to do
if (boost::python::len(t) != 1)
throw invalid_argument(("Collider optionally takes exactly one list of BoundFunctor's as non-keyword argument for constructor ("
+ boost::lexical_cast<string>(boost::python::len(t)) + " non-keyword ards given instead)")
.c_str());
typedef std::vector<shared_ptr<BoundFunctor>> vecBound;
vecBound vb = boost::python::extract<vecBound>(t[0])();
FOREACH(shared_ptr<BoundFunctor> bf, vb) this->boundDispatcher->add(bf);
t = boost::python::tuple(); // empty the args
}
#ifdef YADE_MPI
bool Collider::handleFluidDomainCollision(const Body* b1, const Body* b2) const
{
if (b1 == NULL && b2 == NULL) {
return false;
}
// skip interactions between the fluid solvers mesh bounds (openfoam coupling) and yade subdomains.
else if (b1->getIsFluidDomainBbox() && b2->getIsFluidDomainBbox()) {
return false;
} else {
return true;
}
}
#endif
} // namespace yade
|
