How do I do this with the arrayop macro in Symbolics.jl?
using Tullio, Symbolics
@variables x₀
q = x₀*rand(4, 4)
t = x₀*rand(4, 4, 4)
@tullio r[a, b, c] := q[a, i]*t[i, b, c]
I tried this
r = Symbolics.@arrayop (a, b, c) q[a, i]*t[i, b, c]
r = collect(r)
but it does not work.
Thank you