A procedure analogous to Gibbs sampling, with the exception that the mode of each conditional posterior distribution is determined at each update, rather than sampling a value from these conditional distributions.