A Program for Vectorized Finite-Difference Migration
, by Stewart A. Levin
In SEP-38 and SEP-41 I described a method for finite-difference migration
on a parallel computer. Unfortunately, parallel computers do not abound
yet. But pipelined vector computers do. Here I show by example
that my parallel ideas also vectorize. I provide a Fortran program
for 15 degree time-domain finite-difference migration that
1) vectorizes every inner loop,
2) handles both seismic data and migration velocities in natural trace ordering,
3) is short, simple and clear, and
4) easily converts to higher order finite-differences.
This program is the fastest and most sensible adaptation of 15 degree time-domain migration to a large, vector computer I have ever seen and I predict that it will soon be used routinely within the geophysical industry.