As with the Helmholtz equation, the complex scale-factor, , means is symmetric, but not Hermitian, so the standard Kolmogorov factorization will fail. As discussed in Chapter , however, the method can be extended to factor any cross-spectrum into a pair of minimum phase wavelets and a delay Claerbout (1998c). The algorithm follows the standard Kolmogorov factorization; however, negative lags are kept separately rather than being discarded.
The Kolmogorov factorization is not exact because the filters are factored in the frequency domain, assuming circular boundary conditions; while the polynomial division is performed in the time domain with transient boundary conditions. As a result the filters must be padded in the time-domain before spectral factorization. Padding does not significantly effect the overall cost of the migration, as the computational expense lies in the polynomial division, not in the factorization.
Alternative methods for cross-spectral factorization may avoid the circular boundary condition problem. For example, the Wilson-Burg algorithm Sava et al. (1998); Wilson (1969), based on Newton's recursive linearization, can efficiently factor polynomials, and is especially suited to the helical coordinate system.