Global gf2e_t16_init (const gf2e *ff, const word a): : this is a bit of overkill, we could do better
Global mzd_poly_randomize (mzd_poly_t *A): Allow the user to provide a RNG callback.
Global mzd_slice_add_elem (mzd_slice_t *A, const rci_t row, const rci_t col, word elem): This function is considerably slower than it needs to be.
Global mzd_slice_randomize (mzd_slice_t *A): Allow the user to provide a RNG callback.
Global mzd_slice_read_elem (const mzd_slice_t *A, const rci_t row, const rci_t col): This function is considerably slower than it needs to be.
Global mzd_slice_write_elem (mzd_slice_t *A, const rci_t row, const rci_t col, word elem): This function is considerably slower than it needs to be.
Global mzed_echelonize_newton_john (mzed_t *A, int full): we don't really compute the upper triangular form yet, we need to implement _mzed_gauss_submatrix() and a better table creation for that.
Global mzed_randomize (mzed_t *A): Allow the user to provide a RNG callback.