This paper concerns the systolic array computation of the generalized singular value decomposition. Numerical algorithms for both one and two-dimensional systolic architectures are discussed