flag-algebra 0.4.0

An implementation of Razborov's flag algebras
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

extern crate flag_algebra;

use crate::flags::*;
use flag_algebra::*;

fn identity1<F>(t: Type<F>, n: usize)
where
    F: Flag,
{
    let b = Basis::new(n).with_type(t);
    let x: QFlag<i64, _> = b.random();
    let y = b.random();
    assert_eq!(&(&x - &y) * &(x), (&(&x * &x) - &(&y * &x)));
}

fn commute<F>(t: Type<F>, n1: usize, n2: usize)
where
    F: Flag,
{
    let b1 = Basis::new(n1).with_type(t);
    let b2 = Basis::new(n2).with_type(t);
    let x: QFlag<i64, _> = b1.random();
    let y: QFlag<i64, _> = b2.random();
    assert_eq!(&x * &y, &y * &x)
}

fn assoc<F>(t: Type<F>, n1: usize, n2: usize, n3: usize)
where
    F: Flag,
{
    let b1 = Basis::new(n1).with_type(t);
    let b2 = Basis::new(n2).with_type(t);
    let b3 = Basis::new(n3).with_type(t);
    let x: QFlag<i64, _> = b1.random();
    let y: QFlag<i64, _> = b2.random();
    let z: QFlag<i64, _> = b3.random();
    assert_eq!(&(&x * &y) * &z, &x * &(&y * &z))
}

fn transitive<F>(t: Type<F>, n1: usize, n2: usize, n3: usize)
where
    F: Flag,
{
    let b1 = Basis::new(n1).with_type(t);
    let b2 = Basis::new(n2).with_type(t);
    let b3 = Basis::new(n3).with_type(t);
    let x: QFlag<i64, _> = b1.random();
    assert_eq!(x.expand(b2).expand(b3), x.expand(b3))
}

#[test]
pub fn identities() {
    let t1 = Type::new(1, 0);
    identity1::<Graph>(Type::empty(), 2);
    identity1::<Graph>(t1, 2);
}

#[test]
pub fn mul_commutativity() {
    commute::<Graph>(Type::new(1, 0), 3, 2);
    commute::<Graph>(Type::empty(), 2, 3);
    commute::<OrientedGraph>(Type::new(1, 0), 3, 2);
}

#[test]
pub fn mul_associativity() {
    assoc::<Graph>(Type::new(1, 0), 2, 2, 2);
    assoc::<Graph>(Type::empty(), 2, 1, 2);
    assoc::<OrientedGraph>(Type::new(1, 0), 2, 1, 2);
}

#[test]
pub fn untype_transitivity() {
    transitive::<Graph>(Type::new(1, 0), 3, 4, 5);
    transitive::<Graph>(Type::empty(), 2, 3, 5);
    transitive::<OrientedGraph>(Type::new(1, 0), 3, 4, 5);
}