qmcpack
qmcpack copied to clipboard
Default value for cutoff/radius in open or semi-open BC
Describe the bug It was spotted that a few features set cutoff/radius to magic numbers when user input doesn't exist. In the case periodic boundary condition (ppp), wigner seitz radius is a safe default. However, with a open or semi-open boundary condition, the current magic number 10 is generally considered harmful. I feel it better to stop the code and demand user input.
It seems distance based jastrow factors and gofr estimator all carry a similar issue.
Is 10 B used as rcut for all Jastrow factors in the open BC case?