Mailman 3 December 2015 - NumPy-Discussion

Numpy 1.10.2rc2 released
by Charles R Harris Dec. 9, 2015

Dec. 9, 2015

Hi All, I'm pleased to announce the release of Numpy 1.10.2rc2. After two months of stomping bugs I think the house is clean and we are almost ready to put it up for sale. However, bugs are persistent and may show up at anytime, so please inspect and test thoroughly. Windows binaries and source releases can be found at the usual place on Sourceforge <https://sourceforge.net/projects/numpy/files/NumPy/1.10.2rc2/>. If there are no reports of problems in the next week I plan to release the … [View More]

3 2

Memory mapping and NPZ files
by Mathieu Dubois Dec. 9, 2015

Dec. 9, 2015

Dear all, If I am correct, using mmap_mode with Npz files has no effect i.e.: f = np.load("data.npz", mmap_mode="r") X = f['X'] will load all the data in memory. Can somebody confirm that? If I'm correct, the mmap_mode argument could be passed to the NpzFile class which could in turn perform the correct operation. One way to handle that would be to use the ZipFile.extract method to write the Npy file on disk and then load it with numpy.load with the mmap_mode argument. Note that the user … [View More]

1 0

array of random numbers fails to construct
by DAVID SAROFF (RIT Student) Dec. 8, 2015

Dec. 8, 2015

This works. A big array of eight bit random numbers is constructed: import numpy as np spectrumArray = np.random.randint(0,255, (2**20,2**12)).astype(np.uint8) This fails. It eats up all 64GBy of RAM: spectrumArray = np.random.randint(0,255, (2**21,2**12)).astype(np.uint8) The difference is a factor of two, 2**21 rather than 2**20, for the extent of the first axis. -- David P. Saroff Rochester Institute of Technology 54 Lomb Memorial Dr, Rochester, NY 14623 david.saroff(a)mail.rit.edu | (434) 227-6242

8 12

Re: [Numpy-discussion] Q: Use of scipy.signal.bilinear
by R Schumacher Dec. 8, 2015

Dec. 8, 2015

Sorry - I'll join there. - Ray At 10:00 AM 12/8/2015, you wrote: >On Tue, Dec 8, 2015 at 9:30 AM, R Schumacher ><<mailto:rays@blue-cove.com>rays(a)blue-cove.com> wrote: >We have a function which describes a frequency >response correction to piezo devices we use. To >flatten the FFT, it is similar to: >Cdis_t = .5 >N = 8192 >for n in range(8192): >Â B3 = n * 2560 / N >Â Fc(n) = 1 / >((B3/((1/(Cdis_t*2*pi))**2+B3**2)**0.5)*(-0.01*log(B3) + 1.… [View More]

1 0

Q: Use of scipy.signal.bilinear
by R Schumacher Dec. 8, 2015

Dec. 8, 2015

We have a function which describes a frequency response correction to piezo devices we use. To flatten the FFT, it is similar to: Cdis_t = .5 N = 8192 for n in range(8192): B3 = n * 2560 / N Fc(n) = 1 / ((B3/((1/(Cdis_t*2*pi))**2+B3**2)**0.5)*(-0.01*log(B3) + 1.04145)) In practice it really only matters for low frequencies. I suggested that we might be able to do a time domain correction as a forward-reverse FFT filter using the function, but another said it can also be applied in … [View More]

2 1

Re: [Numpy-discussion] Where is Jaime?
by Peter Creasey Dec. 8, 2015

Dec. 8, 2015

>> > >> > Is the interp fix in the google pipeline or do we need a workaround? >> > >> >> Oooh, if someone is looking at changing interp, is there any chance >> that fp could be extended to take complex128 rather than just float >> values? I.e. so that I could write: >> >> >>> y = interp(mu, theta, m) >> rather than >> >>> y = interp(mu, theta, m.real) + 1.0j*interp(mu, theta, m.imag) >> >> … [View More]

2 1

Where is Jaime?
by Charles R Harris Dec. 7, 2015

Dec. 7, 2015

Anyone hear from Jaime lately? He seems to have gone missing. Chuck

3 6

Re: [Numpy-discussion] NumPy-Discussion Digest, Vol 111, Issue 9
by Sydney Shall Dec. 7, 2015

Dec. 7, 2015

On 07/12/2015 09:38, numpy-discussion-request(a)scipy.org wrote: > Message: 4 > Date: Sun, 06 Dec 2015 22:01:40 -0500 > From: "DAVID SAROFF (RIT Student)"<dps7802(a)rit.edu> > To: Discussion of Numerical Python<numpy-discussion(a)scipy.org> > Cc: Stefi Baum<stefibaumodea(a)gmail.com> > Subject: Re: [Numpy-discussion] array of random numbers fails to > construct > Message-ID: > <CABMftcSkX9kL=G2sepOn80YCKQygFGHa7WBdmAxy3_1iVAtJcw(a)mail.gmail.com… [View More]

2 1

Re: [Numpy-discussion] Where is Jaime?
by Peter Creasey Dec. 6, 2015

Dec. 6, 2015

> > Is the interp fix in the google pipeline or do we need a workaround? > Oooh, if someone is looking at changing interp, is there any chance that fp could be extended to take complex128 rather than just float values? I.e. so that I could write: >>> y = interp(mu, theta, m) rather than >>> y = interp(mu, theta, m.real) + 1.0j*interp(mu, theta, m.imag) which *sounds* like it might be simple and more (Num)pythonic. Peter

2 1

Travis is busted.
by Charles R Harris Dec. 5, 2015

Dec. 5, 2015

Just to note that travis is broken at the moment, it can't find the whitelisted gfortran package.

1 0