lapack-inject 0.1.2

LAPACK compatible interface backed by runtime-registered Fortran LAPACK function pointers
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
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

//! Test compatibility with lapack-sys.
//!
//! This test verifies that lapack-inject provides all the symbols
//! that lapack-sys expects.

// Link lapack-inject to provide the symbols
extern crate lapack_inject;

// Now lapack-sys functions should resolve to lapack-inject's exports
use lapack_sys;
use libloading::Library;

fn load_built_lapack_inject() -> Library {
    let profile = if cfg!(debug_assertions) {
        "debug"
    } else {
        "release"
    };
    let lib_name = if cfg!(target_os = "macos") {
        "liblapack_inject.dylib"
    } else if cfg!(target_os = "windows") {
        "lapack_inject.dll"
    } else {
        "liblapack_inject.so"
    };
    let path = std::path::PathBuf::from(env!("CARGO_MANIFEST_DIR"))
        .join("target")
        .join(profile)
        .join(lib_name);
    unsafe { Library::new(&path) }
        .unwrap_or_else(|err| panic!("failed to load {}: {err}", path.display()))
}

/// Test that we can take the address of lapack-sys functions
/// (they resolve to lapack-inject's symbols)
#[test]
fn symbols_are_linkable() {
    // These function pointer assignments verify that the symbols exist
    // and have the correct type signature.

    // GESV
    let _ = lapack_sys::sgesv_ as *const ();
    let _ = lapack_sys::dgesv_ as *const ();
    let _ = lapack_sys::cgesv_ as *const ();
    let _ = lapack_sys::zgesv_ as *const ();

    // GETRF
    let _ = lapack_sys::sgetrf_ as *const ();
    let _ = lapack_sys::dgetrf_ as *const ();
    let _ = lapack_sys::cgetrf_ as *const ();
    let _ = lapack_sys::zgetrf_ as *const ();

    // GETRI
    let _ = lapack_sys::sgetri_ as *const ();
    let _ = lapack_sys::dgetri_ as *const ();
    let _ = lapack_sys::cgetri_ as *const ();
    let _ = lapack_sys::zgetri_ as *const ();

    // POTRF
    let _ = lapack_sys::spotrf_ as *const ();
    let _ = lapack_sys::dpotrf_ as *const ();
    let _ = lapack_sys::cpotrf_ as *const ();
    let _ = lapack_sys::zpotrf_ as *const ();

    // GEQRF
    let _ = lapack_sys::sgeqrf_ as *const ();
    let _ = lapack_sys::dgeqrf_ as *const ();
    let _ = lapack_sys::cgeqrf_ as *const ();
    let _ = lapack_sys::zgeqrf_ as *const ();

    // GESVD
    let _ = lapack_sys::sgesvd_ as *const ();
    let _ = lapack_sys::dgesvd_ as *const ();
    let _ = lapack_sys::cgesvd_ as *const ();
    let _ = lapack_sys::zgesvd_ as *const ();

    // GESDD
    let _ = lapack_sys::sgesdd_ as *const ();
    let _ = lapack_sys::dgesdd_ as *const ();
    let _ = lapack_sys::cgesdd_ as *const ();
    let _ = lapack_sys::zgesdd_ as *const ();

    // GEEV
    let _ = lapack_sys::sgeev_ as *const ();
    let _ = lapack_sys::dgeev_ as *const ();
    let _ = lapack_sys::cgeev_ as *const ();
    let _ = lapack_sys::zgeev_ as *const ();

    // SYEV
    let _ = lapack_sys::ssyev_ as *const ();
    let _ = lapack_sys::dsyev_ as *const ();

    // HEEV
    let _ = lapack_sys::cheev_ as *const ();
    let _ = lapack_sys::zheev_ as *const ();

    // GELS
    let _ = lapack_sys::sgels_ as *const ();
    let _ = lapack_sys::dgels_ as *const ();
    let _ = lapack_sys::cgels_ as *const ();
    let _ = lapack_sys::zgels_ as *const ();

    println!("All tested symbols are linkable!");
}

#[test]
fn tenferro_lapack_symbols_are_exported_by_cdylib() {
    let lib = load_built_lapack_inject();
    for symbol in [
        b"dgesvd_\0".as_slice(),
        b"zgesvd_\0".as_slice(),
        b"dgeqrf_\0".as_slice(),
        b"zgeqrf_\0".as_slice(),
        b"dorgqr_\0".as_slice(),
        b"zungqr_\0".as_slice(),
        b"dtrtrs_\0".as_slice(),
        b"ztrtrs_\0".as_slice(),
        b"dpotrf_\0".as_slice(),
        b"zpotrf_\0".as_slice(),
        b"dgetrf_\0".as_slice(),
        b"zgetrf_\0".as_slice(),
        b"dsyev_\0".as_slice(),
        b"zheev_\0".as_slice(),
        b"dgeev_\0".as_slice(),
        b"zgeev_\0".as_slice(),
        b"dgetc2_\0".as_slice(),
        b"dgesc2_\0".as_slice(),
        b"zgetc2_\0".as_slice(),
        b"zgesc2_\0".as_slice(),
    ] {
        unsafe {
            lib.get::<*const ()>(symbol)
                .unwrap_or_else(|err| panic!("missing {}: {err}", String::from_utf8_lossy(symbol)));
        }
    }
}