diff --git a/src/DIELECTRIC/pppm_dielectric.cpp b/src/DIELECTRIC/pppm_dielectric.cpp
index e5b0c5abb1..40ebbe4794 100644
--- a/src/DIELECTRIC/pppm_dielectric.cpp
+++ b/src/DIELECTRIC/pppm_dielectric.cpp
@@ -282,7 +282,7 @@ void PPPMDielectric::qsum_qsq(int warning_flag)
   double qsqsume;
 
 #if defined(_OPENMP)
-#pragma omp parallel for default(shared) reduction(+:qsum_local,qsqsum_local)
+#pragma omp parallel for default(shared) reduction(+:qsum_local,qsqsum_local,qsqsume_local)
 #endif
   for (int i = 0; i < nlocal; i++) {
     qsum_local += q[i];
diff --git a/src/DIELECTRIC/pppm_disp_dielectric.cpp b/src/DIELECTRIC/pppm_disp_dielectric.cpp
index 052b617351..e0652195b7 100644
--- a/src/DIELECTRIC/pppm_disp_dielectric.cpp
+++ b/src/DIELECTRIC/pppm_disp_dielectric.cpp
@@ -556,7 +556,7 @@ void PPPMDispDielectric::qsum_qsq(int warning_flag)
   double qsqsume;
 
 #if defined(_OPENMP)
-#pragma omp parallel for default(shared) reduction(+:qsum_local,qsqsum_local)
+#pragma omp parallel for default(shared) reduction(+:qsum_local,qsqsum_local,qsqsume_local)
 #endif
   for (int i = 0; i < nlocal; i++) {
     qsum_local += q[i];