arthroprod 0.3.0

Clifford Algebra based computation for the theory of Absolute Relativity.
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

#[macro_use]
extern crate arthroprod;

use arthroprod::algebra::{full, Form, MultiVector, Sign, ALLOWED_ALPHA_FORMS, AR};
use arthroprod::prelude::*;

// negate everything but the bivectors
fn double_dagger<T: AR>(m: &T) -> MultiVector {
    MultiVector::from_terms(
        m.as_terms()
            .iter()
            .map(|t| match t.alpha().form() {
                Form::Bivector(_, _) => t.clone(),
                _ => -t.clone(),
            })
            .collect(),
    )
}

fn main() {
    let m = mvec![term!(0 1 2 3), term!(2 3), term!(3 1), term!(1 2)];
    let mut mres: MultiVector = full(&m, &double_dagger(&m));
    mres.simplify();

    println!("m = {}", m);
    println!("m ^ m_ddagger = {}\n", mres);

    println!("Dmu G = {}", Dmu().left_apply(G()));
    let res = Dmu().left_apply(G());

    for form in ALLOWED_ALPHA_FORMS.iter() {
        if let Some(terms) = res.get(form) {
            let signs = terms
                .iter()
                .map(|t| if t.sign() == Sign::Pos { "+" } else { "-" })
                .collect::<Vec<&str>>()
                .join("");
            println!("{} {}", signs, terms[0].form());
        }
    }
}