Abelian sandpile model/Identity: Difference between revisions

Content added Content deleted

Inline

@@ Line 990: / Line 990: @@
     '''The stable final state of a sand-pile.
     '''
-    def p(xs):
-        return all([4 > x for x in xs])
     xs = list(rows)
     w = len(xs)
+    def p(xs):
+        return all([x <= w for x in xs])
     return list(chunksOf(w)(
         until(p)(
@@ Line 1,007: / Line 1,008: @@
     '''
     def go(xs):
-        mbi = findIndex(lambda x: 3 < x)(xs)
+        mbi = findIndex(lambda x: w < x)(xs)
+        surplus = 1 + w
         if None is mbi:
             return xs
@@ Line 1,014: / Line 1,016: @@
             return [
 + k if i in neighbours else (
-                    k - 4 if mbi == i else k
+                    k - surplus if mbi == i else k
                 ) for (i, k) in enumerate(xs)
             ]