95 |
|
indexArray.push_back(make_pair(i, atoms[i]->getGlobalIndex())); |
96 |
|
|
97 |
|
sort(indexArray.begin(), indexArray.end(), indexSortingCriterion); |
98 |
– |
|
99 |
– |
//for (int i = 0; i < mpiSim->getMyNlocal(); i++) { |
100 |
– |
// printf("node %d has global %d at local %d\n", worldRank, indexArray[i].second, indexArray[i].first); |
101 |
– |
//} |
102 |
– |
|
98 |
|
} |
99 |
|
|
100 |
|
#endif |
390 |
|
haveError = 0; |
391 |
|
which_atom = i; |
392 |
|
|
398 |
– |
//local_index = -1; |
399 |
– |
|
400 |
– |
//for (j=0; (j<mpiSim->getMyNlocal()) && (local_index < 0); j++) { |
401 |
– |
// if (atoms[j]->getGlobalIndex() == which_atom) local_index = j; |
402 |
– |
//} |
403 |
– |
|
404 |
– |
//if (local_index != -1) { |
405 |
– |
|
393 |
|
local_index = indexArray[currentIndex].first; |
394 |
|
|
395 |
|
if (which_atom == indexArray[currentIndex].second) { |
518 |
|
} |
519 |
|
which_atom = i; |
520 |
|
|
534 |
– |
//local_index = -1; |
535 |
– |
|
536 |
– |
//for (j=0; (j<mpiSim->getMyNlocal()) && (local_index < 0); j++) { |
537 |
– |
// if (atoms[j]->getGlobalIndex() == which_atom) local_index = j; |
538 |
– |
//} |
539 |
– |
|
540 |
– |
//if (local_index != -1) { |
541 |
– |
|
521 |
|
local_index = indexArray[currentIndex].first; |
522 |
|
|
523 |
|
if (which_atom == indexArray[currentIndex].second) { |