Bug when using Tullio + LoopVectorization about tullio.jl HOT 5 CLOSED

julia> A = ones(2, 2, 2, 2);  # also wrong for e.g. A = ones(8,10,8,10)

julia> s_for = sum(A[a,b,a,b] for a in axes(A,1), b in axes(A,2))
4.0

julia> using LoopVectorization

julia> function sum169(A, 𝒶𝓍a=axes(A,1), 𝒶𝓍b=axes(A,2)) 
           𝒜𝒸𝒸 = 0.0
           LoopVectorization.@avx unroll = 0 for b in 𝒶𝓍b
              for a = 𝒶𝓍a
                  𝒜𝒸𝒸 = 𝒜𝒸𝒸 + A[a, b, a, b]
              end
           end
           𝒜𝒸𝒸
       end
sum169 (generic function with 3 methods)

julia> sum169(A)
4.677771489e-314

julia> @macroexpand LoopVectorization.@avx unroll = 0 for b in 𝒶𝓍b
              for a = 𝒶𝓍a
                  𝒜𝒸𝒸 = 𝒜𝒸𝒸 + A[a, b, a, b]
              end
           end
quote
    begin
        #= REPL[15]:5 =#
        nothing
    end
    begin
        #= REPL[15]:4 =#
        nothing
    end
    begin
        #= REPL[15]:3 =#
        nothing
    end
    begin
        #= REPL[15]:2 =#
        nothing
    end
    begin
        #= REPL[15]:1 =#
        nothing
    end
    var"###looprangeb###1###" = LoopVectorization.canonicalize_range(begin
                $(Expr(:inbounds, true))
                local var"#68#val" = 𝒶𝓍b
                $(Expr(:inbounds, :pop))
                var"#68#val"
            end)
    var"###looplenb###2###" = StaticArrayInterface.static_length(var"###looprangeb###1###")
    var"###b_loop_lower_bound###3###" = LoopVectorization.maybestaticfirst(var"###looprangeb###1###")
    var"###b_loop_upper_bound###4###" = LoopVectorization.maybestaticlast(var"###looprangeb###1###")
    var"###b_loop_step###5###" = LoopVectorization.static_step(var"###looprangeb###1###")
    var"###looprangea###6###" = LoopVectorization.canonicalize_range(begin
                $(Expr(:inbounds, true))
                local var"#69#val" = 𝒶𝓍a
                $(Expr(:inbounds, :pop))
                var"#69#val"
            end)
    var"###looplena###7###" = StaticArrayInterface.static_length(var"###looprangea###6###")
    var"###a_loop_lower_bound###8###" = LoopVectorization.maybestaticfirst(var"###looprangea###6###")
    var"###a_loop_upper_bound###9###" = LoopVectorization.maybestaticlast(var"###looprangea###6###")
    var"###a_loop_step###10###" = LoopVectorization.static_step(var"###looprangea###6###")
    if LoopVectorization.check_args(A, typeof(𝒜𝒸𝒸)) && (true && $(Expr(:&&, :((LoopVectorization.can_turbo)(LoopVectorization.add_fast, Val{2}())))))
        (var"##vptr##_A", var"#A#preserve#buffer#") = LoopVectorization.stridedpointer_preserve(A)
        var"##vptr##_A##ind##1##repeated##3##" = LoopVectorization.VectorizationBase.double_index(var"##vptr##_A", Val{0}(), Val{2}())
        var"##vptr##_A##ind##1##repeated##3####ind##3##repeated##3##" = LoopVectorization.VectorizationBase.double_index(var"##vptr##_A##ind##1##repeated##3##", Val{2}(), Val{2}())
        var"####grouped#strided#pointer####12###" = (getfield)(LoopVectorization.grouped_strided_pointer((LoopVectorization.densewrapper(LoopVectorization.gespf1(var"##vptr##_A##ind##1##repeated##3####ind##3##repeated##3##", (var"###a_loop_lower_bound###8###", var"###b_loop_lower_bound###3###")), A),), Val{()}()), 1)
        begin
            $(Expr(:gc_preserve, :(var"##𝒜𝒸𝒸_##onevec##" = begin
          var"##vargsym#257" = ((LoopVectorization.zerorangestart(var"###looprangeb###1###"), LoopVectorization.zerorangestart(var"###looprangea###6###")), (var"####grouped#strided#pointer####12###", Base.eltype(𝒜𝒸𝒸)))
          var"##Tloopeltype##" = LoopVectorization.promote_type(LoopVectorization.eltype(A), Base.eltype(𝒜𝒸𝒸))
          var"##Wvecwidth##" = LoopVectorization.pick_vector_width(var"##Tloopeltype##")
          LoopVectorization._turbo_!(LoopVectorization.avx_config_val(Val{(false, 0, 0, 0, false, 0x0000000000000001, 1, true)}(), var"##Wvecwidth##"), Val{(Symbol("##DROPPED#CONSTANT##"), Symbol("##DROPPED#CONSTANT##"), LoopVectorization.OperationStruct(0x00000000000000000000000000000000, 0x00000000000000000000000000000000, 0x00000000000000000000000000000021, 0x00000000000000000000000000000000, 0x00000000000000000000000000000000, 0x00000000000000000000000000000000, 0x00000000000000000000000000000000, LoopVectorization.constant, 0x0001, 0x0000), :LoopVectorization, :getindex, LoopVectorization.OperationStruct(0x00000000000000000000000000000021, 0x00000000000000000000000000000000, 0x00000000000000000000000000000000, 0x00000000000000000000000000000000, 0x00000000000000000000000000000000, 0x00000000000000000000000000000000, 0x00000000000000000000000000000000, LoopVectorization.memload, 0x0002, 0x0001), :LoopVectorization, :add_fast, LoopVectorization.OperationStruct(0x00000000000000000000000000000021, 0x00000000000000000000000000000021, 0x00000000000000000000000000000000, 0x00000000000000000000000000010002, 0x00000000000000000000000000000000, 0x00000000000000000000000000000000, 0x00000000000000000000000000000000, LoopVectorization.compute, 0x0001, 0x0000))}(), Val{(LoopVectorization.ArrayRefStruct{:A, Symbol("##vptr##_A##ind##1##repeated##3####ind##3##repeated##3##")}(0x00000000000000000000000000000101, 0x00000000000000000000000000000201, 0x00000000000000000000000000000000, 0x00000000000000000000000000000101),)}(), Val{(0, (3,), (1,), (), (), (), ())}(), Val{(:b, :a)}(), Base.Val(Base.typeof(var"##vargsym#257")), LoopVectorization.flatten_to_tuple(var"##vargsym#257")...)
      end), Symbol("#A#preserve#buffer#")))

        (var"##vptr##_A", var"#A#preserve#buffer#") = LoopVectorization.stridedpointer_preserve(A)
        var"##vptr##_A##ind##1##repeated##3##" = LoopVectorization.VectorizationBase.double_index(var"##vptr##_A", Val{0}(), Val{2}())
        var"##vptr##_A##ind##1##repeated##3####ind##3##repeated##3##" = LoopVectorization.VectorizationBase.double_index(var"##vptr##_A##ind##1##repeated##3##", Val{2}(), Val{2}())