Update SIMDPirates version and correctly calculate loopdependencies of variables marked as possible reductions that are not.

Author: chriselrod

Project.toml

@@ -15,7 +15,7 @@ VectorizationBase = "3d5dd08c-fd9d-11e8-17fa-ed2836048c2f"
 [compat]
 DocStringExtensions = "0.8"
 OffsetArrays = "1"
-SIMDPirates = "0.7.13"
+SIMDPirates = "0.7.14"
 SLEEFPirates = "0.4.4"
 UnPack = "0"
 VectorizationBase = "0.10.4"

benchmark/driver.jl

@@ -45,7 +45,7 @@ logdettriangle_bench = benchmark_logdettriangle(sizes)
 
 v = 1
 filetype = "svg"
-const PICTURES = joinpath(pkgdir("LoopVectorization"), "docs", "src", "assets")
+const PICTURES = joinpath(pkgdir(LoopVectorization), "docs", "src", "assets")
 save(joinpath(PICTURES, "bench_filter2d_dynamic_v$v.$filetype"), plot(filter2d_dynamic_bench));
 save(joinpath(PICTURES, "bench_filter2d_3x3_v$v.$filetype"), plot(filter2d_3x3_bench));
 save(joinpath(PICTURES, "bench_filter2d_unrolled_v$v.$filetype"), plot(filter2d_unrolled_bench));

docs/src/assets/bench_AmulB_v1.svg

@@ -240,6 +240,7 @@ function add_compute!(
         # parent = getop(ls, var, elementbytes)
         if length(reduceddeps) == 0
             insert!(vparents, reduction_ind, parent)
+            mergesetv!(deps, loopdependencies(parent))
             op = Operation(length(operations(ls)), var, elementbytes, instruction(ls,instr), compute, deps, reduceddeps, vparents)
             pushop!(ls, op, var)
         else

docs/src/assets/bench_AmulBt_v1.svg

@@ -153,6 +153,20 @@ using Test
         end
     end
 
+    function setcolumstovectorplus100!(Z::AbstractArray{T}, A) where {T}
+        for i = axes(A,1), j = axes(Z,2)
+            acc = zero(T)
+            acc = acc + A[i] + 100
+            Z[i, j] = acc
+        end
+    end
+    function setcolumstovectorplus100avx!(Z::AbstractArray{T}, A) where {T} 
+        @avx for i = axes(A,1), j = axes(Z,2)
+            acc = zero(T)
+            acc = acc + A[i] + 100
+            Z[i, j] = acc
+        end
+    end
 
     function mvp(P, basis, coeffs::Vector{T}) where {T}
         len_c = length(coeffs)
@@ -557,6 +571,11 @@ using Test
         @test dot3v2avx(x, A, y) ≈ d3
         @test dot3_avx(x, A, y) ≈ d3
 
+        A2 = similar(A);
+        setcolumstovectorplus100!(A, x)
+        setcolumstovectorplus100avx!(A2, x)
+        @test A == A2
+        
         maxdeg = 20; nbasis = 1_000; dim = 15;
         r = T == Float32 ? (Int32(1):Int32(maxdeg+1)) : (1:maxdeg+1)
         basis = rand(r, (dim, nbasis));