Function gl_order 

pub fn gl_order(n: u32) -> u128

|GL(n,2)| via the orbit–stabilizer (ordered-basis) factorization Π_{i=0}^{n-1}(2ⁿ − 2ⁱ). Each factor 2ⁿ − 2ⁱ is the count of vectors outside the i-dimensional span of the basis chosen so far — the symmetry-broken enumeration of the bases on which GL(n,2) acts simply transitively. Valid up to n = 10 before the u128 product overflows.