Midterm Project on Iterative Techniques for Solving Linear Systems.
Due March 23, 1998

MATH 4660-001 and C SC 4666-001: Numerical Analysis II
Spring 1998. University of Colorado at Denver

Students are welcome to communicate with each other and the teacher on-line using the Web page
http://www-math.cudenver.edu/~aknyazev/teaching/98/4660/comments

Computer Project Description:

Write and test the MATLAB program of the Jacobi method for sparse matrices. Use it to solve the system with the matrix NOS6:from set LANPRO, from the Harwell-Boeing Collection. Take the right-hand side vector with all components to be one. Use the zero vector as the initial guess. Use MATLAB build-in functions \ and pcg to solve it, too. Discuss numerical results and perfomance. Write a report and publish it on the Internet.

Hints:

• In MATLAB, always start your program with:
  clear all; format long e
• To write an efficient program, take advantange of the sparse representation of the matrix. See help sparfun to learn about build-in functions.
• I tested the MATLAB conversion program from MATRIXMART format.

• Solve for a family of marices with increasing diagonal dominance. Make a MATLAB movie of convergence history when diagonal dominance increases. Use MATLAB program 'mpgwrite' to convert the MATLAB movie into MPEG format. You do not have to make our own copy of mpgwrite any longer. We just made mpgwrite a part of our distribution of MATLABs; type help mpgwrite to learn how to use it. Feel free to use my sample MATLAB program makemovie.m to make the MPEG movie movie1.mpg. Include the MPEG movie in the report on the Internet.
• Find a matrix in the MatrixMart collection such that the Jacobi method works well.

Comments: Projects reports should be put on student's Web pages. Students are responsible for learning the Internet and MATLAB on their own. I will try to avoid discussing such issues in class.

Sample Report 1 and Report 2.