bug fixes, at spin symmetric limit, result is right, but AF limit seems still wrong