On 13 Nov 2011, at 23:02, c. wrote: > > On 13 Nov 2011, at 22:20, Alexander Barth wrote: > >> Dear Michele, >> >> I would like to try the suggestion of Filippone, but I sumbled upon a >> a segmentation fault for the following matrix multiplication: >> >> load test.mat % http://modb.oce.ulg.ac.be/mediawiki/upload/Alex/test.matbut >> [i,j,s]=find(A); >> A2 = sparsersb(i,j,s); >> [i,j,s]=find(extx); >> extx2 = sparsersb(i,j,s); >> extx2*A2; >> >> panic: Segmentation fault -- stopping myself... >> attempting to save variables to `octave-core'... >> save to `octave-core' complete >> Segmentation fault >> >> Any help would be appreciated, >> Cheers, >> Alex > > Alex, > FYI, your example is working fine for me: > >>> load test.mat >>> A2 = sparsersb (A); >>> extx2 = sparsersb (extx); >>> x2 = A2 * extx2; >>> [ia, ja, va] = find (A); >>> [ie, je, ve] = find (extx); >>> A3 = sparsersb (ia, ja, va); >>> extx3 = sparsersb (ie, je, ve); >>> x3 = A3 * extx3; > > What version of Octave are you using? how did you install sparsersb? > c.
OTH, if I run the example repeatedly I also get a crash, eventually: $ RSB_USER_SET_MEM_HIERARCHY_INFO="L2:4/64/3M,L1:8/64/32K" octave -q >> load test.mat >> A2 = sparsersb (A); >> extx2 = sparsersb (extx); >> x2 = A2 * extx2; >> x2 = A2 * extx2; >> x2 = A2 * extx2; c. ------------------------------------------------------------------------------ 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
