Add Eigh_

Author: termoshtt

src/impl2/eigh.rs

@@ -0,0 +1,31 @@
+
+use lapack::c;
+use num_traits::Zero;
+
+use types::*;
+use error::*;
+use layout::Layout;
+
+use super::{into_result, UPLO};
+
+pub trait Eigh_: AssociatedReal {
+    fn eigh(calc_eigenvec: bool, Layout, UPLO, a: &mut [Self]) -> Result<Vec<Self::Real>>;
+}
+
+macro_rules! impl_eigh {
+    ($scalar:ty, $ev:path) => {
+impl Eigh_ for $scalar {
+    fn eigh(calc_v: bool, l: Layout, uplo: UPLO, mut a: &mut [Self]) -> Result<Vec<Self::Real>> {
+        let (n, _) = l.size();
+        let jobz = if calc_v { b'V' } else { b'N' };
+        let mut w = vec![Self::Real::zero(); n as usize];
+        let info = $ev(l.lapacke_layout(), jobz, uplo as u8, n, &mut a, n, &mut w);
+        into_result(info, w)
+    }
+}
+}} // impl_eigh!
+
+impl_eigh!(f64, c::dsyev);
+impl_eigh!(f32, c::ssyev);
+impl_eigh!(c64, c::zheev);
+impl_eigh!(c32, c::cheev);

src/impl2/mod.rs

@@ -4,17 +4,19 @@ pub mod qr;
 pub mod svd;
 pub mod solve;
 pub mod cholesky;
+pub mod eigh;
 
 pub use self::opnorm::*;
 pub use self::qr::*;
 pub use self::svd::*;
 pub use self::solve::*;
 pub use self::cholesky::*;
+pub use self::eigh::*;
 
 use super::error::*;
 
-pub trait LapackScalar: OperatorNorm_ + QR_ + SVD_ + Solve_ + Cholesky_ {}
-impl<A> LapackScalar for A where A: OperatorNorm_ + QR_ + SVD_ + Solve_ + Cholesky_ {}
+pub trait LapackScalar: OperatorNorm_ + QR_ + SVD_ + Solve_ + Cholesky_ + Eigh_ {}
+impl<A> LapackScalar for A where A: OperatorNorm_ + QR_ + SVD_ + Solve_ + Cholesky_ + Eigh_ {}
 
 pub fn into_result<T>(info: i32, val: T) -> Result<T> {
     if info == 0 {