so it turns out that we need int32 indices

Author: junikimm717

include/binsparse/tensor.h

@@ -23,7 +23,11 @@ typedef struct {
 // corresponds to BSP_TENSOR_DENSE
 typedef struct {
   int rank;
+  // pointers_to, while it will only ever point to one bsp_array_t, must be kept
+  // as a pointer (rather than a struct) because there are cases where it MUST
+  // be null.
   bsp_array_t* pointers_to;
+  // indices is supposed to be an array of bsp_array_t's.
   bsp_array_t* indices;
   bsp_level_t* child;
 } bsp_sparse_t;

test/julia/tensor_test.jl

@@ -57,3 +57,35 @@ tensortest(
   joinpath(test_files, "input4.bsp.h5"),
   joinpath(test_files, "output4.bsp.h5")
 )
+
+tensortest(
+  Tensor(
+    Dense(SparseList{Int32}(Element{0.0,Float64,Int32}())),
+    [
+      0 1 0 3;
+      1 0 0 4;
+    ]
+  ),
+  joinpath(test_files, "inputcsr.bsp.h5"),
+  joinpath(test_files, "outputcsr.bsp.h5")
+)
+tensortest(
+  Tensor(
+    Dense(SparseList{Int32}(Element{0.0,Float64,Int32}())),
+    fsprand(10, 10, 0.1)
+  ),
+  joinpath(test_files, "inputcsr2.bsp.h5"),
+  joinpath(test_files, "outputcsr2.bsp.h5")
+)
+
+tensortest(
+  Tensor(
+    Dense(Dense(Element(0.0))),
+    [
+      1 1 2 3;
+      6 1 5 4;
+    ]
+  ),
+  joinpath(test_files, "inputdense.bsp.h5"),
+  joinpath(test_files, "outputdense.bsp.h5")
+)