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

use crate::{
    traits::{Dot, Outer},
    Matrix, Vector,
};
use core::ops::{Add, Mul};

impl<T, const M: usize, const N: usize> Outer<Vector<T, N>> for Vector<T, M>
where
    T: Mul<Output = T> + Copy,
{
    type Output = Matrix<T, M, N>;
    fn outer(self, other: Vector<T, N>) -> Self::Output {
        Matrix::indices().map(|(i, j)| self[i] * other[j])
    }
}

impl<T, const M: usize, const N: usize> Dot<Vector<T, N>> for Matrix<T, M, N>
where
    T: Mul<Output = T> + Add<Output = T> + Copy,
{
    type Output = Vector<T, M>;
    fn dot(self, vec: Vector<T, N>) -> Self::Output {
        Vector::indices().map(|i| self.row(i).dot(vec))
    }
}

impl<T, const M: usize, const N: usize> Dot<Matrix<T, M, N>> for Vector<T, M>
where
    T: Mul<Output = T> + Add<Output = T> + Copy,
{
    type Output = Vector<T, N>;
    fn dot(self, mat: Matrix<T, M, N>) -> Self::Output {
        Vector::indices().map(|j| self.dot(mat.col(j)))
    }
}

impl<T, const L: usize, const M: usize, const N: usize> Dot<Matrix<T, M, N>> for Matrix<T, L, M>
where
    T: Mul<Output = T> + Add<Output = T> + Copy,
{
    type Output = Matrix<T, L, N>;
    fn dot(self, mat: Matrix<T, M, N>) -> Self::Output {
        Matrix::indices().map(|(i, j)| self.row(i).dot(mat.col(j)))
    }
}