}
/* Avoid integer overflows */
-static float comm_vol(int npme, int a, int b, int c)
+static float comm_pme_cost_vol(int npme, int a, int b, int c)
{
float comm_vol;
temp *= ir->nky;
temp *= ir->nkz;
temp /= nk;
- comm += temp;
+ comm_pme += temp;
/* Old line comm_pme += npme[i]*overlap*ir->nkx*ir->nky*ir->nkz/nk; */
}
}
* are similar and therefore these formulas also prefer load balance
* in the FFT and pme_solve calculation.
*/
- comm_pme += comm_vol(npme[YY], ir->nky, ir->nkz, ir->nkx);
- comm_pme += comm_vol(npme[XX], ir->nkx, ir->nky, ir->nkz);
+ comm_pme += comm_pme_cost_vol(npme[YY], ir->nky, ir->nkz, ir->nkx);
+ comm_pme += comm_pme_cost_vol(npme[XX], ir->nkx, ir->nky, ir->nkz);
/* Add cost of pbc_dx for bondeds */
cost_pbcdx = 0;