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