[Pkg-octave-devel] Bug#528078: octave3.0: Segfault when doing a big calcul.

[David Kremer]

Package: octave3.0
Version: 1:3.0.1-6lenny3
Severity: important


Ok, I'm trying to do a big computation with a exponential time 
differencing method, using lot of fft, ifft and vector
element-by-element multiplication. The matter is that , when
I'm trying to run my script, octave3.0 crash with a segfault.

I have run it in gdb :

# START OF BACKTRACE

(gdb) run equa_diff_SSM_ETD1.m 
Starting program: /usr/bin/octave equa_diff_SSM_ETD1.m
(no debugging symbols found)
(no debugging symbols found)
[Thread debugging using libthread_db enabled]
[New Thread 0xb5b4cab0 (LWP 19081)]
GNU Octave, version 3.0.1
Copyright (C) 2008 John W. Eaton and others.
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'.

Octave was configured for "i486-pc-linux-gnu".

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 <bug at octave.org> (but first, please read
http://www.octave.org/bugs.html to learn how to write a helpful report).

For information about changes from previous versions, type `news'.


Program received signal SIGSEGV, Segmentation fault.
[Switching to Thread 0xb5b4cab0 (LWP 19081)]
0xb794fc42 in symbol_record::clear (this=0xaad5fd8) at ov.h:392
392	ov.h: No such file or directory.
	in ov.h
(gdb) 

# END OF BACKTRACE

Step to reproduce : run this script with octave interpreter.

%%%%%%%%%%%%%%%%%%%%%START OF SCRIPT%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all ;

delta_z = 1e-7 ; %
fenetre_t = 7 ;

N_z = 1000 ; 
N_t = 1024 ;%256 ;
N_nu = N_t ;

delta_t = fenetre_t ./ ( N_t - 1) ;

delta_nu = 1 ./ ( ( N_nu - 1 ) .* delta_t ) ;

vecteur_t = [ -fenetre_t/2: delta_t : fenetre_t/2 ] ;
vecteur_z = [ 0 : delta_z : ( N_z - 1 ) * delta_z ] ;
vecteur_nu = [ 0 :delta_nu : (N_nu/2)*delta_nu ,  -(N_nu/2-1)*delta_nu : delta_nu : -delta_nu ] ;

p = 5 ;
gamma = 0 ;

eta  = vecteur_t * p + (p.^3) * vecteur_z(1) - gamma ; 
u = zeros ( N_z , N_t ) ;
u(1,:) = ( sqrt(6)*p ) ./ cosh(eta) ;
c = ( 2 .* pi .* i .* vecteur_nu ).^3 ; 

g1=exp( c .* delta_z ) ;

g2=zeros(1,N_t);
g2(2:N_t)=(g1(2:N_t)-1)./c(2:N_t);
g2(1) = delta_z ;


for k = [ 1:N_z-1 ] 
    u_chapeau_0 = fft( u( k , : )  ) ; %  
    vecteur_derive = real(ifft( i .* 2 .*  pi .* vecteur_nu .* u_chapeau_0 ));
    terme_non_lineaire = fft( ( u( k , : ).^2 ) .* vecteur_derive  )  ; 
    u_chapeau_1 = g1 .* u_chapeau_0 +  terme_non_lineaire .*g2 ;
    u(k+1,:) = real(ifft(u_chapeau_1)) ;
end ;

imagesc(vecteur_t,vecteur_z,real(u)) ;
%%%%%%%%%%%%%%END OF SCRIPT%%%%%%%%%%%%%%%%%%%%%

-- System Information:
Debian Release: 5.0.1
  APT prefers stable
  APT policy: (500, 'stable')
Architecture: i386 (i686)

Kernel: Linux 2.6.26-2-686 (SMP w/2 CPU cores)
Locale: LANG=fr_FR.UTF-8, LC_CTYPE=fr_FR.UTF-8 (charmap=UTF-8)
Shell: /bin/sh linked to /bin/bash

Versions of packages octave3.0 depends on:
ii  libatlas3gf-base [libl 3.6.0-22          Automatically Tuned Linear Algebra
ii  libblas3gf [libblas.so 1.2-2             Basic Linear Algebra Subroutines 3
ii  libc6                  2.7-18            GNU C Library: Shared libraries
ii  libcurl3-gnutls        7.18.2-8lenny2    Multi-protocol file transfer libra
ii  libfftw3-3             3.1.2-3.1         library for computing Fast Fourier
ii  libgcc1                1:4.3.2-1.1       GCC support library
ii  libgfortran3           4.3.2-1.1         Runtime library for GNU Fortran ap
ii  libglpk0               4.29-2            linear programming kit with intege
ii  libhdf5-serial-1.6.6-0 1.6.6-4           Hierarchical Data Format 5 (HDF5) 
ii  liblapack3gf [liblapac 3.1.1-1           library of linear algebra routines
ii  libncurses5            5.7+20081213-1    shared libraries for terminal hand
ii  libpcre3               7.6-2.1           Perl 5 Compatible Regular Expressi
ii  libqhull5              2003.1-9lenny1    calculate convex hulls and related
ii  libreadline5           5.2-3.1           GNU readline and history libraries
ii  libstdc++6             4.3.2-1.1         The GNU Standard C++ Library v3
ii  libsuitesparse-3.1.0   3.1.0-3           collection of libraries for comput
ii  texinfo                4.11.dfsg.1-4     Documentation system for on-line i
ii  zlib1g                 1:1.2.3.3.dfsg-12 compression library - runtime

Versions of packages octave3.0 recommends:
ii  gnuplot                       4.2.2-1.2  A command-line driven interactive 
ii  libatlas3gf-base              3.6.0-22   Automatically Tuned Linear Algebra

Versions of packages octave3.0 suggests:
ii  octave3.0-doc            1:3.0.1-6lenny3 PDF documentation on the GNU Octav
ii  octave3.0-emacsen        1:3.0.1-6lenny3 Emacs support for the GNU Octave l
ii  octave3.0-headers        1:3.0.1-6lenny3 header files for the GNU Octave la
ii  octave3.0-htmldoc        1:3.0.1-6lenny3 HTML documentation on the GNU Octa
ii  octave3.0-info           1:3.0.1-6lenny3 GNU Info documentation on the GNU 

-- no debconf information