Something weird is happening with dot
, simpler examples:
julia> using LinearAlgebra
julia> dot(ones(10), view(ones(10), 1:10))
10.0
julia> dot(ones(10), view(ones(10), 10:-1:1))
1.0
julia> [dot(ones(10), view(ones(10), 10:-1:1)) for _ in 1:5, _ in 1:10]
5×10 Array{Float64,2}:
1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
1.0 1.0 1.0 1.0 1.0 1.0 1.16189e165 1.0 1.0 1.0
1.0 1.0 1.0 1.0 1.0 1.0 1.0 2.52324e180 1.0 1.0
1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
1.0 1.0 1.0 1.0 2.52324e180 1.0 1.0 1.0 1.0 1.0
julia> view(ones(10), 10:-1:1) isa StridedArray
true