spenso 0.6.0

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

use std::ops::Div;

use symbolica::atom::{Atom, AtomView};

use crate::tensors::parametric::{to_param::ToAtom, ConcreteOrParam, ConcreteOrParamRef};

impl<T, U, Out> Div<ConcreteOrParamRef<'_, U>> for ConcreteOrParam<T>
where
    Atom: Div<U, Output = Atom> + for<'a> Div<AtomView<'a>, Output = Atom>,
    T: ToAtom + Div<U, Output = Out>,
{
    type Output = ConcreteOrParam<Out>;
    fn div(self, rhs: ConcreteOrParamRef<'_, U>) -> ConcreteOrParam<Out> {
        match (self, rhs) {
            (ConcreteOrParam::Concrete(c1), ConcreteOrParamRef::Concrete(c2)) => {
                ConcreteOrParam::Concrete(c1 / c2)
            }
            (ConcreteOrParam::Param(a), ConcreteOrParamRef::Concrete(c)) => {
                ConcreteOrParam::Param(a / c)
            }
            (ConcreteOrParam::Param(a1), ConcreteOrParamRef::Param(a2)) => {
                ConcreteOrParam::Param(a1 / a2)
            }
            (mut c, a) => {
                c.to_param();
                c / a
            }
        }
    }
}

// impl<T: Into<Coefficient>> Div<Complex<T>> for Atom {
//     type Output = Atom;
//     fn div(self, rhs: Complex<T>) -> Atom {
//         self / (Atom::num(rhs.re) + Atom::i() * rhs.im)
//     }
// }

// impl<T> Div<&Complex<T>> for Atom
// where
//     for<'a> &'a T: Into<Coefficient>,
// {
//     type Output = Atom;
//     fn div(self, rhs: &Complex<T>) -> Atom {
//         self / (Atom::num(&rhs.re) + Atom::i() * &rhs.im)
//     }
// }

impl<T, U, Out> Div<ConcreteOrParam<T>> for ConcreteOrParam<U>
where
    Atom: Div<T, Output = Atom> + Div<Atom, Output = Atom>,
    U: Div<T, Output = Out> + ToAtom,
{
    type Output = ConcreteOrParam<Out>;
    fn div(self, rhs: ConcreteOrParam<T>) -> Self::Output {
        match (self, rhs) {
            (ConcreteOrParam::Concrete(c1), ConcreteOrParam::Concrete(c2)) => {
                ConcreteOrParam::Concrete(c1 / c2)
            }
            (ConcreteOrParam::Param(a), ConcreteOrParam::Concrete(c)) => {
                ConcreteOrParam::Param(a / c)
            }
            (ConcreteOrParam::Param(a1), ConcreteOrParam::Param(a2)) => {
                ConcreteOrParam::Param(a1 / a2)
            }
            (mut c, a) => {
                c.to_param();
                c / a
            }
        }
    }
}