| 98 |
|
AtomType* atype2; |
| 99 |
|
int atid1, atid2; |
| 100 |
|
|
| 101 |
< |
// We only need to worry about the types that are actually in the simulation: |
| 102 |
< |
|
| 101 |
> |
// We only need to worry about the types that are actually in the |
| 102 |
> |
// simulation: |
| 103 |
> |
|
| 104 |
|
set<AtomType*> atypes = info_->getSimulatedAtomTypes(); |
| 105 |
|
|
| 106 |
|
lj_->setSimulatedAtomTypes(atypes); |
| 426 |
|
|
| 427 |
|
int& sHash = sHash_[sdat.atid]; |
| 428 |
|
|
| 429 |
< |
if ((sHash & EAM_INTERACTION) != 0) eam_->calcFunctional(sdat); |
| 429 |
> |
if ((sHash & EAM_INTERACTION) != 0) eam_->calcFunctional(sdat); |
| 430 |
|
if ((sHash & SC_INTERACTION) != 0) sc_->calcFunctional(sdat); |
| 431 |
|
|
| 432 |
|
// set<NonBondedInteraction*>::iterator it; |
| 513 |
|
for (it = interactions_[*atid][*atid].begin(); |
| 514 |
|
it != interactions_[*atid][*atid].end(); |
| 515 |
|
++it) |
| 516 |
< |
cutoff = max(cutoff, (*it)->getSuggestedCutoffRadius(make_pair(atype, atype))); |
| 516 |
> |
cutoff = max(cutoff, (*it)->getSuggestedCutoffRadius(make_pair(atype, |
| 517 |
> |
atype))); |
| 518 |
|
return cutoff; |
| 519 |
|
} |
| 520 |
|
|
| 528 |
|
|
| 529 |
|
for (it = interactions_[atid][atid].begin(); |
| 530 |
|
it != interactions_[atid][atid].end(); ++it) |
| 531 |
< |
cutoff = max(cutoff, (*it)->getSuggestedCutoffRadius(make_pair(atype, atype))); |
| 531 |
> |
cutoff = max(cutoff, (*it)->getSuggestedCutoffRadius(make_pair(atype, |
| 532 |
> |
atype))); |
| 533 |
|
return cutoff; |
| 534 |
|
} |
| 535 |
|
} //end namespace OpenMD |