Tim Leslie wrote:
>  On 11/23/06, Nils Wagner <nwagner@[...].de> wrote:
>    
> > Hi all,
> >
> > Can someone reproduce the segfault by running the following test
> >     
> 
>  I can confirm a segfault using:
> 
>  Python 2.4.4c1 (#2, Oct 11 2006, 20:00:03)
>  [GCC 4.1.2 20060928 (prerelease) (Ubuntu 4.1.1-13ubuntu5)] on linux2
>  Type "help", "copyright", "credits" or "license" for more information.
>    
> >>> import scipy, numpy
> >>>         
>  s>>> scipy.__version__, numpy.__version__
>  ('0.5.2.dev2322', '1.0.1.dev3435')
> 
>  HTH
> 
>  Tim
> 
>    
Thank you Tim. Meanwhile I have added a line A = A.tocsc() see below

from scipy import *

n = 15
A = sparse.lil_matrix((n,n))
for i in arange(0,n):
  A[i,:n] = random.rand(n)
A = A.tocsc()
B = 2.*sparse.speye(n,n)
sigma =  1.0
sigma_solve = linsolve.splu(A - sigma*B).solve
print sigma_solve

If I run the script in interactive mode I get

python -i test_linsolve.py
Use minimum degree ordering on A'+A.
*** glibc detected *** free(): invalid next size (fast):
0x00000000008d3c30 ***
Abort

If I run the script without -i I get

python test_linsolve.py
Use minimum degree ordering on A'+A.
<built-in method solve of factored_lu object at 0x2aaab36d8c60> 

Any idea ?

Nils

> > from scipy import *
> >
> > n = 15
> > A = sparse.lil_matrix((n,n))
> > for i in arange(0,n):
> >   A[i,:n] = random.rand(n)
> > B = 2.*sparse.speye(n,n)
> > sigma =  1.0
> > sigma_solve = linsolve.splu(A - sigma*B).solve
> >
> >
> > Nils
> >
> > I am using the latest svn versions of numpy/scipy.
> > This is the output of gdb
> >
> > GNU gdb 6.3
> > Copyright 2004 Free Software Foundation, Inc.
> > GDB is free software, covered by the GNU General Public License, and you are
> > welcome to change it and/or distribute copies of it under certain
> > conditions.
> > Type "show copying" to see the conditions.
> > There is absolutely no warranty for GDB.  Type "show warranty" for details.
> > This GDB was configured as "x86_64-suse-linux"...(no debugging symbols
> > found)
> > Using host libthread_db library "/lib64/tls/libthread_db.so.1".
> >
> > (gdb) run test_linsolve.py
> > Starting program: /usr/bin/python test_linsolve.py
> > (no debugging symbols found)
> > (no debugging symbols found)
> > [Thread debugging using libthread_db enabled]
> > [New Thread 46912509653888 (LWP 6613)]
> > (no debugging symbols found)
> > (no debugging symbols found)
> > (no debugging symbols found)
> > (no debugging symbols found)
> > (no debugging symbols found)
> > (no debugging symbols found)
> > (no debugging symbols found)
> > (no debugging symbols found)
> > (no debugging symbols found)
> > Use minimum degree ordering on A'+A.
> >
> > Program received signal SIGSEGV, Segmentation fault.
> > [Switching to Thread 46912509653888 (LWP 6613)]
> > 0x00002aaab0e6e9a1 in genmmd_ (neqns=0x7fffffa847ac, xadj=0x8d6d60,
> > adjncy=0x998550, invp=0x8e19ec, perm=0x88cb20, delta=0x7fffffa847a4,
> >     dhead=0x8e196c, qsize=0x8e6040, llist=0x92bcbc, marker=0x94c78c,
> > maxint=0x7fffffa847a0, nofsub=0x7fffffa8479c) at mmd.c:162
> > 162             perm[nextmd] = -mdeg;
> >
> > _______________________________________________
> > Scipy-dev mailing list
> > Scipy-dev@[...].org
> > http://projects.scipy.org/mailman/listinfo/scipy-dev
> >
> >     
>  _______________________________________________
>  Scipy-dev mailing list
>  Scipy-dev@[...].org
>  http://projects.scipy.org/mailman/listinfo/scipy-dev
>    

 
_______________________________________________
Scipy-dev mailing list
Scipy-dev@[...].org
http://projects.scipy.org/mailman/listinfo/scipy-dev