Flow chart depicting the multi-step iterative block sampling procedure. In the first step, a three-dimensional density field will be realized conditional on the galaxy observations (top left corner). In subsequent steps observer velocity, bias parameters 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 generated by the BORG algorithm.