A random-sampling method as an efficient alternative to variational Monte Carlo for solving Gutzwiller wavefunctions. (7th December 2021)