Expression of the data consistency layer

[zaccharieramzi]

I have a question regarding the expression of the output of the data consistency layer in the noisy case.

We have out = (x + v * x_sampled) / (1 + v) both here and here.

But this is done without consideration of whether the sampled point is in sampled region or not. I.e., the formal expression represented by this code would be, instead of eq. 6 of your paper:

x_rec(k) = x_cnn(k) / (1 + v) if k \notin Omega; (x_cnn(k) + v * x_u(k)) / (1 + v) else

I removed the \hat to make it more readable and replaced lambda by v.

I think the correct implementation of eq. 6 should be:

out = x + (v * (x_sampled - mask * x)) / (1 + v)

What do you think?