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 * R/KhatriRao.R: Efficient KhatriRao() by Michael Cysouw 2012-11-10 Martin Maechler * DESCRIPTION (Version): 1.0-11 * 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 * R/dsyMatrix.R: add BunchKaufman() methods. * R/dspMatrix.R: * src/init.c: add dsyMatrix_trf 2008-03-23 Douglas Bates * DESCRIPTION (Version): release 0.999375-7 * src/CHMfactor.c (CHMfactor_update): fix 2008-03-22 Martin Maechler * src/dsCMatrix.c (dsCMatrix_LDL_D): cleanup, also using internal_chm_factor(). * R/AllGeneric.R: do *not* define a "chol" generic (but rather use the implicit one *without* making pivot part of the signature) * R/*.R: drop the 'pivot' from chol's signature and make 'pivot=FALSE' a default argument of method definitions. * .Rbuildignore: add 'wrld_1deg': I.e. do not put it into released version of Matrix 2008-03-18 Martin Maechler * R/Tsparse.R (.T.2.n): drop 0's before coercion to "nMatrix" * R/sparseMatrix.R (is.na): new simple method * R/denseMatrix.R (is.na): ditto. * R/diagMatrix.R (.symDiagonal): newly exported utility. * R/diagMatrix.R (Ops): * should not become dgeMatrix! * src/UFsparse_download.sh: --> * src/CHOLMOD/: update to CHOLMOD version 1.6 (Nov.2007) 2008-03-17 Martin Maechler * src/dsCMatrix.c (dsCMatrix_LDL_D): even faster utility with same functionality; barely tested in * tests/factorizing.R * src/Csparse.c (diag_tC): new functionality to enable faster determinant(, .) in * R/dsCMatrix.R 2008-03-17 18:53 Douglas Bates * R/CHMfactor.R, inst/include/Matrix.h, inst/include/Matrix_stubs.c, src/CHMfactor.c, src/CHMfactor.h, man/CHMfactor-class.Rd, src/init.c, tests/factorizing.R: Log-determinant of the parent matrix from a CHMfactor object as chm_factor_ldetL2; documentation and support 2008-03-15 Martin Maechler * R/dsCMatrix.R: enable determinant(, .) via chol(.) 2008-03-14 Martin Maechler * R/dsCMatrix.R: setAs(., "dsCMatrix") but with a deprecation warning. 2008-03-13 Martin Maechler * DESCRIPTION (Version, Date): bug-fix release 0.999375-6 * R/diagMatrix.R (diag2tT.u): utility to be smarter in o setAs(., "nMatrix") added. * R/diagMatrix.R (Diagonal): Diagonal(4, x=3) now works too * R/Auxiliaries.R (.diagU2N): more careful coercion in 2 steps new argument 'checkDense = FALSE'. 2008-03-07 Martin Maechler * src/dgeMatrix.c (dgeMatrix_exp): fix the octave-origin bug in the back-permutation of the matrix exponential. * tests/matr-exp.R: test the fix (for an example where expm() was wrong). * DESCRIPTION (Date,Version): ready to release 0.999375-5 * tests/simple.R: testing diagN2U 2008-03-06 Martin Maechler * R/Auxiliaries.R (diagN2U): interface to * src/Csparse.c (Csparse_diagN2U): .Call()able SEXP version of chm_diagN2U() * tests/matprod.R: test for triangularity preserving %*% 2008-03-05 Martin Maechler