implementation with unsigned long matrix
wip-bitmatrix-gauss