[R] 'solve' algorithm

[Prof Brian Ripley]

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
> 
> ______________________________________________
> R-help at stat.math.ethz.ch 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,                  ripley at stats.ox.ac.uk
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