move solve_general_linear to matrices and name gaussj_solve

This method is now used by default in solve instead of inversion
methods.