blob: f000b93709ac8c3c8e0ab113634c9c594e3217ce (
plain)
1
2
3
4
5
6
7
8
9
10
11
|
--- 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)):
|
