R help - Aug 2007 - invert 160000x160000 matrix

Can R invert a 160000x160000 matrix with all positive numbers?  Thanks a lot!

You would need 200GB to store a since image, so if you have about 1TB
of physical memory on your computer, it might be possible.

On 8/13/07, Jiao Yang <yj316 at gwu.edu> wrote:
> Can R invert a 160000x160000 matrix with all positive numbers?  Thanks a
lot!
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

-- 
Jim Holtman
Cincinnati, OH
+1 513 646 9390

What is the problem you are trying to solve?

I don't think you can define a matrix this large in R, even if you have 
the memory. Then, of course, inverting it there may be other programs 
that have limitations.

Paul

Jiao Yang wrote:
> Can R invert a 160000x160000 matrix with all positive numbers?  Thanks a
lot!
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
===================================================================================
La version fran?aise suit le texte anglais.

------------------------------------------------------------------------------------

This email may contain privileged and/or confidential inform...{{dropped}}

A variety of tricks would need to be used to invert a matrix of this size. If
there are any other properties of the matrix that you know (symmetric, positive
definite, etc, sparse) then they could be useful too. You could partition the
matrix first, then use an in-place inverse technique for each block to
individually calculate the blocks-inverses, then combine to get the inverse of
the initial matrix. Again, if the implementation is actually solving an Ax-B = 0
system of equations, then there are specific methods for these too, like an LU
decomp, for instance. You might also want to check out some texts for this, like
the Numerical Recipes.
How's the matrix stored right now?

Best,
-Tir



Tirthankar Patnaik
India Strategy
Citigroup Investment Research
+91-22-6631 9887

For important disclosures regarding Citigroup Investment Research, including
with respect to any issuers mentioned herein, please refer to the Citigroup
Investment Research disclosure website at
https://www.citigroupgeo.com/geopublic/Disclosures/disclosure.html

> -----Original Message-----
> From: r-help-bounces at stat.math.ethz.ch 
> [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Moshe Olshansky
> Sent: Tuesday, August 14, 2007 6:40 AM
> To: Paul Gilbert; Jiao Yang
> Cc: r-help at stat.math.ethz.ch
> Subject: Re: [R] invert 160000x160000 matrix
> 
> While inverting the matrix may be a problem, if you need to 
> solve an equation A*x = b you do not need to invert A, there 
> exist iterative methods which do need A or inv(A) - all you 
> need to provide is a function that computes A*x for an 
> arbitrary vector x.
> For such a large matrix this may be slow but possible.
> 
> --- Paul Gilbert <pgilbert at bank-banque-canada.ca>
> wrote:
> 
> > I don't think you can define a matrix this large in R, even if you
> > have the memory. Then, of course, inverting it there may be other 
> > programs that have limitations.
> > 
> > Paul
> > 
> > Jiao Yang wrote:
> > > Can R invert a 160000x160000 matrix with all
> > positive numbers?  Thanks a lot!
> > > 
> > > ______________________________________________
> > > R-help at stat.math.ethz.ch mailing list 
> > > https://stat.ethz.ch/mailman/listinfo/r-help
> > > PLEASE do read the posting guide
> > http://www.R-project.org/posting-guide.html
> > > and provide commented, minimal, self-contained,
> > reproducible code.
> >
> =============================================================>
=====================> >
> > La version fran?aise suit le texte anglais.
> > 
> >
> --------------------------------------------------------------
> ----------------------
> > 
> > This email may contain privileged and/or confidential 
> > inform...{{dropped}}
> > 
> > ______________________________________________
> > R-help at stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide
> > http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
> >
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide 
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>