1 | initial version |
SVD method has very high numerical robustness, but this robustness has its price both in space and in time. If you don't expect a degenerate case you can you much lighter methods that should be able to run on x86. For example:
solve(A.t()*A, A.t()*b, x, DECOMP_CHOLESKY);