sparse_matrix 0.6.3

A library to do linear algebra with sparse matrix
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

use std::ops:: {Add, Sub, Mul};
use serde::{Serialize, Deserialize};

#[derive(Debug, PartialEq, Clone, Serialize, Deserialize)]
pub struct Vector {
    pub values : Vec<f64>
}

impl Vector {
    pub fn len(&self) -> usize {
        self.values.len()
    }

    pub fn is_empty(&self) -> bool {
        self.values.is_empty()
    }
    
    pub fn norm(&self) -> f64 {
        let mut sum = 0.;
        for x in self.values.iter() {
            sum += x*x;
        }
        sum.sqrt()
    }

    pub fn null(n : usize) -> Vector {
        Vector { values : vec![0.;n] }
    }
}

impl Add for &Vector {
    type Output = Result<Vector, &'static str>;

    fn add(self, other: Self) -> Self::Output {
        if self.values.len() == other.values.len() {
            let mut result = vec![];
            for i in 0..self.values.len() {
                result.push(self.values[i] + other.values[i])
            }
            Ok(Vector { values: result })
        } else {
            Err("Wrong shapes")
        }
    }
}

impl Sub for &Vector {
    type Output = Result<Vector, &'static str>;

    fn sub(self, other: Self) -> Self::Output {
        if self.values.len() == other.values.len() {
            let mut result = vec![];
            for i in 0..self.values.len() {
                result.push(self.values[i] - other.values[i])
            }
            Ok(Vector { values: result })
        } else {
            Err("Wrong shapes")
        }
    }
}

impl Mul for &Vector {
    type Output = Result<f64, &'static str>;

    fn mul(self, other: Self) -> Self::Output {
        if self.values.len() == other.values.len() {
            let mut result = 0.;
            for i in 0..self.values.len() {
                result += self.values[i] * other.values[i];
            }
            Ok(result)
        } else {
            Err("Wrong shapes")
        }
    }
}

impl Mul<f64> for &Vector {
    type Output = Vector;

    fn mul(self, other: f64) -> Self::Output {
        let mut result = self.clone();
        for i in 0..self.len() {
            result.values[i] *= other;
        }
        result
    }
}