opensrdk-linear-algebra 0.7.0

Standard linear algebra library using blas and lapack for OpenSRDK toolchain.
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

use super::SymmetricTridiagonalMatrix;
use crate::matrix::*;
use lapack::dstev;

impl SymmetricTridiagonalMatrix<f64> {
  /// # Eigen decomposition
  /// return (lambda, pt)
  pub fn stev(self) -> Result<(Vec<f64>, Matrix), MatrixError> {
    let (mut d, mut e) = self.eject();
    let n = d.len();
    let mut z = Matrix::new(n, n);
    let mut work = vec![0.0; 1.max(2 * (n - 2))];
    let mut info = 0;

    let n = n as i32;

    unsafe {
      dstev(
        'V' as u8,
        n,
        &mut d,
        &mut e,
        &mut z.elems,
        n,
        &mut work,
        &mut info,
      )
    }

    match info {
      0 => Ok((d, z)),
      _ => Err(MatrixError::LapackRoutineError {
        routine: "dstev".to_owned(),
        info,
      }),
    }
  }
}