Symbolics Einstein notation

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


1 Like