extended_matrix 0.9.9

A matrix calculation module.
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

// external imports
use extended_matrix_float::MyFloatTrait;

use crate::FloatTrait;
use crate::Position;
use crate::{BasicOperationsTrait, IntoMatrixTrait};

pub trait VectorTrait: IntoMatrixTrait {
    fn vector_shape_conformity_check(&self) -> Result<(), String> {
        if self.get_shape().0 != 1 && self.get_shape().1 != 1 {
            return Err("Not a vector".to_string());
        }
        Ok(())
    }

    fn norm(&self) -> Result<<Self as BasicOperationsTrait>::Value, String>
    where
        <Self as BasicOperationsTrait>::Value:
            FloatTrait<Output = <Self as BasicOperationsTrait>::Value>,
    {
        self.vector_shape_conformity_check()?;
        Ok(self
            .get_elements()
            .values()
            .fold(
                <<Self as BasicOperationsTrait>::Value>::from(0f32),
                |acc, x| acc + *x * *x,
            )
            .my_sqrt())
    }

    fn dot_product(&self, other: &Self) -> Result<<Self as BasicOperationsTrait>::Value, String>
    where
        <Self as BasicOperationsTrait>::Value:
            FloatTrait<Output = <Self as BasicOperationsTrait>::Value>,
        Self: Clone,
    {
        self.vector_shape_conformity_check()?;
        other.vector_shape_conformity_check()?;
        let mut lhs = self.clone();
        let mut rhs = other.clone();
        if lhs.get_shape().1 == 1 {
            lhs = lhs.transpose();
        }
        if rhs.get_shape().0 == 1 {
            rhs = rhs.transpose();
        }
        Ok(*lhs.multiply(&rhs)?.get_element_value(&Position(0, 0))?)
    }
}