JuliaSIMD / JuliaSIMD/VectorizedRNG.jl

Support sampling into non-contiguous arrays

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

Description

This is currently silently broken, which is dangerous:
```julia
julia> randn!(local_rng(), view(X, :,2,:))
3×3 view(::Array{Float64, 3}, :, 2, :) with eltype Float64:
-0.601185 0.0 0.0
0.529745 0.0 0.0
-1.30729 0.0 0.0

julia> X
3×3×3 Array{Float64, 3}:
[:, :, 1] =
0.0 -0.601185 -0.662154
0.0 0.529745 1.66186
0.0 -1.30729 0.848472

[:, :, 2] =
0.891294 0.0 0.0
0.928313 0.0 0.0
-1.01618 0.0 0.0

[:, :, 3] =
0.0 0.0 0.0
0.0 0.0 0.0
0.0 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.