ENH: linalg: add the function 'solve_circ' for solving a circulant system