Correct output for cubical complex.

1 files changed, 3 insertions, 2 deletions

diff --git a/phstuff/diphawrapper.py b/phstuff/diphawrapper.py
index a40aba1..02f577c 100644
--- a/phstuff/diphawrapper.py
+++ b/phstuff/diphawrapper.py
@@ -221,9 +221,10 @@ def save_cubical(fname, array):
         for n in array.shape:
             f.write(struct.pack("<q", n))
 
-        for w in array.flat:
+        flat = array.flattened(order="F")
+        for w in flat:
             f.write(struct.pack("<d", float(w)))
-
+            
 def save_simplicial(fname, complex):
     """*EXPERIMENTAL*. Save a filtered simplicial complex.