Implement sparse * diagonal and diagonal * sparse operations, double-prec only. Date: Sun, 8 Mar 2009 16:28:18 -0400 These preserve sparsity, so eye(5) * sprand (5, 5, .2) is *sparse* and not dense. This may affect people who use multiplication by eye() rather than full(). The liboctave routines do *not* check if arguments are scalars in disguise. There is a type problem with checking at that level. I suspect we want diag * "sparse scalar" to stay diagonal, but we have to return a sparse matrix at the liboctave. Rather than worrying about that in liboctave, we cope with it when binding to Octave and return the correct higher-level type. The implementation is in Sparse-diag-op-defs.h rather than Sparse-op-defs.h to limit recompilation. And the implementations are templates rather than macros to produce better compiler errors and debugging information.

If you have trouble transferring Octave from this site, you might try
one of those listed below.  They mirror the /pub/octave directory on
ftp.octave.org

FTP:

  site                          directory
  ----                          ---------
  ftp.eos.hokudai.ac.jp         /pub/GNU/misc/octave
  ftp.u-aizu.ac.jp              /pub/SciEng/numanal/Octave
  neacm.fe.up.pt                /pub/octave

HTTP:

  url
  ---
  http://neacm.fe.up.pt/pub/octave
  http://octave.oss-mirror.org
  http://octave.linux-mirror.org


John W. Eaton
jwe@bevo.che.wisc.edu
University of Wisconsin-Madison
Department of Chemical & Biological Engineering

Wed, 31 Oct 2007 16:32:44 EDT