Function gl_order_p 

pub fn gl_order_p(n: u32, p: u64) -> u128

|GL(n,p)| over GF(p) via the orbit–stabilizer (ordered-basis) factorization Π_{i=0}^{n-1}(pⁿ − pⁱ). GL(n,p) acts simply transitively on ordered bases of GF(p)ⁿ, so each factor pⁿ − pⁱ counts the vectors outside the i-dimensional span built so far — the same symmetry break as over GF(2), now across the field. The p = 2 case is gf2::gl_order.