Add ImplOuter

Author: termoshtt

Cargo.toml

@@ -10,7 +10,8 @@ keywords = ["ndarray", "lapack", "matrix"]
 license = "MIT"
 
 [dependencies]
-lapack = "0.11.1"
+lapack = "0.11"
+blas = "0.15"
 num-traits = "0.1.36"
 
 [dependencies.ndarray]

src/lib.rs

@@ -32,6 +32,7 @@
 //! - [symmetric square root](hermite/trait.HermiteMatrix.html#tymethod.ssqrt)
 //! - [Cholesky factorization](hermite/trait.HermiteMatrix.html#tymethod.cholesky)
 
+extern crate blas;
 extern crate lapack;
 extern crate num_traits;
 #[macro_use(s)]
@@ -45,6 +46,7 @@ pub mod square;
 pub mod hermite;
 pub mod triangular;
 
+pub mod outer;
 pub mod qr;
 pub mod svd;
 pub mod eigh;

src/outer.rs

@@ -0,0 +1,18 @@
+
+use blas::c::{Layout, dger, sger};
+
+pub trait ImplOuter: Sized {
+    fn outer(m: usize, n: usize, a: &[Self], b: &[Self], ab: &mut [Self]);
+}
+
+macro_rules! impl_cholesky {
+    ($scalar:ty, $ger:path) => {
+impl ImplOuter for $scalar {
+    fn outer(m: usize, n: usize, a: &[Self], b: &[Self], mut ab: &mut [Self]) {
+        $ger(Layout::ColumnMajor, m as i32, n as i32, 1.0, a, 1, b, 1, ab, m as i32);
+    }
+}
+}} // end macro_rules
+
+impl_cholesky!(f64, dger);
+impl_cholesky!(f32, sger);