Hi,

Sorry if this has been answered already but I am trying to do some post-processing with the MPS Tensor components and it would be really helpful to just convert them into regular matrices. Specifically, I have some MPS tensor "T" on N sites on a spin-1/2 chain. I can grab the two local components of the MPS using commands like -

```
using ITensors
N =10
sites = siteinds("S=1/2", N)
T = randomMPS(sites, 20)
T5_1=T[5]*setelt(sites[5]=>1)
T5_2=T[5]*setelt(sites[5]=>2)
```

However, I would like to be able to perform matrix multiplications like (T5*1).\rho.(T5*1)^{\dagger} (where \rho is a suitable density matrix) and for this it would be best to have them converted into matrices.

PS: I know that the other way of doing this is to encode the density matrix into an ITensor object and change the link indices suitably to get the multiplication to work out.