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
> stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
R-project.org/posting-guide.html
>
>
--
Brian D. Ripley, ripley at stats.ox.ac.uk
Professor of Applied Statistics, 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