--- trunk/OOPSE/libmdtools/SimInfo.cpp 2004/06/01 18:42:58 1214 +++ trunk/OOPSE/libmdtools/SimInfo.cpp 2004/06/02 14:21:54 1218 @@ -191,25 +191,27 @@ void SimInfo::calcHmatInv( void ) { if( oldOrtho != orthoRhombic ){ - if( orthoRhombic ){ + if( orthoRhombic ) { sprintf( painCave.errMsg, - "OOPSE is switching from the default Non-Orthorhombic\n" + "\n\tOOPSE is switching from the default Non-Orthorhombic\n" "\tto the faster Orthorhombic periodic boundary computations.\n" "\tThis is usually a good thing, but if you wan't the\n" "\tNon-Orthorhombic computations, make the orthoBoxTolerance\n" "\tvariable ( currently set to %G ) smaller.\n", orthoTolerance); + painCave.severity = OOPSE_INFO; simError(); } else { sprintf( painCave.errMsg, - "OOPSE is switching from the faster Orthorhombic to the more\n" + "\n\tOOPSE is switching from the faster Orthorhombic to the more\n" "\tflexible Non-Orthorhombic periodic boundary computations.\n" "\tThis is usually because the box has deformed under\n" "\tNPTf integration. If you wan't to live on the edge with\n" "\tthe Orthorhombic computations, make the orthoBoxTolerance\n" "\tvariable ( currently set to %G ) larger.\n", orthoTolerance); + painCave.severity = OOPSE_WARNING; simError(); } } @@ -459,8 +461,8 @@ void SimInfo::refreshSim(){ //Most of the compilers will organize the memory of vector continuously setFsimulation( &fInfo, &n_global, &n_atoms, identArray, &n_exclude, excl, &nGlobalExcludes, globalExcludes, molMembershipArray, - &mfact[0], &ngroup, &FglobalGroupMembership[0], &isError); - + &mfact[0], &ngroup, &FglobalGroupMembership[0], &isError); + if( isError ){ sprintf( painCave.errMsg, @@ -504,7 +506,7 @@ void SimInfo::checkCutOffs( void ){ if( rCut > maxCutoff ){ sprintf( painCave.errMsg, - "cutoffRadius is too large for the current periodic box.\n" + "\n\tcutoffRadius is too large for the current periodic box.\n" "\tCurrent Value of cutoffRadius = %G at time %G\n " "\tThis is larger than half of at least one of the\n" "\tperiodic box vectors. Right now, the Box matrix is:\n" @@ -516,14 +518,16 @@ void SimInfo::checkCutOffs( void ){ Hmat[0][0], Hmat[0][1], Hmat[0][2], Hmat[1][0], Hmat[1][1], Hmat[1][2], Hmat[2][0], Hmat[2][1], Hmat[2][2]); + painCave.severity = OOPSE_ERROR; painCave.isFatal = 1; simError(); } } else { // initialize this stuff before using it, OK? sprintf( painCave.errMsg, - "Trying to check cutoffs without a box.\n" + "\n\tTrying to check cutoffs without a box.\n" "\tOOPSE should have better programmers than that.\n" ); + painCave.severity = OOPSE_ERROR; painCave.isFatal = 1; simError(); } @@ -596,6 +600,7 @@ void SimInfo::getFortranGroupArrays(SimInfo* info, #endif for (int i = 0; i < numAtom; i++) FglobalGroupMembership.push_back(globalGroupMembership[i] + 1); + myMols = info->molecules; numMol = info->n_mol;