520 |
|
//return the position of joint atom apears in rigidbody's definition |
521 |
|
//for the time being, we will use the most inefficient algorithm, the complexity is O(N2) |
522 |
|
//actually we could improve the complexity to O(NlgN) by sorting the atom index in rigid body first |
523 |
< |
vector<pair<int, int> > MoleculeStamp::getJointAtoms(int rb1, int rb2){ |
523 |
> |
std::vector<pair<int, int> > MoleculeStamp::getJointAtoms(int rb1, int rb2){ |
524 |
|
RigidBodyStamp* rbStamp1; |
525 |
|
RigidBodyStamp* rbStamp2; |
526 |
|
int natomInRb1; |
527 |
|
int natomInRb2; |
528 |
|
int atomIndex1; |
529 |
|
int atomIndex2; |
530 |
< |
vector<pair<int, int> > jointAtomIndexPair; |
530 |
> |
std::vector<pair<int, int> > jointAtomIndexPair; |
531 |
|
|
532 |
|
rbStamp1 = this->getRigidBody(rb1); |
533 |
|
natomInRb1 =rbStamp1->getNMembers(); |