SCM

SCM Repository

[matrix] Diff of /pkg/ChangeLog
ViewVC logotype

Diff of /pkg/ChangeLog

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1654, Fri Oct 27 16:58:15 2006 UTC revision 1845, Wed May 16 21:24:41 2007 UTC
# Line 1  Line 1 
1    2007-05-16  Douglas Bates  <bates@stat.wisc.edu>
2    
3            * R/CHMfactor.R: Added solve methods for a CHMfactor object.
4    
5    2007-05-16  Martin Maechler  <maechler@stat.math.ethz.ch>
6    
7            * R/Auxiliaries.R (sparsapply): new utility, much faster than
8            tapply1() for large sparse matrices.
9    
10    2007-05-15  Martin Maechler  <maechler@stat.math.ethz.ch>
11    
12            * R/Matrix.R (dim<-): reshape now via sparseVector.
13    
14            * R/sparseVector.R: methods and function for
15            * R/AllClass.R: new "sparseVector" class and daughters.
16            * NAMESPACE: export new classes
17    
18    2007-05-14  Martin Maechler  <maechler@stat.math.ethz.ch>
19    
20            * DESCRIPTION (Version): 0.99875-1
21            * src/Makefile.win: also remove Lapack code from here (cf. 04-25).
22    
23    2007-05-11  Martin Maechler  <maechler@stat.math.ethz.ch>
24    
25            * R/Tsparse.R ([, Tsparse): fix last case: *duplicated*, symmetric
26            indexing
27            * tests/indexing.R: test set for that.
28    
29    2007-05-08  Martin Maechler  <maechler@stat.math.ethz.ch>
30    
31            * R/Tsparse.R (replTmat): fix the case of *duplicated* index
32            entries.
33            * tests/indexing.R(out): add regression test for it
34    
35    2007-04-30  Martin Maechler  <maechler@stat.math.ethz.ch>
36    
37            * R/(l(dense|sparse))?Matrix.R (!): use 'x', not 'e1' as argument
38            name for "!" method definitions.
39    
40    2007-04-26  Martin Maechler  <maechler@stat.math.ethz.ch>
41    
42            * R/Tsparse.R (intI): new utility, used for "[" :
43                    Cleanup up there, and fixes for duplicated indices - more TODO!
44    
45            * tests/indexing.R(out): more tests
46    
47    2007-04-25  Douglas Bates  <bates@stat.wisc.edu>
48    
49            * DESCRIPTION,src/Makefile: require R>= 2.5.0 and remove Lapack
50            code that is now part of the R Lapack library.
51    
52            * src/init.c,inst/include/{Matrix_stubs.c,cholmod.h}:export
53            cholmod_factorize_p (used in lme4 for GLMMs and NLMMs).
54    
55    2007-04-21  Martin Maechler  <maechler@stat.math.ethz.ch>
56    
57            * R/Matrix.R (image): method for all Matrices, not just sparse ones.
58    
59    2007-04-17  Martin Maechler  <maechler@stat.math.ethz.ch>
60    
61            * R/Auxiliaries.R (tapply1): unname(.) -> colSums() etc don't end
62            up with extraneous names '0'...'<n-1>'
63    
64    2007-04-12  Martin Maechler  <maechler@stat.math.ethz.ch>
65    
66            * R/dgTMatrix.R (mat2dgT): care about NAs
67    
68    2007-04-11  Martin Maechler  <maechler@stat.math.ethz.ch>
69    
70            * R/kronecker.R: triangularity preserving methods
71    
72    2007-03-27  Martin Maechler  <maechler@stat.math.ethz.ch>
73    
74            * R/kronecker.R: new file collecting kronecker() methods in one
75                    place.  Goal: become much faster!
76    
77    2007-03-23  Martin Maechler  <maechler@stat.math.ethz.ch>
78    
79            * src/dtCMatrix.c (dtCMatrix_solve): use the new code from Kasper
80            Kristensen based cs_spsolve() instead of _lsolve & _usolve which
81            can be much faster.
82    
83            * tests/matprod.R: add regression tests for these (upper & lower).
84    
85    2007-03-19  Martin Maechler  <maechler@stat.math.ethz.ch>
86    
87            * R/Matrix.R (diff): method for our Matrices.
88    
89            * R/sparseMatrix.R (isDiagonal): check dim()!
90    
91    2007-03-17  Martin Maechler  <maechler@stat.math.ethz.ch>
92    
93            * R/Matrix.R (dim<-): new method for "reshape()" built on a
94            proposal from Tamas Papp.
95    
96    2007-03-16  Martin Maechler  <maechler@stat.math.ethz.ch>
97    
98            * R/AllGeneric.R: remove all  if(!isGeneric(.))  clauses
99    
100            * R/zzz.R (.onLoad, .onUnload): do *not* leave bind_activation(TRUE);
101            rather define and export  cBind() and rBind() only.
102            --> useRs *must* change code that used to have cbind()/rbind() !!
103    
104            * R/bind.R: change tests from cbind() to cBind() and similarly to rBind()
105            * R/bind.Rout.save: ditto
106    
107    2007-02-16  Douglas Bates  <bates@r-project.org>
108    
109            * DESCRIPTION (Date, Version): 0.9975-11 with new date
110            * src/dgCMatrix.c (R_to_CMatrix, compressed_to_TMatrix): remove
111            const modifier on declaration of the array 'valid'
112    
113    2007-02-12  Douglas Bates  <bates@stat.wisc.edu>
114    
115            * R/CHMfactor.R: Add image method (coercion to sparseMatrix).
116    
117    2007-02-05  Martin Maechler  <maechler@stat.math.ethz.ch>
118    
119            * DESCRIPTION (Date, Version): 0.9975-10 with new date.
120    
121            * R/Ops.R (Arith): make sure Csparse o Csparse also works for e.g. ntCMatrix
122            * tests/simple.R: test the above and some of these coercions:
123    
124            * R/nsparseMatrix.R: coercing "Csparse" to "lsparseMatrix";
125                    be careful to avoid infinite recursion, using new coercions in
126            * R/ngCMatrix.R and nsC... and ntC...
127            * R/lsparseMatrix.R: ditto
128    
129            * R/SparseM-conv.R: more conversion, notably for triplet matrices.
130    
131            * src/dgCMatrix.c (R_to_C_Matrix): port Doug's fix and
132            * R/Rsparse.R: reactivate .Call()s
133    
134            * tests/Class+Meth.R: a bit more on actual classes
135    
136    2007-02-04  Douglas Bates  <bates@stat.wisc.edu>
137    
138            * src/dgCMatrix.c (compressed_to_TMatrix): fix memory bug using strdup()
139    
140    2007-02-03  Martin Maechler  <maechler@stat.math.ethz.ch>
141    
142            * DESCRIPTION (Version): 0.9975-10 to upload
143    
144            * tests/Class+Meth.R (tstMatrixClass): require coercions to
145            specific classes less unconditionally.
146    
147            * R/Auxiliaries.R: get rid of as_Tsparse() and as_Rsparse()
148    
149            * R/Tsparse.R (triu): etc, use as(*, "TsparseMatrix") instead of as_Tsparse()
150    
151            * R/Rsparse.R (.R.2.T): R-level workaround using compressed_to_TMatrix.
152            * R/Rsparse.R (.R.2.C): R-level workaround since C-level
153            R_to_CMatrix segfaults on one platform.
154            Eliminate most coercion method to *specific* classes, and replace
155            with virtual classes coercions.
156    
157    
158    2007-02-01  Martin Maechler <maechler@stat.math.ethz.ch>
159    
160            * src/init.c: export the CHM...._validate() placeholders, since
161            they *are* called.
162    
163            * tests/Class+Meth.R (classCanCoerce): and starting to test
164            all as(<from>, <to>)
165    
166    2007-01-30  Martin Maechler  <maechler@stat.math.ethz.ch>
167    
168            * R/Tsparse.R ([): more care when subsetting triangular Tsparse
169            * tests/indexing.R: tested now
170            * tests/indexing.Rout.save: updated
171    
172            * src/Csparse.c (Csparse_to_dense): use Rkind = -1 for PATTERN to
173            * src/chm_common.c (chm_dense_to_SEXP): return "ngeMatrix" when
174            appropriate.
175    
176            * NAMESPACE: export a trivial
177            * R/Matrix.R: drop() Matrix-method
178    
179            * R/AllClass.R: moved all prototypes to virtual super classes.
180    
181            * R/Rsparse.R: many more coercions to have less exceptions in
182            * tests/Class+Meth.R:
183    
184            * R/Ops.R (Compare): tweak for case with NA
185            * tests/simpl.R: hence another 'FIXME' eliminated
186    
187    2007-01-29  Martin Maechler  <maechler@stat.math.ethz.ch>
188    
189            * R/diagMatrix.R (solve): the obvious methods for diagonalMatrix
190            objects.
191    
192            * tests/Class+Meth.R (tstMatrixClass): now testing diag(), nnzero(),
193              and more of "!", "&", "|", all, any; coercions
194    
195            * R/Rsparse.R: many coercions (which enable quite a few other
196            methods), thanks to enhancements in
197            * src/dgCMatrix.c (R_to_CMatrix): new, and
198            * src/dgCMatrix.c (compressed_to_TMatrix): now for (d,l,n) ,
199            symmetric & triangular and ..RMatrix objects.
200    
201            * src/TMatrix_as.c (Matrix_T_as_DENSE,Matrix_T_as_GENERAL):
202            renamed file from src/dsTMatrix.c;
203            now dealing with symmetric and triangular Tsparse coercions, both
204            to dense and general.
205    
206    2007-01-27  Martin Maechler  <maechler@stat.math.ethz.ch>
207    
208            * src/dsTMatrix.c: has now "l" and "n" methods besides the "d" ones.
209    
210            * R/Ops.R (Arith): <dgCMatrix> o <numeric> now remains sparse
211            where sensible when the <numeric> is of length > 1.
212    
213    2007-01-26  Martin Maechler  <maechler@stat.math.ethz.ch>
214    
215            * R/Matrix.R ([<-): for M[] <- value: fix length
216    
217    2007-01-25  Martin Maechler  <maechler@stat.math.ethz.ch>
218    
219            * R/Auxiliaries.R (n2l_Matrix): new, to be used in
220            * R/ndenseMatrix.R: new coercions n* -> l*
221    
222    2007-01-22  Martin Maechler  <maechler@stat.math.ethz.ch>
223    
224            * R/triangularMatrix.R: new file; simple triu() and tril() methods.
225    
226            * R/Ops.R ("Logic"): and other "Ops", many updates
227    
228    2007-01-18  Martin Maechler  <maechler@stat.math.ethz.ch>
229    
230            * src/Mutils.h (SET_DimNames): new utility
231    
232            * R/Auxiliaries.R (nnzero): improved and now exported via
233    
234            * NAMESPACE: + nnzero(); length() == prod(dim(.)) method for all "Matrix" objects
235    
236    2007-01-17  Martin Maechler  <maechler@stat.math.ethz.ch>
237    
238            * R/diagMatrix.R (!): fix typo.
239    
240    2007-01-16  Martin Maechler  <maechler@stat.math.ethz.ch>
241    
242            * R/Auxiliaries.R (as_Csparse): and quite a few others:
243                    allow to pass class definition --> speedup
244            * R/sparseMatrix.R: apply the above
245    
246            * R/Csparse.R: coercion Csparse* to dense* now preserves shape properties.
247    
248            * src/Mutils.h (mMatrix_as_geMatrix): new, based on
249            * src/Mutils.c (dup_mMatrix_as_geMatrix): new; generalization of
250            old dup_mMatrix_as_dgeMatrix), eliminating a long-standing "FIXME".
251    
252            * src/dense.c (dense_to_Csparse): use new mMatrix_as_geMatrix()
253    
254            * R/denseMatrix.R (.dense2C): based on dense_to_Csparse: name it,
255            and use it for "sparse*" as well, since it's faster than the
256            as_Csparse(.) way.
257    
258    2007-01-15  Martin Maechler  <maechler@stat.math.ethz.ch>
259    
260            * R/Ops.R ("Logic"): more methods, notably an  <lgC> o <lgC> one.
261    
262    2007-01-12  Martin Maechler  <maechler@stat.math.ethz.ch>
263    
264            * R/Tsparse.R (.TM.repl.i.2col): new internal function to be used
265            as method for   M[ ij ] <- v
266            * R/Csparse.R:: go via Tsparse for "M[ij] <- v"
267    
268            * R/Ops.R: "Compare" for (C|R)sparse: need pointer slot for all
269            FALSE answer
270    
271            * R/Csparse.R (replCmat): fix the "all non-zero" case with reordering
272            * tests/indexing.R: test it, and some of the above
273    
274    2007-01-05  Martin Maechler  <maechler@stat.math.ethz.ch>
275    
276            * R/Auxiliaries.R (is_duplicatedT): new utility
277    
278    2007-01-05  Douglas Bates  <bates@stat.wisc.edu>
279    
280            * src/init.c (R_init_Matrix): export cholmod_scale
281    
282    2006-12-30  Martin Maechler  <maechler@stat.math.ethz.ch>
283    
284            * R/zzz.R (tmp): for R >= 2.5.0, extend formals of our
285            base::as.matrix to (x, ...)
286    
287    2006-12-28  Martin Maechler  <maechler@stat.math.ethz.ch>
288    
289            * R/Ops.R ("Arith" etc): move almost all "Ops" methods to new R
290            file; start using "Logic", hence
291    
292            * DESCRIPTION (Depends): R >= 2.4.1 (since we want "Logic")
293            * NAMESPACE: import and export "Logic"
294    
295    2006-12-27  Martin Maechler  <maechler@stat.math.ethz.ch>
296    
297            * src/zpotfr.f and dependencies: use LAPACK 3.1 version
298            only needed previously to R version 2.5.0.
299    
300    2006-12-26  Martin Maechler  <maechler@stat.math.ethz.ch>
301    
302            * DESCRIPTION (Date, Version): 0.9975-8, ready for release
303    
304            * R/Tsparse.R (replTmat): fix subassignment of triangular
305            * R/Csparse.R (replCmat): ditto
306            * tests/indexing.R: more tests, incl the above fix
307    
308    2006-12-23  Martin Maechler  <maechler@stat.math.ethz.ch>
309    
310            * R/Auxiliaries.R (drop0): extend for non CSparse
311            * R/Auxiliaries.R (diagU2N): should work for all sparseMatrix
312    
313            * src/Csparse.c (Csparse_to_Tsparse, Csparse_general_to_symmetric):
314                    use uplo correctly (!); other places: use uplo_P() macro
315    
316            * R/Csparse.R (replCmat): call diagU2N() when needed
317            * R/Tsparse.R (replTmat): ditto
318    
319            * src/dtCMatrix.c (tCMatrix_validate): new
320            * src/dtTMatrix.c (tTMatrix_validate): new, used in
321            * R/AllClass.R: for validity of dtC, dsC, and dtT, dsT.
322    
323            * R/diagMatrix.R (replDiag): to use in [<-
324    
325    2006-12-22  Martin Maechler  <maechler@stat.math.ethz.ch>
326    
327            * R/Auxiliaries.R (as_Csparse2, as_geSimpl): new functions;
328                    also more general diagU2N().
329    
330    2006-12-21  Martin Maechler  <maechler@stat.math.ethz.ch>
331    
332            * R/bind2.R: new file for all cbind2(), rbind() methods moved here
333            from R/<foo>Matrix.R files. Better diagonal & improved sparse methods.
334    
335    2006-12-20  Martin Maechler  <maechler@stat.math.ethz.ch>
336    
337            * tests/bind.R: a few more cases
338            * R/Auxiliaries.R (.M.kind): also work for atomic vectors
339    
340            * R/denseMatrix.R (cbind2/rbind2): moved here (and generalized) from
341            * R/ddenseMatrix.R (cbind2/rbind2)
342            * R/Tsparse.R (replTmat): final(?!) fix for "[<-" ..
343            * tests/indexing.R
344            * tests/indexing.Rout.save: updated
345    
346    2006-12-18  Martin Maechler  <maechler@stat.math.ethz.ch>
347    
348            * R/Tsparse.R (replTmat): fixed a remaining "[<-" bug in
349            * tests/indexing.R
350    
351    2006-12-15  Martin Maechler  <maechler@stat.math.ethz.ch>
352    
353            * R/sparseMatrix.R (prSpMatrix): "." alignment much improved:
354            align with proper position of "0", i.e., right for integers.
355            argument 'align' by default is "fancy".
356    
357    2006-12-14  Martin Maechler  <maechler@stat.math.ethz.ch>
358    
359            * R/sparseMatrix.R: delegate "Compare" to "Csparse.."
360            * R/Csparse.R: and fix "Compare" for more cases.
361            * tests/Class+Meth.R: test some of these (m == m, m != m)
362    
363    2006-12-13  Martin Maechler  <maechler@stat.math.ethz.ch>
364    
365            * R/lsparseMatrix.R: all() and any() methods
366            * R/ldenseMatrix.R:      ditto
367            * NAMESPACE, R/Matrix.R: ditto
368    
369            * man/all-methods.Rd: document them minimally
370    
371            * tests/simple.R: add a few examples for these
372    
373    2006-12-11  Martin Maechler  <maechler@stat.math.ethz.ch>
374    
375            * R/Tsparse.R ([): fix long standing typo in symmetric case
376            * man/dsCMatrix-class.Rd: add example exhibiting the above case
377    
378    2006-12-10  Douglas Bates  <bates@stat.wisc.edu>
379    
380            * src/CHMfactor.c (CHMfactor_to_sparse): change LDL factorization
381            to LL before converting to a sparse matrix. (The LDL form can be
382            converted to a sparse matrix but it is implicitly a unit
383            triangular matrix and a diagonal matrix overwritten on the diagonal.)
384    
385    2006-12-09  Douglas Bates  <bates@stat.wisc.edu>
386    
387            * src/chm_common.c (chm_factor_to_SEXP): allocate and fill the Dim slot.
388    
389    2006-12-08  Douglas Bates  <bates@stat.wisc.edu>
390    
391            * DESCRIPTION (Version): updated -> release 0.9975-7
392            * src/{init.c,chm_common.c}, inst/include/*: export cholmod_analyze_p
393    
394    2006-11-30  Martin Maechler  <maechler@stat.math.ethz.ch>
395    
396            * R/diagMatrix.R (%*%): write a direct [diag o Csparse] method
397    
398    2006-11-29  Douglas Bates  <bates@stat.wisc.edu>
399    
400            * src/dgeMatrix.c (dgeMatrix_solve): Check error code from dgetri.
401            * tests/dg_Matrix.R: Add Barry Rowlingson's test of a matrix that
402            is exactly singular.
403    
404    2006-11-07  Martin Maechler  <maechler@stat.math.ethz.ch>
405    
406            * DESCRIPTION (Date): updated -> release 0.9975-6
407    
408    2006-11-06  Martin Maechler  <maechler@stat.math.ethz.ch>
409    
410            * R/Csparse.R (replCmat): symmetric indexing of symmetric matrix
411            now returns symmetric.
412    
413            * R/zzz.R ("diag<-"): replace "diag<-" in base for R <= 2.4.x
414    
415            * R/Matrix.R (.M.sub.i.2col):  new, for  M[ cbind(i,j) ] indexing.
416            * R/Matrix.R (.M.repl.i.2col): new, for  M[ cbind(i,j) ] <- value
417    
418            * R/Auxiliaries.R (.type.kind): added
419    
420    2006-11-04  Martin Maechler  <maechler@stat.math.ethz.ch>
421    
422            * src/cs.[ch]: updated to CSparse Version 2.0.3 by simply
423            running src/CSparse_install.sh
424    
425            * R/denseMatrix.R: "[": keep symmetric on symmetric indexing.
426    
427    2006-11-03  Martin Maechler  <maechler@stat.math.ethz.ch>
428    
429            * src/dsCMatrix.c (dsCMatrix_Csparse_solve): new
430            * R/dsCMatrix.R (solve): "fully-sparse" using the above.
431    
432            * R/AllClass.R: "pMatrix" now also inherits from "generalMatrix"
433    
434            * tests/Class+Meth.R (tstMatrixClass): now assure
435            the (important in method programming) property :
436    
437                    ###>>  Every "Matrix" is either
438                    ###>>  "general*", "symmetric*", "triangular*" or "diagonal*"
439    
440                    (where "*" stands for "Matrix")
441    
442            * R/Auxiliaries.R (diagU2N): now .Call()s Csparse_diagU2N for <Csparse>
443    
444            * R/dMatrix.R (Compare(<dMatrix>,<dMatrix>): update and
445            * tests/validObj.R: checks for "comparison"
446    
447            * R/sparseMatrix.R ([): improved indexing for sparse;
448                    trying to keep <symmetric>[ n, n] symmmetric
449    
450            * tests/indexing.R: indexing for logical sparse now ok
451    
452    
453    2006-11-02  Martin Maechler  <maechler@stat.math.ethz.ch>
454    
455            * src/Tsparse.c: use xTsparse_validate() , and hence remove
456    
457            * src/{ltC,lsC,lgT}Matrix.[ch]: removed
458    
459    2006-11-02  Martin Maechler  <maechler@stat.math.ethz.ch>
460    
461            * R/AllClass.R (Matrix-class): check length of dimnames in validity.
462    
463            * tests/simple.R: validObject() checking the above.
464    
465            * src/dgCMatrix.c (xCMatrix_validate): new, small and simple,
466            replacing both dgCMatrix_validate and lgCM*.
467    
468            * src/Csparse.c (Csparse_dense_prod, etc): do not lose dimnames;
469            fix dimnames setting in other places.
470            * src/chm_common.c (chm_dense_to_SEXP): now can pass dimnames
471    
472    2006-11-01  Martin Maechler  <maechler@stat.math.ethz.ch>
473    
474            * R/Csparse.R,src/Csparse.c, etc: tcrossprod(<CsparseM>,<CsparseM>)
475    
476            * R/sparseMatrix.R (isSymmetric): drop 'factors' slot for
477            symmetry test, via
478            * R/Auxiliaries.R (.as.dgC.0.factors):
479    
480    2006-11-01  Douglas Bates  <bates@r-project.org>
481    
482            * R/Csparse.R,src/Csparse.c,tests/matprod.R,
483            man/CsparseMatrix-class.Rd: crossprod(<CsparseMatrix>,
484            <CsparseMatrix>) added
485    
486    2006-10-30  Martin Maechler  <maechler@stat.math.ethz.ch>
487    
488            * tests/matprod.R: add a variation of Harri's example
489    
490            * R/dsparseMatrix.R: fix crossprod(<dsparse>, <dge>) to *not*
491            recursive infinitely.
492    
493            * R/dgCMatrix.R: + solve(<sparse>, <sparse>)
494    
495            * tests/indexing.R: add test for the "<" bug fixed 10-27 in R/dMatrix.R
496    
497    2006-10-28  Martin Maechler  <maechler@stat.math.ethz.ch>
498    
499            * tests/Class+Meth.R (tstMatrixClass): more: use non-trivial
500            matrix if possible; test m+m == 2*m; now test dgRMatrix.
501            * R/dgRMatrix.R (.to.dgR): a few more coercions, in order to
502            satisfy the above test.
503    
504  2006-10-27  Martin Maechler  <maechler@stat.math.ethz.ch>  2006-10-27  Martin Maechler  <maechler@stat.math.ethz.ch>
505    
506            * R/Matrix.R (Ops):  <Matrix> o <matrix>  method added
507    
508          * R/dgCMatrix.R: solve(a, b="missing") based on          * R/dgCMatrix.R: solve(a, b="missing") based on
509          * src/dgCMatrix.c (dgCMatrix_matrix_solve): extend to work with          * src/dgCMatrix.c (dgCMatrix_matrix_solve): extend to work with
510          RHS = NULL.          RHS = NULL.

Legend:
Removed from v.1654  
changed lines
  Added in v.1845

root@r-forge.r-project.org
ViewVC Help
Powered by ViewVC 1.0.0  
Thanks to:
Vienna University of Economics and Business Powered By FusionForge