| 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(); |
