BLAS эквивалент функции LAPACK для графических процессоров

Question

В LAPACK есть это function для диагональных

SUBROUTINE DSPGVX(ITYPE, JOBZ, RANGE, UPLO, N, AP, BP, VL, VU, 
$     IL, IU, ABSTOL, M, W, Z, LDZ, WORK, IWORK, 
$     IFAIL, INFO) 

*

Я ищу ее реализацию на GPU. Я пытаюсь найти, была ли эта функция уже реализована в CUDA (или OpenCL), но только нашла CULA, которая не является открытым исходным кодом.

Поэтому и боковая CUBLAS существует, интересно, как я могу узнать, доступен ли эквивалент BLAS или CUBLAS этой подпрограммы.

Answer 1

В BLAS не включены решатели на основе собственных значений, и CUBLAS в этом отношении не отличается. Разработанная UTK библиотека Magma включает в себя несколько программных задач, связанных с ускорением собственных значений GPU. Я не думаю, что xSPGV реализован, но несколько других. В зависимости от характеристик вашей матрицы есть что-то, что вы могли бы использовать.

Answer 2

MAGMA содержит исходный код. CUDA + только CPU (пока?)

«Проект МАГМЫ направлен на развитие плотной линейной библиотеки алгебры аналогичную LAPACK, но и для гетерогенных/гибридных архитектур, начиная с текущей„Multicore + GPU“системой.»

http://icl.cs.utk.edu/magma/index.html

http://web.eecs.utk.edu/~luszczek/pubs/parcocudaopencl.pdf (МАГМА для OpenCL продолжается)

EDIT: Я наткнулся на эту библиотеку, когда я искал OCL вещи ... http://viennacl.sourceforge.net/ Я собираюсь дать ему ближе смотреть. (Вы можете сказать, что это хорошо по стилю сайта)