spenso 0.5.5

A tensor (n-dim array) network, iterating, and contraction (using automatic abstract index matching) library.
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
113
114
115

use super::ScalarMul;
use crate::{
    algebra::complex::{Complex, RealOrComplex},
    algebra::upgrading_arithmetic::FallibleMul,
    iterators::IteratableTensor,
    structure::{HasStructure, TensorStructure},
    tensors::complex::RealOrComplexTensor,
    tensors::data::{DataTensor, DenseTensor, SetTensorData, SparseTensor},
};

impl<T, U, I, Out> ScalarMul<T> for DenseTensor<U, I>
where
    U: FallibleMul<T, Output = Out>,
    I: TensorStructure + Clone,
{
    type Output = DenseTensor<Out, I>;
    fn scalar_mul(&self, rhs: &T) -> Option<Self::Output> {
        let data: Option<Vec<Out>> = self.iter_flat().map(|(_, u)| u.mul_fallible(rhs)).collect();

        data.map(|data| DenseTensor {
            structure: self.structure().clone(),
            data,
        })
    }
}

impl<T, U, I, Out> ScalarMul<T> for SparseTensor<U, I>
where
    U: FallibleMul<T, Output = Out>,
    I: TensorStructure + Clone,
{
    type Output = SparseTensor<Out, I>;
    fn scalar_mul(&self, rhs: &T) -> Option<Self::Output> {
        let mut data = SparseTensor::empty(self.structure().clone(), self.zero.mul_fallible(rhs)?);
        for (indices, u) in self.iter_flat() {
            data.set_flat(indices, u.mul_fallible(rhs)?).unwrap();
        }
        Some(data)
    }
}

impl<T, U, I, Out> ScalarMul<T> for DataTensor<U, I>
where
    DenseTensor<U, I>: ScalarMul<T, Output = DenseTensor<Out, I>>,
    SparseTensor<U, I>: ScalarMul<T, Output = SparseTensor<Out, I>>,
    I: TensorStructure + Clone,
{
    type Output = DataTensor<Out, I>;
    fn scalar_mul(&self, rhs: &T) -> Option<Self::Output> {
        match self {
            DataTensor::Dense(a) => Some(DataTensor::Dense(a.scalar_mul(rhs)?)),
            DataTensor::Sparse(a) => Some(DataTensor::Sparse(a.scalar_mul(rhs)?)),
        }
    }
}

impl<T, U, S, Out> ScalarMul<RealOrComplex<T>> for RealOrComplexTensor<U, S>
where
    DataTensor<U, S>: ScalarMul<T, Output = DataTensor<Out, S>>
        + ScalarMul<Complex<T>, Output = DataTensor<Complex<Out>, S>>,
    DataTensor<Complex<U>, S>: ScalarMul<T, Output = DataTensor<Complex<Out>, S>>
        + ScalarMul<Complex<T>, Output = DataTensor<Complex<Out>, S>>,
    S: TensorStructure + Clone,
{
    type Output = RealOrComplexTensor<Out, S>;

    fn scalar_mul(&self, rhs: &RealOrComplex<T>) -> Option<Self::Output> {
        Some(match (self, rhs) {
            (RealOrComplexTensor::Real(rt), RealOrComplex::Real(rs)) => {
                RealOrComplexTensor::Real(rt.scalar_mul(rs)?)
            }
            (RealOrComplexTensor::Complex(rt), RealOrComplex::Real(rs)) => {
                RealOrComplexTensor::Complex(rt.scalar_mul(rs)?)
            }
            (RealOrComplexTensor::Real(rt), RealOrComplex::Complex(rs)) => {
                RealOrComplexTensor::Complex(rt.scalar_mul(rs)?)
            }
            (RealOrComplexTensor::Complex(rt), RealOrComplex::Complex(rs)) => {
                RealOrComplexTensor::Complex(rt.scalar_mul(rs)?)
            }
        })
    }
}

impl<U, S> ScalarMul<U> for RealOrComplexTensor<U, S>
where
    DataTensor<U, S>: ScalarMul<U, Output = DataTensor<U, S>>,
    DataTensor<Complex<U>, S>: ScalarMul<U, Output = DataTensor<Complex<U>, S>>,
    S: TensorStructure + Clone,
{
    type Output = RealOrComplexTensor<U, S>;

    fn scalar_mul(&self, rhs: &U) -> Option<Self::Output> {
        Some(match self {
            RealOrComplexTensor::Real(rt) => RealOrComplexTensor::Real(rt.scalar_mul(rhs)?),
            RealOrComplexTensor::Complex(rt) => RealOrComplexTensor::Complex(rt.scalar_mul(rhs)?),
        })
    }
}

impl<U, S> ScalarMul<Complex<U>> for RealOrComplexTensor<U, S>
where
    DataTensor<U, S>: ScalarMul<Complex<U>, Output = DataTensor<Complex<U>, S>>,
    DataTensor<Complex<U>, S>: ScalarMul<Complex<U>, Output = DataTensor<Complex<U>, S>>,
    S: TensorStructure + Clone,
{
    type Output = RealOrComplexTensor<U, S>;

    fn scalar_mul(&self, rhs: &Complex<U>) -> Option<Self::Output> {
        Some(match self {
            RealOrComplexTensor::Real(rt) => RealOrComplexTensor::Complex(rt.scalar_mul(rhs)?),
            RealOrComplexTensor::Complex(rt) => RealOrComplexTensor::Complex(rt.scalar_mul(rhs)?),
        })
    }
}