Implement test for COO

Author: BenBrock

include/binsparse/binsparse.hpp

@@ -294,5 +294,3 @@ inline auto inspect(std::string fname) {
 }
 
 } // end binsparse
-
-#include <binsparse/c_bindings/bc_read_matrix.hpp>

include/binsparse/matrix_market/matrix_market_inspector.hpp

@@ -14,7 +14,7 @@ namespace binsparse {
 // 2 - number of values in matrix
 // 3 - type of the matrix (real / integer / complex / pattern)
 // 4 - comments
-auto mmread_metadata(std::string file_path) {
+inline auto mmread_metadata(std::string file_path) {
   std::string type;
 
   std::ifstream f;

test/gtest/CMakeLists.txt

@@ -20,6 +20,7 @@ enable_testing()
 add_executable(
   binsparse-tests
   csr_test.cpp
+  coo_test.cpp
 )
 
 target_link_libraries(binsparse-tests binsparse fmt GTest::gtest_main)

test/gtest/coo_test.cpp

@@ -0,0 +1,46 @@
+#include <gtest/gtest.h>
+
+#include <fmt/core.h>
+
+#include <binsparse/binsparse.hpp>
+
+inline std::vector file_paths({"1138_bus/1138_bus.mtx",
+                               "chesapeake/chesapeake.mtx",
+                               "mouse_gene/mouse_gene.mtx"});
+
+TEST(BinsparseReadWrite, COOFormat) {
+  using T = float;
+  using I = std::size_t;
+
+  std::string binsparse_file = "out.bsp.hdf5";
+
+  for (auto&& file_path : file_paths) {
+    auto x = binsparse::__detail::mmread<T, I, binsparse::__detail::coo_matrix_owning<T, I>>(file_path);
+
+    auto&& [num_rows, num_columns] = x.shape();
+    binsparse::coo_matrix<T, I> matrix{x.values().data(), x.rowind().data(), x.colind().data(), num_rows, num_columns, I(x.size())};
+    binsparse::write_coo_matrix(binsparse_file, matrix);
+
+    auto matrix_ = binsparse::read_coo_matrix<T, I>(binsparse_file);
+
+    EXPECT_EQ(matrix.nnz, matrix_.nnz);
+    EXPECT_EQ(matrix.m, matrix_.m);
+    EXPECT_EQ(matrix.n, matrix_.n);
+
+    for (I i = 0; i < matrix.nnz; i++) {
+      EXPECT_EQ(matrix.values[i], matrix_.values[i]);
+    }
+
+    for (I i = 0; i < matrix.nnz; i++) {
+      EXPECT_EQ(matrix.rowind[i], matrix_.rowind[i]);
+    }
+
+    for (I i = 0; i < matrix.nnz; i++) {
+      EXPECT_EQ(matrix.colind[i], matrix_.colind[i]);
+    }
+
+    delete matrix_.values;
+    delete matrix_.rowind;
+    delete matrix_.colind;
+  }
+}