JuliaSIMD / JuliaSIMD/VectorizedRNG.jl

Support `view`?

Open
#16 2 comments 0 reactions 0 assignees View on GitHub
Dominant language
Julia
Stars
36
Forks
7
PR merge metrics
No merged PRs in 30d

Description

It seems the reordering of matrix leads to wrong result.

```julia
julia> using VectorizedRNG

julia> m = zeros(Float32, 4,4)
4×4 Matrix{Float32}:
0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0

julia> n = view(m, 2:3,2:3)
2×2 view(::Matrix{Float32}, 2:3, 2:3) with eltype Float32:
0.0 0.0
0.0 0.0

julia> rand!(local_rng(), n)
2×2 view(::Matrix{Float32}, 2:3, 2:3) with eltype Float32:
0.383252 0.0
0.110364 0.0

julia> m
4×4 Matrix{Float32}:
0.0 0.0 0.943732 0.0
0.0 0.383252 0.0 0.0
0.0 0.110364 0.0 0.0
0.0 0.603976 0.0 0.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.