Re: [OctDev] Panic in sparse matrix code

SourceForge Headquarters 1320 Columbia Street Suite 310 San Diego, CA 92101 +1 (858) 422-6466

Is 2.9.x octave available for Fedora?  It has replacement code for the  
octave-forge sparse functions and the octave-forge functions are no  
longer supported [though you are welcome to provide patches if you  
want].

- Paul

On Sep 21, 2005, at 5:19 PM, Orion Poplawski wrote:

> Running octave 2.1.71 with octave-forge-2005-06.13 on Fedora Core 4.  
> With the attached files:
>
> octave test.m
>
> runs ok.
>
> octave test2.m
>
> aborts with the following.  test is basically test2 with the function  
> call folded into the main code.  The code is matlab code I'm trying to  
> run in octave that I haven't explored very much.
>
> GNU Octave, version 2.1.71 (i686-pc-linux-gnu).
> Copyright (C) 2005 John W. Eaton.
> This is free software; see the source code for copying conditions.
> There is ABSOLUTELY NO WARRANTY; not even for MERCHANTIBILITY or
> FITNESS FOR A PARTICULAR PURPOSE.  For details, type `warranty'.
>
> Additional information about Octave is available at  
> http://www.octave.org.
>
> Please contribute if you find this software useful.
> For more information, visit http://www.octave.org/help-wanted.html
>
> Report bugs to <bu...@oc...> (but first, please read
> http://www.octave.org/bugs.html to learn how to write a helpful  
> report).
>
> *** glibc detected *** octave: double free or corruption (!prev):  
> 0x09ce7748 ***
> ======= Backtrace: =========
> /lib/libc.so.6[0x289124]
> /lib/libc.so.6(__libc_free+0x77)[0x28965f]
> /usr/libexec/octave/2.1.71/site/oct/i686-pc-linux-gnu/octave-forge/ 
> sparse.oct(oct_sparse_free+0x21)[0xd894fb]
> /usr/libexec/octave/2.1.71/site/oct/i686-pc-linux-gnu/octave-forge/ 
> sparse.oct(Destroy_CompCol_Matrix+0x32)[0xdbd9ff]
> /usr/libexec/octave/2.1.71/site/oct/i686-pc-linux-gnu/octave-forge/ 
> sparse.oct(_Z30oct_sparse_Destroy_SuperMatrix11SuperMatrix+0x52)[0xd896 
> 16]
> /usr/libexec/octave/2.1.71/site/oct/i686-pc-linux-gnu/octave-forge/ 
> sparse.oct(_ZN13octave_sparseD0Ev+0x40)[0xd7df18]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN12octave_valueD1Ev+0x44)[0x417712a]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZSt8_DestroyIN9__gnu_cxx17__normal_iteratorIP12octave_ 
> valueSt6vectorIS2_SaIS2_EEEES5_EvT_S8_T0_+0x1a)[0x3f8439e]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZNSt6vectorI12octave_valueSaIS0_EEaSERKS2_+0x1cf)[0x3f 
> 8581f]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN19tree_statement_list4evalEbi+0x14c)[0x4299d54]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN20octave_user_function17do_multi_index_opEiRK17octav 
> e_value_list+0x553)[0x41a6d39]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN20octave_user_function7subsrefERKSsRKSt4listI17octav 
> e_value_listSaIS3_EEi+0x23c)[0x41a7fb4]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN12octave_value7subsrefERKSsRKSt4listI17octave_value_ 
> listSaIS3_EEi+0x9c)[0x4178d9a]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN21tree_index_expression6rvalueEi+0x415)[0x4272a33]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN21tree_index_expression6rvalueEv+0x3f)[0x4271549]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN22tree_simple_assignment6rvalueEv+0x4f)[0x4256755]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN22tree_simple_assignment6rvalueEi+0x603)[0x4259269]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN14tree_statement4evalEbib+0x116)[0x42992ba]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_ZN19tree_statement_list4evalEbi+0x12e)[0x4299d36]
> /usr/lib/octave-2.1.71/ 
> liboctinterp.so(_Z9main_loopRKSs+0xec)[0x40e1c5e]
> /usr/lib/octave-2.1.71/liboctinterp.so(octave_main+0x11e7)[0x443de9d]
> octave(main+0x19)[0x804864d]
> /lib/libc.so.6(__libc_start_main+0xdf)[0x23ad5f]
> octave[0x80485b1]
> ======= Memory map: ========
> 00111000-001a1000 r-xp 00000000 03:02 109433      
> /usr/lib/libfftw3.so.3.0.1
> 001a1000-001a5000 rwxp 00090000 03:02 109433      
> /usr/lib/libfftw3.so.3.0.1
> 001a5000-001ae000 r-xp 00000000 03:02 864932      
> /lib/libnss_files-2.3.5.so
> 001ae000-001af000 r-xp 00008000 03:02 864932      
> /lib/libnss_files-2.3.5.so
> 001af000-001b0000 rwxp 00009000 03:02 864932      
> /lib/libnss_files-2.3.5.so
> 00208000-00222000 r-xp 00000000 03:02 867425     /lib/ld-2.3.5.so
> 00222000-00223000 r-xp 00019000 03:02 867425     /lib/ld-2.3.5.so
> 00223000-00224000 rwxp 0001a000 03:02 867425     /lib/ld-2.3.5.so
> 00226000-00349000 r-xp 00000000 03:02 867426     /lib/libc-2.3.5.so
> 00349000-0034b000 r-xp 00123000 03:02 867426     /lib/libc-2.3.5.so
> 0034b000-0034d000 rwxp 00125000 03:02 867426     /lib/libc-2.3.5.so
> 0034d000-0034f000 rwxp 0034d000 00:00 0
> 00351000-00374000 r-xp 00000000 03:02 867427     /lib/libm-2.3.5.so
> 00374000-00375000 r-xp 00022000 03:02 867427     /lib/libm-2.3.5.so
> 00375000-00376000 rwxp 00023000 03:02 867427     /lib/libm-2.3.5.so
> 00378000-0037a000 r-xp 00000000 03:02 867428     /lib/libdl-2.3.5.so
> 0037a000-0037b000 r-xp 00001000 03:02 867428     /lib/libdl-2.3.5.so
> 0037b000-0037c000 rwxp 00002000 03:02 867428     /lib/libdl-2.3.5.so
> 0037e000-00390000 r-xp 00000000 03:02 111368      
> /usr/lib/libz.so.1.2.2.2
> 00390000-00391000 rwxp 00011000 03:02 111368      
> /usr/lib/libz.so.1.2.2.2
> 00398000-00399000 r-xp 00398000 00:00 0
> 00469000-00477000 r-xp 00000000 03:02 867435      
> /lib/libpthread-2.3.5.so
> 00477000-00478000 r-xp 0000d000 03:02 867435      
> /lib/libpthread-2.3.5.so
> 00478000-00479000 rwxp 0000e000 03:02 867435      
> /lib/libpthread-2.3.5.so
> 00479000-0047b000 rwxp 00479000 00:00 0
> 0048e000-004b5000 r-xp 00000000 03:02 110766      
> /usr/lib/libreadline.so.5.0
> 004b5000-004b9000 rwxp 00027000 03:02 110766      
> /usr/lib/libreadline.so.5.0
> 004b9000-004ba000 rwxp 004b9000 00:00 0
> 004bc000-004fe000 r-xp 00000000 03:02 128642  
> /usr/lib/octave-2.1.71/libcruft.so.2.1.71
> 004fe000-00502000 rwxp 00041000 03:02 128642  
> /usr/lib/octave-2.1.71/libcruft.so.2.1.71
> 00502000-00503000 rwxp 00502000 00:00 0
> 00505000-0055c000 r-xp 00000000 03:02 103091  
> /usr/lib/libgfortran.so.0.0.0
> 0055c000-0055e000 rwxp 00056000 03:02 103091  
> /usr/lib/libgfortran.so.0.0.0
> 0055e000-00562000 rwxp 0055e000 00:00 0
> 0057b000-0058a000 r-xp 00000000 03:02 867429      
> /lib/libresolv-2.3.5.so
> 0058a000-0058b000 r-xp 0000e000 03:02 867429      
> /lib/libresolv-2.3.5.so
> 0058b000-0058c000 rwxp 0000f000 03:02 867429      
> /lib/libresolv-2.3.5.so
> 0058c000-0058e000 rwxp 0058c000 00:00 0
> 00590000-00599000 r-xp 00000000 03:02 867440  
> /lib/libgcc_s-4.0.1-20050727.so.1
> 00599000-0059a000 rwxp 00009000 03:02 867440  
> /lib/libgcc_s-4.0.1-20050727.so.1
> 0059c000-0067c000 r-xp 00000000 03:02 111392      
> /usr/lib/libstdc++.so.6.0.5
> 0067c000-00681000 rwxp 000df000 03:02 111392      
> /usr/lib/libstdc++.so.6.0.5
> 00681000-00687000 rwxp 00681000 00:00 0
> 00689000-006df000 r-xp 00000000 03:02 99344       
> /usr/lib/libblas.so.3.0.3
> 006df000-006e0000 rwxp 00055000 03:02 99344       
> /usr/lib/libblas.so.3.0.3
> 0070a000-0070c000 r-xp 00000000 03:02 864925     /lib/libcom_err.so.2.1
> 0070c000-0070d000 rwxp 00001000 03:02 864925     /lib/libcom_err.so.2.1
> 00787000-0079d000 r-xp 00000000 03:02 201709  
> /usr/libexec/octave/2.1.71/site/oct/i686-pc-linux-gnu/octave-forge/ 
> dispatch.oct
> 0079d000-0079e000 rwxp 00016000 03:02 201709  
> /usr/libexec/octave/2.1.71/site/oct/i686-pc-linux-gnu/octave-forge/ 
> dispatch.oct
> 0097c000-009eb000 r-xp 00000000 03:02 105049      
> /usr/lib/libkrb5.so.3.2
> 009eb000-009ee000 rwxp 0006e000 03:02 105049      
> /usr/lib/libkrb5.so.3.2
> 009f0000-00a25000 r-xp 00000000 03:02 864979     /lib/libssl.so.0.9.7f
> 00a25000-00a28000 rwxp 00035000 03:02 864979     /lib/libssl.so.0.9.7f
> 00a2a000-00a41000 r-xp 00000000 03:02 105050  
> /usr/lib/libgssapi_krb5.so.2.2
> 00a41000-00a42000 rwxp 00017000 03:02 105050  
> /usr/lib/libgssapi_krb5.so.2.2
> 00cee000-00d11000 r-xp 00000000 03:02 111383      
> /usr/lib/libk5crypto.so.3.0
> 00d11000-00d12000 rwxp 00023000 03:02 111383      
> /usr/lib/libk5crypto.so.3.0
> 00d14000-00d16000 r-xp 00000000 03:02 106698  
> /usr/lib/libkrb5support.so.0.0
> 00d16000-00d17000 rwxp 00001000 03:02 106698  
> /usr/lib/libkrb5support.so.0.0
> 00d64000-00ddb000 r-xp 00000000 03:02 201832  
> /usr/libexec/octave/2.1.71/site/oct/i686-pc-linux-gnu/octave-forge/ 
> make_sparse.oct
> 00ddb000-00ddd000 rwxp 00077000 03:02 201832  
> /usr/libexec/octave/2.1.71/site/oct/i686-pc-linux-gnu/octave-forge/ 
> make_sparse.oct
> 02bdb000-02f92000 r-xp 00000000 03:02 99339       
> /usr/lib/liblapack.so.3.0.3
> 02f92000-02f95000 rwxp 003b7000 03:02 99339       
> /usr/lib/liblapack.so.3.0.3
> 02f95000-03099000 rwxp 02f95000 00:00 0
> 0309b000-03675000 r-xp 00000000 03:02 128643  
> /usr/lib/octave-2.1.71/liboctave.so.2.1.71
> 03675000-03681000 rwxp 005da000 03:02 128643  
> /usr/lib/octave-2.1.71/liboctave.so.2.1.71
> 03681000-03682000 rwxp 03681000 00:00 0
> 03dc0000-03dfe000 r-xp 00000000 03:02 109746      
> /usr/lib/libncurses.so.5.4
> 03dfe000-03e07000 rwxp 0003d000 03:02 109746      
> /usr/lib/libncurses.so.5.4
> 03e09000-0457d000 r-xp 00000000 03:02 128644  
> /usr/lib/octave-2.1.71/liboctinterp.so.2.1.71
> 0457d000-04599000 rwxp 00773000 03:02 128644  
> /usr/lib/octave-2.1.71/liboctinterp.so.2.1.71
> 04599000-0459b000 rwxp 04599000 00:00 0
> 0460f000-04707000 r-xp 00000000 03:02 867431      
> /lib/libcrypto.so.0.9.7f
> 04707000-04719000 rwxp 000f8000 03:02 867431      
> /lib/libcrypto.so.0.9.7f
> 04719000-0471c000 rwxp 04719000 00:00 0
> 064f8000-0662c000 r-xp 00000000 03:02 107290      
> /usr/lib/libhdf5.so.0.0.0
> 0662c000-0662f000 rwxp 00134000 03:02 107290      
> /usr/lib/libhdf5.so.0.0.0
> 08048000-08049000 r-xp 00000000 03:02 99342      /usr/bin/octave-2.1.71
> 08049000-0804d000 rw-p 00000000 03:02 99342      /usr/bin/octave-2.1.71
> 09152000-09d2a000 rw-p 09152000 00:00 0          [heap]
> b6c00000-b6c21000 rw-p b6c00000 00:00 0
> b6c21000-b6d00000 ---p b6c21000 00:00 0
> b6de6000-b7c4c000 rw-p b6de6000 00:00 0
> b7c4c000-b7c4d000 r--p 00edd000 03:02 107632  
> /usr/lib/locale/locale-archive
> b7c4d000-b7c7e000 r--p 00e5c000 03:02 107632  
> /usr/lib/locale/locale-archive
> b7c7e000-b7e7e000 r--p 00000000 03:02 107632  
> /usr/lib/locale/locale-archive
> b7e7e000-b7fe1000 rw-p b7e7e000 00:00 0
> b7ffb000-b7ffc000 rw-p b7ffb000 00:00 0
> bfee7000-bfefc000 rw-p bfee7000 00:00 0          [stack]
> panic: Aborted -- stopping myself...
> attempting to save variables to `octave-core'...
> save to `octave-core' complete
> Abort
>
> -- 
> Orion Poplawski
> System Administrator                   303-415-9701 x222
> Colorado Research Associates/NWRA      FAX: 303-415-9702
> 3380 Mitchell Lane, Boulder CO 80301   http://www.co-ra.com
> function [grid,xg,yg]=make_grid(minx,maxx,miny,maxy,dx,dy)
>
> xg=minx:dx:maxx;
> yg=maxy:-dy:miny;
>
> k=0;
> for x=xg
>    for y=yg
>       k=k+1;
>       grid(k,1)=x;
>       grid(k,2)=y;
>    end
> end
>
> %%% function MedBhm_preproc.m
>
> %**************************************************
> %**** Make grids for atm
> %**************************************************
>
> [grdxyA,xgA,ygA]=make_grid(14.000,37.000,30.000,46.0000,.25,.25);
> nxA=length(xgA)-2;
> nyA=length(ygA)-2;
>
> MtilB=make_MtilB_sparse(nxA,nyA);
> %%% function MedBhm_preproc.m
>
> %**************************************************
> %**** Make grids for atm
> %**************************************************
>
> [grdxyA,xgA,ygA]=make_grid(14.000,37.000,30.000,46.0000,.25,.25);
> nx=length(xgA)-2;
> ny=length(ygA)-2;
>
> Nx=nx+2;
> Ny=ny+2;
> N=Nx*Ny;
> n=nx*ny;
> NB=2*Nx+2*Ny-4;
>
> MtilB=sparse(NB,N);
>
> for k=1:Ny
>    MtilB(k,k)=1;
>    MtilB(Ny+Nx-2+k,N+1-k)=1;
> end
>
> for i=1:nx
>    MtilB(Ny+i,Ny*(i+1))=1;
>    MtilB(NB+1-i,(Ny*i)+1)=1;
> end
>
> function MtilB=make_MtilB_sparse(nx,ny)
>
> Nx=nx+2;
> Ny=ny+2;
> N=Nx*Ny;
> n=nx*ny;
> NB=2*Nx+2*Ny-4;
>
> MtilB=sparse(NB,N);
>
> for k=1:Ny
>    MtilB(k,k)=1;
>    MtilB(Ny+Nx-2+k,N+1-k)=1;
> end
>
> for i=1:nx
>    MtilB(Ny+i,Ny*(i+1))=1;
>    MtilB(NB+1-i,(Ny*i)+1)=1;
> end
>

Re: [OctDev] Panic in sparse matrix code

A collection of packages providing extra functionality for GNU Octave

Re: [OctDev] Panic in sparse matrix code