diff --git a/openfpm_numerics b/openfpm_numerics
index cddf4cadcaab4c6c9bd425b67fd144d3c8644143..71fba44b3fbffce4b56adbcea3c715f16da074be 160000
--- a/openfpm_numerics
+++ b/openfpm_numerics
@@ -1 +1 @@
-Subproject commit cddf4cadcaab4c6c9bd425b67fd144d3c8644143
+Subproject commit 71fba44b3fbffce4b56adbcea3c715f16da074be
diff --git a/src/Grid/grid_dist_id.hpp b/src/Grid/grid_dist_id.hpp
index 6007fec2f4d631b9f0a002cee2cce8a9c04bf5ff..1a0e885f9f6d373d33954d10ab7f1b3ddfba2783 100644
--- a/src/Grid/grid_dist_id.hpp
+++ b/src/Grid/grid_dist_id.hpp
@@ -1671,7 +1671,7 @@ public:
 	 * \return The information about the local grids
 	 *
 	 */
-	const openfpm::vector<GBoxes<device_grid::dims>> & getLocalGridsInfo()
+	const openfpm::vector<GBoxes<device_grid::dims>> & getLocalGridsInfo() const
 	{
 #ifdef SE_CLASS2
 		check_valid(this,8);
@@ -2957,6 +2957,18 @@ public:
 		return loc_grid.get(i);
 	}
 
+	/*! \brief Get the i sub-domain grid
+	 *
+	 * \param i sub-domain
+	 *
+	 * \return local grid
+	 *
+	 */
+	const device_grid & get_loc_grid(size_t i) const
+	{
+		return loc_grid.get(i);
+	}
+
 	/*! \brief Get the i sub-domain grid
 	 *
 	 * \param i sub-domain
@@ -2993,7 +3005,7 @@ public:
 	 * \return the number of local grid
 	 *
 	 */
-	size_t getN_loc_grid()
+	size_t getN_loc_grid() const
 	{
 		return loc_grid.size();
 	}