diff -Nru rmatrix-1.6-1.1/ChangeLog rmatrix-1.6-5/ChangeLog --- rmatrix-1.6-1.1/ChangeLog 2015-04-29 16:14:56.000000000 +0000 +++ rmatrix-1.6-5/ChangeLog 1970-01-01 00:00:00.000000000 +0000 @@ -1,5087 +0,0 @@ -2015-04-01 Martin Maechler - - * inst/NEWS.Rd: Moving ChangeLog to new (markup) NEWS file - -2015-03-16 Martin Maechler - - * R/Auxiliaries.R (anyDuplicatedT): renamed from is_duplicatedT(), - and exported. uniqTsparse(): exported too. - -2014-06-15 Martin Maechler - - * DESCRIPTION (Version): 1.1-4, released to CRAN on 2014-06-14 - * src/dsyMatrix.c (dsyMatrix_matrix_mm): fix crossprod(, ) bug - -2014-04-26 Martin Maechler - - * new rsparsematrix() - -2014-03-30 Martin Maechler - - * DESCRIPTION (Version): 1.1-3, released to CRAN on 2014-03-30 - -2014-03-12 Martin Maechler - - * R/dgTMatrix.R (image): fix bug in default ylim computation. - -2014-01-28 Martin Maechler - - * R/products.R: matrix products overhauled; should work with sparseVectors; - speedup of crossprod(v, ), thanks to nudge by Niels Richard Hansen. - * man/matrix-products.Rd: all matrix products documented in one file. - * tests/matprod.R: more extensive testing - -2014-01-20 Martin Maechler - - * DESCRIPTION (Version): 1.1-2, released to CRAN on 2014-01-28 - - * NAMESPACE: export fast power-user coercion utilities - .dsy2mat(), .dxC2mat(), .T2Cmat(), ..2dge(). - -2013-12-23 Martin Maechler - - * R/dgTMatrix.R (image): (xlim, ylim) get a slightly changed - default, plus ylim := sort(ylim, "decreasing"). - This is strictly not back-compatible but should never harm. - -2013-09-26 Martin Maechler - - * R/spModels.R (fac2sparse, fac2Sparse): newly exported; plus - 'giveCsparse' option. - -2013-09-16 Martin Maechler - - * src/scripts/0get-SuiteSparse.sh: new download script, for - * src/CHOLMOD/*, src/AMD/*, ...: getting SuiteSparse version 4.2.1 - - * R/zzz.R (.SuiteSparse_version): new function - -2013-09-13 Martin Maechler - - * R/dsCMatrix.R (solve.dsC.*): finally fix the long-lasting - undetected solve() bug (only in case Cholmod fails) [r2908]. - - * DESCRIPTION (Version): 1.0-15, CRAN-*non*-released 2013-09-26 - -2013-09-12 Martin Maechler - - * DESCRIPTION (Version): 1.0-14, CRAN-released 2013-09-12 - * R/dgCMatrix.R: "stop gap fix" for R 3.0.2: partly revert - solve(,*) changes in 1.0-13. - -2013-08-27 Fabian Scheipl - - * man/indMatrix-class.Rd: new "indMatrix" class, a natural - superclass of "pMatrix". Many methods moved from "pMatrix" to - "indMatrix". - -2013-05-09 Martin Maechler - - * DESCRIPTION (Version): 1.0-13, CRAN-released 2013-09-10 - - * R/KhatriRao.R: Efficient KhatriRao() by Michael Cysouw - -2013-03-26 Martin Maechler - - * DESCRIPTION (Version): 1.0-12, CRAN-released 2013-03-26 - -2012-11-10 Martin Maechler - - * DESCRIPTION (Version): 1.0-11, CRAN-released 2013-02-02 - - * R/SparseM-conv.R: as(, "dgCMatrix") now works again - * tests/other-pkgs.R: test that. - - * src/Mutils.c: do *not* use '#if R_VERSION < ..' so this also - *runs* in older R when installed in R >= 2.15.2. - -2012-10-15 Martin Maechler - - * src/Mutils.c (Mmatrix): new, to be used by .External() - in order to replace .Internal(matrix(....)) - -2012-10-05 Martin Maechler - - * R/diagMatrix.R (.sparseDiagonal): new 'unitri' argument; - more flexibility; new: solve(, ) - -2012-09-10 Martin Maechler - - * DESCRIPTION (Version): 1.0-10, CRAN-released: 2012-10-16, r2845 - -2012-09-01 Martin Maechler - - * R/sparseVector.R (sparseVector): new constructor - - * inst/test-tools-Matrix.R (rspMat): smarter; also useful for large - dimensions. - -2012-07-23 Martin Maechler - - * tests/group-methods.R: now do much more testing, notably of - pairs of matrices ... however only when 'doExtras' is true, which - it is not by default, e.g., on CRAN. - -2012-07-21 Martin Maechler - - * R/Ops.R, R/diagMatrix.R: many fixes, notably for rare operations - that were not triggered before. - - * R/Auxiliaries.R (allTrueMat): new utility. - -2012-07-20 Martin Maechler - - * R/dsparseMatrix.R, R/sparseVector.R, ...: Newly defined - is.finite() and is.infinite() methods for all our *Matrix and *sparseVector. - -2012-07-14 Martin Maechler - - * R/d??Matrix.R (diag<-): many "diag<-" methods, which preserve - symmetricity, triangularity (even uni-triangularity sometimes); - Partly by also making A[cbind(i,i)] preserve such properties. - - * src/Mutils.c (SET_packed_setDiag) {and similar}: - *_setDiag() C functions implementing "diag<-" R methods. - -2012-06-30 Martin Maechler - - * R/Ops.R (Ops.x.x): now, dense symmetric and triangular matrices - are preserved for many arithmetic and logic binary operators. - - * src/ldense.c (lsyMatrix_as_lspMatrix, ..): more coercions keep - the dimnames(). - - * R/symmetricMatrix.R (pack, unpack): new pack() {"inverse" of - unpack(), including pack()}; new unpack() methods. - -2012-06-20 Douglas Bates - - * src/scripts/DEPS.mkf, ... - * src/*.c: Update to version 4.0.0 of SuiteSparse - -2012-06-19 Martin Maechler - - * DESCRIPTION (Version): 1.0-8, CRAN-released: 2012-06-20, r2789 - - * R/CHMfactor.R (update): *DO* allow non-symmetric parent. - * man/CHMfactor-class.Rd: be more clear about the two ways. - * tests/factorizing.R: more update() testing - -2012-06-12 Martin Maechler - - * tests/matprod.R (chkDnProd): new testing function - -2012-06-08 Martin Maechler - - * R/CHMfactor.R (update): now *warn* when parent is not formally - symmetric and coerce it. - -2012-06-05 Martin Maechler - - * R/Auxiliaries.R (chk.s): "check dots" - new utility -- possibly - for base R ? - -2012-04-16 Martin Maechler - - * R/sparseMatrix.R (sparseMatrix): now also works for - 'symmetric=TRUE' and lower triangular indices. - -2012-04-15 Martin Maechler - - * R/CHMfactor.R (updown): new generic and methods, - * man/updown.Rd: provided by Nicholas Nagle. - -2012-03-30 Martin Maechler - - * DESCRIPTION (Version): 1.0-7, CRAN-released: for a few days only. - -2012-03-16 Martin Maechler - - * DESCRIPTION (Version): 1.0-6, CRAN-released: 2012-03-30, r2775 - * DESCRIPTION (Depends): R >= 2.15.0 - -2012-03-15 Martin Maechler - - * R/spModels.R (sparseInt.r): recursion free (which does not help much). - -2012-03-05 Martin Maechler - - * src/dtCMatrix.c (dtCMatrix_sparse_solve): no longer use Alloca() - here. - * tests/factorizing.R (checkSchur): check against segfault example. - - * R/Matrix.R (chol2inv()) new method. - -2012-03-01 Martin Maechler - - * R/spModels.R (sparse.model.matrix, model.spmatrix): add 'verbose' - argument in order to show to the user what's going on. - * man/sparse.model.matrix.Rd: ditto - -2012-02-27 Martin Maechler - - * R/Ops.R (A.M.n, A.n.M): o now correct, newly - via sparseVector. - -2012-02-25 Martin Maechler - - * DESCRIPTION (Version): 1.0-5, CRAN-released: 2012-03-15, r2773 - * src/chm_common.c (chm_factor_to_SEXP): in case of failure, - * src/dsCMatrix.c (internal_chm_factor): ensure memory cleanup; - memory leak reported by Kasper Kristensen at dtu.dk. - -2012-02-17 Martin Maechler - - * DESCRIPTION (Version): 1.0-4, CRAN-released: 2012-02-21, r2765 - * R/Ops.R: Fix "-" method for diagonalMatrix types. - -2012-01-12 Martin Maechler - - * DESCRIPTION (Suggests): MASS, as we use a data set in a vignette - -2011-12-09 Martin Maechler - - * DESCRIPTION (Version): 1.0-3, CRAN-released: 2012-01-12, r2749 - * R/diagMatrix.R (.bdiag): now works correctly when all blocks are "lMatrix" - -2011-11-02 Martin Maechler - - * R/zzz.R (.onLoad): when R version >= 2.15.0, no longer - need to assignInNamespace( ns = "base"), methods provides S4 generic. - -2011-10-30 Martin Maechler - - * DESCRIPTION (Version): 1.0-2, CRAN-released: 2011-11-19, r2739 - - * tests/other-pkgs.R: print more, don't run SparseM on Solaris for now. - * tests/Simple.R: encoding warning should not be promoted to error. - -2011-10-22 Martin Maechler - - * R/rankMatrix.R (rankMatrix): 'method = "qrLINPACK"' now also - works for sparse matrices, but using \code{sparseQR()}. - - * man/sparseQR-class.Rd: document options "Matrix.quiet" (old) and - new "Matrix.quiet.qr.R" for suppressiong the permutation warning. - * R/sparseQR.R: - -2011-10-17 Douglas Bates - - * src/Csparse.c (Csparse_submatrix): plugging memory leak - -2011-10-08 Martin Maechler - - * R/bind2.R (cbind2Sparse, rbind2Sparse): auxiliaries, - used, also in new methods for . - * tests/bind.R: testing some of these. - - * R/Matrix.R, man/Matrix.Rd: optional argument 'doDiag'. - - * DESCRIPTION (Version): 1.0-1, CRAN-released: 2011-10-18, r2732 - -2011-09-27 Martin Maechler - - * DESCRIPTION (Version): 1.0-0 -- to be released for R 2.14.0 - -2011-09-16 Martin Maechler - - * R/dsCMatrix.R (solve.dsC.mat): new utility which calls lu() if - CHOLMOD'S Cholesky() errors (when matrix is not pos.def.). - -2011-09-15 Martin Maechler - - * R/bandSparse.R (bandSparse): and - * R/sparseMatrix.R (sparseMatrix): add 'giveCsparse = TRUE' - argument and allow returning Tsparse*, useful e.g. when used in - bdiag(). - -2011-08-17 Martin Maechler - - * NAMESPACE: export diagN2U() & diagU2N(). They were "missing" for - at least one user (GG). - * man/diagU2N.Rd: docu + example. - - * DESCRIPTION (Version): 0.9996875-4 (*not* yet released) - -2011-08-12 Martin Maechler - - * DESCRIPTION (Version): 0.9996875-3, for CRAN - -2011-08-05 Martin Maechler - - * R/sparseVector.R (head): method; used in a few cases, eliminating - two FIXMEs. - - * DESCRIPTION (Version): 0.9996875-1 - - * R/ngTMatrix.R: stop() instead of warning() when NA's are coerced. - * R/Csparse.R, R/Tsparse.R - * src/t_Csparse_subassign.c, src/Csparse.c: [..] <- val - now works via .Call(*Csparse_subassign, ...) and no longer suffers - from unnecessary memory-blowup. - -2011-07-29 Martin Maechler - - * R/dgTMatrix.R (image): add 'useRaster = FALSE' argument, - providing the possibility of using raster (instead of rectangle - drawing vector) images. - -2011-07-27 Martin Maechler - - * R/nearPD.R: allow 'ensureSymmetry' argument for speedup. - -2011-06-10 Martin Maechler - - * R/diagMatrix.R (Cspdiagprod, diagCspprod): fixup for symmetric - sparse, and non constant-diagonal. - -2011-05-20 Martin Maechler - - * R/dsCMatrix.R (determinant(): fix for Matrix(0, 1) case. - -2011-05-18 Martin Maechler - - * R/sparseMatrix.R (sparseMatrix): add 'symmetric' argument. - -2011-04-04 Martin Maechler - - * src/Csparse.c (Csparse_subassign): unfinished prototype - - * src/....: Finally no longer hack "UFlong := int", but rather - * src/UFconfig/UFconfig.h: use standard CHOLMOD headers - - * DESCRIPTION (Version): 0.9996875-0 - -2011-03-31 Martin Maechler - - * DESCRIPTION (Version): 0.999375-49 - -2011-03-30 Douglas Bates - - * Matrix/src/chm_common.c: [r2658] Install symbols first - - {preventing seg.fault under -gct} - -2011-03-17 Martin Maechler - - * DESCRIPTION (Version): 0.999375-48 - *only* difference to CRAN released *-47, is the work around Sweave bug - in inst/doc/sparseModels.Rnw. - -2011-02-23 Martin Maechler - - * src/factorizations.c (LU_expand): now also works for non-square (m x n) - * tests/factorizing.R: testing that - -2011-02-22 Martin Maechler - - * R/Auxiliaries.R (t_geMatrix): drop 'factors', as they can be - wrong. - -2011-02-18 Martin Maechler - - * R/Tsparse.R (replTmat): fix bug for M[i,j] <- v, - when j had duplicated entries. - * tests/indexing.R (chkAssign): new function; testing the above. - -2011-02-17 Martin Maechler - - * R/AllClass.R, R/sparseVector.R, man/sparseVector-class.Rd: now - require explicitly that i-slot must be sorted for sparseVectors. - -2011-02-16 Martin Maechler - - * R/sparseMatrix.R (formatSparseM): align="right" accidentally did - not use zero.print at all. - print/format sparse Matrix: fix align="right" and improve docu - -2011-02-17 Douglas Bates - - * DESCRIPTION: Remove Encoding: directive. - -2011-02-10 Martin Maechler - - * inst/doc/sparseModels.Rnw: use png for large graphics (suggestion from - Brian) - -2011-02-05 Martin Maechler - - * man/symmpart.Rd: update, thanks to Spencer Graves' prompting. - -2011-01-07 Martin Maechler - - * R/CHMfactor.R (determinant()): no longer warn about - incompatible change of 2009-09-01. - -2011-01-04 Martin Maechler - - * R/nearPD.R (nearPD): better error message when all eigenvalues - are (at least close to) negative. - -2010-12-18 Martin Maechler - - * DESCRIPTION (Version): 0.999375-47, CRAN-released: 2011-02-23, r2653 - - * R/spModels.R, NAMESPACE: remove model.Matrix(); - we had deprecated it for about four months now. - -2010-12-12 Martin Maechler - - * R/eigen.R,...: use full argument names; - * R/*.R: get rid of more from - checkUsagePackage("Matrix", suppressPartialMatchArgs = FALSE) - -2010-12-11 Martin Maechler - - * DESCRIPTION (Version): 0.999375-46, CRAN-released: 2010-12-14, r2633 - - * R/products.R: dimension fixing in some Matrix o vector [t]crossprod()s. - - * src/Csparse.c (nz2Csparse, nz_pattern_to_Csparse): new utilities, - callable from C and R. - (Csparse_dense_prod): check pattern matrix and coerce to "d..". - * tests/Simple.R: testing %*% - - * R/ngCMatrix.R, R/nsCMatrix.R, .. : use the new fast coercions. - -2010-10-08 Martin Maechler - - * DESCRIPTION (Version): 0.999375-45, CRAN-released: 2010-11-10, r2624 - - * R/sparseMatrix.R (graph.wgtMatrix): add 'graph::' in a few - places; as 'graph' package is not imported and may well be loaded - only. - -2010-09-09 Martin Maechler - - * R/sparseMatrix.R (setAs): graph |-> Matrix: via CsparseMatrix - -2010-08-21 Martin Maechler - - * R/spModels.R (sparse.model.matrix): argument 'drop.unused.levels = FALSE' - NB: was *true* implicitly, before. Compatibility with model.matrix(). - - * R/sparseMatrix.R (formatSpMatrix, formatSparseM): factored out of - printSpMatrix(); export as potentially useful for standard - matrices, even. - -2010-08-11 Martin Maechler - - * R/eigen.R (Schur): correct setMethod() such that - Schur() works. - -2010-08-10 Martin Maechler - - * R/diagMatrix.R, man/bdiag.Rd, NAMESPACE: export .bdiag() as well. - -2010-08-09 Martin Maechler - - * DESCRIPTION (Version): 0.999375-44, CRAN-released: 2010-09-11, r2618 - - * R/diagMatrix.R (diagCspprod, Cspdiagprod): drop (possibly wrong) @factors - * R/Ops.R (.Arith.CM.atom, .Arith.atom.CM, A.M.n, A.n.M): ditto - - * tests/factorizing.R: check some of the above. - - -2010-08-04 Douglas Bates - - * R/spModels.R (fac2sparse): Fix name resolution problem - (R-SIG-Mixed-Models post by Florent Duyme). - -2010-07-25 Martin Maechler - - * DESCRIPTION (Depends): require R >= 2.10.0 --> can clean up - * R/spModels.R: prepare to move most parts to new package - MatrixModels - -2010-07-23 Martin Maechler - - * R/Auxiliaries.R (prMatrix): add " (unitriangular)" as we already - have for sparse matrices. - -2010-07-22 Martin Maechler - - * R/Auxiliaries.R (.diagU2N): implement for "dtpMatrix" = old 'FIXME'; - (.dense.diagU2N): new utility, called from .diagU2N() - -2010-07-19 Martin Maechler - - * src/dtrMatrix.c (dtrMatrix_dtrMatrix_mm): new for tri %*% tri - * R/products.R (%*%): ditto - * tests/matprod.R: test it - -2010-07-16 Martin Maechler - - * R/spModels.R (do.defaults): add 'nonMatched.action' with default - ensuring that typos are caught. - -2010-07-16 Douglas Bates - - * R/spModels.R (do.defaults): utility function; TODO: move to R - -2010-07-16 Martin Maechler - - * DESCRIPTION (Version): 0.999375-43, CRAN-released: 2010-08-05, r 2599 - - * R/AllClass.R (Model): as mother class (of "glpModel") - - * R/spModels.R (IRLS): more options() - (updateModel): update() - -2010-07-13 Martin Maechler - - * DESCRIPTION (Version): 0.999375-42, CRAN-released: 2010-07-15, r 2566 - - * R/spModels.R (glm4, IRLS): glm4 [was 'glm1']; tweaks. - -2010-07-12 Martin Maechler - - * NAMESPACE: rename, export and - * man/glpModel-class.Rd: document Doug's new "glpModel" class. - -2010-07-08 Douglas Bates - - * R/AllClass.R: new "lpMod" class (-> later =: "glpModel"), and - working function: - * R/spModels.R (glm1): using linear pred.Model class, and - Bates-Watts convergence criterion. - -2010-07-06 Martin Maechler - - * R/lMatrix.R (whichDense): use arrayInd() - * R/zzz.R (arrayInd): provide for older R versions - -2010-07-05 Martin Maechler - - * src/chm_common.c (chm_triplet_to_SEXP): deal more carefully with - NAs, needed e.g., on Solaris; thanks to Ruth and Brian Ripley. - - * R/Ops.R (Compare ): fix bug uncovered by "alongside" the - above. - -2010-07-02 Martin Maechler - - * R/sparseMatrix.R (x[] <- 0): fix shortcut code. - * tests/indexing.R: and test - - * R/nearPD.R (nearPD): 'doDykstra = TRUE' and *do* use Dykstra's - correction which was *not* used in Jens' code; thanks to Bernhard - Spangl for a report - -2010-06-26 Martin Maechler - - * R/Matrix.R: fix mean() method; add sparse one, remaining sparse - at least for trim = 0. - -2010-06-08 Martin Maechler - - * DESCRIPTION (Version): 0.999375-41, CRAN-released: 2010-07-03, r 2555 - - * R/sparseVector.R (spV2M): enable sparseVector -> (sparse)Matrix - as(*,.) coercion. - -2010-06-07 Martin Maechler - - * R/Tsparse.R (.TM.repl.i.mat): renamed from .TM.repl.i.2col(). - * R/Matrix.R (.repl.i.lSMat): implement logical sparse - sub-assignment: M[] <- v; ditto for dense & nsparse. - - * R/Csparse.R (.CM.repl.i.lSMat, ..): direct logical sparse - sub-assignment for "Csparse". - -2010-06-04 Martin Maechler - - * R/sparseMatrix.R (sparseMatrix): re-enable 'dimnames' argument. - -2010-06-03 Martin Maechler - - * R/spModels.R (model.Matrix): tweak for NULL contrasts in dense case. - -2010-06-02 Martin Maechler - - * tests/spModel.matrix.R (Sparse.model.matrix): adapt to the fact, that - sparse.model.matrix() returns not just a dgCMatrix. - -2010-05-29 Martin Maechler - - * DESCRIPTION (Version): 0.999375-40, CRAN-released: 2010-06-04, r 2546 - - * R/AllClass.R: new classes "ModelMatrix", "sparseModelMatrix", etc. - - * R/spModels.R (sparse.model.matrix): now return - "dsparseModelMatrix" object, notably with 'assign' slot. - * R/spModels.R (model.spmatrix): faster, using lower level - cbind2/rbind2 directly. - - * R/spModels.R (model.Matrix): new function, - returning "ddenseModelMatrix". - - * NAMESPACE: export new classes. - -2010-05-18 Martin Maechler - - * src/Csparse.c (Csparse_horzcat, Csparse_vertcat): ensure that - rBind()/cBind() i.e., rbind2()/cbind2() return logical sparse - matrices when the components are. - * tests/bind.R: test - -2010-05-15 Martin Maechler - - * R/sparseMatrix.R: A[] <- v ; differentiate dense & sparse - * R/pMatrix.R: disallow [] <- v more consequently - - * tests/indexing.R: test above - - -2010-05-08 Martin Maechler - - * R/spModels.R (model.spmatrix): deal with "AsIs" components - * tests/spModel.matrix.R: test that - -2010-05-01 Martin Maechler - - * R/condest.R (onenormest, condest): allow to estimate condition - number for large sparse matrices. - * condest.Rd: docu - - * R/pMatrix.R (.inv.perm): utility; add [t]crossprod() methods - - * man/sparseLU-class.Rd: A = P'LUQ; add examples, with "identities" - - * R/Auxiliaries.R (mmultCheck): new arg. 'kind' - -2010-04-28 Martin Maechler - - * DESCRIPTION (Version): 0.999375-39, CRAN-released: 2010-05-19, r 2540 - - * R/spModels.R (fac2sparse): using names(formals(new))[[1]] to - adapt to a future change in new()'s first argument *name*. - -2010-03-31 Martin Maechler - - * R/spModels.R (lm.fit.sparse): update, allowing weights; also - return residuals, notably for "cholesky" case. - - * man/lm.fit.sparse.Rd: examples; comparing with dense case. - - * src/dgCMatrix.c (dgCMatrix_cholsol): comments; also compute - residuals. - -2010-03-30 Martin Maechler - - * R/spModels.R (sparse.model.matrix, model.spmatrix): border case - '~ 1' should also work. Add 'transpose = FALSE' argument. - - * tests/spModel.matrix.R: test that. - -2010-03-27 Martin Maechler - - * R/sparseMatrix.R (printSpMatrix): ensure returning original argument - -2010-03-26 Martin Maechler - - * R/sparseVector.R (coercion from TsparseMatrix): diagU2N() when - needed. - - * inst/test-tools.R (checkMatrix): explicit which() test for - "l" and "nMatrix". New sparseVector (coercion "and back") check. - -2010-03-25 Martin Maechler - - * R/lMatrix.R (which): define methods for which(<[ln]Matrix>). - - * inst/test-tools.R (Q.eq.symmpart): new utility, now called in - checkMatrix(). - - * R/nearPD.R: use symmpart() for non-symmetric x - * man/nearPD.Rd: improve title - -2010-03-24 Martin Maechler - - * R/colSums.R (.diag.Mean): define methods for "diagonalMatrix" - - * src/Mutils.c (m_encodeInd, do_ii_FILL): - coerce ij to integer if necessary; check that ij are within "dim[]" - values. Parametrize do_ii_FILL() to be used in m_encodeInd2() as well: - * src/Mutils.c (m_encodeInd2): also check bounds (if desired). - - * tests/indexing.R: test the above. - -2010-03-19 Martin Maechler - - * src/dgeMatrix.c (dgeMatrix_solve): compute the recip.cond.number - and also bail out for *computational* singularity {as "base R" - does}, from (code) suggestion by Daniel Sabanés Bové. - - * tests/dg_Matrix.R: "test" the above. - -2010-03-01 Martin Maechler - - * man/rep2abI.Rd: rep2abI() utility is exported now. - - * R/Csparse.R (subCsp_cols, subCsp_rows, subCsp_ij): dimnames() <- - fix for character subsetting. - - * tests/indexing.R: testing it. - -2010-02-26 Martin Maechler - - * R/spModels.R (model.spmatrix): warn and coerce to sparseMatrix if - result would end up dense (e.g., in case of *no* factors in formula). - - * tests/spModel.matrix.R: test the above. - -2010-02-12 Martin Maechler - - * R/dtrMatrix.R: add solve(, ): e.g., for - solve(, ) in lme4. - -2010-02-09 Martin Maechler - - * DESCRIPTION (Version): 0.999375-38, CRAN-released: 2010-03-31, r 2529 - - * NAMESPACE, R/AllGeneric.R, R/zzz.R: change det() into a regularly exported - function (masking base::det) instead of load-time hack. - * man/Matrix-class.Rd: \alias, docu - -2010-02-05 Martin Maechler - - * DESCRIPTION (Version): 0.999375-37, CRAN-released: 2010-02-05 - - * inst/test-tools.R (Qidentical.DN): new - (Qidentical): all Qidentical.DN() - - * R/Csparse.R (subCsp_ij, subCsp_cols, subCsp_rows): use - CHOLMOD's submatrix C code, instead of matrix multiplication; - now *do* keep dimnames of result, wherever classical matrix - subsetting does. - -2010-02-04 Martin Maechler - - * DESCRIPTION (Version): 0.999375-36, CRAN-released: 2010-02-04 - - * R/Csparse.R (subCsp_ij): Fix [0,0] bug - -2010-02-03 Martin Maechler - - * R/Tsparse.R (.TM.repl.i.2col): [ ] <- FALSE fix - * tests/indexing.R, *.Rout.save: test that - -2010-01-28 Martin Maechler - - * src/Csparse.c (Csparse_crossprod): PROTECT() Tsparse_diagU2N() - result, from a suggestion by Andrew Runnalls. - -2010-01-22 Martin Maechler - - * R/SparseM-conv.R (setAs(., "matrix.csc")): fix typo in method. - -2010-01-20 Martin Maechler - - * R/AllGeneric.R: nnzero() is now generic, - * R/nnzero.R: newly containing all nnzero() methods. - - * R/zzz.R (det): assign base::det such that it uses S4 generic kronecker. - -2010-01-18 Martin Maechler - - * R/spModels.R (contr.poly): [the back-compatible-only version]: - do not use a default for 'scores'; rather rely on stats::contr.poly. - - * tests/spModel.matrix.R: test that case - -2009-12-28 Douglas Bates - - * DESCRIPTION (Version): 0.999375-35, CRAN-released: 2010-02-03 - - * src/init.c, inst/include/Matrix_stubs.c: cholmod_band_inplace() - exported. - -2009-12-23 Martin Maechler - - * tests/indexing.R: slightly adapt to the very slight [] changes. - * inst/test-tools.R (Q.C.identical): + checkClass - - * R/Tsparse.R ([-methods): for now go via "Csparse" and drop - all the sophisticated code dealing with the many cases. - - * R/Csparse.R (subCsp_cols, etc): faster [i,j] via matrix - multiplication, thanks to suggestions by Greg Jorstad. - - * R/Auxiliaries.R (paste0): more use of paste0() - -2009-12-22 Martin Maechler - - * R/diagMatrix.R (.sparseDiagonal): made more general, allowing - to directly build the equivalent of Diagonal(n)[, j] - - * man/Diagonal.Rd: document .sparseDiagonal() as well. - -2009-12-21 Martin Maechler - - * R/AllClass.R: abIndex@x and rleDiff@first are now "numLike", - such that also logical can be converted to "abIndex". - - * R/abIndex.R (.diff): new utility, used in num2abI() - -2009-12-19 Martin Maechler - - * src/abIndex.c: include new - * src/t_Matrix_rle.c (Matrix_RLE_): is template for - Matrix_rle_i() and Matrix_rle_d(); now obeys a 'force' argument. - - * R/abIndex.R: implement methods, at least with scalars. - (all.equal.abI): add also all.equal() methods. - - * tests/abIndex-tsts.R: testing , using all.equal(). - - * R/AllClass.R: classUnion "numLike" := {"numeric", "logical"} - -2009-12-18 Martin Maechler - - * src/abIndex.c (Matrix_int_rle): UNPROTECT() needed in trivial case. - - * R/abIndex.R (abIseq1, abIseq): new functions for building - "abIndex" vectors. - - * tests/abIndex-tsts.R (tst.c.abI): test new c("") method. - - * DESCRIPTION (Version): 0.999375-34, CRAN-released: -never- - - * R/Ops.R: use prod(d) instead of d[1]*d[2], as the latter may - integer overflow; fixes o - - * tests/Simple.R: test that. - -2009-12-11 Martin Maechler - - * R/sparseVector.R (TsparseM* -> sparseV): symmetricMatrix needs to - be expanded. This fixes a bug reported by Yu-Sung Su. - * tests/indexing.R: testing the fix. - - * inst/test-tools.R (all.equalX): new util - -2009-12-09 Martin Maechler - - * R/Ops.R (A.n.M, A.n.M): o : remain - sparse also when *majority* (instead of all) of 0 o v is 0. - - * tests/group-methods.R: test one such case. - -2009-12-06 Martin Maechler - - * DESCRIPTION (Version): 0.999375-33, CRAN-released: 2009-12-11 - - * R/Ops.R (Compare): fix case with NA x. - - * R/not.R: fix "typo" in ! - - * R/Ops.R (Ops.spV.spV): fix thinko - -2009-12-05 Martin Maechler - - * R/sparseVector.R: setAs(nsparseV*, lsparseV*) etc - - * R/Ops.R (Ops.spM.spV, Ops.spV.spM): sparseVec. o sparseMat. - -2009-11-20 Martin Maechler - - * R/Ops.R (Ops.spV.spV): enable sparseVector operations in more cases. - - * R/sparseVector.R (is.na): methods defined. - * R/sparseVector.R (intIv): also accept "sparseVector"s - - * tests/Simple.R: check the above - -2009-11-19 Martin Maechler - - * R/sparseVector.R (newSpV, newSpVec): new utility, dropping 0 "on the fly". - * R/sparseVector.R (atomic -> sparse*): fix for NA case. - - * R/Ops.R (): using newSpVec() - - * R/not.R: fix thinko in ! - -2009-11-17 Martin Maechler - - * tests/other-pkgs.R: detach(*, unload) Rgraphviz too - -2009-11-14 Martin Maechler - - * R/AllClass.R: "abIndex" (and "rleDiff") class - * R/abIndex.R: with some methods; commit these finally, even if - it's mostly unfinished. - * src/abIndex.[ch]: new: currently only for .Call(Matrix_int_rle,*) - * tests/abIndex-tsts.R: basic consistency checks for "abIndex". - - * R/diagMatrix.R (diagOdiag): "exploding" Matrix.msg() only level 2; - * tests/indexing.Rout.save: update - -2009-11-11 Martin Maechler - - * DESCRIPTION (Version): 0.999375-32, CRAN-released: 2009-11-20 - - * src/Csparse.c (Csparse_Csparse_prod, Csparse_Csparse_crossprod): - PROTECT(.) the dimnames; thanks to Kaspar Kristensen - -2009-10-24 Martin Maechler - - * R/Ops.R (Logic.lCMat): to be used for lsC* and ltC* as well, - effectively replacing previous suboptimal methods. - - * src/chm_common.c (chm2Ralloc): Fix unidiagonal ntC segfault: - assign 'x' only when non-pattern. - * src/chm_common.c (as_cholmod_triplet): reallocate now in a way - that works; fix documentation about return value in diagU2N case; - ditto for - * src/chm_common.c (as_cholmod_sparse): - - - * R/sparseMatrix.R (printSpMatrix): add 'cld' argument, typically - passed from printSpMatrix2; and indicate "unit-diagonal" - -2009-10-22 Martin Maechler - - * R/lsparseMatrix.R (C2l): fix for case with NA. - - * R/Csparse.R (replCmat): drop "stale" cached @factors - factorizations after sub-assignments. - * R/Tsparse.R (replTmat, .TM.repl.i.2col): ditto - -2009-10-19 Martin Maechler - - * src/dgCMatrix.c (dgCMatrix_LU): new boolean argument - 'error_on_sing' to allow no error in case of singularity; - needed for determinant(), e.g. - - * R/Auxiliaries.R (detSparseLU): using lu(*, errSing=FALSE) - - * R/dgCMatrix.R, R/dsparseMatrix.R: lu() methods, using 'errSing' - - * R/sparseMatrix.R (printSpMatrix): fix bug introduced on *-09-10 - * tests/Simple.R: test for that. - -2009-10-18 Martin Maechler - - * src/dgeMatrix.c (dgeMatrix_crossprod): do not fail in 0-column case. - - * inst/test-tools.R (Q.eq): new utility - (checkMatrix): minimally check %*%, crossprod() and tcrossprod() - - * R/products.R: more '%*%' methods for "[ln]?Matrix", "missing" - -2009-10-06 Martin Maechler - - * DESCRIPTION (Version): 0.999375-31, CRAN-released: 2009-10-06 - - * inst/include/Matrix_stubs.c (M_R_cholmod_error): revert (2009-09-18), - i.e., no longer restore cholmod_common. {{M_cholmod_defaults() - still seems not usable from lme4's init.c}} - -2009-10-05 Martin Maechler - - * src/dtrMatrix.c (dtrMatrix_chol2inv): use "dpoMatrix" instead of "dsy" - - * R/dtrMatrix.R: make use of implicit generic for chol2inv() in - newer R versions. - -2009-09-30 Martin Maechler - - * R/CHMfactor.R (solve): fix methods for "ddiMatrix" and "missing" RHS. - * tests/factorizing.R: test these - - * R/Matrix.R (image): fix Matrix method to work for "ddiMatrix" - * R/diagMatrix.R: coercion to "dsparse*" - * tests/Simple.R: test image() - -2009-09-29 Martin Maechler - - * R/AllGeneric.R: rcond is implicit generic in latest versions of R - -2009-09-22 Martin Maechler - - * R/Ops.R (A.M.n,A.n.M): replace "Ops" methods by explicit "Arith", - "Logic", etc, getting rid of ambiguity (notes), and of infinite recursions. - - * tests/group-methods.R: test these systematically. - -2009-09-18 Martin Maechler - - * inst/include/Matrix_stubs.c (M_R_cholmod_start): print_function = NULL - as in src/chm_common.c (2009-07-20) - (M_R_cholmod_error): ditto, using new M_cholmod_default(), declared in - * inst/include/cholmod.h - - * R/Tsparse.R (intI): do *not* allow logical subscript (row or - column) to be too long, compatibly with traditional matrix - indexing. - * tests/indexing.R: and assert the error. - -2009-09-17 Martin Maechler - - * R/pMatrix.R: as(sparseMatrix, pMatrix) - - * R/CHMfactor.R (solve): method for (CHMfactor, missing) - - * inst/test-tools.R (assertError): use tryCatch() - - * R/diagMatrix.R (.sparseDiagonal): fix shape "g" case. - - * R/Auxiliaries.R (isTriC): do not wrongly return TRUE for a *sCMatrix. - - * man/chol2inv-methods.Rd: document & example - -2009-09-16 Douglas Bates - - * NAMESPACE, R/dtrMatrix.R, src/dtrMatrix.[ch]: add chol2inv() - method for dtrMatrix. - -2009-09-12 Martin Maechler - - * R/sparseVector.R ([): allow *indexing* with "lsparseVector" - -2009-09-11 Martin Maechler - - * R/sparseVector.R (prSpVector): using ":" as in printSpMatrix() - (Summary): add "Summary" group methods - * man/sparseVector-class.Rd: and test a bit - -2009-09-10 Martin Maechler - - * R/sparseMatrix.R (printSpMatrix): visually differentiate - non-structural zeros in *logical* sparse matrices, using ":" - - * R/Auxiliaries.R (setparts): new utility, for - * R/Ops.R (Ops.spV.spV): start implementing sparseVector arithmetic etc - -2009-09-08 Martin Maechler - - * R/dgCMatrix.R (qr): for sparseMatrix must coerce to "dgCMatrix". - * tests/factorizing.R: test qr() - -2009-09-01 Martin Maechler - - * R/CHMfactor.R (determinant): divide previous log(det(.)) by 2; - now returning det(L); and modify the "CHMfactor.warn" message. - - * man/CHMfactor-class.Rd: modify documentation accordingly. - -2009-08-21 Martin Maechler - - * R/spModels.R (`contrasts<-`): a version that can also work with a - "sparseMatrix" value. This is put into R 2.10.0 (devel) as well. - - * src/Mutils.h: rename any_NA to any_NA_in_x - * src/Mutils.c (check_scalar_string): add doxygen doc - -2009-08-15 Martin Maechler - - * R/spModels.R (fac2Sparse): make also work for 'contrasts.arg' = - matrix; - -2009-07-28 Martin Maechler - - * R/spModels.R (contr.sum): need also to define contr.*() as long - as we document (./man/) them. - -2009-07-27 Martin Maechler - - * DESCRIPTION (Version): 0.999375-30, CRAN-released: 2009-07-28 - - * R/Matrix.R (all.equal_Mat): add factorsCheck=FALSE argument - * R/Auxiliaries.R (attr.all_Mat, attrSlots): ditto - -2009-07-25 Martin Maechler - - * R/Auxiliaries.R (attr.all_Mat): fix checking of non-content slots. - * R/Matrix.R (all.equal_Mat): thinko needed s/&&/&/ - * R/sparseMatrix.R (all.equal(.) methods): ditto - ===> Note: all.equal() is more stringent for "Matrix" arguments now! - -2009-07-23 Martin Maechler - - * R/spModels.R (model.spmatrix, sparse2int): "complete" re-write - - * tests/spModel.matrix.R: many tests added - -2009-07-20 Martin Maechler - - * src/chm_common.c (R_cholmod_l_start): set print_function to NULL, - as we have long suggested ==> get rid of random strings seen in - some cholmod warnings. - (R_cholmod_error): call cholmod_l_defaults() before error(), so - we restore the Cholmod global in case of errors. - - * R/ldenseMatrix.R (.rcond_via_d): fix thinko - -2009-07-18 Martin Maechler - - * R/CHMfactor.R (isLDL): need a "!" as 'type' is "is_ll" - * src/dsCMatrix.c (dsCMatrix_Cholesky): update, notably when caching. - - * tests/indexing.R: test col.names printing of sparseMatrix - -2009-07-16 Martin Maechler - - * inst/test-tools.R (allCholesky): new testing function - -2009-07-15 Martin Maechler - - * src/dsCMatrix.c (dsCMatrix_Cholesky): add possibility to set each - of 'perm', 'LDL', and 'super' to NA (in addition to TRUE / FALSE). - in these case, a CHOLMOD-heuristic choses the option "sensibly". - * man/Cholesky.Rd: document the new possibility. - -2009-07-14 Martin Maechler - - * R/rankMatrix.R (rankMatrix): diff(sval) <= 0 - - * R/spModels.R (model.spmatrix): fix case of missing main effect - * tests/spModel.matrix.R: new file - -2009-07-11 Martin Maechler - - * R/sparseMatrix.R (show, printSpMatrix2): both print() and show() - now use printSpMatrix2(), and that now already prints - " x sparse Matrix of class ...". - - * R/CHMfactor.R (isLDL): fix and - * NAMESPACE: export isLDL() - -2009-07-10 Martin Maechler - - * R/spModels.R (model.spmatrix): mf may be simple data.frame - -2009-07-09 Martin Maechler - - * NAMESPACE: export sparse.model.matrix() and - * man/sparse.model.matrix.Rd: document it - -2009-07-08 Martin Maechler - - * R/Tsparse.R (intI): also work for integer dimnames (well ..) - -2009-07-07 Martin Maechler - - * R/sparseMatrix.R: "factor out" sparse model things into - - * R/spModels.R (sparse.model.matrix): new model matrix functions - -2009-06-20 Douglas Bates - - * src/CHMfactor.c: Ensure updated LL stays LL - -2009-06-10 Martin Maechler - - * DESCRIPTION (Version): 0.999375-29, CRAN-released: 2009-06-11 - -2009-06-10 Douglas Bates - - * [r2404] src/Mutils.c: Change value of set_factor to be the cached factor - * [r2403] src/dgCMatrix.c, src/dgCMatrix.h, src/init.c: Comment out - unused dgCMatrix_lusol function - * [r2402] R/dgCMatrix.R: R-level implementation of solve("dgCMatrix", - "missing") - * [r2401] src/dgCMatrix.c: Re-arrange LU factorization to always use the - cached value. - -2009-06-09 Douglas Bates - - * [r2399] src/dgCMatrix.c: PROTECT the result from dgCMatrix_LU in - dgCMatrix_matrix_solve - -2009-06-06 Martin Maechler - - * R/Tsparse.R: add numeric -> Tsparse* coercion - -2009-06-05 Martin Maechler - - * src/Mutils.h, src/dgeMatrix.c: using dngettext(.) only - on future versions of R. - -2009-06-04 Martin Maechler - - * DESCRIPTION (Version): 0.999375-28, CRAN-released: 2009-06-08 - - * po/de.po, po/R-de.po: German translations from Chris Leick. - * inst/po/de/: ditto, via po/update-me.sh - -2009-05-28 Martin Maechler - - * src/chm_common.c, src/cs_utils.c, etc: internationalize more - messages; fix some; thanks to feedback from Chris Leick. - -2009-05-27 Martin Maechler - - * man/denseMatrix-class.Rd, etc: 'factors' is *not* a slot in this - class; found by the upcoming R 2.10.0 codocClasses(). - -2009-05-25 Martin Maechler - - * po/update-me.sh, Matrix.pot, etc: updated *.pot / *.po files - - * DESCRIPTION (Version, Date): 0.999375-27, CRAN-released today. - - * R/sparseVector.R: add as(, "dsparseVector") - (spV2M): now works (again!?) for "isparseVector" (-> "dgTMatrix"). - - * tests/matprod.R: tcrossprod(), sparseVector multiplications, i.e., - features of next paragraph. - -2009-05-23 Martin Maechler - - * R/products.R: move almost all %*%, crossprod(), tcrossprod() - methods to a new file. tcrossprod() compatibility with *fixed* - base-R; enable operations with "sparseVector"s; some extra methods - to avoid ambiguity messages. - -2009-05-20 Martin Maechler - - * R/Auxiliaries.R (.M.v, .v.M): slight speedup, and *use* them in - * R/Matrix.R (crossprod) - -2009-05-18 Martin Maechler - - * R/sparseVector.R (dim<-): prod(as.integer(.)) may overflow! - - * R/Matrix.R (Matrix): Matrix(sv, ..) now works for a sparseVector - - * R/sparseVector.R (spV2M): allow zero nrow or ncol. - -2009-05-16 Martin Maechler - - * R/sparseMatrix.R (dim<-): should also work for diagonalMatrix. - - * inst/test-tools.R (vec): add test for dim(x) <- c(n, 1) - - -2009-05-13 Martin Maechler - - * R/rankMatrix.R (rankMatrix): new function for matrix rank; - "to be discussed". - -2009-05-07 Doug Bates and Martin Maechler - - * src/Mutils.c (Matrix_check_class_etc): diverse smallish - improvements, stuffing a small leak, adding PROTECT()ion, - but, alas, not enough. - -2009-05-06 Martin Maechler - - * R/HBMM.R (readMM): add checkIJ() to produce nicer error messages, - on malformed input. - * tests/write-read.R: and test that.. - -2009-04-18 Martin Maechler - - * DESCRIPTION (Version, Date): 0.999375-26, CRAN-released on 2009-04-29 - - * src/Mutils.h (Matrix_check_class_etc): try to ensure it works - also when Matrix is loaded but not attached. - - * src/init.c (R_init_Matrix): get Matrix namespace in C. - - * R/zzz.R (.onLoad): a *working* fix aka "hack" - - -2009-04-15 Martin Maechler - - * DESCRIPTION (Version): 0.999375-25 packaged --> CRAN - -2009-04-09 Martin Maechler - - * R/Auxiliaries.R (Matrix.msg): new utility - (.M.vectorSub): ditto, for vector-indexing; in - * R/Matrix.R ([): M[i, drop=] should do *vector*-indexing. - * R/Tsparse.R ([): ditto; - * R/diagMatrix.R ([, subDiag): ditto. - - * R/Tsparse.R ([): more careful indexing of (triangular) - TsparseMatrix. - - * tests/indexing.R: testing the above - - * R/Auxiliaries.R (gT2tT, check.gT2tT): consistency and efficiency - improvement, using 'do.n' argument. - -2009-04-08 Martin Maechler - - * R/Matrix.R: add as(., "vector") etc - * man/Matrix-class.Rd: ditto - * inst/test-tools.R (checkMatrix): and check them - -2009-04-07 Martin Maechler - - * DESCRIPTION (Version, Date): 0.999375-24 ... released to CRAN - - * R/sparseVector.R: fix coercion from xsparse* - * tests/Simple.R: and check it. - - * man/lu.Rd: document 'warnSing' argument - * tests/dg_Matrix.R: and test it - - * src/dgeMatrix.c (dgeMatrix_LU): missing 'return' - -2009-04-06 Martin Maechler - - * DESCRIPTION (Version, Date): 0.999375-24 only for 2.9.0 (and newer) - -2009-03-30 Martin Maechler - - * src/Mutils.h (Matrix_check_class_etc): new version that *computes* - the environment to be used. - * src/Tsparse.c, src/cs_utils.c, src/chm_common.c, src/dgCMatrix.c: - use the above in place of Matrix_check_class_and_super(). - -2009-03-26 Martin Maechler - - * R/Auxiliaries.R (MatrixClass): use cld@package, not packageSlot(.) ! - -2009-03-25 Martin Maechler - - * DESCRIPTION (Version, Date): 0.999375-23 - - * tests/Class+Meth.R (dotestMat): use getClassDef() for speed; - adaptions to also work with "Matrix-extensions". - * inst/test-tools.R (checkMatrix): similar adaptions. - -2009-03-24 Martin Maechler - - * R/Auxiliaries.R (MatrixClass, l2d_meth, ...): make use of - MatrixClass(): functionality for setClass("foo", contains="dgCMatrix") - - * src/Mutils.h (Matrix_check_class_and_super): new utility, to be - used in lieu of Matrix_check_class() in most cases. - * src/Tsparse.c, src/cs_utils.c, src/chm_common.c, src/dgCMatrix.c: use it, - currently only with R_GlobalEnv {otherwise: API change} - -2009-03-12 Martin Maechler - - * man/band.Rd: note and example about band() |-> dense - - * R/ddenseMatrix.R (.bandDense): fix typo in argument check. - * R/Csparse.R: ditto - * src/dense.c (dense_band): limit index range; thanks to Bill Dunlap. - * tests/Simple.R (band): check it - -2009-03-11 Martin Maechler - - * R/dMatrix.R (Summary): (dsparseMatrix): new logic, fixing - prod(). - - * inst/test-tools.R (eqDeterminant): improve after suggestions from Doug - * inst/test-tools.R (checkMatrix): message(), not warning(), - for differing prod(). - - * src/dgeMatrix.h, src/init.c: dgeMatrix_LU(x, warn_singularity) - * src/dgeMatrix.c (dgeMatrix_LU): allow to suppress singularity warning. - (dgeMatrix_determinant, dgeMatrix_rcond): do *not* warn on exact singularity - -2009-03-10 Martin Maechler - - * R/Matrix.R (Summary.l, Summary.np): no conversion to "dMatrix" for all()/any() - - * tests/Simple.R: do not use memory-expensive all(I.. == Diagonal(n)) - for R >= 2.9.0, keep option 'warn = 2', i.e. warnings - producing errors ==> - - * R/Ops.R: small fixes, getting rid of warnings in tests/Simple.R - - * TODO: think about providing something like allCompare(.,., "==") - -2009-03-02 Martin Maechler - - * DESCRIPTION (Version): 0.999375-22, ready but not submitted - -2009-03-01 Martin Maechler - - * inst/test-tools.R (eqDeterminant): fix for NAs - -2009-02-28 Martin Maechler - - * R/diagMatrix.R (prod, all): fix methods (cut & paste remnant); - for NA. - - * R/Auxiliaries.R (detSparseLU): determinant(<..NA..>) now gives NaN - - * R/sparseMatrix.R (printSpMatrix): workaround format.info() R bug - * tests/Simple.R: test that - -2009-02-27 Martin Maechler - - * R/Matrix.R (Matrix): Matrix(x,*) returns x unaltered when - is(x, "diagonalMatrix"); OTOH, Matrix(x,1,1) should typically - *not* return a diagonalMatrix. - - * R/diagMatrix.R (setAs(matrix, *)): fix for NAs. - * tests/Simple.R: test things like Matrix(NA, 1,1) - -2009-02-25 Martin Maechler - - * NAMESPACE: add bandSparse() "properly" - * man/bandSparse.Rd: doc including examples - -2009-02-20 Martin Maechler - - * R/bandSparse.R (bandSparse): constructor for band(diagonal) - sparse matrices. - -2009-02-13 Martin Maechler - - * DESCRIPTION (Version): 0.999375-21, released to CRAN - - * src/Mutils.h (Matrix_with_SPQR): #define (or #undef), in order to - allow easy building "SPQR-free version" of Matrix. - -2009-02-11 Martin Maechler - - * R/Csparse.R (replCmat): another check for 'has.x' - * tests/indexing.R: very large (very sparse) sub-indexing. - -2009-01-30 Martin Maechler - - * DESCRIPTION (Version): 0.999375-20 - * man/spqr.Rd: disable example on Windows for now - - * inst/test-tools.R (checkMatrix): simpler for "logical" pMatrix. - - * R/ngTMatrix.R: fix ngT -> nge / lge coercion, and - * tests/Simple.R: adapt test. - - * R/pMatrix.R: coercion to "matrix": more sensical to coerce to logical - instead of 0/1 integer. This a mild back-compatibility breach. - * man/pMatrix-class.Rd: adapt, including example - * R/sparseMatrix.R (printSpMatrix): print as logical, too. - -2009-01-29 Martin Maechler - - * R/Auxiliaries.R (geClass): define for "pMatrix" - * R/pMatrix.R: pMatrix -> ngeMatrix coercion - * man/pMatrix-class.Rd: - - * DESCRIPTION (Version): 0.999375-19 considering release. - - * R/ngTMatrix.R: coercing correctly to ngeMatrix via lge*. - - * man/dgTMatrix-class.Rd: remove no-more-existing coercion method; - mention a bit less, and note. - - * R/diagMatrix.R: ensure that * is diagonal - even with new method dispatch - - * R/Matrix.R (.M.sub.i.logical): change comments only, and extend - error message which is never called currently. - - * R/Ops.R (Ops.x.x): fix for new "inherited method" dispatch - -2009-01-28 Martin Maechler - - * R/spqr.R: put spqr() / class "SPQR" related code in one file for now. - * NAMESPACE: export "SPQR" class - * man/SPQR-class.Rd: document it - -2009-01-27 Martin Maechler - - * R/Auxiliaries.R (is.na_nsp): should produce "nsCMatrix" only when - dimnames are symmetric too. - - * R/sparseQR.R (solve): method for (, ) - -2009-01-26 Douglas Bates - - * src/cs_utils.c (Matrix_as_cs): in diagU2N case: Force sorted - columns after expanding unit diagonal. - -2009-01-21 Martin Maechler - - * R/Tsparse.R (intI): for now give intelligible error on NA indices. - - * R/Matrix.R (subset.ij): should deal correctly with NA indices in - the 2-column-matrix index case. - - * src/Mutils.c (m_encodeInd, m_encodeInd2): prepare to deal better - with NA in indices. - -2009-01-20 Martin Maechler - - * inst/doc/Intro2Matrix.Rnw: mention sparseMatrix - * man/sparseMatrix.Rd: fix typos found by parse_Rd - -2009-01-19 Martin Maechler - - * DESCRIPTION (Version,Date): release 0.999375-18 (as of yesterday) - -2009-01-18 Douglas Bates - - * [r2319] src/Csparse.c, src/Csparse.h: Added an untested create_Csparse - function - -2009-01-17 Martin Maechler - - * R/sparseMatrix.R (sparseMatrix): recycle 'x' if necessary; - use 'index1' rather than 'index0' {as proposed by Doug}. - - * R/dgeMatrix.R: drop two unnecessary (maybe very very slightly - faster) methods for %*%. - -2009-01-12 Douglas Bates - - * [r2313] R/sparseMatrix.R, man/sparseMatrix.rd: Use intermediate - triplet rep in sparseMatrix. - -2009-01-07 Martin Maechler - - * R/sparseMatrix.R (sparseMatrix): new function to be used in - place of new(), notably for CsparseMatrix objects. - * man/sparseMatrix.Rd: its doc; plus examples - - * man/CsparseMatrix-class.Rd: fix long-standing thinko; update the - \note{} to reflect the fact that the validity method no longer - modifies its argument. - - * R/Csparse.R (.validateCsparse): new utility accessing - * src/Csparse.c (Csparse_validate_, Csparse_validate2): new - utilities allowing the "sort-in-place" that used to be part of - former Csparse_validate(). - -2008-12-10 Douglas Bates - - * [r2309] DESCRIPTION: Release 0.99375-17 - -2008-12-05 Douglas Bates - - * [r2308] inst/include/Matrix_stubs.c, inst/include/cholmod.h: Adding - const for picky compilers; Using CHM_FR typedef - -2008-11-17 Martin Maechler - - * [r2307] ChangeLog, tests/validObj.R: update (as of Oct.26) - -2008-10-31 Douglas Bates - - * [r2306] inst/include/Matrix_stubs.c: Consistency with cholmod.h re - const qualifier - -2008-10-26 Douglas Bates - - * src/Csparse.c (Csparse_validate): do not sort, but report error - on unsorted row-indices within column. - - * tests/validObj.R: adapt the regression test. - -2008-10-17 Douglas Bates - - * [r2300] inst/include/Matrix_stubs.c: Consistency with SparseSuite names - * [r2299] src/AMD/Source/Makefile, src/CHMfactor.c, - src/CHOLMOD/Lib/Makefile, src/COLAMD/Source/Makefile, src/Csparse.c, - src/Tsparse.c, src/UFconfig/UFconfig.h, src/chm_common.c, - src/chm_common.h, src/dense.c, src/dgCMatrix.c, src/dsCMatrix.c, - src/init.c, src/t_gCMatrix_colSums.c: Compile only the UF_long version - of SparseSuite but setting UF_long to be int - -2008-10-17 Martin Maechler - - * src/Mutils.h: include also for AIX - -2008-10-16 Martin Maechler - - * src/scripts/ : move *.sh, *.mkf and *.mk files from src/ to src/scripts/ - - * src/*/{Lib|Source}/Makefile: change to non-GNU-make style, and - other improvements from Brian Ripley. - -2008-10-15 Martin Maechler - - * src/Makevars, src/Makevars.win: replacing - * src/Makefile, src/Makefile.win which are no longer needed - -2008-10-14 Martin Maechler - - * src/Makefile, src/*/Makefile: changes from Brian Ripley enabling - parallel make. - -2008-10-13 Douglas Bates - - * [r2285] src/SPQR/Include/spqr.hpp: Include the C string declarations. - -2008-10-06 Douglas Bates - - * [r2284] src/Makefile, src/SPQR/Lib/Makefile_SPQR, src/SPQR/Makefile, - src/Win.mk: Modifications for compilation under Windows - * [r2283] tests/validObj.R: remove reference to undefined object - * [r2282] inst/doc/UFsparse/SPQR.txt: SPQR license information - -2008-10-04 Douglas Bates - - * [r2281] src/Makefile, src/chm_common.c, src/chm_common.h, - src/dgCMatrix.c, src/init.c: SparseSuiteQR interface - -2008-10-03 Martin Maechler - - * src/Csparse.c (isValid_Csparse): new utility - * src/chm_common.c (as_cholmod_sparse): add validity check early - * tests/validObj.R: test the above - -2008-10-02 Douglas Bates - - * [r2277] TODO: Think of a better way of structuring include files - * [r2276] src/AMD/Source/Makefile, src/CHOLMOD/Lib/Makefile, - src/COLAMD/Source/Makefile, src/Makefile, - src/SPQR/Include/SuiteSparseQR_C.h, src/chm_common.h, src/dgCMatrix.c, - src/dgCMatrix.h, src/init.c, src/sparseQR.h: Add dgCMatrix_SPQR and - modify other code to allow its compilation and linkage - -2008-10-02 Martin Maechler - - * R/sparseMatrix.R (printSpMatrix2): another validObject(.) call - * src/Makefile: getting rid of SUBSTAMP etc, thanks to Simon U. - -2008-10-01 Douglas Bates - - * src/Makefile, inst/doc/UFsparse/, inst/include/: Adjustments for - SuiteSparse configuration. - * src/CHOLMOD, src/UFconfig: update versions of SuiteSparse - libraries. - * src/SPQR/* add "SPQR", from Tim Davis' "SuiteSparse" collection, - not yet with an R interface. - -2008-09-25 Martin Maechler - - * Release 0.999375-15 to CRAN. - -2008-09-23 Martin Maechler - - * src/dsyMatrix.c (dsyMatrix_as_dspMatrix): copy 'factors' slot - * tests/dpo-test.R: test for it - - * R/Tsparse.R (intI): improve one error message. - -2008-09-22 Martin Maechler - - * DESCRIPTION (Version): 0.999375-15 to be released - - * R/diagMatrix.R (diag o ): explicit setMethods, - in order to keep result diagonal in many more cases. - (coerce) to denseMatrix now *does* coerce. - - * man/diagonalMatrix-class.Rd: \alias{} for these. - - * R/Auxiliaries.R (.dense.prefixes): ".diMatrix" has not been dense - anymore! - - * R/dMatrix.R: as(, ): enable Cholesky/BunchKaufman - -2008-09-20 Martin Maechler - - * R/lsparseMatrix.R (C2l): as(, "lMatrix") should - preserve NAs. - - * R/Ops.R (Arith, Compare): keep diag="U" in more cases when sensible - -2008-09-19 Martin Maechler - - * R/Ops.R (Arith): better o which - preserves symmetry / triangularity in "obvious" cases. - - * R/dpoMatrix.R: setAs(., "lMatrix") and setAs(., "nMatrix") - * R/dppMatrix.R: ditto - * man/dpoMatrix-class.Rd: - - * inst/test-tools.R (checkMatrix): add [dln]Matrix <-> [dln]Matrix - coercion checks. - * tests/indexing.Rout.save: more ambiguity warnings from the new checks. - - * R/dMatrix.R: dMatrix -> nMatrix: fix dense case. - -2008-09-18 Martin Maechler - - * tests/factorizing.R: test expand() - - * src/dgCMatrix.c (dgCMatrix_LU): fill @Dim slot correctly. - - * R/Tsparse.R (replTmat): also optionally warn when sub-assignment - loses symmetry of Matrix. - -2008-09-17 Martin Maechler - - * R/Tsparse.R (.TM.repl.i.2col): fix sub-assignment of "dsCMatrix"; - bug reported by Jacob van Etten. - - * tests/indexing.R, tests/indexing.Rout.save: testing it - -2008-09-10 Douglas Bates - - * [r2260] src/Mutils.h, src/chm_common.c, src/cs_utils.c, - src/dsCMatrix.c: Update doxygen comments - -2008-09-10 Martin Maechler - - * man/ddiMatrix-class.Rd: docu update : sparse, not dense (see 2008-07-28). - * man/ldiMatrix-class.Rd , man/diagonalMatrix-class.Rd: ditto. - -2008-09-08 Martin Maechler - - * DESCRIPTION (Version): 0.999375-14 to be released to CRAN, - in order to comply to to pending changes in R-devel (2.8.0). - -2008-09-05 Martin Maechler - - * R/Matrix.R (solve,..): use (Matrix,ANY)... instead of (Matrix,Matrix) - * R/Ops.R: ditto; most of these from John Chambers. - * man/Matrix-class.Rd - - * R/Auxiliaries.R (.diagU2N): fix for Rsparse* - * tests/Simple.R: test it - -2008-09-02 Martin Maechler - - * man/drop0.Rd: - * R/Auxiliaries.R (drop0): new argument 'tol = 0' (and - 'is.Csparse'); entails updates in - * R/Csparse.R, R/lMatrix.R, R/dMatrix.R - -2008-08-30 Martin Maechler - - * DESCRIPTION (Version): 0.999375-13 released to CRAN - -2008-08-29 Martin Maechler - - * R/Auxiliaries.R (attrSlotNames, attrSlots, attr.all_Mat): new - utilities, for now to be used in all.equal() methods. - - * R/Matrix.R (all.equal_Mat): be more careful (less tolerant) in - all.equal() methods. - * R/sparseMatrix.R: ditto - -2008-08-28 Martin Maechler - - * DESCRIPTION (Version): 0.999375-12 released to CRAN - - * R/Ops.R (Compare(,): fix for nsparseMatrix - -2008-08-27 Douglas Bates - - * R/sparseMatrix.R: fac2sparse() for NA's - -2008-08-26 Martin Maechler - - * R/sparseVector.R (all.equal.sparseV): make work for nsparseVector's - -2008-08-25 Martin Maechler - - * src/dgCMatrix.c (dgCMatrix_LU): partially revert change - r2175 (2008-04-23) and do give an error for a singular matrix. - * man/lu.Rd, R/Auxiliaries.R (detSparseLU): adapt to that. - - * R/LU.R: expand( ) - - * NAMESPACE, man/all.equal-methods.Rd: - * R/Matrix.R (all.equal): methods for Matrices, - * R/sparseMatrix.R, R/sparseVector.R: sparseMatrices and -Vectors - - -2008-08-23 Douglas Bates - - * [r2243] R/CHMfactor.R, man/CHMfactor-class.Rd: Added expand method and - documentation for CHMfactor class - * [r2241] R/CHMfactor.R: Added trivial coercion method for CHMfactor -> - pMatrix - -2008-08-19 Martin Maechler - - * R/nsCMatrix.R (setAs(., dgTMatrix))): remove unused method - -2008-08-18 Martin Maechler - - * R/Ops.R (.Ops2dge.via.x, and many others): eliminate never-used - variable assignments {from "next version" codetools}. - * R/Csparse.R, R/Tsparse.R: ditto - -2008-08-17 Martin Maechler - - * R/sparseVector.R (sp2vec, coerce to sparseVector): - make sure no integer overflows happen, and use double precision - (n, index) where appropriate. - - * tests/Simple.R: test "large" sparse vectors. - -2008-08-14 Martin Maechler - - * R/AllClass.R (sparseVector): 'length' and 'i' are "numeric", not - just integer (allowing much longer length). - -2008-07-28 Martin Maechler - - * R/AllClass.R (diagonalMatrix): extend "sparseMatrix" instead of - "denseMatrix". This renders "scarceMatrix" dispensable and - invalidates part of MM's presentations on "space of Matrix - classes", but seems cleaner overall. - - * R/diagMatrix.R, etc: eliminate "scarceMatrix", replacing it by - "sparseMatrix" in method signatures; - further, instead of coercing to "sparseMatrix", coerce to - "TsparseMatrix" now. - -2008-07-26 Martin Maechler - - * src/dgCMatrix.c (dgCMatrix_qrsol): allow third argument 'order' in - .Call(.) - * R/sparseMatrix.R (lm.fit.sparse), src/dgCMatrix.h, src/init.c: ditto - -2008-07-24 Martin Maechler - - * R/dgeMatrix.R: need solve(, ) against infinite - recursion in new test in - * tests/matprod.R: testing the above and another solve() case, below - - * R/sparseMatrix.R (lm.fit.sparse): slightly more efficient for - e.g. triangular 'x'. - - * src/dgCMatrix.c (dgCMatrix_qrsol): use AS_CSP() hence allowing dtC* - * src/dgCMatrix.c (dgCMatrix_cholsol): use AS_CHM_SP() to work - correctly with unit-triangular x - - * src/dsCMatrix.c (dsCMatrix_Csparse_solve): use AS_CHM_SP() - instead of not checking diagU. - - * R/diagMatrix.R, R/Auxiliaries.R: tweaks to allow later - experiments where diagonalMatrix would extend sparse*. - - -2008-07-23 Martin Maechler - - * src/dgCMatrix.c (compressed_non_0_ij): for nnzero(), use - "x@p[ncol(x)]" instead of length(x@i). - * R/Auxiliaries.R (nnzero): ditto - - * src/chm_common.c (as_cholmod_sparse): Do not use LENGTH(islot) - since that fails for over-allocated i. - - * tests/validObj.R: more testing of over-allocated (i,x) slots - -2008-07-23 Martin Maechler - - * src/chm_common.c (as_cholmod_sparse): Add 4th argument - 'sort_in_place' and set it to TRUE in call from - * src/Csparse.c (Csparse_validate): sort in place, making sure that - validObject() continues sort the columns if needed. - - * tests/validObj.R: now tests that more thoroughly, and - * man/CsparseMatrix-class.Rd: documents it. - -2008-07-22 Douglas Bates - - * [r2233] src/chm_common.c: sort columns more carefully in - as_cholmod_sparse - * [r2231] src/chm_common.c: check for sorted columns in - as_cholmod_sparse instead of assuming it - * [r2228] src/cs_utils.c: Use R_alloc for semi-permanent storage when - expanding diagonal - - -2008-07-21 Martin Maechler - - * src/cs_utils.c (Matrix_as_cs): add check_Udiag switch - * src/cs_utils.h (AS_CSP & AS_CSP__): two versions - * src/dtCMatrix.c (dtCMatrix_sparse_solve): no longer needs - diagU2N(). - - * R/diagMatrix.R (.sparseDiagonal): new utility, used in "old" - .spDiagonal() and new .trDiagonal(). - - * R/dtCMatrix.R (solve): make use of .trDiagonal() - -2008-07-19 Martin Maechler - - * R/Auxiliaries.R (dimNamesCheck): fix very long standing buglet, - stumbled upon by Michael H. - - * tests/validObj.R: testing the bug-fix - - * src/chm_common.h (AS_CHM_SP__, AS_CHM_TR__): the non-diagU2N - versions. - - * src/chm_common.c (as_cholmod_sparse, as_cholmod_triplet): new - boolean argument 'check_Udiag' (formerly implicitly was FALSE). - - * src/Csparse.c (Csparse_Csparse_prod, ...): etc: use the new - AS_CHM_SP() which includes diagU2N. - - * inst/include/Matrix_stubs.c (M_as_cholmod_sparse): similarly - adapt to 3 arguments. - -2008-07-17 Douglas Bates - - * [r2220] src/dtCMatrix.c: Correct number of columns for rhs in - dtCMatrix_sparse_solve - * [r2219] src/cs_utils.c: eye->nz should be -1 for compressed column - storage - * [r2217] R/dtCMatrix.R, src/dtCMatrix.c, src/dtCMatrix.h, src/init.c: - Replace dtCMatrix_solve by more general dtCMatrix_sparse_solve and add - new solve method - * [r2216] src/cs_utils.c: Utility csp_eye in case we want to fold - Csparse_diagU2N functionality into Matrix_as_cs - -2008-07-17 Martin Maechler - - * inst/doc/factor2sparse.Rnw: revive a year-old unfinished vignette - -2008-07-16 Douglas Bates - - * [r2212] R/sparseMatrix.R: fac2sparse gets a drop.unused.levels - argument with default TRUE - * [r2211] inst/include/Matrix.h, inst/include/Matrix_stubs.c: Export - Csparse_diagU2N - * [r2210] R/dtCMatrix.R, man/dtCMatrix-class.Rd, src/dtCMatrix.c, - src/dtCMatrix.h, src/init.c: Remove vestigial methods based on the - parent array; new solve method for signature(a = "dtCMatrix", b = - "sparseMatrix") - -2008-07-14 Martin Maechler - - * R/diagMatrix.R (.bdiag): also keep result triangular - -2008-07-13 Douglas Bates - - * [r2208] src/dtCMatrix.c, src/dtTMatrix.c: Revise wording in error - messages - -2008-07-12 Martin Maechler - - * R/diagMatrix.R (.bdiag): make more sophisticated, keeping - symmetric- or nsparse-Matrix, also fixing the bug introduced with - the new version (.999375-10). - - * tests/Simple.R: regression test for the bugs. - -2008-07-07 Martin Maechler - - * R/sparseVector.R (spV2M): fix for "isparseVector". - Further, fix -> coercion - (repSpV): add rep(, .) method - - * R/sparseMatrix.R ([<-): when the RHS is a "scarce"Matrix, do *not* - use as.vector(.). - - * R/Tsparse.R ([<-): & [CR]sparse.R: methods for value = "sparseVector" - - * R/AllClass.R: new (hidden) class union "scarceMatrix" - * R/AllClass.R: sparseVector: add prototype - - * src/dsCMatrix.c (chk_nm): fix "Cholesky" check, thanks to - Kasper Kristensen - -2008-06-28 Martin Maechler - - * tests/other-pkgs.R: add more interesting Matrix -> graph - -2008-06-27 Martin Maechler - - * R/sparseMatrix.R (Tsp2grNEL): add 'need.uniq' argument for speed - improvement in "[CR]sparseMatrix" -> "graphNEL" coercion. - -2008-06-25 Martin Maechler - - * DESCRIPTION: release 0.999375-10 to CRAN - - * R/diagMatrix.R (.bdiag): use more efficient Bates-algorithm for - .bdiag(). - * man/bdiag.Rd: update, including test. - -2008-06-24 Martin Maechler - - * tests/Simple.R: add minimal bdiag() checks - - * R/diagMatrix.R (bdiag): fix bdiag(diag(4:5)) case; - factor out the Tsparse - internal computation into new .bdiag() - -2008-06-14 Martin Maechler - - * R/nsparseMatrix.R: specific fast all() and any() methods. - - * src/dgTMatrix.c (MAKE_gTMatrix_to_geMatrix): better error message - when trying to produce too large result. - - * inst/test-tools.R (checkMatrix): add 'do.matrix' with a check for - very large matrices. - - * R/sparseMatrix.R (printSpMatrix2): new function, used by - show(). - - * R/diagMatrix.R (show): print only diagonal entries when nrow(.) >= 50. - -2008-06-13 Martin Maechler - - * src/Mutils.c (m_encodeInd, m_encodeInd2): - * R/*.R: .Call(m_encodeInd2?, ..) instead of encodeInd2?(..). - - * R/Auxiliaries.R (encodeInd2?): care against integer overflow (for - big matrices): -> 2nd argument 'di' instead of 'nr' ! - -2008-06-09 Martin Maechler - - * R/dtCMatrix.R: faster dgC* -> dtC* method - * tests/Simple.R: (hence removed "FIXME") - - * R/Auxiliaries.R (copyClass): add 'check = .copyClass.check' - which can be be turned off for possible speed gains. - -2008-06-02 Martin Maechler - - * R/dMatrix.R, R/Tsparse.R: get rid of some disambiguation warnings - - * R/not.R: fix ! - implement two old "FIXME"s: ! & ! - - * R/diagMatrix.R (Ops): fix one-off thinko in - o methods - - * inst/test-tools.R: checkMatrix(): drop0(m12) unconditionally - -2008-05-30 Martin Maechler - - * R/diagMatrix.R (Ops): more o methods - * tests/Simple.R: testing the above - -2008-05-07 Martin Maechler - - * NAMESPACE: also import "Ops" {found via new tools:::get_S4_generics_with_methods()} - -2008-05-03 Martin Maechler - - * R/Matrix.R (dimnames<-): dimnames(.) <- NULL works with a message; - * NAMESPACE (unname): finally exported - -2008-04-28 Martin Maechler - - * R/Auxiliaries.R (cholMat): possibly keep packed - - * inst/test-tools.R (checkMatrix): add fixed determinant.matrix() - for R < 2.7.0. - - * R/Tsparse.R ([): for triangularMatrix, check if result may be - triangular. - * tests/indexing.R: regression test - -2008-04-26 Martin Maechler - - * inst/test-tools.R (checkMatrix): test not only for dMatrix - * R/Ops.R: more methods, for lMatrix/nMatrix - -2008-04-24 Martin Maechler - - * R/expm.R: new file for all expm() methods; notably new ones for - diagonal*, triangular*, symmetric* and also simple "matrix". - -2008-04-23 Martin Maechler - - * R/dsyMatrix.R: setAs(.) rather than setIs(.) "to dpoMatrix" - - * inst/test-tools.R (assert.EQ.mat): better message when not equal - - * src/dgeMatrix.c (dgeMatrix_matrix_crossprod): allow integer RHS. - - * R/AllClass.R (ddiMatrix,ldiMatrix): extend [dl]Matrix instead of - [dl]denseMatrix {identical slots; but more logical method dispatch} - -2008-04-23 Martin Maechler - - * R/sparseMatrix.R (summary): use length() instead of nnzero() - - * R/diagMatrix.R (determinant): method - -2008-04-22 Martin Maechler - - * src/dsCMatrix.c (dsCMatrix_Cholesky): suppress CHOLMOD printing too - - * R/Rsparse.R (.viaC.to.R): uplo fix for symmetric & triangular - - * R/dsCMatrix.R (determinant): switch to use Cholesky( 'LDL' ) and - its diagonal - -2008-04-21 Martin Maechler - - * R/dMatrix.R (Summary): short-cut for prod() - - * R/not.R: fix !<"U"-diag-ltrMatrix> ; drop {R < 2.6.0}-branch - - * R/Auxiliaries.R (tT2gT): fix for "U"-diag-ltTMatrix - - * R/AllClass.R: ltTMatrix validity: call (more stringent) tTMatrix_validate - -2008-04-19 Martin Maechler - - * R/Ops.R (Ops.x.x): aux.function, not just for (dMatrix,dMatrix) - - * R/Ops.R (.do.Logic.lsparse): fix NA case - - * R/Tsparse.R (replTmat): fix r.sym case, using indTri() - - * R/Auxiliaries.R (nnzero): fix special cases, using - * R/Auxiliaries.R (indDiag): new utility - * R/Auxiliaries.R (indTri): new argument 'diag' - - * R/dMatrix.R: Summmary(): different branch for dsy* - -2008-04-18 Martin Maechler - - * R/diagMatrix.R: "Summary" methods, and more "Arith" / "Ops" - - * src/Csparse.c (Csparse_drop): preserve (uplo,diag) for - ".tCMatrix" triangular matrix. - - * R/Auxiliaries.R (drop0): use diagU2N(); according to helpfile, - 'clx' can be class *or* class representation. - (nnzSparse): new. - (non0ind): new 'xtendSymm' argument; used in - * R/Ops.R: - - * R/pMatrix.R: more efficient determinant() method - -2008-04-17 Martin Maechler - - * [deactivated] R/Matrix.R (det): det() as base::det(), but with Matrix environment - * [deactivated] NAMESPACE: export(det) - - * R/diagMatrix.R (mkDiag): new substitute for diag() - - * R/Auxiliaries.R (nnzero): fix for "diagU2N"-case - (as0, as1): new utility. - - * R/Csparse.R (diag, band): need diagU2N() - - * src/dgeMatrix.c (get_norm): return NA when 'x' slot has NA/NaN. - - * R/lMatrix.R: coerce(nMatrix |-> lMatrix) fix - - * R/Ops.R (Compare): fix o case - - * R/nsparseMatrix.R, R/pMatrix.R: is.na() method - -2008-04-16 Martin Maechler - - * R/Auxiliaries.R (signPerm): new utility for - * R/Auxiliaries.R (detSparseLU): determinant() via sparseLU - - * src/dsCMatrix.c (dsCMatrix_LDL_D): suppress 'CHOLMOD warning'; - since we have our own error handler (and can use tryCatch() from R). - -2008-04-15 Martin Maechler - - * R/dgTMatrix.R (image): ha! set col=NA if lwd is small -- very nice! - - * R/dsCMatrix.R (determinant): use tryCatch() and switch to lu() - when not positive definite. - - * R/Auxiliaries.R (asTri): new auxiliary - (non0.i): make *much* faster for Tsparse with many non-zero - entries; add 'uniqT = TRUE' argument to be used when sorting is undesired. - - * tests/Class+Meth.R (dotestMat): now build on checkMatrix() from - test-tools.R , see below - - * R/dMatrix.R: fix "Summary" method: + diagU2N() - - * NAMESPACE, R/Matrix.R: add mean() method - -2008-04-12 Martin Maechler - - * R/dgTMatrix.R (image): allow to *not* take abs(), and use - cold-hot colors; try *changed* default useAbs = FALSE - * man/image-methods.Rd: nice examples - -2008-04-11 Martin Maechler - - * inst/test-tools.R (checkMatrix): new function for modularizing part - of tstMatrixClass() in tests/Class+Meth.R - - * R/dsCMatrix.R: coercion from "matrix" - * R/ltTMatrix.R, R/ntTMatrix.R: ditto (change it). - - * tests/Class+Meth.R (tstMatrixClass): some cleanup; - add test for as(matrix(,0,0), ) - -2008-04-09 Martin Maechler - - * src/dgeMatrix.c (dgeMatrix_determinant): for n == 0, - work as determinant(matrix(,0,0)). - - * src/dgCMatrix.c (dgCMatrix_LU): return {L,U} as dtCMatrix - * man/sparseLU-class.Rd: ditto - - * R/dgCMatrix.R (determinant): method building on lu() - * R/sparseMatrix.R, Matrix.R, ...: ditto - - * R/Auxiliaries.R (mkDet): auxiliary for determinant() - -2008-04-07 Martin Maechler - - * R/sparseMatrix.R (summary): no 'x' column for - pattern matrices. - -2008-04-02 Martin Maechler - - * src/dense.c (dense_to_Csparse): all declarations *before* - R_CheckStack(); needed e.g. for ancient gcc 2.96. - - -2008-03-29 Martin Maechler - - * DESCRIPTION (Version): 0.999375-9 --- need to release - for R-devel (and R 2.7.0 alpha soon). - - * R/AllClass.R: drop "double" from "atomicVector" class union - - * R/AllGeneric.R (rcond): check for base::rcond instead of R version - * R/dgeMatrix.R: ditto - - * R/sparseMatrix.R (summary): count NAs - - * inst/doc/Intro2Matrix.Rnw: changes, aim for *one* introduction. - -2008-03-28 Martin Maechler - - * R/AllGeneric.R: rcond(., norm) instead of rcond(., type), in - order to become compatible to new R 2.7.0 base::rcond(). - -2008-03-25 Martin Maechler - - * DESCRIPTION (Version,Date): 0.999375-8 --- released to CRAN - - * R/diagMatrix.R (Ops): fix newly introduce bug in <.di> o - - * inst/test-tools.R (isValid): new utility, used much in - * tests/simple.R: - - * man/BunchKaufman-methods.Rd: added too - -2008-03-24 Martin Maechler