lapack-traits 0.4.0

Generic real/complex scalar trait wrappers for BLAS and LAPACK routines
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

use crate::Scalar;
use lapacke::{Layout, sgesv, dgesv, cgesv, zgesv};
use num_complex::Complex32 as c32;
use num_complex::Complex64 as c64;

pub trait Tgesv : Scalar {
    unsafe fn gesv(layout: Layout,
            n: i32,
            nrhs: i32,
            a: &mut [Self],
            lda: i32,
            ipiv: &mut [i32],
            b: &mut [Self],
            ldb: i32,) -> i32;
}

macro_rules! impl_tgesv(
    ($N: ty, $tgesv: path) => (
        impl Tgesv for $N{
            #[inline]
            unsafe fn gesv(
                layout: Layout,
                n: i32,
                nrhs: i32,
                a: &mut [Self],
                lda: i32,
                ipiv: &mut [i32],
                b: &mut [Self],
                ldb: i32,) -> i32
            {
                $tgesv(layout, n, nrhs, a, lda, ipiv, b, ldb)
            }
        }
    )
);

impl_tgesv!(f32, sgesv);
impl_tgesv!(f64, dgesv);
impl_tgesv!(c32, cgesv);
impl_tgesv!(c64, zgesv);