pbc_dx(pbc, fr.x[C], fr.x[D], vDC.as_vec());
pbc_dx(pbc, fr.x[C], fr.x[B], vBC.as_vec());
- for(auto &i : {vBA, vDC, vBC}){
- for(std::size_t j {XX}; j <= ZZ; ++j){
- i[j] *= gmx::c_nm2A;
- }
+ for(std::size_t j {XX}; j <= ZZ; ++j){
+ vBA[j] *= gmx::c_nm2A;
+ vDC[j] *= gmx::c_nm2A;
+ vBC[j] *= gmx::c_nm2A;
}
-
for(std::size_t i{XX}, j{j + 1}, k{i + 2}; i <= ZZ; ++i, ++j, ++k){
if (j > 2){
j -= 3;