Re: Memory exhausted when using \

[Mathieu Dubois]

Le 19/07/2012 16:31, Martin Helm a écrit :

> Am 19.07.2012 15:55, schrieb Mathieu Dubois:
> > Le 17/07/2012 13:16, Martin Helm a écrit :
> > > Am 17.07.2012 13:07, schrieb Mathieu Dubois:
> > > >      Hi Martin,
> > > >
> > > > Sorry for the delay (I have a lot of other problems :).
> > > >
> > > > I have installed lsqrSOL and used your snippet (except that I set atol
> > > > = btol = 1e6). It seems to work (but it's rather slow compared to the
> > > > matlab implementation which use \).
> > > >
> > > > I will try to make sure that the results are similar.
> > > >
> > > > Thanks for your help,
> > > > Mathieu
> > > Do not forget what Jordi posted, maybe you simply ran into a bug in your
> > > old octave version and the \ operator works for you with that problem
> > > you try to solve when you update to the current version 3.6
> > Hello,
> >
> > I have tried today with octave 3.6.1 (on a 64 bits machine but without
> > 64 bits adressing) and ran into the same problem.
> >
> > I have tried to recompile octave with --enable-64 but without success.
> >
> > Mathieu
> Could you upload the matrix and the vector somewhere, it is otherwise
> difficult to see why you get that problem. Since it is a sparse matrix
> and we do so far not know anything about its structure or how sparse or
> full it is?
> If that is no problem for you to make its content public

Sorry for the dealy (it took a long time to rebuild octave :).

I have no problem to make the data available but as I said the code is not mine.

You can download the dump at:  
http://www.ensta-paristech.fr/~mdubois/memory-exhausted

The command I used was:
save memory-exhausted A b
right before the line:
tic; X = A \ b; toc;
(executing this line raises the error)

Mathieu