Is there a direct solver for sparse A^TAx=A^Tb, with b having shape n,m?

[Richard Roberts]

Hi there,

I have a sparse matrix A and a dense matrix b. I want to solve for the unknown x in the linear equation A^TAx=A^Tb (A^T is the transpose of A). I would prefer to use a direct, rather than iterative method. Does one exist in SciPy's sparse linear algebra library? (all I can find are methods that assume b has the shape n,1).

Thanks for any help!