/* Update potential energies */
kernel_data->energygrp_elec[ggid] += velecsum;
kernel_data->energygrp_polarization[ggid] += vgbsum;
- dvda[nri] = dvda[nri] + dvdasum*isai0*isai0;
+ dvda[inr] = dvda[inr] + dvdasum*isai0*isai0;
/* Increment number of inner iterations */
inneriter += j_index_end - j_index_start;
fshift[i_shift_offset+YY] += ty;
fshift[i_shift_offset+ZZ] += tz;
- dvda[nri] = dvda[nri] + dvdasum*isai0*isai0;
+ dvda[inr] = dvda[inr] + dvdasum*isai0*isai0;
/* Increment number of inner iterations */
inneriter += j_index_end - j_index_start;