diff --git a/tests/lbpm_disc_pp.cpp b/tests/lbpm_disc_pp.cpp
index 2bb35abe..12ed4f7f 100644
--- a/tests/lbpm_disc_pp.cpp
+++ b/tests/lbpm_disc_pp.cpp
@@ -92,6 +92,7 @@ inline void SignedDistanceDiscPack(double *Distance, int ndiscs, double *List_cx
 		imin = int ((cx-2*r)/hx);
 		imax = int ((cx+2*r)/hx)+2;
 		jmin = int ((cy-2*r)/hy);
+		jmax = int ((cy+2*r)/hy)+2;
 
 		// Obviously we have to do something at the edges
 		if (imin<0)		imin = 0;