Scalable techniques for computing band linear recurrences on massively parallel and vector supercomputers

Presents a new scalable algorithm, called the `regular schedule´, for parallel evaluation of band linear recurrences (BLRs), i.e. mth-order linear recurrences for m>1. Its scalability and simplicity make it well suited for vector supercomputers and massively parallel computers. We describe our implementation of the regular schedule on two types of machines: the Convex C240 and the MasPar MP-2. The scalability of our scheduling techniques is demonstrated on the two machines. Significant improvements in CPU performance for a range of programs containing BLRs implemented using the regular schedule in C over the same programs implemented using the highly-optimized coded-in-assembly BLAS routines are demonstrated on the Convex C240. We also demonstrate the scalability of this schedule on the MasPar MP-2 with up to 2000 processors. Our approach can be used both at the user level, in parallel programming code containing BLRs, and in compiler parallelization of such programs combined with recurrence recognition techniques for massively parallel and vector supercomputers