Using Polynomials.jl with CUDA.jl

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])

Is there any simple way to solve this?

Open an issue. This should be possible.

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

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