Keep dimensionality of slices

General Usage
5

Sounds like you want the selectdim function:

julia> x = reshape(1:24, (2,3,4))
2×3×4 reshape(::UnitRange{Int64}, 2, 3, 4) with eltype Int64:
[:, :, 1] =
 1  3  5
 2  4  6

[:, :, 2] =
 7   9  11
 8  10  12

[:, :, 3] =
 13  15  17
 14  16  18

[:, :, 4] =
 19  21  23
 20  22  24

julia> y = (selectdim(x, 2, i:i) for i in axes(x,2))
Base.Generator{Base.OneTo{Int64}, var"#3#4"}(var"#3#4"(), Base.OneTo(3))

julia> first(y)
2×1×4 view(reshape(::UnitRange{Int64}, 2, 3, 4), :, 1:1, :) with eltype Int64:
[:, :, 1] =
 1
 2

[:, :, 2] =
 7
 8

[:, :, 3] =
 13
 14

[:, :, 4] =
 19
 20

julia> last(y)
2×1×4 view(reshape(::UnitRange{Int64}, 2, 3, 4), :, 3:3, :) with eltype Int64:
[:, :, 1] =
 5
 6

[:, :, 2] =
 11
 12

[:, :, 3] =
 17
 18

[:, :, 4] =
 23
 24

with the obvious generalization:

julia> sliceiter(x, m) = (selectdim(x, m, i:i) for i in axes(x,m))
3 Likes