seems synchronization is a real issue. Also not sure why I declared matrices in Davidson to be localized. At least it is better as diagonalization is always sync'ed in this way