Initialization of density matrix via NDM

[Sukesh01]

Hi, I couldn't figure out the method for preparing the variational mixed state in arbitrary ( prior) defined density matrix, for example, if |S>= singlet state between two spins, rho = |S><S|. Though the documentation explains about finding the steady state, I cannot find any code snippets, which helps me. Can someone clarify and provide code/example?

[PhilipVinc]

Nothing is built into netket.
In general, you'd want to maximise the overlap with a given state of which you know the entries.
But you have to write it yourself.

[Sukesh01]

Can you guide me how to do that? Ive tried to generate a random NDM variational state and try to optimise params to overlap with known state. But it is raising few effors. Some guidance would be really appreciated, professor!