[Pkg-octave-devel] Bug#706376: Bug#706376: Bug#706376: octave: sparse matrix n*2^16

Jordi Gutiérrez Hermoso jordigh at octave.org
Mon Apr 29 21:07:16 UTC 2013

On 29 April 2013 16:46, Michael D. Godfrey <michaeldgodfrey at gmail.com>
> > On 04/29/2013 03:12 PM, Jordi Gutiérrez Hermoso wrote:
> >
> > On 29 April 2013 14:50, Ed Meyer <eem2314 at gmail.com> wrote:
> >
> >> I'm not proposing using anything but octave_idx_type for indexing
> >> or changing the return type of numel() - I just question why
> >> numel() is used for sparse matrices. It should be irrelevant for
> >> anything but ccs2full().
> >
> > The numel function is just one example where this sparse matrix
> > with big dimensions broke. Sparse matrices this large can still
> > break in other ways. Furthermore, what do you propose to do if
> > numel is called for sparse matrices, despite your suggestion to
> > not call it for sparse matrices? A lot of code out there already
> > assumes that you can treat a sparse matrix like any other matrix.
> > I don't think numel is to blame.

> Will 64bit Octave "automatically" fix this?

Yes, if you mean compiling Octave with 64-bit indexing. Or at least it
will delay the problem until someone tries to create a matrix of
size 2^64 by 2^64.

- Jordi G. H.

More information about the Pkg-octave-devel mailing list