worked on mean-field example
merge
for consistency, use LR/SR with sp.linalg.eigs