Goertzel algorithm produces incorrect phase
The short answer is that when Goertzel algorithm is used for omega, which doesn't correspond to of the finite set of frequencies computed by DFT (which correspond to integer k), then an additional correction factor should be applied. I've got an answer from Peter K after describing my question more clearly here: What&#39;s wrong with my Goertzel algorithm implementation?