pub fn solve_space(
equations: &[ModpEquation],
num_vars: usize,
p: u64,
) -> Option<SolutionSpaceP>Expand description
Solve a GF(p) system for its entire solution space via Gaussian elimination to reduced row
echelon form, returning the symmetry-broken SolutionSpaceP (particular solution + kernel basis),
or None iff the system is inconsistent. Generalizes solve, which returns just one witness, to the
full coset — the GF(p) analogue of crate::gf2::solve_gf2. p must be prime.