Update blas and lapack

Author: IvanUkhov

Cargo.toml

@@ -21,5 +21,5 @@ enum-error-derive = "0.1"
 num-traits  = "0.1"
 num-complex = "0.1"
 ndarray = { version = "0.9",  default-features = false, features = ["blas"] }
-lapack  = { version = "0.11", default-features = false }
-blas    = { version = "0.15", default-features = false }
+lapack  = { version = "0.13", default-features = false }
+blas    = { version = "0.16", default-features = false }

src/lapack_traits/cholesky.rs

@@ -17,7 +17,7 @@ macro_rules! impl_cholesky {
 impl Cholesky_ for $scalar {
     fn cholesky(l: MatrixLayout, uplo: UPLO, mut a: &mut [Self]) -> Result<()> {
         let (n, _) = l.size();
-        let info = $potrf(l.lapacke_layout(), uplo as u8, n, &mut a, n);
+        let info = unsafe { $potrf(l.lapacke_layout(), uplo as u8, n, &mut a, n) };
         into_result(info, ())
     }
 }

src/lapack_traits/eigh.rs

@@ -21,7 +21,7 @@ impl Eigh_ for $scalar {
         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);
+        let info = unsafe { $ev(l.lapacke_layout(), jobz, uplo as u8, n, &mut a, n, &mut w) };
         into_result(info, w)
     }
 }

src/lapack_traits/opnorm.rs

@@ -32,8 +32,8 @@ macro_rules! impl_opnorm {
 impl OperatorNorm_ for $scalar {
     fn opnorm(t: NormType, l: MatrixLayout, a: &[Self]) -> Self::Real {
         match l {
-            MatrixLayout::F((col, lda)) => $lange(cm, t as u8, lda, col, a, lda),
-            MatrixLayout::C((row, lda)) => $lange(cm, t.transpose() as u8, lda, row, a, lda),
+            MatrixLayout::F((col, lda)) => unsafe { $lange(cm, t as u8, lda, col, a, lda) },
+            MatrixLayout::C((row, lda)) => unsafe { $lange(cm, t.transpose() as u8, lda, row, a, lda) },
         }
     }
 }

src/lapack_traits/qr.rs

@@ -24,14 +24,14 @@ impl QR_ for $scalar {
         let (row, col) = l.size();
         let k = min(row, col);
         let mut tau = vec![Self::zero(); k as usize];
-        let info = $qrf(l.lapacke_layout(), row, col, &mut a, l.lda(), &mut tau);
+        let info = unsafe { $qrf(l.lapacke_layout(), row, col, &mut a, l.lda(), &mut tau) };
         into_result(info, tau)
     }
 
     fn q(l: MatrixLayout, mut a: &mut [Self], tau: &[Self]) -> Result<()> {
         let (row, col) = l.size();
         let k = min(row, col);
-        let info = $gqr(l.lapacke_layout(), row, k, k, &mut a, l.lda(), &tau);
+        let info = unsafe { $gqr(l.lapacke_layout(), row, k, k, &mut a, l.lda(), &tau) };
         into_result(info, ())
     }
 

src/lapack_traits/solve.rs

@@ -25,21 +25,21 @@ impl Solve_ for $scalar {
         let (row, col) = l.size();
         let k = ::std::cmp::min(row, col);
         let mut ipiv = vec![0; k as usize];
-        let info = $getrf(l.lapacke_layout(), row, col, a, l.lda(), &mut ipiv);
+        let info = unsafe { $getrf(l.lapacke_layout(), row, col, a, l.lda(), &mut ipiv) };
         into_result(info, ipiv)
     }
 
     fn inv(l: MatrixLayout, a: &mut [Self], ipiv: &Pivot) -> Result<()> {
         let (n, _) = l.size();
-        let info = $getri(l.lapacke_layout(), n, a, l.lda(), ipiv);
+        let info = unsafe { $getri(l.lapacke_layout(), n, a, l.lda(), ipiv) };
         into_result(info, ())
     }
 
     fn solve(l: MatrixLayout, t: Transpose, a: &[Self], ipiv: &Pivot, b: &mut [Self]) -> Result<()> {
         let (n, _) = l.size();
         let nrhs = 1;
         let ldb = 1;
-        let info = $getrs(l.lapacke_layout(), t as u8, n, nrhs, a, l.lda(), ipiv, b, ldb);
+        let info = unsafe { $getrs(l.lapacke_layout(), t as u8, n, nrhs, a, l.lda(), ipiv, b, ldb) };
         into_result(info, ())
     }
 }

src/lapack_traits/svd.rs

@@ -52,7 +52,7 @@ impl SVD_ for $scalar {
         };
         let mut s = vec![Self::Real::zero(); k as usize];
         let mut superb = vec![Self::Real::zero(); (k-2) as usize];
-        let info = $gesvd(l.lapacke_layout(), ju as u8, jvt as u8, m, n, &mut a, lda, &mut s, &mut u, ldu, &mut vt, ldvt, &mut superb);
+        let info = unsafe { $gesvd(l.lapacke_layout(), ju as u8, jvt as u8, m, n, &mut a, lda, &mut s, &mut u, ldu, &mut vt, ldvt, &mut superb) };
         into_result(info, SVDOutput {
             s: s,
             u: if ldu > 0 { Some(u) } else { None },

src/lapack_traits/triangular.rs

@@ -27,7 +27,7 @@ impl Triangular_ for $scalar {
     fn inv_triangular(l: MatrixLayout, uplo: UPLO, diag: Diag, a: &mut [Self]) -> Result<()> {
         let (n, _) = l.size();
         let lda = l.lda();
-        let info = $trtri(l.lapacke_layout(), uplo as u8, diag as u8, n, a, lda);
+        let info = unsafe { $trtri(l.lapacke_layout(), uplo as u8, diag as u8, n, a, lda) };
         into_result(info, ())
     }
 
@@ -42,7 +42,7 @@ impl Triangular_ for $scalar {
         println!("lda = {}", lda);
         println!("nrhs = {}", nrhs);
         println!("ldb = {}", ldb);
-        let info = $trtrs(al.lapacke_layout(), uplo as u8, Transpose::No as u8, diag as u8, n, nrhs, a, lda, &mut b, ldb);
+        let info = unsafe { $trtrs(al.lapacke_layout(), uplo as u8, Transpose::No as u8, diag as u8, n, nrhs, a, lda, &mut b, ldb) };
         into_result(info, ())
     }
 }