@cgw I tried to check everything. For example the preconditioner choice did not affect the results or savage=2 in older versionis the same as the default value in 2021 version. I just want to let you know that we can get really better results using the code by some modifications. I will try my best to implement those modifications in the new version to make everything easier. However, some subroutines have evolved completely.