Numerical Linear Algebra in arbitrary precision

[Ken]

Brand new Python user and a bit overwhelmed with the variety of
packages available.  Any recommendation for performing numerical
linear algebra (specifically least squares and generalized least
squares using QR or SVD) in arbitrary precision?  I've been looking at
mpmath but can't seem to find much info on built in functions except
for LU decomposition/solve.

Appreciate any comments.

Ken