SCM

SCM Repository

[matrix] Annotation of /pkg/TODO
ViewVC logotype

Annotation of /pkg/TODO

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1226 - (view) (download)

1 : bates 984 ------
2 : bates 918 - Sparse matrix methods can now be based on the CHOLMOD package. We
3 :     will need to migrate from the current code to CHOLMOD-based code
4 :     using #ifdef USE_CHOLMOD. Some of the things to be done
5 :    
6 :     - Move documentation from subdirectories of src to inst/doc
7 :     - Write utilities to create a cholmod_sparse pointer from a
8 :     dgCMatrix or lgCMatrix (or zgCMatrix) object without copying and
9 :     allocating.
10 :     - Start adding simple S4 methods (rcond, %*%, +, cbind, t).
11 :    
12 : bates 344 - Report the problem in the Linux ldexp manual page. The second and
13 :     third calls in the Synopsis should be to ldexpf and ldexpl.
14 : maechler 472
15 : maechler 868 - [,] indexing: for sparse "works", but not yet for negative indices!
16 : maechler 472
17 : bates 536 - consider moving alloc3Darray from ./src/Mutils.c to
18 : maechler 538 $(RSRC)/src/base/array.c
19 : bates 536
20 : maechler 517 -------
21 :    
22 : maechler 634 - provide methods for "dspMatrix" and "dppMatrix"!
23 :    
24 : maechler 538 - implement (more) methods for supporting "packed" (symmetric / triangular)
25 : maechler 634 matrices; particularly something like pack() and unpack() [to/from our
26 :     classes from/to "numeric"] --- have already man/unpack.Rd but no method yet!
27 : maechler 538
28 : maechler 834 (have some dtr* <-> dtp*)
29 :    
30 : bates 645 ---
31 :    
32 :     - combine the C functions for multiplication by special forms and
33 :     solution wrt special forms by using a 'right' argument and a
34 :     'classed' argument.
35 : maechler 675 [done with dgeMatrix_matrix_mm(); not yet for other classes;
36 :     and for _crossprod()]
37 :    
38 :     - add more comprehensive examples / tests for Schur decomposition
39 : maechler 696
40 :     - arithmetic for sparse matrices:
41 : maechler 935 <sparseMatrix> o <same-dim-sparseMatrix>
42 : maechler 696 should return a sparse matrix for at least "+" and "*" , also %%,
43 :     and "/" and "%/%" at least when the RHS is non-zero a scalar.
44 : maechler 935 Challenge: nice implementation (``common non-0''; but Tsparse* is not uniq).
45 : maechler 696
46 : bates 712 ---
47 :    
48 : bates 881 - Create a Harwell-Boeing version of the matrix mm and the response
49 :     vector y in inst/external and remove them from the data directory.
50 :     Modify any examples that use them and modify the Comparisons vignette.
51 : maechler 935
52 :     - "Math2" , "Math", "Arith":
53 :     keep triangular and symmetric Matrices when appropriate:
54 :     particularly desirable for "Math2": round(), signif()
55 : maechler 956
56 :     - "d" <-> "l" coercion for all "[TCR]" sparse matrices is really trivial:
57 :     "d" -> "l" : drops the 'x' slot
58 :     "l" -> "d" : construct an 'x' slot of all '1'
59 :     We currently have many of these conversions explicitly, e.g.
60 :     setAs("dsTMatrix", "lsTMatrix",
61 :     function(from) new("lsTMatrix", i = from@i, j = from@j, uplo = from@uplo,
62 :     Dim = from@Dim, Dimnames = from@Dimnames))
63 :     but I would rather want to automatically construct all these coercion
64 :     methods at once by a ``method constructor'', i.e.,
65 :     for all "dsparse*" -> "lsparse*" and vice versa.
66 :     How can one do this {in a documented way} ?
67 : maechler 1087
68 : maechler 1097 - tcrossprod(x, y) : do provide methods for y != NULL
69 :     calling Lapack's DGEMM for "dense"
70 : maechler 1109 [2005-12-xx: done for dgeMatrix at least]
71 :    
72 :     - "ltTmatrix" :
73 :     2) new("ltTMatrix", diag = "U", Dim = c(2:2, 2:2)) shows lack of ltTM*_validate
74 : maechler 1201
75 :     - BUGlet: Shouldn't lose factorization here:
76 :     h6 <- Hilbert(6); chol(h6) ; str(h6) # has factor
77 :     str(H6 <- as(h6, "dspMatrix")) # has lost factor
78 :     ## and the same in a similar situation involving "dpo", "dpp"
79 : maechler 1226
80 :     - "Compare" for dgC & dgT --- want to be able to do M[M > 10] etc

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