454 |
|
char* ffPath_env = "FORCE_PARAM_PATH"; |
455 |
|
char* ffPath; |
456 |
|
char temp[200]; |
457 |
– |
char errMsg[1000]; |
457 |
|
|
458 |
|
headAtomType = NULL; |
459 |
|
currentAtomType = NULL; |
475 |
|
// Init the atomStruct mpi type |
476 |
|
|
477 |
|
atomStruct atomProto; // mpiPrototype |
478 |
< |
int atomBC[3] = {15,6,4}; // block counts |
478 |
> |
int atomBC[3] = {15,11,4}; // block counts |
479 |
|
MPI_Aint atomDspls[3]; // displacements |
480 |
|
MPI_Datatype atomMbrTypes[3]; // member mpi types |
481 |
|
|
657 |
|
|
658 |
|
void DUFF::readParams( void ){ |
659 |
|
|
661 |
– |
int i, a, b, c, d; |
660 |
|
int identNum; |
663 |
– |
char* atomA; |
664 |
– |
char* atomB; |
665 |
– |
char* atomC; |
666 |
– |
char* atomD; |
661 |
|
|
662 |
|
atomStruct atomInfo; |
663 |
|
bondStruct bondInfo; |
779 |
|
int isGB = 0; |
780 |
|
int isLJ = 1; |
781 |
|
int isEAM =0; |
782 |
< |
double GB_dummy = 0.0; |
789 |
< |
|
790 |
< |
|
782 |
> |
|
783 |
|
currentAtomType = headAtomType->next;; |
784 |
|
while( currentAtomType != NULL ){ |
785 |
|
|
1302 |
|
} |
1303 |
|
|
1304 |
|
gBend = new GhostBend( *the_atoms[a], |
1305 |
< |
*the_atoms[b] ); |
1305 |
> |
*the_atoms[b]); |
1306 |
> |
|
1307 |
|
gBend->setConstants( currentBendType->k1, |
1308 |
|
currentBendType->k2, |
1309 |
|
currentBendType->k3, |
1319 |
|
currentBendType->k3, |
1320 |
|
currentBendType->t0 ); |
1321 |
|
bendArray[i] = qBend; |
1322 |
< |
} |
1322 |
> |
} |
1323 |
|
} |
1324 |
|
} |
1325 |
|
} |
1631 |
|
else{ |
1632 |
|
sprintf( painCave.errMsg, |
1633 |
|
"Unknown DUFF bond type \"%s\" at line %d\n", |
1634 |
< |
info.type, |
1634 |
> |
bondType, |
1635 |
|
lineNum ); |
1636 |
|
painCave.isFatal = 1; |
1637 |
|
simError(); |