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
|
/////////////////////////////////////////////////////////////
// //
// Copyright (c) 2003-2014 by The University of Queensland //
// Centre for Geoscience Computing //
// http://earth.uq.edu.au/centre-geoscience-computing //
// //
// Primary Business: Brisbane, Queensland, Australia //
// Licensed under the Open Software License version 3.0 //
// http://www.apache.org/licenses/LICENSE-2.0 //
// //
/////////////////////////////////////////////////////////////
#include <mpi.h>
#include "Model/RotThermPairInteraction.h"
ARotThermPairInteraction::~ARotThermPairInteraction()
{
}
ARotThermPairInteraction::ARotThermPairInteraction():AInteraction()
{
m_p1=NULL;
m_p2=NULL;
m_id.push_back(-1);
m_id.push_back(-1);
}
ARotThermPairInteraction::ARotThermPairInteraction(CRotThermParticle* p1,CRotThermParticle* p2) : AInteraction()
{
setPP(p1,p2);
m_init=true;
}
void ARotThermPairInteraction::setPP(CRotThermParticle *p1,CRotThermParticle *p2)
{
m_p1=p1;
m_p2=p2;
m_id.clear();
m_id.push_back(m_p1->getID());
m_id.push_back(m_p2->getID());
}
void ARotThermPairInteraction::setPP(const vector<CRotThermParticle*> pp)
{
setPP(pp[0], pp[1]);
}
void ARotThermPairInteraction::checkIDs()
{
if ((m_id[0]!=m_p1->getID())||(m_id[1]!=m_p2->getID())){
std::cout
<< "inconsistent IDs : " << m_id[0] << "-" << m_id[1]
<< " vs. "
<< m_p1->getID() << "-" << m_p2->getID() << std::endl << std::flush;
}
}
/*!
check if any of the particles in the interaction fits tag & mask
\param tag the tag
\param mask the mask
*/
bool ARotThermPairInteraction::hasTag(int tag ,int mask) const
{
int tag1=m_p1->getTag();
int tag2=m_p2->getTag();
return (((tag1 & mask)==(tag & mask))||((tag2 & mask)==(tag & mask)));
}
|
