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

2013-04-29 Thread Ed Meyer
On Mon, Apr 29, 2013 at 12:12 PM, Jordi Gutiérrez Hermoso < jord...@octave.org> wrote: > On 29 April 2013 14:50, Ed Meyer 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 ma

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

2013-04-29 Thread Michael D. Godfrey
On 04/29/2013 03:12 PM, Jordi Gutiérrez Hermoso wrote: On 29 April 2013 14:50, Ed Meyer 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 anyth

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

2013-04-29 Thread Jordi Gutiérrez Hermoso
On 29 April 2013 14:50, Ed Meyer 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 e

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

2013-04-29 Thread Ed Meyer
On Mon, Apr 29, 2013 at 11:07 AM, Jordi Gutiérrez Hermoso < jord...@octave.org> wrote: > On 29 April 2013 14:00, Jordi Gutiérrez Hermoso > wrote: > > And yes, sparse matrices can be indexed by a single index instead of > > two, like any other matrix. Internally in Octave's source, the > > assumpt

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

2013-04-29 Thread Jordi Gutiérrez Hermoso
On 29 April 2013 14:00, Jordi Gutiérrez Hermoso wrote: > And yes, sparse matrices can be indexed by a single index instead of > two, like any other matrix. Internally in Octave's source, the > assumption that a single index of octave_idx_type is available is used > throughout. The specific case w

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

2013-04-29 Thread Jordi Gutiérrez Hermoso
On 29 April 2013 13:21, Ed Meyer wrote: > > > On Mon, Apr 29, 2013 at 9:47 AM, Jordi Gutiérrez Hermoso > wrote: >> >> >> All matrices need to be linearly indexable, and of course, this is how >> they are actually stored in memory, as a single long array indexed by >> a single index. Thus, the tot

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

2013-04-29 Thread Ed Meyer
On Mon, Apr 29, 2013 at 9:47 AM, Jordi Gutiérrez Hermoso wrote: > > All matrices need to be linearly indexable, and of course, this is how > they are actually stored in memory, as a single long array indexed by > a single index. Thus, the total number of indexable elements of a > matrix can't be

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

2013-04-29 Thread Jordi Gutiérrez Hermoso
On 29 April 2013 12:40, Ed Meyer wrote: > > > On Mon, Apr 29, 2013 at 6:10 AM, Jordi Gutiérrez Hermoso > wrote: >> >> On 29 April 2013 06:25, Miroslaw Kwasniak >> wrote: >> > it's something wrong whith sparse matrices A(n,n) when n is a multiple >> > of 65536=2^16. >> > >> > Demonstration code =

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

2013-04-29 Thread Ed Meyer
On Mon, Apr 29, 2013 at 6:10 AM, Jordi Gutiérrez Hermoso wrote: > On 29 April 2013 06:25, Miroslaw Kwasniak > wrote: > > it's something wrong whith sparse matrices A(n,n) when n is a multiple > > of 65536=2^16. > > > > Demonstration code == > > > > for i=1:3;

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

2013-04-29 Thread Jordi Gutiérrez Hermoso
On 29 April 2013 06:25, Miroslaw Kwasniak wrote: > it's something wrong whith sparse matrices A(n,n) when n is a multiple > of 65536=2^16. > > Demonstration code == > > for i=1:3; > for n=i*2^16+(-1:1); > A=spdiags(ones(n,1),0,n,n); > t=trace(A); >