numrs2 0.3.3

A Rust implementation inspired by NumPy for numerical computing (NumRS2)
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

#![allow(deprecated)]
#![allow(clippy::result_large_err)]

use numrs2::prelude::*;
use numrs2::random;

fn main() -> Result<()> {
    println!("NumRS2 v0.1.1 Basic Usage Example");
    println!("==========================================");

    // Create arrays
    let a = Array::from_vec(vec![1.0, 2.0, 3.0, 4.0, 5.0, 6.0]).reshape(&[2, 3]);
    let b = Array::from_vec(vec![7.0, 8.0, 9.0, 10.0, 11.0, 12.0]).reshape(&[2, 3]);

    println!("Array a:");
    println!("{}", a);

    println!("\nArray b:");
    println!("{}", b);

    // Element-wise operations
    let c = a.add(&b);
    println!("\na + b = {}", c);

    let d = a.subtract(&b);
    println!("\na - b = {}", d);

    let e = a.multiply(&b);
    println!("\na * b = {}", e);

    let f = a.divide(&b);
    println!("\na / b = {}", f);

    // Matrix multiplication
    let g = a.matmul(&b.transpose())?;
    println!("\na @ b.T = {}", g);

    // Create arrays with built-in methods
    let zeros = Array::<f64>::zeros(&[2, 3]);
    println!("\nzeros(2, 3) = {}", zeros);

    let ones = Array::<f64>::ones(&[2, 3]);
    println!("\nones(2, 3) = {}", ones);

    let full = Array::<f64>::full(&[2, 3], 5.0);
    println!("\nfull(2, 3, 5.0) = {}", full);

    // Using array methods
    println!("\nArray information:");
    println!("shape: {:?}", a.shape());
    println!("ndim: {}", a.ndim());
    println!("size: {}", a.size());

    // Slicing
    let slice = a.slice(0, 1)?;
    println!("\nSlice a[1, :] = {}", slice);

    // Using map operations
    let squared = a.map(|x| x * x);
    println!("\nSquared elements: {}", squared);

    // Parallel operations
    let parallel_sqrt = a.par_map(|x: f64| x.sqrt());
    println!("\nParallel sqrt: {}", parallel_sqrt);

    // Random arrays
    let rand_arr = random::rand::<f64>(&[2, 3])?;
    println!("\nrandom(2, 3) = {}", rand_arr);

    // Vector operations
    let v1 = Array::from_vec(vec![1.0, 2.0, 3.0]);
    let v2 = Array::from_vec(vec![4.0, 5.0, 6.0]);
    let dot_product = v1.dot(&v2)?;
    println!("\nDot product: {}", dot_product);

    // Matrix condition number
    let well_conditioned = Array::from_vec(vec![3.0, 1.0, 1.0, 2.0]).reshape(&[2, 2]);
    let moderately_ill_conditioned =
        Array::from_vec(vec![1.0, 0.9999, 0.9999, 1.0]).reshape(&[2, 2]);
    let very_ill_conditioned = Array::from_vec(vec![1.0, 0.999999, 0.999999, 1.0]).reshape(&[2, 2]);

    #[cfg(feature = "lapack")]
    {
        println!("\nCondition Number Examples (requires lapack feature):");
        println!("Well-conditioned matrix:");
        println!("{}", well_conditioned);
        println!(
            "Moderately ill-conditioned matrix:\n{}",
            moderately_ill_conditioned
        );
        println!("Very ill-conditioned matrix:\n{}", very_ill_conditioned);
        // Note: condition number functionality requires LAPACK feature
        println!("Condition number calculation requires 'lapack' feature");
        println!("Enable with: cargo run --example basic_usage --features lapack");
    }

    #[cfg(not(feature = "lapack"))]
    {
        println!("\nCondition Number Examples:");
        println!("Condition number functionality requires the 'lapack' feature.");
        println!("Enable with: cargo run --example basic_usage --features lapack");
        println!("For now, showing matrix contents:");
        println!("Well-conditioned matrix:\n{}", well_conditioned);
        println!(
            "Moderately ill-conditioned matrix:\n{}",
            moderately_ill_conditioned
        );
        println!("Very ill-conditioned matrix:\n{}", very_ill_conditioned);
    }

    Ok(())
}