8 |
|
|
9 |
|
namespace oopse { |
10 |
|
|
11 |
– |
static IntegratorBuilder<NPTf>* NPTfCreator = new IntegratorBuilder<NPTf>("NPTf"); |
12 |
– |
|
11 |
|
// Basic non-isotropic thermostating and barostating via the Melchionna |
12 |
|
// modification of the Hoover algorithm: |
13 |
|
// |
248 |
|
void NPTf::loadEta() { |
249 |
|
eta= currentSnapshot_->getEta(); |
250 |
|
|
251 |
< |
if (!eta.isDiagonal()) { |
252 |
< |
sprintf( painCave.errMsg, |
253 |
< |
"NPTi error: the diagonal elements of are eta matrix is not same or etaMat is not a diagonal matrix"); |
254 |
< |
painCave.isFatal = 1; |
255 |
< |
simError(); |
256 |
< |
} |
251 |
> |
//if (!eta.isDiagonal()) { |
252 |
> |
// sprintf( painCave.errMsg, |
253 |
> |
// "NPTf error: the diagonal elements of eta matrix are not the same or etaMat is not a diagonal matrix"); |
254 |
> |
// painCave.isFatal = 1; |
255 |
> |
// simError(); |
256 |
> |
//} |
257 |
|
} |
258 |
|
|
259 |
|
void NPTf::saveEta() { |