Using Polynomials.jl with CUDA.jl

Specific Domains GPU
1

I’m on my very first trial with CUDA.jl, and I’m trying to use it together with Polynomials.jl to calculate the roots of thousands of polynomials with my NVidia GPU. However, when I try to create a polynomial using CuArrays (or CuVector) the REPL freezes:

using CUDA, Polynomials
a = CuVector([1.0f0, 1.0f0, 1.0f0])
Polynomial(a)

Is there any simple way to solve this?

2

Open an issue. This should be possible.

3

Thanks! I opened an issue at CUDA.jl github page.

4

Doesn’t hang here. Polynomials is GPU ‘incompatible’ though, it performs scalar indexing, so it won’t just work out of the box.