Tullio broadcasting issue

New to Julia
, ,
2

Hi @swanchristmas, can you provide a minimal reproducable example? i.e. Please read: make it easier to help you

For what it’s worth, if I try to construct such an example, I see no difference here:

julia> using Tullio, TensorOperations

julia> let 
           N, P = 2,3
           A = rand(N)
           B = rand(P)
           C = rand(N)
           @tullio out1[n, p] := A[n] * B[p] + C[n]
           @tullio out2[n, p] := A[n] * B[p] + C[n] * ones(P)[p]
           @tensor out3[n, p] := A[n] * B[p] + C[n] * ones(P)[p]
          
           @info "" out1 out2 out3
       end
┌ Info: 
│   out1 =
│    2×3 Matrix{Float64}:
│     0.670981  0.719456  0.856811
│     0.577745  0.641696  0.822899
│   out2 =
│    2×3 Matrix{Float64}:
│     0.670981  0.719456  0.856811
│     0.577745  0.641696  0.822899
│   out3 =
│    2×3 Matrix{Float64}:
│     0.670981  0.719456  0.856811
└     0.577745  0.641696  0.822899

Is it possible there’s an error or typo somewhere in your code that’s causing the discrepancy you’re seeing?