diff --git a/src/USER-DPD/pair_multi_lucy_rx.cpp b/src/USER-DPD/pair_multi_lucy_rx.cpp
index 0a15796666ed3c76476a15b1697dbfc55990cdbc..43d41147416bc6ffd28fcfea29d48b02cae0bb29 100644
--- a/src/USER-DPD/pair_multi_lucy_rx.cpp
+++ b/src/USER-DPD/pair_multi_lucy_rx.cpp
@@ -923,15 +923,15 @@ void PairMultiLucyRX::computeLocalDensity()
           rho_i += factor;
           if (newton_pair || j < nlocal)
             rho[j] += factor;
-        } else if (rsq < cutsq[itype][jtype]) {
-          const double rcut = sqrt(cutsq[itype][jtype]);
-          const double tmpFactor = 1.0-sqrt(rsq)/rcut;
-          const double tmpFactor4 = tmpFactor*tmpFactor*tmpFactor*tmpFactor;
-          const double factor = (84.0/(5.0*pi*rcut*rcut*rcut))*(1.0+3.0*sqrt(rsq)/(2.0*rcut))*tmpFactor4;
-          rho_i += factor;
-          if (newton_pair || j < nlocal)
-            rho[j] += factor;
         }
+      } else if (rsq < cutsq[itype][jtype]) {
+        const double rcut = sqrt(cutsq[itype][jtype]);
+        const double tmpFactor = 1.0-sqrt(rsq)/rcut;
+        const double tmpFactor4 = tmpFactor*tmpFactor*tmpFactor*tmpFactor;
+        const double factor = (84.0/(5.0*pi*rcut*rcut*rcut))*(1.0+3.0*sqrt(rsq)/(2.0*rcut))*tmpFactor4;
+        rho_i += factor;
+        if (newton_pair || j < nlocal)
+          rho[j] += factor;
       }
     }