[R] 'solve' algorithm
Prof Brian Ripley
ripley at stats.ox.ac.uk
Wed Jan 14 10:52:33 CET 2004
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
More information about the R-help
mailing list