[Scipy-tickets] [SciPy] #376: sandbox.delaunay crash on eps-perturbed data

SciPy scipy-tickets@scipy....
Tue Feb 20 05:52:01 CST 2007


#376: sandbox.delaunay crash on eps-perturbed data
--------------------+-------------------------------------------------------
 Reporter:  pv      |       Owner:  somebody
     Type:  defect  |      Status:  new     
 Priority:  normal  |   Milestone:          
Component:  Other   |     Version:  0.5.2   
 Severity:  normal  |    Keywords:          
--------------------+-------------------------------------------------------
 sandbox.delaunay.Triangulation appears to crash if two data points
 are machine-epsilon close to each other.

 On scipy.sandbox.delaunay shipped with Scipy 0.5.2:
 {{{
 from scipy.sandbox.delaunay import Triangulation
 import scipy as S

 x  = S.array([7, 8, 9, 7, 8, 9], dtype=S.float64)
 y  = S.array([0, 0, 0, 1, 1, 1], dtype=S.float64)

 xp = S.r_[x, x + 1e-16] # xp differs from x by epsilon
 yp = S.r_[y, y]

 print "%.20g %.20g" % (x[-1], xp[-1])

 def test_1():
     tri = Triangulation(x, y)
     print "Test 1 OK"

 def test_2():
     tri = Triangulation(xp, yp)
     print "Test 2 OK"

 test_1()
 test_2()
 }}}
 crashes on test_2 with
 {{{
 Test 1 OK
 *** glibc detected *** python: corrupted double-linked list: 0x081fa6d8
 ***
 }}}

 Traceback attached.

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


More information about the Scipy-tickets mailing list