C++ library with an implementation of the Conjugate Gradient method

This is not directly related to CUDA, but Im asking it here because I think you guys may know…I’m looking for a decent implementation of the Conjugate Gradient method in C/C++, CPU code, GPU only in the future perhaps. Google failed at helping me [or perhaps I just don’t know how to use it :rolleyes:]

Thanks.

http://trilinos.sandia.gov/CapabilityWebpa…s/overview.html

Try MtxAlg at www.gh-numerics.com. It’s a sparse matrix library with an optimized conjugate gradient solver that runs in parallel on the CPU and GPU. It also comes with sample source code and sample projects to get you going.