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, 6 insertions, 6 deletions

diff --git a/misc/benchmarks/various/cspan_bench.c b/misc/benchmarks/various/cspan_bench.c
index 02ae3237..6ca7425d 100644
--- a/misc/benchmarks/various/cspan_bench.c
+++ b/misc/benchmarks/various/cspan_bench.c
@@ -28,8 +28,8 @@ static void MDRanges_setup(intptr_t state)
 
     for (intptr_t s = 0; s < state; ++s)
     {
-        MD3 r_in = cspan_md(Vin, nx, ny, nz);
-        MD3 r_out = cspan_md(Vout, nx, ny, nz);
+        MD3 r_in = cspan_md('C', Vin, nx, ny, nz);
+        MD3 r_out = cspan_md('C', Vout, nx, ny, nz);
 
         r_in = cspan_slice(MD3, &r_in, {lx, hx}, {ly, hy}, {lz, hz});
         r_out = cspan_slice(MD3, &r_out, {lx, hx}, {ly, hy}, {lz, hz});
@@ -64,8 +64,8 @@ static void TraditionalForLoop(intptr_t state)
 
 static void MDRanges_nested_loop(intptr_t state)
 {
-    MD3 r_in = cspan_md(Vin, nx, ny, nz);
-    MD3 r_out = cspan_md(Vout, nx, ny, nz);
+    MD3 r_in = cspan_md('C', Vin, nx, ny, nz);
+    MD3 r_out = cspan_md('C', Vout, nx, ny, nz);
     r_in = cspan_slice(MD3, &r_in, {lx, hx}, {ly, hy}, {lz, hz});
     r_out = cspan_slice(MD3, &r_out, {lx, hx}, {ly, hy}, {lz, hz});
 
@@ -91,8 +91,8 @@ static void MDRanges_nested_loop(intptr_t state)
 
 static void MDRanges_loop_over_joined(intptr_t state)
 {
-    MD3 r_in = cspan_md(Vin, nx, ny, nz);
-    MD3 r_out = cspan_md(Vout, nx, ny, nz);
+    MD3 r_in = cspan_md('C', Vin, nx, ny, nz);
+    MD3 r_out = cspan_md('C', Vout, nx, ny, nz);
     r_in = cspan_slice(MD3, &r_in, {lx, hx}, {ly, hy}, {lz, hz});
     r_out = cspan_slice(MD3, &r_out, {lx, hx}, {ly, hy}, {lz, hz});