JuliaGPU / JuliaGPU/CUDA.jl

Support for matrix exp

Open
#168 1 comment 4 reactions 0 assignees View on GitHub
cuda array enhancement good first issue
Dominant language
Julia
Stars
1.4k
Forks
281
Avg merge
1d 7h
Merged PRs (30d)
30

Description

```julia
using CuArrays
n = 4
A = randn(n, n) |> cu
exp(A)
```
```
ERROR: ReadOnlyMemoryError()
Stacktrace:
[1] exp!(::CuArray{Float32,2}) at /buildworker/worker/package_linux64/build/usr/share/julia/stdlib/v0.7/LinearAlgebra/src/lapack.jl:209
[2] exp(::CuArray{Float32,2}) at /buildworker/worker/package_linux64/build/usr/share/julia/stdlib/v0.7/LinearAlgebra/src/dense.jl:508
[3] top-level scope at none:0
```

I get the same error without a stack trace for `A \ cu(randn(n))`

(This is using CuArrays#vc/0.7 (a55b12a180273c80c34c1d) depending on GPUArrays#master, since CuArrays#master doesn't seem to support indexing `A[i,j]` that `exp(m::AbstractMatrix)` uses)

How much work would it take to support LinearAlgebra functions like these? If reasonable for a beginner like myself, any pointers to where to start?

Thanks!

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.