Basic HTML version of Foils prepared 15 March 1996

Foil 23 Basis of Estimates for Communication Costs for ADI

From Spatial Differencing and ADI Solution of the NAS Benchmarks CPSP713 Case studies in Computational Science -- Spring Semester 1996. by Geoffrey C. Fox


Communication is needed when one sweeps in a given index, say i, and that index location changes
Note Algorithms for sequential tridiagonal solvers are given in the appendix of volume 1 of Hirsch
For a scalar tridiagonal matrix with N rows, the solution takes time 5N
For a block tridiagonal matrix, one increases this by a factor between 25 and 125 for 5 by 5 blocks.
Note communication proportional to "5" i.e. size of vector being calculated as one changes block size
However calculation proportional to 52 to 53 as one changes block size i.e. it increases faster with block size than communication
  • This illustrates that extra CFD components (1 goes to 5) increases efficiency of parallel CFD codes compared to comparable scalar equations



Northeast Parallel Architectures Center, Syracuse University, npac@npac.syr.edu

If you have any comments about this server, send e-mail to webmaster@npac.syr.edu.

Page produced by wwwfoil on Sun Feb 22 1998