void ParrinelloRahmanBarostat::integrateBoxVelocityEquations(Step step)
{
- auto box = statePropagatorData_->constBox();
+ const auto* box = statePropagatorData_->constBox();
parrinellorahman_pcoupl(fplog_,
step,
inputrec_,
void ParrinelloRahmanBarostat::scaleBoxAndPositions()
{
// Propagate the box by the box velocities
- auto box = statePropagatorData_->box();
+ auto* box = statePropagatorData_->box();
for (int i = 0; i < DIM; i++)
{
for (int m = 0; m <= i; m++)
// Scale the coordinates
const int start = 0;
const int homenr = mdAtoms_->mdatoms()->homenr;
- auto x = as_rvec_array(statePropagatorData_->positionsView().paddedArrayRef().data());
+ auto* x = as_rvec_array(statePropagatorData_->positionsView().paddedArrayRef().data());
for (int n = start; n < start + homenr; n++)
{
tmvmul_ur0(mu_, x[n], x[n]);
if (inputrecPreserveShape(inputrec_))
{
- auto box = statePropagatorData_->box();
+ auto* box = statePropagatorData_->box();
const int ndim = inputrec_->epct == PressureCouplingType::SemiIsotropic ? 2 : 3;
do_box_rel(ndim, inputrec_->deform, boxRel_, box, true);
}
// the scaling matrix is calculated, without updating the box velocities.
// The call to parrinellorahman_pcoupl is using nullptr for fplog (since we don't expect any
// output here) and for the pressure (since it might not be calculated yet, and we don't need it).
- auto box = statePropagatorData_->constBox();
+ const auto* box = statePropagatorData_->constBox();
parrinellorahman_pcoupl(nullptr,
initStep_,
inputrec_,
{
wallcycle_start(wcycle_, WallCycleCounter::Update);
- auto xp = as_rvec_array(statePropagatorData_->positionsView().paddedArrayRef().data());
- auto x = as_rvec_array(statePropagatorData_->constPositionsView().paddedArrayRef().data());
- auto v = as_rvec_array(statePropagatorData_->constVelocitiesView().paddedArrayRef().data());
+ auto* xp = as_rvec_array(statePropagatorData_->positionsView().paddedArrayRef().data());
+ const auto* x = as_rvec_array(statePropagatorData_->constPositionsView().paddedArrayRef().data());
+ const auto* v = as_rvec_array(statePropagatorData_->constVelocitiesView().paddedArrayRef().data());
int nth = gmx_omp_nthreads_get(ModuleMultiThread::Update);
int homenr = mdAtoms_->mdatoms()->homenr;
{
wallcycle_start(wcycle_, WallCycleCounter::Update);
- auto v = as_rvec_array(statePropagatorData_->velocitiesView().paddedArrayRef().data());
- auto f = as_rvec_array(statePropagatorData_->constForcesView().force().data());
- auto invMassPerDim = mdAtoms_->mdatoms()->invMassPerDim;
+ auto* v = as_rvec_array(statePropagatorData_->velocitiesView().paddedArrayRef().data());
+ const auto* f = as_rvec_array(statePropagatorData_->constForcesView().force().data());
+ const auto* invMassPerDim = mdAtoms_->mdatoms()->invMassPerDim;
const real lambdaStart = (numStartVelocityScalingValues == NumVelocityScalingValues::Single)
? startVelocityScaling_[0]
{
wallcycle_start(wcycle_, WallCycleCounter::Update);
- auto xp = as_rvec_array(statePropagatorData_->positionsView().paddedArrayRef().data());
- auto x = as_rvec_array(statePropagatorData_->constPositionsView().paddedArrayRef().data());
- auto v = as_rvec_array(statePropagatorData_->velocitiesView().paddedArrayRef().data());
- auto f = as_rvec_array(statePropagatorData_->constForcesView().force().data());
- auto invMassPerDim = mdAtoms_->mdatoms()->invMassPerDim;
+ auto* xp = as_rvec_array(statePropagatorData_->positionsView().paddedArrayRef().data());
+ const auto* x = as_rvec_array(statePropagatorData_->constPositionsView().paddedArrayRef().data());
+ auto* v = as_rvec_array(statePropagatorData_->velocitiesView().paddedArrayRef().data());
+ const auto* f = as_rvec_array(statePropagatorData_->constForcesView().force().data());
+ const auto* invMassPerDim = mdAtoms_->mdatoms()->invMassPerDim;
const real lambdaStart = (numStartVelocityScalingValues == NumVelocityScalingValues::Single)
? startVelocityScaling_[0]
{
wallcycle_start(wcycle_, WallCycleCounter::Update);
- auto xp = as_rvec_array(statePropagatorData_->positionsView().paddedArrayRef().data());
- auto x = as_rvec_array(statePropagatorData_->constPositionsView().paddedArrayRef().data());
- auto v = as_rvec_array(statePropagatorData_->velocitiesView().paddedArrayRef().data());
- auto f = as_rvec_array(statePropagatorData_->constForcesView().force().data());
- auto invMassPerDim = mdAtoms_->mdatoms()->invMassPerDim;
+ auto* xp = as_rvec_array(statePropagatorData_->positionsView().paddedArrayRef().data());
+ const auto* x = as_rvec_array(statePropagatorData_->constPositionsView().paddedArrayRef().data());
+ auto* v = as_rvec_array(statePropagatorData_->velocitiesView().paddedArrayRef().data());
+ const auto* f = as_rvec_array(statePropagatorData_->constForcesView().force().data());
+ const auto* invMassPerDim = mdAtoms_->mdatoms()->invMassPerDim;
const real lambdaStart = (numStartVelocityScalingValues == NumVelocityScalingValues::Single)
? startVelocityScaling_[0]