On 20111115@01:42, Michele Martone wrote: > On 20111114@10:42, Michele Martone wrote: > > On 20111113@23:09, Carlo de Falco wrote: > > > ... > > Dear Alex, Carlo, > > > > Thanks, you catched another bug(s)! > > > > Indeed, sparse-sparse product is not a well tested feature: I wrote it > > more for the challenge/fun of having it, to be honest. > > It's parallel, however I wonder how much efficient. > > > > By the way, it seems like even: > > A=[1];B=sparsersb(A);C=sparsersb(A);D=B*C; > > is enough to make it crash. > > I'll definitely have to investigate --- please do not report about it > > until I propose some update about it. > > > > I wonder how much critical this operation is in linear system solving > > applications. > > This information would be useful in order to settle priorities in > > bugfixing; do you have an idea about this ? > > Alex, Carlo, > > I uploaded an updated tarball http://www.ipp.mpg.de/~mima/librsb.tar.gz > where it seems safe to call: > > load test.mat % http://modb.oce.ulg.ac.be/mediawiki/upload/Alex/test.mat > [i,j,s]=find(A); > A2 = sparsersb(i,j,s); > [i,j,s]=find(extx); > extx2 = sparsersb(i,j,s); > extx2*A2; > > A=[1];B=sparsersb(A);C=sparsersb(A);D=B*C; > > However, I notice that handling of empty matrices (e.g.:[0]) is > not yet stable.. > ...
Update: now in sparsersb, error() is used to handle critical errors. Therefore now, when issuing 'sparsersb([0])' (for now, librsb does not support empty matrices), with librsb API functions failing in allocating a matrix, an exception (I'm guessing it is such) is thrown by error(), and it is safe to continue octave execution. p.s.: now one may ./configure sparsersb with the --with-verbosity-level=1 option, allowing the tracking all of sparsersb's methods/functions executions, whose name is dumped on screen (this is a debug option).
pgpfovTp7YgsR.pgp
Description: PGP signature
------------------------------------------------------------------------------ RSA(R) Conference 2012 Save $700 by Nov 18 Register now http://p.sf.net/sfu/rsa-sfdev2dev1
_______________________________________________ Octave-dev mailing list Octave-dev@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/octave-dev
