[SciPy-user] extracting elements of a matrix using arrays as indices (SEGFAULT!)

Stéfan van der Walt stefan@sun.ac...
Wed May 21 10:29:23 CDT 2008


Hi Mike

2008/5/21 Michael Hearne <mhearne@usgs.gov>:
> If I try that on a larger example, similar to the second one I made
>
> yesterday:
>
> nrows = 648
>
> ncols = 690
>
> data = rand(nrows,ncols)
>
> i,j = (data < 0.14).nonzero()
>
> data[i,:][:,j] = data[i,:][:,j]*0
>
> I get another segmentation fault.

I can confirm that this bug exists in r5120.

Regards
Stéfan


More information about the SciPy-user mailing list