SCM

SCM Repository

[matrix] Diff of /pkg/Matrix/TODO
ViewVC logotype

Diff of /pkg/Matrix/TODO

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2687, Sat Aug 6 13:24:21 2011 UTC revision 2688, Sat Aug 6 15:33:33 2011 UTC
# Line 97  Line 97 
97    be for triangularMatrix, and probably also "dense*general*Matrix" since the    be for triangularMatrix, and probably also "dense*general*Matrix" since the
98    above currently goes via "matrix" and back instead of using the 'x' slot    above currently goes via "matrix" and back instead of using the 'x' slot
99    directly; in particular, the triangular* "class property" is lost!    directly; in particular, the triangular* "class property" is lost!
100      [current ??]
101    
102    - The "[<-" now uses  src/t_Csparse_subassign.c and no longer explodes
103      memory.  *However* it is still too slow when the replacment region is large.
104    
   Note that 'diag(M[,-1]) <- val' is deadly slow (*) for large sparse M,  
   but that's because of the 2nd line assignment in the equivalent  
     tmpM <- `diag<-`(M[,-1], val)  
     M[,-1] <- tmpM  
   (*): gives *error* {about negative integer} when  
        prod(dim(M)) > .Machine$integer.max, e.g. for square (n x n) M when  
        n >= 46341 == ceiling(2^15.5)  
   This is "the same" as Ashley Ford's report (25 Feb 2010),  
   MM @ ~/R/MM/Pkg-ex/Matrix/nsp-2col-index-bug.R  
   --> This no longer happens for CsparseMat.. (but still for Tsparse*)  
105    
106  - examples for solve( Cholesky(.), b, system = c("A", "LDLt"....))  - examples for solve( Cholesky(.), b, system = c("A", "LDLt"....))
107    probably rather in man/CHMfactor-class.Rd than man/Cholesky.Rd    probably rather in man/CHMfactor-class.Rd than man/Cholesky.Rd
# Line 221  Line 215 
215  - (systematically check that LAPACK-calling functions check for  - (systematically check that LAPACK-calling functions check for
216     0-dimensional input themselves; LAPACK gives an integer error code)     0-dimensional input themselves; LAPACK gives an integer error code)
217    
218  - the  f[,5762] <- thisCol # now ...  - the  f[,5762] <- thisCol  now go via Csparse_subassign() call ...
219    line in tests/indexing.R uses very large objects unnecessarily;    [ in tests/indexing.R ].
220    Improve  replTmat()  in R/Tsparse.R, making use of new "abIndex"    Still would be nice to be able to use abIndex (see replTmat in R/Tsparse.R)
   vectors.  
221    
222  - {IS THIS CURRENT?}  - {IS THIS CURRENT?}
223    Sept. 2009:    Sept. 2009:
# Line 252  Line 245 
245    correct and equivalent to  src/CHOLMOD/Include/cholmod_core.h  and  siblings ???    correct and equivalent to  src/CHOLMOD/Include/cholmod_core.h  and  siblings ???
246    {currently need to do this manually (Emacs M-x compare-windows) for the   typedefs}    {currently need to do this manually (Emacs M-x compare-windows) for the   typedefs}
247    
248  - finalize and activate the new *unused* code in src/sparseVector.c  - finalize and activate the new *unused* code in src/t_sparseVector.c

Legend:
Removed from v.2687  
changed lines
  Added in v.2688

root@r-forge.r-project.org
ViewVC Help
Powered by ViewVC 1.0.0  
Thanks to:
Vienna University of Economics and Business Powered By FusionForge