add linear system solver using Gaussian elimination with partial pivoting.
remove unused member from matrix_vector
add conversion operator to value type when matrix_vector is 1x1