7 |
|
|
8 |
|
#ifdef IS_MPI |
9 |
|
#include <mpi.h> |
10 |
– |
#include <mpi++.h> |
10 |
|
#endif //is_mpi |
11 |
|
|
12 |
|
#include "ForceFields.hpp" |
357 |
|
int isDipole = 0; |
358 |
|
int isSSD = 0; |
359 |
|
int isGB = 0; |
361 |
– |
double w0 = 0.0; |
362 |
– |
double v0 = 0.0; |
360 |
|
double dipole = 0.0; |
364 |
– |
double GB_dummy = 0.0; |
365 |
– |
|
361 |
|
|
362 |
|
currentAtomType = headAtomType; |
363 |
|
while( currentAtomType != NULL ){ |
372 |
|
&(currentAtomType->epslon), |
373 |
|
&(currentAtomType->sigma), |
374 |
|
&dipole, |
380 |
– |
&w0, |
381 |
– |
&v0, |
382 |
– |
&GB_dummy, |
383 |
– |
&GB_dummy, |
384 |
– |
&GB_dummy, |
385 |
– |
&GB_dummy, |
386 |
– |
&GB_dummy, |
387 |
– |
&GB_dummy, |
375 |
|
&isError ); |
376 |
|
if( isError ){ |
377 |
|
sprintf( painCave.errMsg, |
423 |
|
|
424 |
|
if( bigSigma < currentAtomType->sigma ) bigSigma = currentAtomType->sigma; |
425 |
|
} |
439 |
– |
|
440 |
– |
|
441 |
– |
|
442 |
– |
#ifdef IS_MPI |
443 |
– |
sprintf( checkPointMsg, "LJ_FF atoms initialized succesfully" ); |
444 |
– |
MPIcheckPoint(); |
445 |
– |
#endif // is_mpi |
446 |
– |
|
426 |
|
} |
427 |
|
|
428 |
|
void LJ_FF::initializeBonds( int nBonds, Bond** BondArray, |
434 |
|
painCave.isFatal = 1; |
435 |
|
simError(); |
436 |
|
} |
458 |
– |
#ifdef IS_MPI |
459 |
– |
MPIcheckPoint(); |
460 |
– |
#endif // is_mpi |
461 |
– |
|
437 |
|
} |
438 |
|
|
439 |
|
void LJ_FF::initializeBends( int nBends, Bend** bendArray, |
445 |
|
painCave.isFatal = 1; |
446 |
|
simError(); |
447 |
|
} |
473 |
– |
#ifdef IS_MPI |
474 |
– |
MPIcheckPoint(); |
475 |
– |
#endif // is_mpi |
476 |
– |
|
448 |
|
} |
449 |
|
|
450 |
|
void LJ_FF::initializeTorsions( int nTorsions, Torsion** torsionArray, |
456 |
|
painCave.isFatal = 1; |
457 |
|
simError(); |
458 |
|
} |
488 |
– |
#ifdef IS_MPI |
489 |
– |
MPIcheckPoint(); |
490 |
– |
#endif // is_mpi |
491 |
– |
|
459 |
|
} |
460 |
|
|
461 |
|
void LJ_FF::fastForward( char* stopText, char* searchOwner ){ |