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