| 14 |
|
|
| 15 |
|
#include "MatVec3.h" |
| 16 |
|
|
| 17 |
+ |
#include "ConstraintManager.hpp" |
| 18 |
+ |
|
| 19 |
|
#ifdef IS_MPI |
| 20 |
|
#include "mpiSimulation.hpp" |
| 21 |
|
#endif |
| 78 |
|
|
| 79 |
|
ngroup = 0; |
| 80 |
|
|
| 81 |
+ |
consMan = NULL; |
| 82 |
+ |
|
| 83 |
|
wrapMeSimInfo( this ); |
| 84 |
|
} |
| 85 |
|
|
| 92 |
|
|
| 93 |
|
for(i = properties.begin(); i != properties.end(); i++) |
| 94 |
|
delete (*i).second; |
| 95 |
< |
|
| 95 |
> |
|
| 96 |
> |
if (!consMan) |
| 97 |
> |
delete consMan; |
| 98 |
|
} |
| 99 |
|
|
| 100 |
|
void SimInfo::setBox(double newBox[3]) { |
| 199 |
|
|
| 200 |
|
if( orthoRhombic ) { |
| 201 |
|
sprintf( painCave.errMsg, |
| 202 |
< |
"\n\tOOPSE is switching from the default Non-Orthorhombic\n" |
| 202 |
> |
"OOPSE is switching from the default Non-Orthorhombic\n" |
| 203 |
|
"\tto the faster Orthorhombic periodic boundary computations.\n" |
| 204 |
|
"\tThis is usually a good thing, but if you wan't the\n" |
| 205 |
|
"\tNon-Orthorhombic computations, make the orthoBoxTolerance\n" |
| 210 |
|
} |
| 211 |
|
else { |
| 212 |
|
sprintf( painCave.errMsg, |
| 213 |
< |
"\n\tOOPSE is switching from the faster Orthorhombic to the more\n" |
| 213 |
> |
"OOPSE is switching from the faster Orthorhombic to the more\n" |
| 214 |
|
"\tflexible Non-Orthorhombic periodic boundary computations.\n" |
| 215 |
|
"\tThis is usually because the box has deformed under\n" |
| 216 |
|
"\tNPTf integration. If you wan't to live on the edge with\n" |
| 474 |
|
sprintf( painCave.errMsg, |
| 475 |
|
"There was an error setting the simulation information in fortran.\n" ); |
| 476 |
|
painCave.isFatal = 1; |
| 477 |
+ |
painCave.severity = OOPSE_ERROR; |
| 478 |
|
simError(); |
| 479 |
|
} |
| 480 |
|
|
| 513 |
|
|
| 514 |
|
if( rCut > maxCutoff ){ |
| 515 |
|
sprintf( painCave.errMsg, |
| 516 |
< |
"\n\tcutoffRadius is too large for the current periodic box.\n" |
| 516 |
> |
"cutoffRadius is too large for the current periodic box.\n" |
| 517 |
|
"\tCurrent Value of cutoffRadius = %G at time %G\n " |
| 518 |
|
"\tThis is larger than half of at least one of the\n" |
| 519 |
|
"\tperiodic box vectors. Right now, the Box matrix is:\n" |
| 532 |
|
} else { |
| 533 |
|
// initialize this stuff before using it, OK? |
| 534 |
|
sprintf( painCave.errMsg, |
| 535 |
< |
"\n\tTrying to check cutoffs without a box.\n" |
| 535 |
> |
"Trying to check cutoffs without a box.\n" |
| 536 |
|
"\tOOPSE should have better programmers than that.\n" ); |
| 537 |
|
painCave.severity = OOPSE_ERROR; |
| 538 |
|
painCave.isFatal = 1; |