rssn 0.2.9

A comprehensive scientific computing library for Rust, aiming for feature parity with NumPy and SymPy.
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

use crate::ffi_apis::common::BincodeBuffer;
use crate::ffi_apis::common::from_bincode_buffer;
use crate::ffi_apis::common::to_bincode_buffer;
use crate::symbolic::core::Expr;
use crate::symbolic::matrix::add_matrices;
use crate::symbolic::matrix::determinant;
use crate::symbolic::matrix::inverse_matrix;
use crate::symbolic::matrix::mul_matrices;
use crate::symbolic::matrix::solve_linear_system;
use crate::symbolic::matrix::transpose_matrix;

/// Performs matrix addition.
/// Takes two bincode-serialized `Expr` representing matrices,
/// and returns a bincode-serialized `Expr` representing their sum.
#[unsafe(no_mangle)]
pub extern "C" fn rssn_bincode_matrix_add(
    m1_buf: BincodeBuffer,
    m2_buf: BincodeBuffer,
) -> BincodeBuffer {
    let m1: Option<Expr> = from_bincode_buffer(&m1_buf);

    let m2: Option<Expr> = from_bincode_buffer(&m2_buf);

    match (m1, m2) {
        | (Some(matrix1), Some(matrix2)) => {
            let result = add_matrices(&matrix1, &matrix2);

            to_bincode_buffer(&result)
        },
        | _ => BincodeBuffer::empty(),
    }
}

/// Performs matrix multiplication.
/// Takes two bincode-serialized `Expr` representing matrices,
/// and returns a bincode-serialized `Expr` representing their product.
#[unsafe(no_mangle)]
pub extern "C" fn rssn_bincode_matrix_mul(
    m1_buf: BincodeBuffer,
    m2_buf: BincodeBuffer,
) -> BincodeBuffer {
    let m1: Option<Expr> = from_bincode_buffer(&m1_buf);

    let m2: Option<Expr> = from_bincode_buffer(&m2_buf);

    match (m1, m2) {
        | (Some(matrix1), Some(matrix2)) => {
            let result = mul_matrices(&matrix1, &matrix2);

            to_bincode_buffer(&result)
        },
        | _ => BincodeBuffer::empty(),
    }
}

/// Performs matrix transposition.
/// Takes a bincode-serialized `Expr` representing a matrix,
/// and returns a bincode-serialized `Expr` representing its transpose.
#[unsafe(no_mangle)]
pub extern "C" fn rssn_bincode_matrix_transpose(matrix_buf: BincodeBuffer) -> BincodeBuffer {
    let matrix: Option<Expr> = from_bincode_buffer(&matrix_buf);

    if let Some(m) = matrix {
        let result = transpose_matrix(&m);

        to_bincode_buffer(&result)
    } else {
        BincodeBuffer::empty()
    }
}

/// Computes the determinant of a matrix.
/// Takes a bincode-serialized `Expr` representing a matrix,
/// and returns a bincode-serialized `Expr` representing its determinant.
#[unsafe(no_mangle)]
pub extern "C" fn rssn_bincode_matrix_determinant(matrix_buf: BincodeBuffer) -> BincodeBuffer {
    let matrix: Option<Expr> = from_bincode_buffer(&matrix_buf);

    if let Some(m) = matrix {
        let result = determinant(&m);

        to_bincode_buffer(&result)
    } else {
        BincodeBuffer::empty()
    }
}

/// Computes the inverse of a matrix.
/// Takes a bincode-serialized `Expr` representing a matrix,
/// and returns a bincode-serialized `Expr` representing its inverse.
#[unsafe(no_mangle)]
pub extern "C" fn rssn_bincode_matrix_inverse(matrix_buf: BincodeBuffer) -> BincodeBuffer {
    let matrix: Option<Expr> = from_bincode_buffer(&matrix_buf);

    if let Some(m) = matrix {
        let result = inverse_matrix(&m);

        to_bincode_buffer(&result)
    } else {
        BincodeBuffer::empty()
    }
}

/// Solves a linear system of equations AX = B.
/// Takes two bincode-serialized `Expr` representing matrix A and vector B,
/// and returns a bincode-serialized `Expr` representing the solution vector X.
#[unsafe(no_mangle)]
pub extern "C" fn rssn_bincode_matrix_solve_linear_system(
    a_buf: BincodeBuffer,
    b_buf: BincodeBuffer,
) -> BincodeBuffer {
    let a: Option<Expr> = from_bincode_buffer(&a_buf);

    let b: Option<Expr> = from_bincode_buffer(&b_buf);

    match (a, b) {
        | (Some(matrix_a), Some(vector_b)) => {
            match solve_linear_system(&matrix_a, &vector_b) {
                | Ok(result) => to_bincode_buffer(&result),
                | Err(_) => BincodeBuffer::empty(),
            }
        },
        | _ => BincodeBuffer::empty(),
    }
}