search for: tauc

Dear R developers, The trace of the hat matrix H~(n,n) is computed as follows: tr(H) = tr(BS^-1B') = tr(S^-1B'B) := tr(X) = sum(diag(X)) with B~(n,p), S~(p,p). Since p is of the order 10^3 but S is sparse I would like to employ Taucs linear solver ( http://www.tau.ac.il/~stoledo/taucs/ ) on SX = B'B. (Further improvement by implying a looping over i=1,...,p, calling taucs_linsolve(S, X[,i], (B'B)[,i]) and saving X[i,i] only is pending.) For this purpose I compiled the C code "hattrace.c" to a shared objec...

...trix_0.6-1.tar.gz, to CRAN where it is in the src/contrib/1.9.0/Other directory. It requires some of the packages that will appear in R-1.9.0. This version marks a major redesign of the Matrix package to use S4 classes and methods and to incorporate sparse matrix manipulations using routines from TAUCS (http://www.tau.ac.il/~stoledo/taucs/), Metis (http://www-users.cs.umn.edu/~karypis/metis/) and UMFPACK (http://www.cise.ufl.edu/research/sparse/umfpack). Although this package is very much a work in progress and the documentation is, shall we say, sketchy, there are two vignettes included with t...

...sure that this is really a bug; but it is annoying, anyway). 2. Installing Matrix failed with a message saying something like "no rule for %_D.o" after compiling a lot of sources and combining them into an archive. I traced this down to a Linuxism. I found the string _D.o in subdirectory taucs of Matrix. Forcing GNU make by adding MAKE=gmake to the environment of R installed the package successfully. I consider this a real non-portability bug.