clarabel 0.11.1

Clarabel Conic Interior Point Solver for Rust / Python
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

#![allow(non_snake_case)]

use crate::algebra::{
    DenseMatrix, FloatT, Matrix, MultiplySYMV, ShapedMatrix, Symmetric,
};

impl<T> MultiplySYMV for Symmetric<'_, Matrix<T>>
where
    T: FloatT,
{
    type T = T;
    // implements y = αA*x + βy
    fn symv(&self, x: &[Self::T], y: &mut [Self::T], α: Self::T, β: Self::T) {
        let (m, n) = self.size();
        assert!(m == n);

        // standard BLAS ?symv arguments for computing matrix-vector product
        let uplo = self.uplo.as_blas_char();
        let n = n.try_into().unwrap();
        let a = self.src.data();
        let lda = n;
        let incx = 1;
        let incy = 1;
        T::xsymv(uplo, n, α, a, lda, x, incx, β, y, incy);
    }
}

macro_rules! generate_test_gsymv {
    ($fxx:ty, $test_name:ident) => {
        #[test]
        fn $test_name() {
            #[rustfmt::skip]
            let A = Matrix::<$fxx>::from(&[
                [ 1.,  2.,   4.], 
                [ 0.,  3.,   5.], 
                [ 0.,  0.,   6.],
            ]);

            let x = vec![1., -2., 3.];
            let mut y = vec![-4., -1., 3.];
            A.sym_up().symv(&x, &mut y, 2.0, 3.0);
            assert_eq!(y, [6.0, 19.0, 33.0]);

            #[rustfmt::skip]
            let A = Matrix::<$fxx>::from(&[
                [ 1.,  0.,   0.], 
                [ 2.,  3.,   0.], 
                [ 4.,  5.,   6.],
            ]);

            let x = vec![1., -2., 3.];
            let mut y = vec![-4., -1., 3.];
            A.sym_lo().symv(&x, &mut y, 2.0, 3.0);
            assert_eq!(y, [6.0, 19.0, 33.0]);
        } 
    };
}

generate_test_gsymv!(f32, test_gsymv_f32);
generate_test_gsymv!(f64, test_gsymv_f64);