SCM Repository

[matrix] View of /pkg/TODO
ViewVC logotype

View of /pkg/TODO

Parent Directory Parent Directory | Revision Log Revision Log

Revision 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

- 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 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

      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).
ViewVC Help
Powered by ViewVC 1.0.0  
Thanks to:
Vienna University of Economics and Business Powered By FusionForge