arpack-ng 0.2.2

Save interface to arpack-ng-sys
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

use arpack_ng::*;
#[cfg(feature = "nalgebra")]
use nalgebra::DMatrix;
#[cfg(feature = "ndarray")]
use ndarray::prelude::*;
#[cfg(feature = "nalgebra")]
use num_complex::Complex64;

#[cfg(feature = "ndarray")]
fn main() -> Result<(), Error> {
    let m = Array2::ones((200, 200));
    let (val, vec) = m.eigenvectors(&Which::LargestRealPart, 2, 50, 100)?;
    println!("{:?} {:?}", val.shape(), vec.shape());
    for i in 0..val.len() {
        println!("{} => {:?}", val[i], vec.slice(s![.., i]));
    }
    Ok(())
}

#[cfg(feature = "nalgebra")]
fn main() -> Result<(), Error> {
    let m = DMatrix::from_element(200, 200, Complex64::new(1., 0.));
    let (val, vec) = m.eigenvectors(&Which::LargestRealPart, 2, 50, 100)?;
    println!("{:?} {:?}", val.shape(), vec.shape());
    for i in 0..val.len() {
        println!("{} => {:.1}", val[i], vec.column(i));
    }
    Ok(())
}