论文信息 - Solving Linear Recurrence Systems Using Level 2 and 3 BLAS Routines

Solving Linear Recurrence Systems Using Level 2 and 3 BLAS Routines

The aim of this paper is to present a new efficient BLAS-based algorithm for solving linear recurrence systems with constant coefficients, which can be easily and efficiently implemented on shared or distributed memory machines and clusters of workstations. The algorithm is based on level 3 and level 2 BLAS routines _GEMM, _GEMV and _TRMV, which are crucial for its efficiency even when the order of a system is relatively high. The results of experiments performed on a dual-processor Pentium III computer are also presented and discussed.

Przemyslaw Stpiczynski | Przemysław Stpiczyński

[1] Przemyslaw Stpiczynski. Parallel Algorithms for Solving Linear Recurrence Systems , 1992, CONPAR.

[2] H. A. van der Vorst,et al. Vectorization of Linear Recurrence Relations , 1989 .

[3] Elizabeth R. Jessup,et al. Toward Memory-Efficient Linear Solvers , 2002, VECPAR.

[4] Josep-Lluís Larriba-Pey,et al. Review of General and Toeplitz Vector Bidiagonal Solvers , 1996, Parallel Comput..

[5] J. Stoer,et al. Introduction to Numerical Analysis , 2002 .

[6] Almerico Murli,et al. Algorithm 682: Talbot's method of the Laplace inversion problems , 1990, TOMS.

[7] Marcin Paprzycki,et al. Parallel solution of linear recurrence systems , 1996 .

[8] Jack J. Dongarra,et al. Solving linear systems on vector and shared memory computers , 1990 .

[9] Roberto Barrio,et al. On the numerical evaluation of linear recurrences , 2003 .

[10] Przemyslaw Stpiczynski,et al. A New Message Passing Algorithm for Solving Linear Recurrence Systems , 2001, PPAM.

[11] D. Sorensen,et al. LAPACK Working Note No. 2: Block reduction of matrices to condensed forms for eigenvalue computations , 1987 .

[12] Begnaud Francis Hildebrand,et al. Introduction to numerical analysis: 2nd edition , 1987 .

[13] Przemyslaw Stpiczynski,et al. Fast parallel algorithms for computing trigonometric sums , 2002, Proceedings. International Conference on Parallel Computing in Electrical Engineering.