Implement TriangularMatrix for RcArray

Author: termoshtt

src/triangular.rs

@@ -1,5 +1,5 @@
 
-use ndarray::{Ix2, Array, NdFloat, ArrayBase, DataMut};
+use ndarray::{Ix2, Array, RcArray, NdFloat, ArrayBase, DataMut};
 
 use matrix::{Matrix, MFloat};
 use square::SquareMatrix;
@@ -32,6 +32,17 @@ impl<A: MFloat> TriangularMatrix for Array<A, Ix2> {
     }
 }
 
+impl<A: MFloat> TriangularMatrix for RcArray<A, Ix2> {
+    fn solve_upper(&self, b: Self::Vector) -> Result<Self::Vector, LinalgError> {
+        let x = self.to_owned().solve_upper(b.to_owned())?;
+        Ok(x.into_shared())
+    }
+    fn solve_lower(&self, b: Self::Vector) -> Result<Self::Vector, LinalgError> {
+        let x = self.to_owned().solve_lower(b.to_owned())?;
+        Ok(x.into_shared())
+    }
+}
+
 pub fn drop_upper<A: NdFloat, S>(mut a: ArrayBase<S, Ix2>) -> ArrayBase<S, Ix2>
     where S: DataMut<Elem = A>
 {