| 12 |
|
#include "UseTheForce/ForceFields.hpp" |
| 13 |
|
#include "primitives/SRI.hpp" |
| 14 |
|
#include "utils/simError.h" |
| 15 |
< |
|
| 15 |
> |
#include "types/AtomType.hpp" |
| 16 |
|
#include "UseTheForce/DarkSide/eam_interface.h" |
| 17 |
< |
#include "UseTheForce/DarkSide/atype_interface.h" |
| 17 |
> |
|
| 18 |
|
//#include "UseTheForce/fortranWrappers.hpp" |
| 19 |
|
|
| 20 |
|
#ifdef IS_MPI |
| 328 |
|
double *eam_rhovals; // rho of r values |
| 329 |
|
double *eam_Frhovals; // F of rho values |
| 330 |
|
char eamPotFile[1000]; |
| 331 |
+ |
|
| 332 |
|
|
| 333 |
|
|
| 334 |
|
bigSigma = 0.0; |
| 467 |
|
int isError; |
| 468 |
|
|
| 469 |
|
// dummy variables |
| 469 |
– |
int isLJ = 0; |
| 470 |
– |
int isDipole = 0; |
| 471 |
– |
int isSSD = 0; |
| 472 |
– |
int isGB = 0; |
| 473 |
– |
int isEAM = 1; |
| 474 |
– |
int isCharge = 0; |
| 475 |
– |
double dipole = 0.0; |
| 476 |
– |
double charge = 0.0; |
| 477 |
– |
double eamSigma = 0.0; |
| 478 |
– |
double eamEpslon = 0.0; |
| 470 |
|
|
| 471 |
|
currentAtomType = headAtomType->next; |
| 472 |
|
while( currentAtomType != NULL ){ |
| 473 |
|
|
| 474 |
|
if( currentAtomType->name[0] != '\0' ){ |
| 475 |
< |
isError = 0; |
| 476 |
< |
makeAtype( &(currentAtomType->ident), |
| 477 |
< |
&isLJ, |
| 478 |
< |
&isSSD, |
| 479 |
< |
&isDipole, |
| 480 |
< |
&isGB, |
| 481 |
< |
&isEAM, |
| 482 |
< |
&isCharge, |
| 483 |
< |
&eamEpslon, |
| 493 |
< |
&eamSigma, |
| 494 |
< |
&charge, |
| 495 |
< |
&dipole, |
| 496 |
< |
&isError ); |
| 497 |
< |
if( isError ){ |
| 498 |
< |
sprintf( painCave.errMsg, |
| 499 |
< |
"Error initializing the \"%s\" atom type in fortran\n", |
| 500 |
< |
currentAtomType->name ); |
| 501 |
< |
painCave.isFatal = 1; |
| 502 |
< |
simError(); |
| 503 |
< |
} |
| 504 |
< |
} |
| 475 |
> |
|
| 476 |
> |
AtomType* at = new AtomType(); |
| 477 |
> |
at->setIdent(currentAtomType->ident); |
| 478 |
> |
at->setName(currentAtomType->name); |
| 479 |
> |
at->setEAM(); |
| 480 |
> |
at->complete(); |
| 481 |
> |
|
| 482 |
> |
} |
| 483 |
> |
|
| 484 |
|
currentAtomType = currentAtomType->next; |
| 485 |
|
} |
| 486 |
|
|
