#include "gromacs/math/vec.h"
#include "gromacs/mdlib/gmx_omp_nthreads.h"
-void calc_mu(int start, int homenr, gmx::ArrayRef<gmx::RVec> x, const real q[], const real qB[],
- int nChargePerturbed,
- dvec mu, dvec mu_B)
+void calc_mu(int start,
+ int homenr,
+ gmx::ArrayRef<gmx::RVec> x,
+ const real q[],
+ const real qB[],
+ int nChargePerturbed,
+ dvec mu,
+ dvec mu_B)
{
int end, m;
double mu_x, mu_y, mu_z;
- end = start + homenr;
+ end = start + homenr;
mu_x = mu_y = mu_z = 0.0;
#pragma omp parallel for reduction(+: mu_x, mu_y, mu_z) schedule(static) \
for (int i = start; i < end; i++)
{
// Trivial OpenMP region that cannot throw
- mu_x += q[i]*x[i][XX];
- mu_y += q[i]*x[i][YY];
- mu_z += q[i]*x[i][ZZ];
+ mu_x += q[i] * x[i][XX];
+ mu_y += q[i] * x[i][YY];
+ mu_z += q[i] * x[i][ZZ];
}
mu[XX] = mu_x;
mu[YY] = mu_y;
for (int i = start; i < end; i++)
{
// Trivial OpenMP region that cannot throw
- mu_x += qB[i]*x[i][XX];
- mu_y += qB[i]*x[i][YY];
- mu_z += qB[i]*x[i][ZZ];
+ mu_x += qB[i] * x[i][XX];
+ mu_y += qB[i] * x[i][YY];
+ mu_z += qB[i] * x[i][ZZ];
}
mu_B[XX] = mu_x * ENM2DEBYE;
mu_B[YY] = mu_y * ENM2DEBYE;
}
}
-gmx_bool read_mu(FILE *fp, rvec mu, real *vol)
+gmx_bool read_mu(FILE* fp, rvec mu, real* vol)
{
/* For backward compatibility */
real mmm[4];
- if (fread(mmm, static_cast<size_t>(4*sizeof(real)), 1, fp) != 1)
+ if (fread(mmm, static_cast<size_t>(4 * sizeof(real)), 1, fp) != 1)
{
return FALSE;
}