Implemented Stefans idea: do the same as Goh but use a sampled version of the binary units where the sampling is expected to be sparse for the latent bias (seems to do something similar, but faster than the Goh method because no sorting is required)