MATLAB Answers

Accelerate eigs with GPU

22 views (last 30 days)
Sven Auschra
Sven Auschra on 29 Jul 2020
Commented: Bruno Luong on 30 Jul 2020
Dear all,
I have implemented a numerical solver (of the Fokker-Planck equation) in MATLAB.
At some point, the algorithm needs to calculate
eigs(L, 1, 0)
of a very large sparse matrix L.
I would like to perform this calculation on a GPU to lower the computational costs. So, I created the array
Lgpu = gpuArray(L);
on the GPU and tried to calculate
eigs(Lgpu, 1, 0)
Unfortunately, I receive the error message: "First argument must be a double matrix or a function."
I am wondering what the cause of this error might be and appreciate any help from you.
Thank you very much.

Answers (2)

Walter Roberson
Walter Roberson on 29 Jul 2020
There is no GPU support for sparse arrays.

Sign in to comment.

Christine Tobler
Christine Tobler on 30 Jul 2020
Edited: Edric Ellis on 30 Jul 2020
The eigs function is not supported on the GPU. There is support for sparse matrices on the GPU, since R2015a: Release notes parallel computing toolbox.
You could pass a function handle to EIGS that would use computation on the GPU, but would need to accept and pass back out vectors on the CPU. I'm not sure how efficient that would be, but it could be worth a try.
Bruno Luong
Bruno Luong on 30 Jul 2020
Yeah I would also expect any speedup using GPU. Mostly EIGS on sparse is mainly an iterative process in double layers. Nothing really leans for GPU computation.

Sign in to comment.





Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!