Closed haakoek closed 4 years ago
When using OATDCC
this is the case as seen here. For the ground state situation I think it makes sense to change basis as this is needed for other observables. However, I think the reason the basis does not change by default is to avoid "silently" mucking up a system. I am inclined to agree that the basis should change.
https://github.com/Schoyen/coupled-cluster/blob/master/coupled_cluster/ccd/oaccd.py
In the call
def compute_ground_state( self, max_iterations=100, tol=1e-4, termination_tol=1e-4, tol_factor=0.1, change_system_basis=False, **mixer_kwargs, )
should change_system_basis=True by default? I think that makes most sense at the point, but there might be some good arguments for why not doing this.