Added support for column-major md cspan.

API change: the create function cspan_md(order, array, d1, d2, ...) has the new first argument order, which must be either 'C' or 'F' (C: row-major or Fortran: column-major). The representation of strides was changed.

1 files changed, 5 insertions, 5 deletions

diff --git a/misc/tests/cspan_test.c b/misc/tests/cspan_test.c
index 9afa63f5..d266fa38 100644
--- a/misc/tests/cspan_test.c
+++ b/misc/tests/cspan_test.c
@@ -8,7 +8,7 @@ using_cspan3(intspan, int);
 
 CTEST(cspan, subdim) {
     int array[] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
-    intspan3 m = cspan_md(array, 2, 2, 3);
+    intspan3 m = cspan_md('C', array, 2, 2, 3);
     
     for (size_t i = 0; i < m.shape[0]; ++i) {
         intspan2 sub_i = cspan_submd3(&m, i);
@@ -23,7 +23,7 @@ CTEST(cspan, subdim) {
 
 CTEST(cspan, slice) {
     int array[] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
-    intspan2 m1 = cspan_md(array, 3, 4);
+    intspan2 m1 = cspan_md('C', array, 3, 4);
 
     size_t sum1 = 0;
     for (size_t i = 0; i < m1.shape[0]; ++i) {
@@ -53,7 +53,7 @@ CTEST(cspan, slice2) {
         c_forrange (i, 10*20*30)
             cstack_int_push(&stack, i);
 
-        intspan3 ms3 = cspan_md(stack.data, 10, 20, 30);
+        intspan3 ms3 = cspan_md('C', stack.data, 10, 20, 30);
         ms3 = cspan_slice(intspan3, &ms3, {1,4}, {3,7}, {20,24});
 
         size_t sum = 0;
@@ -93,7 +93,7 @@ CTEST_SETUP(cspan_cube) {
     c_forrange (i, N)
         cstack_int_push(&_self->stack, i+1);
 
-    intspan3 ms3 = cspan_md(_self->stack.data, CUBE, CUBE, CUBE);
+    intspan3 ms3 = cspan_md('C', _self->stack.data, CUBE, CUBE, CUBE);
 
     c_forrange (i, 0, ms3.shape[0], TSIZE) {
         c_forrange (j, 0, ms3.shape[1], TSIZE) {
@@ -114,7 +114,7 @@ CTEST_TEARDOWN(cspan_cube) {
 
 CTEST_F(cspan_cube, slice3) {
     intptr_t n = cstack_int_size(&_self->stack);
-    //printf("\ntiles: %zi, cells: %zi\n", Tiles_size(&_self->tiles), n);
+    //printf("\ntiles: %d, cells: %d\n", (int)Tiles_size(&_self->tiles), (int)n);
 
     int64_t sum = 0;
     // iterate each 3d tile in sequence