From Chapter 2: modification of problem 14:
Rather than modifying existing LAPACK or CLAPACK
codes as in the book's version of question 2.14,
we will instead do this assignment in Matlab,
as described below.
Your assignment is to implement Gaussian Elimination
with Complete Pivoting (GECP) in Matlab, and evaluate
it as described here. On the class webpage you will
find four programs,
gecp.m is the routine that performs GECP,
but it only includes the calling sequence and
leading comments that describe exactly what its
inputs and outputs should be; your job is to fill
in the implementation. testgecp.m is a program that
tests gecp.m, including a comparison with Matlab's
built-in implmentation of Gaussian Elimination with
partial pivoting (which, ultimately, is the LAPACK
routine dgesv.f, which you can find at
www.netlib.org/lapack if you like). You should run
testgecp.m (which also calls testge.m and randmatgen.m,
and prints out results) and discuss whether each
of test matrix results is what you expect:
You do not need to do timing comparisons between
your version of gecp.m and Matlab's built in gepp,
which will be much faster.