diff --git a/src/gravity/Default/gravity_iact.h b/src/gravity/Default/gravity_iact.h
index b044107b15b72602ff99144569235e0b6c94078c..eca5c2491cbdcf5f0eca01417c8e6b29efc53459 100644
--- a/src/gravity/Default/gravity_iact.h
+++ b/src/gravity/Default/gravity_iact.h
@@ -64,8 +64,6 @@ __attribute__((always_inline)) INLINE static void runner_iact_grav_pp(
 
     kernel_grav_eval(ui, &W);
 
-    if (W < 0) error("W < 0");
-
     /* Get softened gravity */
     fi = mj * hi_inv3 * W * f_lr;
   }
@@ -83,8 +81,6 @@ __attribute__((always_inline)) INLINE static void runner_iact_grav_pp(
 
     kernel_grav_eval(uj, &W);
 
-    if (W < 0) error("W < 0");
-
     /* Get softened gravity */
     fj = mi * hj_inv3 * W * f_lr;
   }
@@ -135,8 +131,6 @@ __attribute__((always_inline)) INLINE static void runner_iact_grav_pp_nonsym(
 
     kernel_grav_eval(ui, &W);
 
-    if (W < 0) error("W < 0");
-
     /* Get softened gravity */
     f = mj * hi_inv3 * W * f_lr;
   }
diff --git a/src/multipole.h b/src/multipole.h
index ff669bdc8f4eb10923b5b860e062888c6dec06f5..b9d49dcf0fc3b605849f7b058aef14843b73517d 100644
--- a/src/multipole.h
+++ b/src/multipole.h
@@ -1626,7 +1626,7 @@ INLINE static void gravity_M2L(struct grav_tensor *l_b,
                   m_a->M_010 * D_021(dx, dy, dz, r_inv) +
                   m_a->M_001 * D_012(dx, dy, dz, r_inv);
 
-    /*  3rd order multipole term (addition to rank 2)*/
+    /*  3rd order multipole term (addition to rank 3)*/
     l_b->F_300 += m_a->M_000 * D_300(dx, dy, dz, r_inv);
     l_b->F_030 += m_a->M_000 * D_030(dx, dy, dz, r_inv);
     l_b->F_003 += m_a->M_000 * D_003(dx, dy, dz, r_inv);
@@ -2151,7 +2151,7 @@ INLINE static void gravity_M2L(struct grav_tensor *l_b,
                   m_a->M_010 * D_soft_021(dx, dy, dz, r, eps_inv) +
                   m_a->M_001 * D_soft_012(dx, dy, dz, r, eps_inv);
 
-    /*  3rd order multipole term (addition to rank 2)*/
+    /*  3rd order multipole term (addition to rank 3)*/
     l_b->F_300 += m_a->M_000 * D_soft_300(dx, dy, dz, r, eps_inv);
     l_b->F_030 += m_a->M_000 * D_soft_030(dx, dy, dz, r, eps_inv);
     l_b->F_003 += m_a->M_000 * D_soft_003(dx, dy, dz, r, eps_inv);