| 327 |
|
bool hasAngularMomentumFluxVector = rnemdParams->haveAngularMomentumFluxVector(); |
| 328 |
|
|
| 329 |
|
if (hasAngularMomentumFluxVector) { |
| 330 |
< |
fluxVector_ = rnemdParams->getAngularMomentumFluxVector(); |
| 330 |
> |
std::vector<RealType> amf = rnemdParams->getAngularMomentumFluxVector(); |
| 331 |
> |
if (amf.size() != 3) { |
| 332 |
> |
sprintf(painCave.errMsg, |
| 333 |
> |
"RNEMDRTheta: Incorrect number of parameters specified for angularMomentumFluxVector.\n" |
| 334 |
> |
"\tthere should be 3 parameters, but %lu were specified.\n", |
| 335 |
> |
amf.size()); |
| 336 |
> |
painCave.isFatal = 1; |
| 337 |
> |
simError(); |
| 338 |
> |
} |
| 339 |
> |
fluxVector_.x() = amf[0]; |
| 340 |
> |
fluxVector_.y() = amf[1]; |
| 341 |
> |
fluxVector_.z() = amf[2]; |
| 342 |
|
} else { |
| 343 |
|
|
| 344 |
|
std::string fluxStr = rnemdParams->getFluxType(); |
