fix ownership

Author: montanalow

Cargo.toml

@@ -42,7 +42,7 @@ approx = "0.4"
 
 ndarray = { version = "0.15", features = ["approx"] }
 ndarray-linalg = { version = "0.16", optional = true }
-sprs = { version = "0.11", default-features = false }
+sprs = { version = "=0.11.1", default-features = false }
 
 thiserror = "1.0"
 

algorithms/linfa-kernel/Cargo.toml

@@ -26,7 +26,7 @@ features = ["std", "derive"]
 [dependencies]
 ndarray = "0.15"
 num-traits = "0.2"
-sprs = { version="0.11", default-features = false }
+sprs = { version="0.11.1", default-features = false }
 
 linfa = { version = "0.7.0", path = "../.." }
 linfa-nn = { version = "0.7.0", path = "../linfa-nn" }

algorithms/linfa-kernel/src/inner.rs

@@ -61,7 +61,19 @@ impl<F: Float> Inner for CsMat<F> {
     type Elem = F;
 
     fn dot(&self, rhs: &ArrayView2<F>) -> Array2<F> {
-        self.mul(&rhs.to_owned())
+        let mut result = Array2::zeros((self.rows(), rhs.ncols()));
+
+        // Handle potential sparse matrices
+        for j in 0..rhs.ncols() {
+            let col = rhs.column(j);
+            let col_result = self.mul(&col.to_owned());
+            // Copy result into appropriate column of output
+            for i in 0..self.rows() {
+                result[[i, j]] = col_result[i];
+            }
+        }
+
+        result
     }
     fn sum(&self) -> Array1<F> {
         let mut sum = Array1::zeros(self.cols());
@@ -106,7 +118,19 @@ impl<'a, F: Float> Inner for CsMatView<'a, F> {
     type Elem = F;
 
     fn dot(&self, rhs: &ArrayView2<F>) -> Array2<F> {
-        self.mul(&rhs.to_owned())
+        let mut result = Array2::zeros((self.rows(), rhs.ncols()));
+
+        // Handle potential sparse matrices
+        for j in 0..rhs.ncols() {
+            let col = rhs.column(j);
+            let col_result = self.mul(&col.to_owned());
+            // Copy result into appropriate column of output
+            for i in 0..self.rows() {
+                result[[i, j]] = col_result[i];
+            }
+        }
+
+        result
     }
     fn sum(&self) -> Array1<F> {
         let mut sum = Array1::zeros(self.cols());

algorithms/linfa-preprocessing/Cargo.toml

@@ -29,7 +29,7 @@ ndarray-rand = { version = "0.14" }
 unicode-normalization = "0.1.8"
 regex = "1.4.5"
 encoding = "0.2"
-sprs = { version = "0.11.0", default-features = false }
+sprs = { version = "0.11.1", default-features = false }
 
 serde_regex = { version = "1.1", optional = true }