I am wondering of one can use ITensors to solve linear equations in the form of

A X = B

where all of the components are tensor networks or high order tensors in some compressed format? Or if not does anybody know a package that does that?

My tensors represent high order polynomials of large number of variables, hence they a symmetric. This needs some compressed format, so tensor networks, tensor train or H-tucker are the candidates. I need some machinery that can substitute polynomials into each other, multiply them and then solve these multi-linear equations.

Many thanksâ€¦