ogdoad 1.0.5

Clifford algebras and quadratic forms over exact, finite, local, transfinite, and game-adjacent scalar backends.
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

//! Row-reduction kernels specialized to the finite nim field.

use crate::scalar::{nim_add, nim_inv, nim_mul};

/// Rank of a rectangular matrix over `F_{2^128}`, with entries stored as raw
/// `u128` nimbers.
pub(crate) fn nim_rank(mut rows: Vec<Vec<u128>>) -> usize {
    let nrows = rows.len();
    if nrows == 0 {
        return 0;
    }
    let ncols = rows[0].len();
    let mut pr = 0usize;
    for col in 0..ncols {
        let Some(p) = (pr..nrows).find(|&r| rows[r][col] != 0) else {
            continue;
        };
        rows.swap(pr, p);
        let inv = nim_inv(rows[pr][col]).expect("nonzero nimber is invertible");
        for c in col..ncols {
            rows[pr][c] = nim_mul(rows[pr][c], inv);
        }
        for r in 0..nrows {
            if r != pr && rows[r][col] != 0 {
                let f = rows[r][col];
                for c in col..ncols {
                    rows[r][c] = nim_add(rows[r][c], nim_mul(f, rows[pr][c]));
                }
            }
        }
        pr += 1;
        if pr == nrows {
            break;
        }
    }
    pr
}