 By introducing a homogeneous reference medium, the heterogeneity of the medium is transformed into a polarization constraint.The Green operator of the reference medium, known explicitly in Fourier space, can be used to iteratively update the polarization field.Several improvements and accelerations have been made to this method, which is now used internationally in dedicated codes.