Now the state specific calculations with the restart option should be working. The normal statespecific calculations is sort of like a full restart, where all the wavefunctions are canonicalized and then a generate block is done and finally the normal dmrg targetting each state in turn is performed. The restart option is used when the normal state specific calculations fails at some stage. The restart flag first canonicalizes all the wavefunctions and then initializes the overlaps between wavefunctions. Finally it starts optimizing the last wavefunction it was working with and begins at the sweep iteration which it had previously died.