| 111 |
|
} |
| 112 |
|
} |
| 113 |
|
|
| 114 |
< |
if (nConstrained){ |
| 115 |
< |
constrainA(); |
| 116 |
< |
} |
| 114 |
> |
consFramework->doConstrainA(); |
| 115 |
|
|
| 116 |
|
// Finally, evolve chi a half step (just like a velocity) using |
| 117 |
|
// temperature at time t, not time t+dt/2 |
| 191 |
|
} |
| 192 |
|
} |
| 193 |
|
|
| 194 |
< |
if (nConstrained){ |
| 197 |
< |
constrainB(); |
| 198 |
< |
} |
| 194 |
> |
consFramework->doConstrainB(); |
| 195 |
|
|
| 196 |
|
if (fabs(prevChi - chi) <= chiTolerance) break; |
| 197 |
|
} |
| 216 |
|
|
| 217 |
|
if (!have_target_temp) { |
| 218 |
|
sprintf( painCave.errMsg, |
| 219 |
< |
"NVT error: You can't use the NVT integrator without a targetTemp!\n" |
| 219 |
> |
"You can't use the NVT integrator without a targetTemp!\n" |
| 220 |
|
); |
| 221 |
|
painCave.isFatal = 1; |
| 222 |
+ |
painCave.severity = OOPSE_ERROR; |
| 223 |
|
simError(); |
| 224 |
|
return -1; |
| 225 |
|
} |
| 228 |
|
|
| 229 |
|
if (!have_tau_thermostat) { |
| 230 |
|
sprintf( painCave.errMsg, |
| 231 |
< |
"NVT error: If you use the constant temperature\n" |
| 232 |
< |
" integrator, you must set tauThermostat.\n"); |
| 231 |
> |
"If you use the constant temperature\n" |
| 232 |
> |
"\tintegrator, you must set tauThermostat.\n"); |
| 233 |
> |
painCave.severity = OOPSE_ERROR; |
| 234 |
|
painCave.isFatal = 1; |
| 235 |
|
simError(); |
| 236 |
|
return -1; |
| 238 |
|
|
| 239 |
|
if (!have_chi_tolerance) { |
| 240 |
|
sprintf( painCave.errMsg, |
| 241 |
< |
"NVT warning: setting chi tolerance to 1e-6\n"); |
| 241 |
> |
"In NVT integrator: setting chi tolerance to 1e-6\n"); |
| 242 |
|
chiTolerance = 1e-6; |
| 243 |
|
have_chi_tolerance = 1; |
| 244 |
+ |
painCave.severity = OOPSE_INFO; |
| 245 |
|
painCave.isFatal = 0; |
| 246 |
|
simError(); |
| 247 |
|
} |