Update.

diff --git a/grids.py b/grids.py
index 3343066..52db2b5 100755 (executable)
--- a/grids.py
+++ b/grids.py
@@ -953,7 +953,7 @@ class Grids(problem.Problem):
             assert i[q] != ai and j[q] != aj
 
             for Z in [X, f_X]:
-                for k in range(1, nb_rec):
+                for k in range(0, nb_rec):
                     Z[i[k] : i[k] + delta, j[k] : j[k] + delta] = c[k]
                 # Z[i[0] + delta // 2 - 1, j[0] + delta // 2 - 1] = c[0]
                 # Z[i[0] + delta // 2 - 1, j[0] + delta // 2 + 1] = c[0]
@@ -962,8 +962,8 @@ class Grids(problem.Problem):
 
             # f_X[i[0] : i[0] + delta, j[0] : j[0] + delta] = c[q]
 
-            # f_X[i[0] + delta // 2, j[0] + delta // 2] = c[q]
-            f_X[i[0] : i[0] + delta, j[0] : j[0] + delta] = c[q]
+            f_X[i[0] + delta // 2, j[0] + delta // 2] = c[q]
+            # f_X[i[0] : i[0] + delta, j[0] : j[0] + delta] = c[q]
 
             ii, jj = (
                 i[0] + delta // 2 + (i[q] - ai).sign().long(),