--- Begin Message ---
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@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
--- End Message ---