test smatrix

acts-project · Nov 29, 2024 · 2a58e77 · 2a58e77
1 parent 4ba6777
commit 2a58e77
Showing 1 changed file with 13 additions and 1 deletion.
diff --git a/storage/smatrix/include/algebra/storage/impl/smatrix_getter.hpp b/storage/smatrix/include/algebra/storage/impl/smatrix_getter.hpp
@@ -137,7 +137,19 @@ struct block_getter {
       const ROOT::Math::SMatrix<scalar_t, ROWS, COLS> &m, unsigned int row,
       unsigned int col) const {
 
-    return m.template SubCol<ROOT::Math::SVector<scalar_t, SIZE>>(col, row);
+    // TODO: SMatrix bug?
+    // return m.template SubCol<ROOT::Math::SVector<scalar_t, SIZE>>(col, row);
+
+    assert(col < COLS);
+    assert(row + SIZE <= ROWS);
+
+    ROOT::Math::SVector<scalar_t, SIZE> ret;
+
+    for (std::size_t irow = row; irow < row + SIZE; ++irow) {
+      ret[irow - row] = m[col][irow];
+    }
+
+    return ret;
   }
 };  // struct block_getter