Mercurial > hg > octave-lyh
diff PROJECTS @ 10527:b4d2080b6df7
Replace nzmax by nnz as needed
author | David Bateman <dbateman@free.fr> |
---|---|
date | Fri, 16 Apr 2010 16:14:45 +0200 |
parents | b9986ce0047c |
children | 83a37dc986f8 |
line wrap: on
line diff
--- a/PROJECTS +++ b/PROJECTS @@ -61,41 +61,21 @@ * Improve QR factorization functions, using idea based on CSPARSE cs_dmsol.m + * Improve QR fqctorization by replace CXSPARSE code with SPQR code, and make + the linear solve return 2-norm solutions for ill-conditioned matrices + based on this new code + * Implement fourth argument to the sprand and sprandn, and addition arguments to sprandsym that the leading brand implements. * Sparse logical indexing in idx_vector class so that something like - "a=sprandn(1e6,1e6,1e-6); a(a<1) = 0" won't cause a memory overflow. - - * Make spalloc (r, c, n) actually create an empty sparse with n - non-zero elements? This allows something like - - sm = spalloc (r, c, n) - for j=1:c - for i=1:r - tmp = foo (i,j); - if (tmp != 0.) - sm (i,j) = tmp; - endif - endfor - endfor - - actually make sense. Otherwise the above will cause massive amounts - of memory reallocation. - - The fact is that this doesn't make sense in any case as the assign - function makes another copy of the sparse matrix. So although spalloc - might easily be made to have the correct behavior, the first assign - will cause the matrix to be resized! There seems to be no simple - way to treat this but a complete rewrite of the sparse assignment - functions... + 'a=sprandn(1e6,1e6,1e-6); a(a<1) = 0' won't cause a memory overflow. * Other missing Functions - symmmd Superseded by symamd - colmmd Superseded by colamd - cholinc - bicg Can this be taken from octave-forge? - - cgs - gmres - lsqr - minres