diff --git a/src/GPU/pppm_gpu.cpp b/src/GPU/pppm_gpu.cpp index fc9318e6ef..c026dca445 100644 --- a/src/GPU/pppm_gpu.cpp +++ b/src/GPU/pppm_gpu.cpp @@ -301,7 +301,7 @@ void PPPMGPU::compute(int eflag, int vflag) energy = energy_all; energy *= 0.5*volume; - energy -= g_ewald*qsqsum/1.772453851 + + energy -= g_ewald*qsqsum/MY_PIS + MY_PI2*qsum*qsum / (g_ewald*g_ewald*volume); energy *= qscale; }