In this paper we introduce fast numerical algorithms for the solution of the model [Total variation regularization cost function for demodulating phase discontinuities, Journal Applied Optics, 53(11):2297-2301, 2014]. For each variable, background illumination, amplitude modulation and phase map, we develop a fixed point method. Then, we write all three algorithms in the same framework and analyze their convergence rates, local smoothing factors by means of Local Fourier Analysis and present experimental evidence of their performance on synthetic and real world problems.