matrix_basis seems to work better