Implement HermiteMatrix for RcArray<A, Ix2>

Author: termoshtt

src/hermite.rs

@@ -1,6 +1,6 @@
 //! Define trait for Hermite matrices
 
-use ndarray::{Ix2, Array};
+use ndarray::{Ix2, Array, RcArray};
 use lapack::c::Layout;
 
 use matrix::{Matrix, MFloat};
@@ -71,3 +71,25 @@ impl<A: HMFloat> HermiteMatrix for Array<A, Ix2> {
         Ok(rt * rt)
     }
 }
+
+impl<A: HMFloat> HermiteMatrix for RcArray<A, Ix2> {
+    fn eigh(self) -> Result<(Self::Vector, Self), LinalgError> {
+        // XXX unnecessray clone (should use into_owned())
+        let (e, v) = self.to_owned().eigh()?;
+        Ok((e.into_shared(), v.into_shared()))
+    }
+    fn ssqrt(self) -> Result<Self, LinalgError> {
+        // XXX unnecessray clone (should use into_owned())
+        let s = self.to_owned().ssqrt()?;
+        Ok(s.into_shared())
+    }
+    fn cholesky(self) -> Result<Self, LinalgError> {
+        // XXX unnecessray clone (should use into_owned())
+        let s = self.to_owned().cholesky()?;
+        Ok(s.into_shared())
+    }
+    fn deth(self) -> Result<Self::Scalar, LinalgError> {
+        // XXX unnecessray clone (should use into_owned())
+        self.to_owned().deth()
+    }
+}