| 347 |
|
return; |
| 348 |
|
} |
| 349 |
|
|
| 350 |
< |
void EAM::calcDensity(DensityData ddat) { |
| 350 |
> |
void EAM::calcDensity(DensityData &ddat) { |
| 351 |
|
|
| 352 |
|
if (!initialized_) initialize(); |
| 353 |
|
|
| 363 |
|
return; |
| 364 |
|
} |
| 365 |
|
|
| 366 |
< |
void EAM::calcFunctional(FunctionalData fdat) { |
| 366 |
> |
void EAM::calcFunctional(FunctionalData &fdat) { |
| 367 |
|
|
| 368 |
|
if (!initialized_) initialize(); |
| 369 |
|
|
| 377 |
|
} |
| 378 |
|
|
| 379 |
|
|
| 380 |
< |
void EAM::calcForce(InteractionData idat) { |
| 380 |
> |
void EAM::calcForce(InteractionData &idat) { |
| 381 |
|
|
| 382 |
|
if (!initialized_) initialize(); |
| 383 |
|
|
| 478 |
|
idat.fshift1 = data1.F->getValueAt( idat.rho1 - rhb ); |
| 479 |
|
idat.fshift2 = data1.F->getValueAt( idat.rho2 - rha ); |
| 480 |
|
|
| 481 |
< |
idat.pot += phab; |
| 481 |
> |
idat.pot[3] += phab; |
| 482 |
|
|
| 483 |
< |
idat.vpair += phab; |
| 483 |
> |
idat.vpair[3] += phab; |
| 484 |
|
} |
| 485 |
|
|
| 486 |
|
return; |