Fix auto-growing in go().
Fix old typo in matmul (simple_diag_matrix). New numpy seems to bring it out.
Fix same bug in eyemat.
Fix an example.
Add block length parameter to heisenberg_af_uniform for S=1/2 case.
Fix pinv tests for multi_uni. Only testing L=1 for now.
Improve sandwich scattering example.