# SCM Repository

# View of /pkg/TODO

Parent Directory | Revision Log

Revision

File size: 3046 byte(s)

**868**- (**download**) (**annotate**)*Fri Aug 19 17:01:11 2005 UTC*(14 years, 3 months ago) by*maechler*File size: 3046 byte(s)

printing of sparse logicals; "[" reorganization (part 1)

- Report the problem in the Linux ldexp manual page. The second and third calls in the Synopsis should be to ldexpf and ldexpl. - [,] indexing: for sparse "works", but not yet for negative indices! - group generics: "Arith" (partly done), but also "Compare", "Math" etc; see ?Math and the examples in ?SetGeneric - methods for rbind and cbind where they make sense (Kurt is particularly interested in dgCMatrix ones, for pkg 'arules') ****> UNFEASIBLE as long cbind is function(..., *) !!!! <******* --> follow John's proposition cbind(x, ...) via Generic cbind2(x,y)? - consider moving alloc3Darray from ./src/Mutils.c to $(RSRC)/src/base/array.c - data/ : names 'mm' and even more 'y' are ``too short''. If we really want to keep them, don't use "LazyData" (such that one needs data(*) explicitly); But MM would rather want something like ex.mm and ex.y - slot "factors" maybe should move up to "Matrix" ------- We have a (at least one) basic problem : Currently the show() method fail sometime after coercion: e.g. 'sy' show()s wrongly, even though it "str()" fine : (po <- crossprod(Matrix(0:3, 2))) # ok (ge <- as(po, "dgeMatrix")) # ok (sy <- as(po, "dsyMatrix")) # BAD str(sy) # looks fine or example(expand) # -> ex$L and ex$U look bad, however as(ex$L, "dgeMatrix") # `works' {Of course, we don't need a workaround but must understand and solve the problem} --- - provide methods for "dspMatrix" and "dppMatrix"! - implement (more) methods for supporting "packed" (symmetric / triangular) matrices; particularly something like pack() and unpack() [to/from our classes from/to "numeric"] --- have already man/unpack.Rd but no method yet! (have some dtr* <-> dtp*) - implement diagonal Matrix class "ddiMatrix" etc using constructor function Diagonal() or Diag(). - FIXME: dtpMatrix(... diag = "U") (i.e., unit-diagonal packed triangular) ----- *does* need 'x' entries for the diagonal but these are never looked at. -> change doc {or code ? - depends on what Lapack expects !} BUG: rcond() of a singular dpoMatrix gives a LaPack error instead of just 0: MM <- crossprod(M <- Matrix(c(1:4,9:6), 2,4)) ; rcond(MM) ##> Error in rcond(MM) : Lapack routine dpotrf returned error code 4 --- - combine the C functions for multiplication by special forms and solution wrt special forms by using a 'right' argument and a 'classed' argument. [done with dgeMatrix_matrix_mm(); not yet for other classes; and for _crossprod()] - add more comprehensive examples / tests for Schur decomposition - arithmetic for sparse matrices: <sparseMatrix> o { <scalar> | <same-dim-sparseMatrix> } should return a sparse matrix for at least "+" and "*" , also %%, and "/" and "%/%" at least when the RHS is non-zero a scalar. --- - in lmer.c check all instances of the use of ZtX and XtX and change them so that having a negative last element of nc means use the response only (but look for it in the right place).

root@r-forge.r-project.org | ViewVC Help |

Powered by ViewVC 1.0.0 |