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 |