optimization can execute colorization tasks better under dual consistency terms in wavelet domain, namely data-consistency
Then, sampling is performed by annealing Langevin dynamics and alternative data-consistency update.