A three-dimensional systolic architecture is described that efficiently performs passive broadband interpolation beam-forming for an array of sensors. The architecture transforms from time domain to beam space without intermediate data storage. The required convolution operations are implemented with linear sets of multiply/accumulate elements.