484 |
|
} |
485 |
|
|
486 |
|
|
487 |
< |
//creat cutoff group for molecule |
487 |
> |
//create cutoff group for molecule |
488 |
|
molInfo.myCutoffGroups.clear(); |
489 |
|
for (j = 0; j < molInfo.nCutoffGroups; j++){ |
490 |
|
|
563 |
|
sprintf(checkPointMsg, "all molecules initialized succesfully"); |
564 |
|
MPIcheckPoint(); |
565 |
|
#endif // is_mpi |
566 |
– |
|
567 |
– |
// clean up the forcefield |
568 |
– |
|
569 |
– |
if (!globals->haveRcut()){ |
566 |
|
|
571 |
– |
the_ff->calcRcut(); |
572 |
– |
|
573 |
– |
} else { |
574 |
– |
|
575 |
– |
the_ff->setRcut( globals->getRcut() ); |
576 |
– |
} |
577 |
– |
|
578 |
– |
the_ff->cleanMe(); |
567 |
|
} |
568 |
|
|
569 |
|
void SimSetup::initFromBass(void){ |
1004 |
|
#endif //is_mpi |
1005 |
|
|
1006 |
|
double theRcut, theRsw; |
1007 |
+ |
|
1008 |
+ |
if (globals->haveRcut()) { |
1009 |
+ |
theRcut = globals->getRcut(); |
1010 |
+ |
|
1011 |
+ |
if (globals->haveRsw()) |
1012 |
+ |
theRsw = globals->getRsw(); |
1013 |
+ |
else |
1014 |
+ |
theRsw = theRcut; |
1015 |
+ |
|
1016 |
+ |
info[i].setDefaultRcut(theRcut, theRsw); |
1017 |
+ |
|
1018 |
+ |
} else { |
1019 |
+ |
|
1020 |
+ |
the_ff->calcRcut(); |
1021 |
+ |
theRcut = info[i].getRcut(); |
1022 |
+ |
|
1023 |
+ |
if (globals->haveRsw()) |
1024 |
+ |
theRsw = globals->getRsw(); |
1025 |
+ |
else |
1026 |
+ |
theRsw = theRcut; |
1027 |
+ |
|
1028 |
+ |
info[i].setDefaultRcut(theRcut, theRsw); |
1029 |
+ |
} |
1030 |
|
|
1031 |
|
if (globals->getUseRF()){ |
1032 |
|
info[i].useReactionField = 1; |
1033 |
< |
|
1033 |
> |
|
1034 |
|
if (!globals->haveRcut()){ |
1035 |
|
sprintf(painCave.errMsg, |
1036 |
|
"SimSetup Warning: No value was set for the cutoffRadius.\n" |
1107 |
|
strcpy(checkPointMsg, "post processing checks out"); |
1108 |
|
MPIcheckPoint(); |
1109 |
|
#endif // is_mpi |
1110 |
+ |
|
1111 |
+ |
// clean up the forcefield |
1112 |
+ |
the_ff->cleanMe(); |
1113 |
|
} |
1114 |
|
|
1115 |
|
void SimSetup::initSystemCoords(void){ |