[Scipy-tickets] [SciPy] #776: Segmentation fault in interpolate.interp2d for too many data points

SciPy scipy-tickets@scipy....
Sun Nov 2 20:30:32 CST 2008


#776: Segmentation fault in interpolate.interp2d for too many data points
-------------------------------+--------------------------------------------
 Reporter:  pv                 |       Owner:  somebody                   
     Type:  defect             |      Status:  new                        
 Priority:  normal             |   Milestone:  0.7.0                      
Component:  scipy.interpolate  |     Version:                             
 Severity:  normal             |    Keywords:  interpolate bisplrep surfit
-------------------------------+--------------------------------------------
 Scipy SVN r4939, 32-bit Ubuntu system.

 interpolate.interp2d segfaults when there are too many data points (but
 not enough to cause immediate lwrk overflow, which would be detected):
 {{{
 import numpy as np
 from scipy.interpolate import interp2d

 x = np.arange(1203, dtype=np.float64)
 y = np.arange(63, dtype=np.float64)
 z = np.zeros(1203*63, dtype=np.float64)
 ip = interp2d(x, y, z)
 }}}
 On another (64-bit Debian) machine, this results to a MemoryError in
 _fitpack._surfit.

 gdb shows
 {{{
 (gdb) bt
 #0  0xb6a7b3d8 in fporde (x=0xb642c008, y=0xb6397008, m=@0xbffda980,
     kx=@0xbffda504, ky=@0xbffda4fc, tx=0xb584a008, nx=@0xbffda4e4,
     ty=0xb584aef0, ny=@0xbffda4e0, nummer=0x4f9f35e0, index=0x4fa3d614,
     nreg=@0xbffda4e8) at scipy/interpolate/fitpack/fporde.f:24
 #1  0xb6a74087 in fpsurf (iopt=@0xbffda984, m=@0xbffda980, x=0xb642c008,
     y=0xb6397008, z=0xb64c1008, w=0xb6302008, xb=@0xbffda928,
 xe=@0xbffda920,
     yb=@0xbffda918, ye=@0xbffda910, kxx=@0xbffda97c, kyy=@0xbffda978,
     s=@0xbffda908, nxest=@0xbffda974, nyest=@0xbffda970, eta=@0xbffda8f8,
     tol=@0xbffda778, maxit=@0xbffda794, nmax=@0xbffda78c, km1=@0xbffda79c,
     km2=@0xbffda798, ib1=@0xbffda7a4, ib3=@0xbffda7a0, nc=@0xbffda790,
     intest=@0xbffda784, nrest=@0xbffda788, nx0=@0xbffda96c, tx=0xb584a008,
     ny0=@0xbffda968, ty=0xb584aef0, c=0xb584bdd8, fp=@0xbffda900,
     fp0=@0xb5a048a0, fpint=0x4f7981d8, coord=0x4f799f78, f=0x4f426c48,
     ff=0x4f5df710, a=0x1bf0d5a0, q=0xb5a048a8, bx=0x4f79dad0,
 by=0x4f7a0788,
     spx=0x4f7a3440, spy=0x4f8cb510, h=0x4f79bd18, index=0x4fa3d614,
     nummer=0x4f9f35e0, wrk=0x4fb18ca4, lwrk=@0xbffda960, ier=@0xbffda958)
     at scipy/interpolate/fitpack/fpsurf.f:164
 #2  0xb6a65fa6 in surfit (iopt=@0xbffda984, m=@0xbffda980, x=0xb642c008,
     y=0xb6397008, z=0xb64c1008, w=0xb6302008, xb=@0xbffda928,
 xe=@0xbffda920,
     yb=@0xbffda918, ye=@0xbffda910, kx=@0xbffda97c, ky=@0xbffda978,
     s=@0xbffda908, nxest=@0xbffda974, nyest=@0xbffda970, nmax=@0xbffda950,
     eps=@0xbffda8f8, nx=@0xbffda96c, tx=0xb584a008, ny=@0xbffda968,
     ty=0xb584aef0, c=0xb584bdd8, fp=@0xbffda900, wrk1=0xb5a048a0,
     lwrk1=@0xbffda964, wrk2=0x4fb18ca4, lwrk2=@0xbffda960,
 iwrk=0x4f9f35e0,
     kwrk=@0xbffda95c, ier=@0xbffda958)
     at scipy/interpolate/fitpack/surfit.f:400
 #3  0xb6a51db2 in fitpack_surfit (dummy=0x0, args=0xa1008b4)
     at scipy/interpolate/src/__fitpack.h:194
 #4  0x080cea39 in PyEval_EvalFrameEx (f=0xa1abf14, throwflag=0)
     at ../Python/ceval.c:3595
 ...
 }}}

-- 
Ticket URL: <http://scipy.org/scipy/scipy/ticket/776>
SciPy <http://www.scipy.org/>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list