A question about your suggestion:
" Here’s a simple one you can try: choose a certain order on the variables
(say row by row in the part of the image you want to complete), and then for each node update all
messages from it to its neighbors). "
Just want to make sure that I understand well.
The messages that are outgoing from observed pixels are always 1 (or some number that we generate randomly, but for simplicity assume that I chose 1), and we are calculating only the outgoing messages of unobserved pixels, and updating only those messages while iterating.
Then the influence of the observed pixels over the unobserved pixels messages will be only by the pairwise function and not by the outgoing messages (because they are always 1).
Correct?