Xref: utzoo comp.lang.c:17773 sci.math.stat:690
Path: utzoo!utgpu!jarvis.csri.toronto.edu!mailrus!csd4.milw.wisc.edu!dogie.macc.wisc.edu!uwvax!tank!paideia!vevea
From: vevea@paideia.uchicago.edu (Jack L. Vevea)
Newsgroups: comp.lang.c,sci.math.stat
Subject: Re: Need matrix inversion C routine.
Keywords: C matrix invert
Message-ID: <2846@tank.uchicago.edu>
Date: 22 Apr 89 01:27:24 GMT
References: <5785@cbnews.ATT.COM> <10087@smoke.BRL.MIL>
Sender: news@tank.uchicago.edu
Reply-To: vevea@paideia.uchicago.edu (Jack L. Vevea)
Organization: University of Chicago, Dept of Education
Lines: 15

In article <10087@smoke.BRL.MIL> gwyn@brl.arpa (Doug Gwyn) writes:
>In article <5785@cbnews.ATT.COM> wkb@cbnews.ATT.COM (Wm. Keith Brummett) writes:
>>    I have need of a small, fast routine in C language that will invert
>>    matrices of order <= 4.


	While we're on this subject, I have need of a long, slow matrix
inversion routine that is at least callable from C.  I say
"long, slow" because I need something I can apply to matrices with dimension
in the hundreds, so it must be sensitive to the possibility of rounding
error, and pick the optimal pivots.  It would be nice if, in addition, it
could track its own precision in some way.  If anyone could at least point me
in the direction of a good algorithm, I would be eternally grateful.

	(Apologies if this doesn't belong in comp.lang.c)