JuliaMath / JuliaMath/GSL.jl

Wrap matrix `copyto!` incorrect

Open
#144 3 comments 0 reactions 0 assignees View on GitHub
Dominant language
Julia
Stars
101
Forks
30
PR merge metrics
No merged PRs in 30d

Description

When using `copyto!` with `wrap_gsl_matrix`, the result looks correct but it is not.
```julia
using GSL
mat = reshape([Float64(10*j+i) for i=1:2 for j=1:4],4,2)
matgsl = GSL.matrix_alloc(4,2)
copyto!(GSL.wrap_gsl_matrix(matgsl), CartesianIndices((4,2)), mat, CartesianIndices((4,2)))
```
Which gives the expected results for `mat` and `GSL.wrap_gsl_matrix(matgsl)`
```
4×2 Matrix{Float64}:
11.0 12.0
21.0 22.0
31.0 32.0
41.0 42.0
```
but e.g. the first column that's actually in the GSL matrix is not correct
```
[GSL.matrix_get(matgsl,i,0) for i=0:3]
4-element Vector{Float64}:
11.0
31.0
12.0
32.0
```

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.