Geometrical tools to map systems of affine recurrence equations on regular arrays

We propose a method based on geometrical tools to map problems onto regular and synchronous processor arrays. The problems we consider are defined by systems of affine recurrence equations (SARE). From such a problem specification we extract the data dependencies in terms of two classes of vectors: the utilization vectors and the dependence vectors. We use these vectors to express constraints on the timing or the allocation functions. We differentiate two classes of constraints. The causal ones are intrinsic timing constraints induced by the system of equations defining the problem. A given choice of target architecture may impose new constraints on the timing or the allocation. We call them the architecture-related constraints. We use these constraints to determine first an affine timing function and next an allocation by projection. We finally illustrate the method with three examples: the matrix multiplication, the recursive convolution and the LLt Cholesky factorization.

[1]  Dan I. Moldovan,et al.  Partitioning and Mapping Algorithms into Fixed Size Systolic Arrays , 1986, IEEE Transactions on Computers.

[2]  H. T. Kung Why systolic architectures? , 1982, Computer.

[3]  P. Clauss Synthèse d'algorithmes systoliques et implantation optimale en place sur réseaux de processeurs synchrones , 1990 .

[4]  Patrice Quinton Automatic synthesis of systolic arrays from uniform recurrent equations , 1984, ISCA '84.

[5]  Philippe Clauss,et al.  A geometrical coding to compile affine recurrence equations on regular arrays , 1991, [1991] Proceedings. The Fifth International Parallel Processing Symposium.

[6]  Thomas Kailath,et al.  Regular iterative algorithms and their implementation on processor arrays , 1988, Proc. IEEE.

[7]  H. T. Kung,et al.  Systolic Arrays for (VLSI). , 1978 .

[8]  Patrice Quinton,et al.  Computability of recurrence equations , 1993, Theor. Comput. Sci..

[9]  C. Mongenet,et al.  Calculus of space-optimal mappings of systolic algorithms on processor arrays , 1990, [1990] Proceedings of the International Conference on Application Specific Array Processors.

[10]  Yves Robert,et al.  Systolic solution of the Algebraic Path Problem , 1987 .

[11]  Guy-René Perrin,et al.  Synthesis of Systolic arrays for Inductive Problems , 1987, PARLE.

[12]  Sanjay V. Rajopadhye Synthesizing systolic arrays with control signals from recurrence equations , 2005, Distributed Computing.

[13]  D.I. Moldovan,et al.  On the design of algorithms for VLSI systolic arrays , 1983, Proceedings of the IEEE.

[14]  Peter R. Cappello,et al.  Converting affine recurrence equations to quasi-uniform recurrence equations , 1988, J. VLSI Signal Process..

[15]  Catherine Mongenet,et al.  Affine Timings for Systems of Affine Recurrence Equations , 1991, PARLE.

[16]  Richard M. Karp,et al.  The Organization of Computations for Uniform Recurrence Equations , 1967, JACM.

[17]  Patrice Quinton,et al.  The mapping of linear recurrence equations on regular arrays , 1989, J. VLSI Signal Process..

[18]  Sun-Yuan Kung,et al.  Optimal Systolic Design for the Transitive Closure and the Shortest Path Problems , 1987, IEEE Transactions on Computers.

[19]  P. Quinton Automatic synthesis of systolic arrays from uniform recurrent equations , 1984, ISCA 1984.

[20]  Ph. Clauss,et al.  Synthesis of process arrays , 1989 .

[21]  W. Blankinship A New Version of the Euclidean Algorith , 1963 .

[22]  Benjamin W. Wah,et al.  Systematic approaches to the design of algorithmically specified systolic arrays , 1985, ICASSP '85. IEEE International Conference on Acoustics, Speech, and Signal Processing.

[23]  Alexander Schrijver,et al.  Theory of linear and integer programming , 1986, Wiley-Interscience series in discrete mathematics and optimization.

[24]  Peter R. Cappello,et al.  Scheduling a system of affine recurrence equations onto a systolic array , 1988, [1988] Proceedings. International Conference on Systolic Arrays.

[25]  Philippe Clauss,et al.  Synthesis of size-optimal toroidal arrays for the Algebraic Path Problem: A new contribution , 1992, Parallel Comput..

[26]  G. Rote A systolic array algorithm for the algebraic path problem (shortest paths; Matrix inversion) , 1985, Computing.