[R] sparse matrix tools

[Thaden, John J]

So far I've used only 'new("dgC.Matrix",...)' and
'new("dgT.Matrix",...)'!  I did not mean to malign
Matrix speed/functionality -- I've not tested these
yet -- nor do I quite know yet what function(s) I
needto perform on my matrices! My questions were  
hypothetical.  Thanks for the additional information
on packages I mentioned. As I become able to reply 
regarding my needs for (Matrix or UMFPACK) functions,
I shall reply offlist.  Thanks and best regards, 
-John Thaden (operating R v.2.3.0 within Windows XP)

> MM = Martin Maechler [maechler at stat.math.ethz.ch]  
> replied on Monday, July 03, 2006 2:45 AM:

>>>>> "JJTh" == Thaden, John J <ThadenJohnJ at uams.edu>
>>>>>     on Sun, 2 Jul 2006 07:38:21 -0500 writes:

    JJTh> Dear R-Help list: I'm using the Matrix library to
    JJTh> operate on 600 X ~5000 element unsymmetrical sparse
    JJTh> arrays. So far, so good, but if I find I need more
    JJTh> speed or functionality, 

MM> Can you be more specific?  
MM> Functionality:
MM>    For asymmetric matrices, in our view, the most 
MM>    important gap to fill is the LU  decomposition
MM>    and hence solve() features.
MM> 
MM> Speed: Are you sure the time your R code spends is
MM>    spent in functions from "Matrix"?  {Did you use
MM> 'Rprof()' ?}? If yes, which ones?

    JJTh> how hard would it be to
    JJTh> utilize other sparse matrix toolsets from within R,
    JJTh> say MUMPS, PARDISO or UMFPACK, that do not have
    JJTh> explicit R interfaces?  More information on these is
    JJTh> available here

    JJTh> www.cise.ufl.edu/research/sparse/umfpack/
    JJTh> www.computational.unibas.ch/cs/scicomp/software/pardiso
    JJTh> www.enseeiht.fr/lima/apo/MUMPS/ 

>From these, only the first one is open source.
Unfortunately, the PARDISO people seem to believe in money
making with scientific software -- a particular shame for since
they only work at most an hour away from me.
MUMPS is said to be "public domain", but then you only get it
after filling out a web form and only for a specific hardware.  Also,
it is about massively parallel computation, very interesting for
sparse matrices, but AFAIK not yet in our main focus.

UMFPACK is different, and even ships with the Matrix
package, because we have planned to interface to it, but haven't
yet got to that.  What parts of UMFPACK functionality would you
be interested in ?

    JJTh> and in these reviews

    JJTh> ftp://ftp.numerical.rl.ac.uk/pub/reports/ghsNAGIR20051r1.pdf
    JJTh> http://www.cise.ufl.edu/research/sparse/codes/ 

    JJTh> neither of which reviewed the R Matrix package, unfortunately.

    JJTh> Thanks, - John Thaden, Ph.D., U. Arkansas for
    JJTh> Med. Sci., Little Rock.


Regards,
Martin Maechler,  ETH Zurich

Confidentiality Notice: This e-mail message, including any a...{{dropped}}