Function guff_matrix::reference_matrix_multiply

pub fn reference_matrix_multiply<S: Simd + Copy, G>(
    xform: &mut impl SimdMatrix<S>, 
    input: &mut impl SimdMatrix<S>, 
    output: &mut impl SimdMatrix<S>, 
    field: &G
) where
    G: GaloisField,
    <S as Simd>::E: From<<G as GaloisField>::E> + Copy,
    <G as GaloisField>::E: From<<S as Simd>::E> + Copy, 

Reference matrix multiply. Doesn’t use SIMD at all, but uses generic Simd types to be compatible with actual Simd implementations. Note that this multiply routine does not check the gcd condition so it can be used to multiply matrices of arbitrary sizes.