What is the best way (in terms of the accuracy of the result) to solve a linear system?

What is the best way (in terms of the accuracy of the result) to solve a linear system?. Is it using iterative methods (like conjugate gradient...) or decomposition methods (like SVD or LU)?

If you can recommend any solver in C++.

Thank you in advance.