ENH: Adding ROOTSMALL to the symmetric matrix in leastSquares

In some geometries the inv of symmetric matrix for leastSquares
seems to become ill-conditioned, as work around we add a
ROOTSMALL value as initilization. This is WIP