On 13 Nov 2011, at 23:09, Carlo de Falco wrote: > > 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.
And here is the backtrace of the crash: #0 0x16e1fb2a in rsb_do_util_csr_csr_sparse_mul_serial () #1 0x1683ad3f in rsb_spgemm_inner.omp_fn.0 () #2 0x1683a9ee in rsb_spgemm_inner () #3 0x1683b272 in rsb_spgemm_tmp () #4 0x1683c8b7 in rsb_do_matrix_mul () #5 0x1682ac6e in rsb_matrix_mul () #6 0x16802f14 in oct_binop_op_spmul () #7 0x0061b3e4 in do_binary_op (op=octave_value::op_mul, v1=@0xbfffc20c, v2=@0xbfffc204) at ov.cc:1902 #8 0x006d6244 in tree_binary_expression::rvalue1 (this=0x2f6fbe0) at pt-binop.cc:132 #9 0x006d36aa in tree_simple_assignment::rvalue1 (this=0x3a83ad0) at pt-assign.cc:205 #10 0x006e6f63 in tree_evaluator::visit_statement (this=0xfa7ce0, stmt=@0x3974d00) at pt-eval.cc:739 #11 0x007188fb in tree_statement::accept (this=0x3974d00, tw=@0xfa7ce0) at pt-stmt.cc:151 #12 0x006e59e5 in tree_evaluator::visit_statement_list (this=0xfa7ce0, lst=@0x3a56540) at pt-eval.cc:775 #13 0x0008cef3 in tree_statement_list::accept (this=0x3a56540, tw=@0xfa7ce0) at pt-stmt.cc:215 #14 0x004718d2 in main_loop () at toplev.cc:574 #15 0x0040573d in octave_main (argc=6, argv=0xbfffc61c, embedded=0) at octave.cc:938 #16 0x00001f80 in main (argc=6, argv=0xbfffc61c) at main.c:35 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
