NAG/IMSL compared to CUDA/CUBLAS

Do you know how these libraries may compare?
I need to do some matrix/statistics calculation and we are evaluating to use NAG libraries or to go ahead with CUDA/CUBLAS custom functions
Do you know if there any standard well know library such as NAG that will be ported to GPU?
Thanks in advance