65 |
|
if( !strcmp( info.name, name ) ){ |
66 |
|
sprintf( painCave.errMsg, |
67 |
|
"Duplicate LJ atom type \"%s\" found in " |
68 |
< |
"the LJ_FF param file./n", |
68 |
> |
"the LJFF param file./n", |
69 |
|
name ); |
70 |
|
painCave.isFatal = 1; |
71 |
|
simError(); |
117 |
|
//**************************************************************** |
118 |
|
|
119 |
|
|
120 |
< |
LJ_FF::LJ_FF(){ |
120 |
> |
LJFF::LJFF(){ |
121 |
|
|
122 |
|
char fileName[200]; |
123 |
|
char* ffPath_env = "FORCE_PARAM_PATH"; |
162 |
|
|
163 |
|
// generate the force file name |
164 |
|
|
165 |
< |
strcpy( fileName, "LJ_FF.frc" ); |
166 |
< |
// fprintf( stderr,"Trying to open %s\n", fileName ); |
165 |
> |
strcpy( fileName, "LJFF.frc" ); |
166 |
> |
fprintf( stderr,"Trying to open %s\n", fileName ); |
167 |
|
|
168 |
|
// attempt to open the file in the current directory first. |
169 |
|
|
201 |
|
#ifdef IS_MPI |
202 |
|
} |
203 |
|
|
204 |
< |
sprintf( checkPointMsg, "LJ_FF file opened sucessfully." ); |
204 |
> |
sprintf( checkPointMsg, "LJFF file opened sucessfully." ); |
205 |
|
MPIcheckPoint(); |
206 |
|
|
207 |
|
#endif // is_mpi |
208 |
|
} |
209 |
|
|
210 |
|
|
211 |
< |
LJ_FF::~LJ_FF(){ |
211 |
> |
LJFF::~LJFF(){ |
212 |
|
|
213 |
|
if( headAtomType != NULL ) delete headAtomType; |
214 |
|
|
223 |
|
#endif // is_mpi |
224 |
|
} |
225 |
|
|
226 |
< |
void LJ_FF::initForceField( int ljMixRule ){ |
226 |
> |
void LJFF::initForceField( int ljMixRule ){ |
227 |
|
|
228 |
|
initFortran( ljMixRule, 0 ); |
229 |
|
} |
230 |
|
|
231 |
< |
void LJ_FF::cleanMe( void ){ |
231 |
> |
void LJFF::cleanMe( void ){ |
232 |
|
|
233 |
|
#ifdef IS_MPI |
234 |
|
|
243 |
|
#endif // is_mpi |
244 |
|
} |
245 |
|
|
246 |
< |
void LJ_FF::readParams( void ){ |
246 |
> |
void LJFF::readParams( void ){ |
247 |
|
|
248 |
|
atomStruct info; |
249 |
|
info.last = 1; // initialize last to have the last set. |
303 |
|
// send out the linked list to all the other processes |
304 |
|
|
305 |
|
sprintf( checkPointMsg, |
306 |
< |
"LJ_FF atom structures read successfully." ); |
306 |
> |
"LJFF atom structures read successfully." ); |
307 |
|
MPIcheckPoint(); |
308 |
|
|
309 |
|
currentAtomType = headAtomType->next; //skip the first element who is a place holder. |
388 |
|
|
389 |
|
#ifdef IS_MPI |
390 |
|
sprintf( checkPointMsg, |
391 |
< |
"LJ_FF atom structures successfully sent to fortran\n" ); |
391 |
> |
"LJFF atom structures successfully sent to fortran\n" ); |
392 |
|
MPIcheckPoint(); |
393 |
|
#endif // is_mpi |
394 |
|
|
395 |
|
} |
396 |
|
|
397 |
|
|
398 |
< |
void LJ_FF::initializeAtoms( int nAtoms, Atom** the_atoms ){ |
398 |
> |
void LJFF::initializeAtoms( int nAtoms, Atom** the_atoms ){ |
399 |
|
|
400 |
|
int i; |
401 |
|
|
425 |
|
} |
426 |
|
} |
427 |
|
|
428 |
< |
void LJ_FF::initializeBonds( int nBonds, Bond** BondArray, |
428 |
> |
void LJFF::initializeBonds( int nBonds, Bond** BondArray, |
429 |
|
bond_pair* the_bonds ){ |
430 |
|
|
431 |
|
if( nBonds ){ |
432 |
|
sprintf( painCave.errMsg, |
433 |
< |
"LJ_FF does not support bonds.\n" ); |
433 |
> |
"LJFF does not support bonds.\n" ); |
434 |
|
painCave.isFatal = 1; |
435 |
|
simError(); |
436 |
|
} |
437 |
|
} |
438 |
|
|
439 |
< |
void LJ_FF::initializeBends( int nBends, Bend** bendArray, |
439 |
> |
void LJFF::initializeBends( int nBends, Bend** bendArray, |
440 |
|
bend_set* the_bends ){ |
441 |
|
|
442 |
|
if( nBends ){ |
443 |
|
sprintf( painCave.errMsg, |
444 |
< |
"LJ_FF does not support bends.\n" ); |
444 |
> |
"LJFF does not support bends.\n" ); |
445 |
|
painCave.isFatal = 1; |
446 |
|
simError(); |
447 |
|
} |
448 |
|
} |
449 |
|
|
450 |
< |
void LJ_FF::initializeTorsions( int nTorsions, Torsion** torsionArray, |
450 |
> |
void LJFF::initializeTorsions( int nTorsions, Torsion** torsionArray, |
451 |
|
torsion_set* the_torsions ){ |
452 |
|
|
453 |
|
if( nTorsions ){ |
454 |
|
sprintf( painCave.errMsg, |
455 |
< |
"LJ_FF does not support torsions.\n" ); |
455 |
> |
"LJFF does not support torsions.\n" ); |
456 |
|
painCave.isFatal = 1; |
457 |
|
simError(); |
458 |
|
} |
459 |
|
} |
460 |
|
|
461 |
< |
void LJ_FF::fastForward( char* stopText, char* searchOwner ){ |
461 |
> |
void LJFF::fastForward( char* stopText, char* searchOwner ){ |
462 |
|
|
463 |
|
int foundText = 0; |
464 |
|
char* the_token; |