Diagonal v. sparse matrices [was Re: deprecated functions]

John W. Eaton jwe at octave.org
Thu Mar 5 01:57:42 CST 2009


On  5-Mar-2009, Jaroslav Hajek wrote:

| On Wed, Mar 4, 2009 at 10:17 PM, John W. Eaton <jwe at octave.org> wrote: 
|
| > What do you mean by "without losing efficiency"?  Does
| >
| >  a = diag (...);
| >  sparse (a);
| >
| > convert A to full first, then sparse?
| 
| ... yes.
| 
| > If so, it seems that should be fairly easy to fix.
| 
| It was really quite straightforward - fix is online.

Thanks.

| > In any case, I think the following change should
| > be made to spdiag.  This is much simpler than having to check
| > arguments in spdiag, and will be efficient if the "sparse (diag)"
| > constructor is fixed.  But the important thing is that for
| > compatibility with previous versions of Octave, spdiag returns a
| > sparse matrix, not a diag matrix.  Do you agree?
| 
| Agreed.

OK, I applied that change.

Thanks,

jwe



More information about the Octave-maintainers mailing list