On Wed, 14 Jan 2004, W. Beldman wrote: > I'm looking for the algorithm of the solve function (base package). As I > understand the help text neither Linpack nor Lapack's ZESV is used (in the
Well, your example is not complex. By default (LINPACK=FALSE), DGESV is used: just read the code. > simple example below) but what algorithm is? Something like Gauss Jordan > Elimination should be applicable to solve a linear system > > In particular I want to be able to translate it into i.e. Delphi, but in the > example below I'm not convinced that I understand what's happening. vD seems > rather small... > > > TIA, > > W. Beldman > > > > vD <- matrix(0.305625, 1, 1) > > vD > [,1] > [1,] 0.305625 > > v <- matrix(c(0.348125, .296875, .305625), 1, 3) > > v > [,1] [,2] [,3] > [1,] 0.348125 0.296875 0.305625 > > A <- solve(vD, v) > > A > [,1] [,2] [,3] > [1,] 1.139059 0.9713701 1 > > ______________________________________________ > [EMAIL PROTECTED] mailing list > https://www.stat.math.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html > > -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595 ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
