--- pyamg/gallery/stencil.py.orig 2011-08-06 20:26:36.000000000 +0100 | |
+++ pyamg/gallery/stencil.py 2011-08-06 20:27:12.000000000 +0100 | |
@@ -85,7 +85,7 @@ | |
# compute index offset of each dof within the stencil | |
strides = numpy.cumprod( [1] + list(reversed(grid)) )[:-1] | |
- indices = S.nonzero() | |
+ indices = tuple(i.copy() for i in S.nonzero()) | |
for i,s in zip(indices,S.shape): | |
i -= s // 2 | |
for stride,coords in zip(strides, reversed(indices)): |