fix issue #407

Author: chriselrod

Project.toml

@@ -1,7 +1,7 @@
 name = "LoopVectorization"
 uuid = "bdcacae8-1622-11e9-2a5c-532679323890"
 authors = ["Chris Elrod <[email protected]>"]
-version = "0.12.113"
+version = "0.12.114"
 
 [deps]
 ArrayInterface = "4fba245c-0d91-5ea0-9b3e-6abc04ee57a9"

src/parse/memory_ops_common.jl

@@ -894,7 +894,11 @@ function array_reference_meta!(
           move_to_last!(offsets, ind_prev_index)
           move_to_last!(strides, ind_prev_index)
           move_to_last!(loopedindex, ind_prev_index)
-          move_to_last!(loopdependencies, ind_prev_index)
+          if length(loopdependencies) == length(indices)
+            move_to_last!(loopdependencies, ind_prev_index)
+          else
+            move_to_last!(loopdependencies, findfirst(Base.Fix2(===, ind), loopdependencies)::Int)
+          end
           vptrarray = repeated_index!(
             ls,
             indices,

test/simplemisc.jl

@@ -1,4 +1,10 @@
-
+import LinearAlgebra
+function issue407!(Hdata::Matrix{Complex{T}}) where {T}
+    Hview=reinterpret(reshape,T,Hdata)
+    @tturbo for i in 1:size(Hdata,1)
+        Hview[1, i, i] = i
+    end 
+end
 function lv_turbo(r::Vector{UInt64}, mask::UInt64)
   @turbo for i = 1:length(r)
     r[i] &= mask
@@ -13,3 +19,10 @@ end
   r2 = r1 .& mask
   @test lv_turbo(r1, mask) == r2
 end
+@testset "issue 407" begin
+  @show @__LINE__
+  A = zeros(Complex{Float64}, 10, 10);
+  issue407!(A)
+  @test real.(A) == Diagonal(1:10)
+  @test all(iszero, imag.(A))
+end