Migrate to lapacke library

Also bumped version to 0.8

Author: alexbool

.gitignore

@@ -10,3 +10,7 @@ Cargo.lock
 
 # cargo fmt
 *.bk
+
+# IntelliJ Rust
+.idea
+*.iml

Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name    = "ndarray-linalg"
-version = "0.7.1"
+version = "0.8.0"
 authors = ["Toshiki Teramura <[email protected]>"]
 
 description   = "Linear algebra package for rust-ndarray using LAPACK"
@@ -11,40 +11,20 @@ license       = "MIT"
 readme        = "README.md"
 categories    = ["algorithms", "science"]
 
-[features]
-default  = ["openblas-static"]
-openblas-shared = ["lapack/openblas"]
-openblas-static = ["lapack/openblas", "openblas-src/static"]
-openblas-system = ["lapack/openblas", "openblas-src/system"]
-netlib-shared   = ["lapack/netlib"]
-netlib-static   = ["lapack/netlib", "netlib-src/static"]
-netlib-system   = ["lapack/netlib", "netlib-src/system"]
-intel-mkl       = ["intel-mkl-src"]
-
 [dependencies]
 rand = "0.3"
 derive-new = "0.5"
 procedurals = "0.2"
 num-traits  = "0.1"
-num-complex = "0.1"
-lapack = { version = "0.13", default-features = false }
-lapack-sys = { version = "0.11", default-features = false }
+num-complex = { version = "0.1", default-features = false }
+lapacke = "0.1.4"
+lapack-src = { git = "https://github.com/alexbool/lapack-src", rev = "da33354", optional = true }
 
 [dependencies.ndarray]
 version = "0.10"
 default-features = false
 
-[dependencies.openblas-src]
-version = "0.5.3"
-default-features = false
-optional = true
-
-[dependencies.netlib-src]
-version = "0.7.0"
-default-features = false
-optional = true
-
-[dependencies.intel-mkl-src]
-version = "0.2.5"
-default-features = false
-optional = true
+[features]
+openblas = ["lapack-src", "lapack-src/openblas"]
+netlib = ["lapack-src", "lapack-src/netlib"]
+intel-mkl = ["lapack-src", "lapack-src/intel-mkl"]

README.md

@@ -10,23 +10,38 @@ Dependencies
 -------------
 
 - [bluss/rust-ndarray](https://github.com/bluss/rust-ndarray)
-- [stainless-steel/lapack](https://github.com/stainless-steel/lapack)
+- [blas-lapack-rs/lapacke](https://github.com/blas-lapack-rs/lapacke)
 
 and more (See Cargo.toml).
 
-Feature flags
---------------
+Choosing LAPACKE implementation
+--------------------------------
+
+For the sake of linking flexibility, you must provide LAPACKE implementation (as an `extern crate`) yourself.
+Currently three LAPACKE implementations are supported and tested:
 
 - [OpenBLAS](https://github.com/cmr/openblas-src)
-  - `openblas-static`: use OpenBLAS with static link (default)
-  - `openblas-shared`: use OpenBLAS with shared link
-  - `openblas-system`: use system OpenBLAS (experimental)
 - [Netlib](https://github.com/cmr/netlib-src)
-  - `netlib-static`: use Netlib with static link (default)
-  - `netlib-shared`: use Netlib with shared link
-  - `netlib-system`: use system Netlib (experimental)
 - [Intel MKL](https://github.com/termoshtt/rust-intel-mkl) (non-free license, see the linked page)
-  - `intel-mkl`: use Intel MKL shared link (experimental)
+
+You should link a LAPACKE implementation to a final crate (like binary executable or dylib) only, not to a Rust library.
+Example:
+
+`Cargo.toml`:
+```toml
+[depencdencies]
+ndarray = "0.10"
+ndarray-linalg = "0.8"
+openblas-src = "0.5" # or another backend of your choice
+
+```
+
+`main.rs`:
+```rust
+extern crate ndarray;
+extern crate ndarray_linalg;
+extern crate openblas_src; // or another backend of your choice
+```
 
 Examples
 ---------

examples/eigh.rs

@@ -1,6 +1,8 @@
 
 extern crate ndarray;
 extern crate ndarray_linalg;
+#[cfg(feature = "lapack-src")]
+extern crate lapack_src;
 
 use ndarray::*;
 use ndarray_linalg::*;

examples/solve.rs

@@ -1,6 +1,8 @@
 
 extern crate ndarray;
 extern crate ndarray_linalg;
+#[cfg(feature = "lapack-src")]
+extern crate lapack_src;
 
 use ndarray::*;
 use ndarray_linalg::*;

examples/solveh.rs

@@ -1,6 +1,8 @@
 
 extern crate ndarray;
 extern crate ndarray_linalg;
+#[cfg(feature = "lapack-src")]
+extern crate lapack_src;
 
 use ndarray::*;
 use ndarray_linalg::*;

src/cholesky.rs

@@ -9,7 +9,7 @@
 //! Using the Cholesky decomposition of `A` for various operations, where `A`
 //! is a Hermitian (or real symmetric) positive definite matrix:
 //!
-//! ```
+//! ```no_run
 //! #[macro_use]
 //! extern crate ndarray;
 //! extern crate ndarray_linalg;

src/lapack_traits/cholesky.rs

@@ -1,6 +1,6 @@
 //! Cholesky decomposition
 
-use lapack::c;
+use lapacke;
 
 use error::*;
 use layout::MatrixLayout;
@@ -46,7 +46,7 @@ impl Cholesky_ for $scalar {
 }
 }} // end macro_rules
 
-impl_cholesky!(f64, c::dpotrf, c::dpotri, c::dpotrs);
-impl_cholesky!(f32, c::spotrf, c::spotri, c::spotrs);
-impl_cholesky!(c64, c::zpotrf, c::zpotri, c::zpotrs);
-impl_cholesky!(c32, c::cpotrf, c::cpotri, c::cpotrs);
+impl_cholesky!(f64, lapacke::dpotrf, lapacke::dpotri, lapacke::dpotrs);
+impl_cholesky!(f32, lapacke::spotrf, lapacke::spotri, lapacke::spotrs);
+impl_cholesky!(c64, lapacke::zpotrf, lapacke::zpotri, lapacke::zpotrs);
+impl_cholesky!(c32, lapacke::cpotrf, lapacke::cpotri, lapacke::cpotrs);

src/lapack_traits/eigh.rs

@@ -1,6 +1,6 @@
 //! Eigenvalue decomposition for Hermite matrices
 
-use lapack::c;
+use lapacke;
 use num_traits::Zero;
 
 use error::*;
@@ -27,7 +27,7 @@ impl Eigh_ for $scalar {
 }
 }} // impl_eigh!
 
-impl_eigh!(f64, c::dsyev);
-impl_eigh!(f32, c::ssyev);
-impl_eigh!(c64, c::zheev);
-impl_eigh!(c32, c::cheev);
+impl_eigh!(f64, lapacke::dsyev);
+impl_eigh!(f32, lapacke::ssyev);
+impl_eigh!(c64, lapacke::zheev);
+impl_eigh!(c32, lapacke::cheev);

src/lapack_traits/opnorm.rs

@@ -1,7 +1,7 @@
 //! Operator norms of matrices
 
-use lapack::c;
-use lapack::c::Layout::ColumnMajor as cm;
+use lapacke;
+use lapacke::Layout::ColumnMajor as cm;
 
 use layout::MatrixLayout;
 use types::*;
@@ -24,7 +24,7 @@ impl OperatorNorm_ for $scalar {
 }
 }} // impl_opnorm!
 
-impl_opnorm!(f64, c::dlange);
-impl_opnorm!(f32, c::slange);
-impl_opnorm!(c64, c::zlange);
-impl_opnorm!(c32, c::clange);
+impl_opnorm!(f64, lapacke::dlange);
+impl_opnorm!(f32, lapacke::slange);
+impl_opnorm!(c64, lapacke::zlange);
+impl_opnorm!(c32, lapacke::clange);