main() { MAT *A, *Acopy; VEC *b, *x; PERM *pivot; A = m_input(MNULL); /* read in matrix from stdin */ Acopy = m_copy(A,MNULL); /* create & copy matrix */ b = v_input(VNULL); /* read in vector from stdin */ pivot = px_get(A->m); /* get pivot permutation -- size = # rows of A */ LUfactor(A,pivot); /* LU factorisation */ x = LUsolve(A,pivot,b,VNULL); /* ... and solve A.x = b; result is in x */ /* .... and print out solution with a comment */ printf("# x =\n"); v_output(x); }For information on the manual, see the bottom of this page.
FAQ.meschach: Meschach's Frequently Asked Questions (with answers)
mesch12a.announce: Announcement notice for version 1.2
mesch12b.index: Netlib index entry for new version
mesch12b.tar.Z: tar file containing complete library (compress'd)
mesch12b.tar.gz ....(gzip'd)
mesch12b.zip ....(zip'd)
BinHex4.0'd Mac setup file for Think C courtesy of Brent Boyer (boyer@jumpjibe.stanford.edu)
meschach0.shar.Z: shar file containing core routines (compress'd)
meschach0.shar.gz ....(gzip'd)
meschach1.shar.Z: shar file containing basic matrix operations (compress'd)
meschach1.shar.gz ....(gzip'd)
meschach2.shar.Z: shar file containing matrix factorisations (compress'd)
meschach2.shar.gz ....(gzip'd)
meschach3.shar.Z: shar file containing sparse and iterative methods (compress'd)
meschach3.shar.gz ....(gzip'd)
meschach4.shar.Z: shar file containing complex operations (compress'd)
meschach4.shar.gz ....(gzip'd)
oldmeschach.shar.Z: shar file containing old iterative routines (compress'd)
oldmeschach.shar.gz ....(gzip'd)
rcs_logs: Change log a la RCS
version1.1b: Directory for version 1.1b
version1.2a: Directory for version 1.2a
The Meschach 1.2 manual (Meschach: Matrix Computations in C)is published by the Centre for Mathematics and its Applications, School of Mathematical Sciences, Australian National University, Canberra, ACT 0200, Australia. To order your copy send a message to the above address by paper mail, or send an order over email to jillian.smith@anu.edu.au with your paper mail address. Cost: about A$30 (that's roughly US$22) + postage and handling. You can pay on delivery by international money order, credit/charge card etc. on delivery.
Note: There is also the tutorial chapter of the manual in ASCII form distributed with Meschach.
David Stewart, Advanced Computation Group, School of Mathematical Sciences, Australian National University
des@thrain.anu.edu.au