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
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137

// Copyright 2014 Michael Yang. All rights reserved.
// Use of this source code is governed by a MIT-style
// license that can be found in the LICENSE file.

//! Bindings for vector functions.

pub mod cblas_s {
    use libc::{c_float, c_int, c_void};

    pub use self::cblas_scopy as copy;
    pub use self::cblas_saxpy as axpy;
    pub use self::cblas_sscal as scal;
    pub use self::cblas_sswap as swap;
    pub use self::cblas_sdsdot as dsdot;
    pub use self::cblas_sdot as dot;
    pub use self::cblas_sasum as asum;
    pub use self::cblas_scasum as casum;
    pub use self::cblas_snrm2 as nrm2;
    pub use self::cblas_scnrm2 as cnrm2;
    pub use self::cblas_srot as rot;
    pub use self::cblas_srotm as rotm;
    pub use self::cblas_srotg as rotg;
    pub use self::cblas_srotmg as rotmg;

    extern {
        pub fn cblas_scopy(n: c_int, x: *const c_float,  inc_x: c_int, y: *mut c_float,  inc_y: c_int);
        pub fn cblas_saxpy(n: c_int, alpha: c_float,       x: *const c_float,  inc_x: c_int, y: *mut c_float,  inc_y: c_int);
        pub fn cblas_sscal(n: c_int, alpha: c_float,       x: *mut c_float,  inc_x: c_int);
        pub fn cblas_sswap(n: c_int, x: *mut c_float,  inc_x: c_int, y: *mut c_float,  inc_y: c_int);
        pub fn cblas_sdsdot(n: c_int, alpha: c_float, x: *const c_float, inc_x: c_int, y: *const c_float, inc_y: c_int) -> c_float;
        pub fn cblas_sdot(n: c_int, x: *const c_float,  inc_x: c_int, y: *const c_float,  inc_y: c_int) -> c_float;
        pub fn cblas_sasum(n: c_int, x: *const c_float,  inc_x: c_int) -> c_float;
        pub fn cblas_scasum(n: c_int, x: *const c_void,   inc_x: c_int) -> c_float;
        pub fn cblas_snrm2(n: c_int, x: *const c_float,  inc_x: c_int) -> c_float;
        pub fn cblas_scnrm2(n: c_int, x: *const c_void,   inc_x: c_int) -> c_float;
        pub fn cblas_srot(n: c_int, x: *mut c_float,  inc_x: c_int, y: *mut c_float,  inc_y: c_int, c: c_float,  s: c_float);
        pub fn cblas_srotm(n: c_int, x: *mut c_float,  inc_x: c_int, y: *mut c_float,  inc_y: c_int, p: *const c_float);
        pub fn cblas_srotg(a: *mut c_float,  b: *mut c_float,  c: *mut c_float,  s: *mut c_float);
        pub fn cblas_srotmg(d1: *mut c_float,  d2: *mut c_float,  b1: *mut c_float,  b2: c_float,  p: *mut c_float);
    }
}

pub mod cblas_d {
    use libc::{c_double, c_float, c_int, c_void};

    pub use self::cblas_dcopy as copy;
    pub use self::cblas_daxpy as axpy;
    pub use self::cblas_dscal as scal;
    pub use self::cblas_dswap as swap;
    pub use self::cblas_dsdot as dsdot;
    pub use self::cblas_ddot as dot;
    pub use self::cblas_dasum as asum;
    pub use self::cblas_dzasum as zasum;
    pub use self::cblas_dnrm2 as nrm2;
    pub use self::cblas_dznrm2 as znrm2;
    pub use self::cblas_drot as rot;
    pub use self::cblas_drotm as rotm;
    pub use self::cblas_drotg as rotg;
    pub use self::cblas_drotmg as rotmg;

    extern {
        pub fn cblas_dcopy(n: c_int, x: *const c_double, inc_x: c_int, y: *mut c_double, inc_y: c_int);
        pub fn cblas_daxpy(n: c_int, alpha: c_double,      x: *const c_double, inc_x: c_int, y: *mut c_double, inc_y: c_int);
        pub fn cblas_dscal (n: c_int, alpha: c_double,      x: *mut c_double, inc_x: c_int);
        pub fn cblas_dswap(n: c_int, x: *mut c_double, inc_x: c_int, y: *mut c_double, inc_y: c_int);
        pub fn cblas_dsdot(n: c_int, x: *const c_float,  inc_x: c_int, y: *const c_float,  inc_y: c_int) -> c_double;
        pub fn cblas_ddot (n: c_int, x: *const c_double, inc_x: c_int, y: *const c_double, inc_y: c_int) -> c_double;
        pub fn cblas_dasum (n: c_int, x: *const c_double, inc_x: c_int) -> c_double;
        pub fn cblas_dzasum(n: c_int, x: *const c_void,   inc_x: c_int) -> c_double;
        pub fn cblas_dnrm2 (n: c_int, x: *const c_double, inc_x: c_int) -> c_double;
        pub fn cblas_dznrm2(n: c_int, x: *const c_void,   inc_x: c_int) -> c_double;
        pub fn cblas_drot(n: c_int, x: *mut c_double, inc_x: c_int, y: *mut c_double, inc_y: c_int, c: c_double, s: c_double);
        pub fn cblas_drotm(n: c_int, x: *mut c_double, inc_x: c_int, y: *mut c_double, inc_y: c_int, p: *const c_double);
        pub fn cblas_drotg(a: *mut c_double, b: *mut c_double, c: *mut c_double, s: *mut c_double);
        pub fn cblas_drotmg(d1: *mut c_double, d2: *mut c_double, b1: *mut c_double, b2: c_double, p: *mut c_double);
    }
}

pub mod cblas_c {
    use libc::{c_float, c_int, c_void};

    pub use self::cblas_ccopy as copy;
    pub use self::cblas_caxpy as axpy;
    pub use self::cblas_cscal as scal;
    pub use self::cblas_csscal as sscal;
    pub use self::cblas_cswap as swap;
    pub use self::cblas_cdotu_sub as dotu_sub;
    pub use self::cblas_cdotc_sub as dotc_sub;

    extern {
        pub fn cblas_ccopy(n: c_int, x: *const c_void,   inc_x: c_int, y: *mut c_void,   inc_y: c_int);
        pub fn cblas_caxpy(n: c_int, alpha: *const c_void, x: *const c_void,   inc_x: c_int, y: *mut c_void,   inc_y: c_int);
        pub fn cblas_cscal (n: c_int, alpha: *const c_void, x: *mut c_void,   inc_x: c_int);
        pub fn cblas_csscal(n: c_int, alpha: c_float,       x: *mut c_void,   inc_x: c_int);
        pub fn cblas_cswap(n: c_int, x: *mut c_void,   inc_x: c_int, y: *mut c_void,   inc_y: c_int);
        pub fn cblas_cdotu_sub(n: c_int, x: *const c_void, inc_x: c_int, y: *const c_void, inc_y: c_int, dotu: *mut c_void);
        pub fn cblas_cdotc_sub(n: c_int, x: *const c_void, inc_x: c_int, y: *const c_void, inc_y: c_int, dotc: *mut c_void);
    }
}

pub mod cblas_z {
    use libc::{c_double, c_int, c_void};

    pub use self::cblas_zcopy as copy;
    pub use self::cblas_zaxpy as axpy;
    pub use self::cblas_zscal as scal;
    pub use self::cblas_zdscal as dscal;
    pub use self::cblas_zswap as swap;
    pub use self::cblas_zdotu_sub as dotu_sub;
    pub use self::cblas_zdotc_sub as dotc_sub;

    extern {
        pub fn cblas_zcopy(n: c_int, x: *const c_void,   inc_x: c_int, y: *mut c_void,   inc_y: c_int);
        pub fn cblas_zaxpy(n: c_int, alpha: *const c_void, x: *const c_void,   inc_x: c_int, y: *mut c_void,   inc_y: c_int);
        pub fn cblas_zscal (n: c_int, alpha: *const c_void, x: *mut c_void,   inc_x: c_int);
        pub fn cblas_zdscal(n: c_int, alpha: c_double,      x: *mut c_void,   inc_x: c_int);
        pub fn cblas_zswap(n: c_int, x: *mut c_void,   inc_x: c_int, y: *mut c_void,   inc_y: c_int);
        pub fn cblas_zdotu_sub(n: c_int, x: *const c_void, inc_x: c_int, y: *const c_void, inc_y: c_int, dotu: *mut c_void);
        pub fn cblas_zdotc_sub(n: c_int, x: *const c_void, inc_x: c_int, y: *const c_void, inc_y: c_int, dotc: *mut c_void);
    }
}

pub mod cblas_i {
    use libc::{c_double, c_float, c_int, c_void, size_t};

    pub use self::cblas_isamax as samax;
    pub use self::cblas_idamax as damax;
    pub use self::cblas_icamax as camax;
    pub use self::cblas_izamax as zamax;

    extern {
        pub fn cblas_isamax(n: c_int, x: *const c_float,  inc_x: c_int) -> size_t;
        pub fn cblas_idamax(n: c_int, x: *const c_double, inc_x: c_int) -> size_t;
        pub fn cblas_icamax(n: c_int, x: *const c_void,   inc_x: c_int) -> size_t;
        pub fn cblas_izamax(n: c_int, x: *const c_void,   inc_x: c_int) -> size_t;
    }
}