Fix a case of wrong aliasing for setindex (#164)

* Fix a case of wrong aliasing for setindex

* Update src/diskarray.jl

Co-authored-by: Rafael Schouten <[email protected]>

* Update src/diskarray.jl

Co-authored-by: Rafael Schouten <[email protected]>

* update version

---------

Co-authored-by: Rafael Schouten <[email protected]>

Author: meggart

Project.toml

@@ -1,7 +1,7 @@
 name = "DiskArrays"
 uuid = "3c3547ce-8d99-4f5e-a174-61eb10b00ae3"
 authors = ["Fabian Gans <[email protected]>"]
-version = "0.4.1"
+version = "0.4.2"
 
 [deps]
 LRUCache = "8ac3fa9e-de4c-5943-b1dc-09c6b5f20637"

src/diskarray.jl

@@ -165,10 +165,10 @@ splitcs(::Tuple{}, csnow, csrem) = (csnow, csrem)
 
 #Determine wether output and temp array can a) be identical b) share memory through reshape or 
 # c) need to be allocated individually
-function output_aliasing(di::DiskIndex)
+function output_aliasing(di::DiskIndex, ndims_dest, ndims_source)
     if all(i->isa(i,Union{Int,AbstractUnitRange,Colon}),di.temparray_indices) && 
         all(i->isa(i,Union{Int,AbstractUnitRange,Colon}),di.output_indices)
-        if di.output_size == di.temparray_size
+        if di.output_size == di.temparray_size && ndims_dest == ndims_source
             return :identical
         else 
             return :reshapeoutput
@@ -229,7 +229,7 @@ function getindex_disk_nobatch!(out,a,i)
     indices = resolve_indices(a, i, NoBatch(allow_steprange(a), 1.0))
     #@debug output_size, temparray_size, output_indices, temparray_indices, data_indices
     outputarray = create_outputarray(out, a, indices.output_size)
-    outalias = output_aliasing(indices)
+    outalias = output_aliasing(indices,ndims(outputarray),ndims(a))
     if outalias === :identical
         readblock!(a, outputarray, indices.data_indices...)
     elseif outalias === :reshapeoutput
@@ -297,7 +297,7 @@ end
 
 function setindex_disk_nobatch!(a,v,i)
     indices = resolve_indices(a, i, NoBatch())
-    outalias = output_aliasing(indices)
+    outalias = output_aliasing(indices,ndims(a),ndims(v))
     if outalias === :identical
         writeblock!(a, v, indices.data_indices...)
     elseif outalias === :reshapeoutput

test/runtests.jl

@@ -75,7 +75,7 @@ function test_setindex(a)
     a[1, 1, 1] = 1
     a[1, 2] = 2
     a[1, 3, 1, 1] = 3
-    a[2, :] = [1, 2, 3, 4, 5]
+    a[2:2, :] = [1, 2, 3, 4, 5]
     a[3, 3:4, 1, 1] = [3, 4]
     # Test bitmask indexing
     m = falses(4, 5, 1)
@@ -468,21 +468,20 @@ end
     a = AccessCountDiskArray(reshape(1:20, 4, 5, 1); chunksize=(4, 1, 1))
     i = (1:3,:,:)
     di = DiskArrays.resolve_indices(a,i,DiskArrays.NoBatch())
-    @test DiskArrays.output_aliasing(di) == :identical
+    @test DiskArrays.output_aliasing(di,3,3) == :identical
+    @test DiskArrays.output_aliasing(di,2,3) == :reshapeoutput
     i = (1,:,:)
     di = DiskArrays.resolve_indices(a,i,DiskArrays.NoBatch())
-    @test DiskArrays.output_aliasing(di) == :reshapeoutput
+    @test DiskArrays.output_aliasing(di,2,2) == :reshapeoutput
     i = ([1,3],:,:)
     di = DiskArrays.resolve_indices(a,i,DiskArrays.NoBatch())
-    @test DiskArrays.output_aliasing(di) == :noalign
+    @test DiskArrays.output_aliasing(di,3,3) == :noalign
     i = (1:3,:)
     di = DiskArrays.resolve_indices(a,i,DiskArrays.NoBatch())
-    @test DiskArrays.output_aliasing(di) == :reshapeoutput
+    @test DiskArrays.output_aliasing(di,2,2) == :reshapeoutput
     i = (1:3,:,:,1,1)
     di = DiskArrays.resolve_indices(a,i,DiskArrays.NoBatch())
-    @test DiskArrays.output_aliasing(di) == :identical
-
-
+    @test DiskArrays.output_aliasing(di,3,3) == :identical
 end