| 30 |
|
double mass; |
| 31 |
|
double lattice_constant; |
| 32 |
|
double eam_drho; // The distance between each of the points indexed by rho. |
| 33 |
+ |
double eam_rcut; // The cutoff radius for eam. |
| 34 |
|
double eam_dr; // The distance between each of the rho points. |
| 35 |
|
int eam_nrho; // Number of points indexed by rho |
| 36 |
|
int eam_nr; // The number of points based on r (Both Phi(r) and Rho(r)). |
| 36 |
– |
double eam_rcut; // The cutoff radius for eam. |
| 37 |
|
int eam_ident; // Atomic number |
| 38 |
|
int ident; |
| 39 |
|
int last; // 0 -> default |
| 76 |
|
void add( atomStruct &info, double *the_eam_rvals, |
| 77 |
|
double *the_eam_rhovals,double *the_eam_Frhovals ){ |
| 78 |
|
|
| 79 |
– |
int i; |
| 80 |
– |
|
| 79 |
|
// check for duplicates |
| 80 |
|
|
| 81 |
|
if( !strcmp( info.name, name ) ){ |
| 162 |
|
char* ffPath_env = "FORCE_PARAM_PATH"; |
| 163 |
|
char* ffPath; |
| 164 |
|
char temp[200]; |
| 167 |
– |
char errMsg[1000]; |
| 165 |
|
|
| 166 |
|
headAtomType = NULL; |
| 167 |
|
currentAtomType = NULL; |
| 179 |
|
// Init the atomStruct mpi type |
| 180 |
|
|
| 181 |
|
atomStruct atomProto; // mpiPrototype |
| 182 |
< |
int atomBC[3] = {15,4,6}; // block counts |
| 182 |
> |
int atomBC[3] = {15,5,5}; // block counts |
| 183 |
|
MPI_Aint atomDspls[3]; // displacements |
| 184 |
|
MPI_Datatype atomMbrTypes[3]; // member mpi types |
| 185 |
|
|
| 302 |
|
info.last = 1; // initialize last to have the last set. |
| 303 |
|
// if things go well, last will be set to 0 |
| 304 |
|
|
| 308 |
– |
int i; |
| 305 |
|
int identNum; |
| 306 |
|
double *eam_rvals; // Z of r values |
| 307 |
|
double *eam_rhovals; // rho of r values |
| 531 |
|
|
| 532 |
|
// initialize the atoms |
| 533 |
|
|
| 538 |
– |
|
| 539 |
– |
Atom* thisAtom; |
| 540 |
– |
|
| 534 |
|
for( i=0; i<nAtoms; i++ ){ |
| 535 |
|
|
| 536 |
|
currentAtomType = headAtomType->find( the_atoms[i]->getType() ); |