Hi Jacopo,

This is a good question, but unfortunately the answer is "no" right now. That's the short answer.

The longer answer is that we do plan to add this feature, so let me tell you our idea for how we want to do it, because maybe you can do it yourself if it's an important thing for you to have. What we have learned is that for making a random MPS of bond dimension m > 1 it's not really right to just randomize each tensor separately. Instead, the right thing to do, we think, is to make a unitary circuit where each 2-site unitary gate is a random unitary. Then by applying D layers of this circuit using a TEBD-like algorithm to a starting product state, one gets a (correctly) random MPS of bond dimension 2^D. This would not be too hard to code, as long as one figures out a correct way to make the random unitaries. We have a sample TEBD code on the website here: http://itensor.org/docs.cgi?vers=cppv3&page=formulas/tevol_trotter

For the quantum number conserving version there are some possibly tricky questions about the right way to make the random unitaries we aren't as sure about. So we'd have to do some thinking and some reading to be sure we know how to do that correctly. But we'd like to develop the non-QN version soon anyway, it's just right now we have some other pressing things to do.

Best regards,

Miles