NVIDIA has been developing a library of high-performance parallel sparse iterative linear solvers, with an emphasis on multilevel and multigrid methods. In this presentation, I will provide an overview of the library''s design and outline many of the challenges we have faced in balancing numerical behavior against parallel scalability. Our library has been integrated into ANSYS Fluent 14.5, and will be released as a fully supported feature in the upcoming Fluent 15. I will describe the collaboration between ANSYS and NVIDIA, and present benchmarking results across a variety of test problems from CFD and other fields. Finally, I will talk about our future plans and discuss some of the open research problems in the area of algebraic multigrid on massively parallel processors.