Flow chart depicting the multi-step iterative block sampling procedure exemplified for two data sets. In the first step, a three-dimensional density field will be realized conditional on the galaxy observations. In subsequent steps, foreground template coefficients, the bias parameters, the power spectrum, and the normalization parameters for the galaxy distribution are sampled conditional on respective previous samples. Iteration of this process yields samples from the full joint posterior distribution.