Proof of Correctness for Sparse Tiling of Gauss-Seidel

Michelle Mills Strout, Larry Carter and Jeanne Ferrante
CS2003-0741
April 1, 2003

Gauss-Seidel is an iterative computation used for solving a set of simultaneous linear equations, A\vec{u}=\vec{f}. If the matrix A uses a sparse matrix representation, storing only nonzeros, then the data dependences in the computation arise from A's nonzero structure. We use this structure to schedule the computation at runtime using a technique called {\em full sparse tiling}. The sparse tiled computation exhibits better data locality and therefore improved performance. This paper gives a complete proof that a serial schedule for full sparse tiled Gauss-Seidel generates results equivalent to those that a typical Gauss-Seidel computation produces. We also provide implementation and correctness details for full sparse tiling with reduced worst-case complexity. [NOTE: Dave Wargo told me this techreport would be assigned #CS2003-0741 and I am assuming it will be April 2003 as the date, please let me know if these assumptions are incorrect.]


How to view this document


The authors of these documents have submitted their reports to this technical report series for the purpose of non-commercial dissemination of scientific work. The reports are copyrighted by the authors, and their existence in electronic format does not imply that the authors have relinquished any rights. You may copy a report for scholarly, non-commercial purposes, such as research or instruction, provided that you agree to respect the author's copyright. For information concerning the use of this document for other than research or instructional purposes, contact the authors. Other information concerning this technical report series can be obtained from the Computer Science and Engineering Department at the University of California at San Diego, techreports@cs.ucsd.edu.


[ Search ]


NCSTRL
This server operates at UCSD Computer Science and Engineering.
Send email to webmaster@cs.ucsd.edu