Mailman 3 April 2010 - NumPy-Discussion

Re: [Numpy-discussion] My GSoC Proposal to Implement a Subset of NumPy for PyPy
by Dan Roberts 21 Apr '10

21 Apr '10

Oops, I intended to dig more through the NumPy source before I sent the final version of that message, so I could be speaking from an informed standpoint. Thanks, Daniel Roberts > > _______________________________________________ > NumPy-Discussion mailing list > NumPy-Discuss... -- Travis Oliphant Enthought Inc. 1-512-536-1057 http://www.enthought.com oliphant(a)enthought.com ...

1 0

reformulating a series of correlations as one fft, ifft pair
by Paul Northug 21 Apr '10

21 Apr '10

I am having trouble reformulating a series of correlations as a single fft, ifft pair. I have a set of kernels phi : (M = channel, N = kernel, T = time) correlated with signal a : (N, P+T-1) yielding x : (M, T). The correlation, for now, is only in the last dimension, with the two other dimensions degenerate. Below are 4 implementations (x == y == z != w). I would like to reformulate the operation as a single pair of 3d ffts to compare how it fares for different problem dimensions but I am not able to do it. Implementation #4 is incorrect. Please note the .conj() to get correlations rather than convolutions. Any tips would be appreciated. Optimally, I would not have to pad with too many zeros as the matrices can be quite large. The challenge is to beat implementation 2 as P grows. M, N, P are typically on the order of 10^2 and T, 10^3. # import numpy as np from scipy.signal import fftn, ifftn, correlate M, N, P, T = 2, 3, 4, 5 phi = np.random.randn(M, N, P) a = np.random.randn(N, P+T-1) x = np.zeros((M, T)) # 1. correlate x = np.array([correlate(phi[i], a, 'valid').squeeze() for i in range(M)]) # 2. dot products y = np.zeros_like(x) for b in range(P): y += np.dot(phi[:,:,b], a[:,b:b+T]) # 3. a series of ffts (code adapted from scipy convolve) s1 = np.array(phi.shape[1:]) s2 = np.array(a.shape) size = s1+s2-1 z = np.empty_like(y) fsize = (2**np.ceil(np.log2(size))).astype(np.int) # power of 2 pad af = fftn(a, fsize) fslice = tuple([slice(0, sz) for sz in size]) for i in range(M): phif = fftn(phi[i], fsize).conj() z[i] = np.real(ifftn(phif * af)[fslice])[0,:T] # 4. <BROKEN> single fft, ifft pair - some rows are correct s1 = np.array(phi.shape) s2 = np.array(a[np.newaxis].shape) size = s1+s2-1 fsize = (2**np.ceil(np.log2(size))).astype(np.int) af = fftn(a[np.newaxis], fsize) phif = fftn(phi, fsize).conj() fslice = tuple([slice(0, sz) for sz in size]) w = np.real(ifftn(phif * af)[fslice])[:,:,:T]

3 3

My GSoC Proposal to Implement a Subset of NumPy for PyPy
by Dan Roberts 20 Apr '10

20 Apr '10

Hello NumPy Users, Hi everybody, my name is Dan Roberts, and my Google Summer of Code proposal was categorized under NumPy rather than PyPy, so it will end up being reviewed by mentors for the NumPy project. I'd like to take this chance to introduce myself and my proposal. I hadn't prepared for review by the NumPy mentors, but this can make my proposal stronger than before. With a bit of help from all of you, I can dedicate my summer to creating more useful code than I would have previously. I realize that from the perspective of NumPy, my proposal might seem lacking, so I'd like to also invite the scrutiny of all of the readers of this list. Why should we bother reimplimenting anything? PyPy, for those who are unfamiliar, has the ability to Just-in-Time compile itself and programs that it's running. One of the major advantages of this is that code operating on NumPy arrays could potentially be written in pure-python, with normal looping constructs, and be nearly as fast as a ufunc painstakingly crafted in C. I'd love to see as much Python and as little C as possible, and I'm sure I'm not alone in that wish. A short introduction: I've been coding in Python for the past few years, and have increasingly become interested in speeding up what has become my favorite language. To that end I've become interested in both the PyPy project and the NumPy projects. I've spent a fair amount of time frustrating the PyPy developers with silly questions, written a bit of code for them, and now my GSoC proposal involves both them, and NumPy. Finally, I'd like to ask all of you: what features are most important to you? It's not practical, wise, or even possible for me to reimpliment more than a small portion of NumPy, but if I can address the most important parts, maybe I can make this project useful enough for some of you to use, and close enough for the rest of you that I can drum up some support for more development in the future. My proposal lives at http://codespeak.net/~dan/gsoc/micronumpy.htmlthanks for making it this far through my long winded introduction! I welcome all constructive criticism and thoughts. Thanks, Daniel Roberts

7 7

scipy error undefined symbol: lsame_
by gerardob 20 Apr '10

20 Apr '10

I installed scipy (and all the required libraries) and the following error appears when i tried run a simple example which uses the optimize package of scipy. I tried also numpy alone and it works ( at least for printing numpy.array([10,20,10])) error: Traceback (most recent call last): File "main_test.py", line 2, in <module> from scipy import optimize File "/home/gberbeglia/python/Python-2.6.5/lib/python2.6/site-packages/scipy/optimize/__init__.py", line 11, in <module> from lbfgsb import fmin_l_bfgs_b File "/home/gberbeglia/python/Python-2.6.5/lib/python2.6/site-packages/scipy/optimize/lbfgsb.py", line 30, in <module> import _lbfgsb ImportError: /home/gberbeglia/python/Python-2.6.5/lib/python2.6/site-packages/scipy/optimize/_lbfgsb.so: undefined symbol: lsame_ gberbeglia@actarus:~/python/mycodes> Any ideas on how to solve this problem? Thanks. PS: the example is below: import numpy from scipy import optimize a = numpy.array([10,20,10]) print a def f_(x): return x*x x,f,d = optimize.fmin_l_bfgs_b(f_,[0.1],fprime=None, approx_grad = True, bounds = [(-10000,10000)], iprint=30, maxfun=150000) -- View this message in context: http://old.nabble.com/scipy-error-undefined-symbol%3A-lsame_-tp28287720p282… Sent from the Numpy-discussion mailing list archive at Nabble.com.

2 1

Building DMG w/ Python 2.6 (could not access...no space left on device)
by Patrick Marsh 20 Apr '10

20 Apr '10

Greetings, I'm trying to build DMG releases of numpy-1.4.1rc3 as part of the release process. I recently got my mac back from repairs and spent the weekend trying to get the build environment configured on it. I thought it was correct, but I've got errors when building with both Python 2.5.4 and Python 2.6.5. (I should add that, ironically, I have no problem building the Windows superpack on OS X [10.6].) Below is the end of the output from running 'paver dmg -p 2.6'. I have no idea why I'm told I'm out of space, nor do I know where to begin in trying to track that down. I'll send a separate email regarding my errors regarding Python 2.5 as it is an entirely different issue. Thanks in advance for any help... ============================================================ creating dist creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages nothing to be installed for scheme purelib nothing to be installed for scheme headers creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources /usr/bin/mkbom build/bdist.macosx-10.3-fat/mpkg/platlib/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/./Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.bom /bin/pax -w -f /Users/pmarsh/git/numpy.release/dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.pax.gz -x cpio -z . creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources/English.lproj copying License.txt -> dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/./Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources nothing to be installed for scheme data creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources /usr/bin/mkbom build/bdist.macosx-10.3-fat/mpkg/scripts/usr/local/bin dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/./Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.bom /bin/pax -w -f /Users/pmarsh/git/numpy.release/dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.pax.gz -x cpio -z . creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources/English.lproj copying License.txt -> dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/./Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Resources creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Resources/English.lproj copying License.txt -> dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Resources removing 'build/bdist.macosx-10.3-fat/mpkg' (and everything under it) Creating disk image... ....................................................................................................................... 2010-04-19 09:17:58.614 copy-helper[10435:903] copy error (canceling): /Volumes/numpy/numpy-1.4.1rc3-py2.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.pax.gz: No space left on device ............................................................................could not access /Volumes/numpy/numpy-1.4.1rc3-py2.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.pax.gz - No space left on device hdiutil: create failed - No space left on device hdiutil: attach failed - No such file or directory mkdir: /Volumes/numpy: No such file or directory cp: /Volumes/numpy/.background/dmgbackground.png: No such file or directory 63:67: execution error: Finder got an error: Can’t get disk "numpy". (-1728) chmod: /Volumes/numpy: No such file or directory hdiutil: detach: no device name specified Usage: hdiutil detach [options] <devname> hdiutil detach -help hdiutil: convert failed - No such file or directory Captured Task Output: --------------------- ---> pavement.dmg ---> pavement.pdf ---> pavement.latex Build pdf doc ---> pavement.clean LDFLAGS='-undefined dynamic_lookup -bundle -arch i386 -arch ppc -Wl,-search_paths_first -Lbuild' /Library/Frameworks/Python.framework/Versions/2.6/bin/python setupegg.py bdist_mpkg Traceback (most recent call last): File "/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/paver/tasks.py", line 170, in _run_task return do_task() File "/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/paver/tasks.py", line 167, in do_task return func(**kw) File "pavement.py", line 470, in dmg shutil.copy(source, target) File "/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/shutil.py", line 88, in copy copyfile(src, dst) File "/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/shutil.py", line 52, in copyfile fsrc = open(src, 'rb') IOError: [Errno 2] No such file or directory: 'tools/numpy-macosx-installer/numpy-1.4.1rc3-py2.6-python.org.dmg' Patrick -- Patrick Marsh Ph.D. Student / NSSL Liaison to the HWT School of Meteorology / University of Oklahoma Cooperative Institute for Mesoscale Meteorological Studies National Severe Storms Laboratory http://www.patricktmarsh.com

2 1

Re: [Numpy-discussion] reformulating a series of correlations as one fft, ifft pair
by Paul Northug 20 Apr '10

20 Apr '10

<josef.pktd <at> gmail.com> writes: > > On Mon, Apr 19, 2010 at 2:36 AM, Paul Northug <pnorthug <at> gmail.com> wrote: <snip> > > > > # > > import numpy as np > > from scipy.signal import fftn, ifftn, correlate > > > > M, N, P, T = 2, 3, 4, 5 > > > > phi = np.random.randn(M, N, P) > > a = np.random.randn(N, P+T-1) > > x = np.zeros((M, T)) > > > > # 1. correlate > > > > x = np.array([correlate(phi[i], a, 'valid').squeeze() for i in range(M)]) > > > > # 2. dot products > > > > y = np.zeros_like(x) > > for b in range(P): y += np.dot(phi[:,:,b], a[:,b:b+T]) > > > > # 3. a series of ffts (code adapted from scipy convolve) > > > > s1 = np.array(phi.shape[1:]) > > s2 = np.array(a.shape) > > size = s1+s2-1 > > z = np.empty_like(y) > > fsize = (2**np.ceil(np.log2(size))).astype(np.int) # power of 2 pad > > af = fftn(a, fsize) > > fslice = tuple([slice(0, sz) for sz in size]) > > for i in range(M): > > phif = fftn(phi[i], fsize).conj() > > z[i] = np.real(ifftn(phif * af)[fslice])[0,:T] > > > > # 4. <BROKEN> single fft, ifft pair - some rows are correct > > > > s1 = np.array(phi.shape) > > s2 = np.array(a[np.newaxis].shape) > > size = s1+s2-1 > > fsize = (2**np.ceil(np.log2(size))).astype(np.int) > > af = fftn(a[np.newaxis], fsize) > > phif = fftn(phi, fsize).conj() > > fslice = tuple([slice(0, sz) for sz in size]) > > w = np.real(ifftn(phif * af)[fslice])[:,:,:T] > > I only checked this in your example > > >>> np.squeeze(signal.fftconvolve(phi[:,::-1,::-1], a[None,...],mode='valid')) > array([[ 0.42695078, 4.3433924 , -2.53084395, -4.22672281, -5.16705981], > [-2.35108765, -2.15060513, -1.23528122, 1.40386599, 1.38189896]]) > > >>> y > array([[ 0.42695078, 4.3433924 , -2.53084395, -4.22672281, -5.16705981], > [-2.35108765, -2.15060513, -1.23528122, 1.40386599, 1.38189896]]) > >>> x > array([[ 0.42695078, 4.3433924 , -2.53084395, -4.22672281, -5.16705981], > [-2.35108765, -2.15060513, -1.23528122, 1.40386599, 1.38189896]]) > > Josef > Thanks Josef. You saved me a lot of work. I was missing the middle ::-1. The fft approach is orders of magnitude slower than method #2 in this case. Using CUDA fft3, it is about 15 times slower for the desired problem size. It is not a fair comparison as the convolution is degenerate along 2 dimensions but it is useful to know.

1 0

patch to pickle np.memmap
by Brent Pedersen 19 Apr '10

19 Apr '10

hi, i posted a patch to allow pickling of np.memmap objects. http://projects.scipy.org/numpy/ticket/1452 currently, it always returns 'r' for the mode. is that the best thing to do there? any other changes? -brent

1 2

Building DMG w/ Python 2.6 (could not access...no space left on device)
by Patrick Marsh 19 Apr '10

19 Apr '10

Greetings, (I apologize if this becomes a duplicate. I originally sent this email from the wrong email address.) I'm trying to build DMG releases of numpy-1.4.1rc3 as part of the release process. I recently got my mac back from repairs and spent the weekend trying to get the build environment configured on it. I thought it was correct, but I've got errors when building with both Python 2.5.4 and Python 2.6.5. (I should add that, ironically, I have no problem building the Windows superpack on OS X [10.6].) Below is the end of the output from running 'paver dmg -p 2.6'. I have no idea why I'm told I'm out of space, nor do I know where to begin in trying to track that down. I'll send a separate email regarding my errors regarding Python 2.5 as it is an entirely different issue. Thanks in advance for any help... ============================================================ creating dist creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages nothing to be installed for scheme purelib nothing to be installed for scheme headers creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources /usr/bin/mkbom build/bdist.macosx-10.3-fat/mpkg/platlib/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/./Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.bom /bin/pax -w -f /Users/pmarsh/git/numpy.release/dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.pax.gz -x cpio -z . creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources/English.lproj copying License.txt -> dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/./Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources nothing to be installed for scheme data creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources /usr/bin/mkbom build/bdist.macosx-10.3-fat/mpkg/scripts/usr/local/bin dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/./Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.bom /bin/pax -w -f /Users/pmarsh/git/numpy.release/dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.pax.gz -x cpio -z . creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources/English.lproj copying License.txt -> dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/./Contents/Packages/numpy-scripts-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Resources creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Resources creating dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Resources/English.lproj copying License.txt -> dist/numpy-1.4.1rc3-py2.6-macosx10.6.mpkg/Contents/Resources removing 'build/bdist.macosx-10.3-fat/mpkg' (and everything under it) Creating disk image... ....................................................................................................................... 2010-04-19 09:17:58.614 copy-helper[10435:903] copy error (canceling): /Volumes/numpy/numpy-1.4.1rc3-py2.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.pax.gz: No space left on device ............................................................................could not access /Volumes/numpy/numpy-1.4.1rc3-py2.6.mpkg/Contents/Packages/numpy-platlib-1.4.1rc3-py2.6-macosx10.6.pkg/Contents/Archive.pax.gz - No space left on device hdiutil: create failed - No space left on device hdiutil: attach failed - No such file or directory mkdir: /Volumes/numpy: No such file or directory cp: /Volumes/numpy/.background/dmgbackground.png: No such file or directory 63:67: execution error: Finder got an error: Can’t get disk "numpy". (-1728) chmod: /Volumes/numpy: No such file or directory hdiutil: detach: no device name specified Usage: hdiutil detach [options] <devname> hdiutil detach -help hdiutil: convert failed - No such file or directory Captured Task Output: --------------------- ---> pavement.dmg ---> pavement.pdf ---> pavement.latex Build pdf doc ---> pavement.clean LDFLAGS='-undefined dynamic_lookup -bundle -arch i386 -arch ppc -Wl,-search_paths_first -Lbuild' /Library/Frameworks/Python.framework/Versions/2.6/bin/python setupegg.py bdist_mpkg Traceback (most recent call last): File "/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/paver/tasks.py", line 170, in _run_task return do_task() File "/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/paver/tasks.py", line 167, in do_task return func(**kw) File "pavement.py", line 470, in dmg shutil.copy(source, target) File "/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/shutil.py", line 88, in copy copyfile(src, dst) File "/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/shutil.py", line 52, in copyfile fsrc = open(src, 'rb') IOError: [Errno 2] No such file or directory: 'tools/numpy-macosx-installer/numpy-1.4.1rc3-py2.6-python.org.dmg' Patrick -- Patrick Marsh Ph.D. Student / NSSL Liaison to the HWT School of Meteorology / University of Oklahoma Cooperative Institute for Mesoscale Meteorological Studies National Severe Storms Laboratory http://www.patricktmarsh.com

1 0

Infix operators for numpy arrays
by John Schulman 19 Apr '10

19 Apr '10

Hi all, I just posted a recipe (http://code.activestate.com/recipes/577201/) for infix operators that work with numpy arrays, so you can write X *dot* Y. It works by making the Infix class a subclass of ndarray. John

1 0

Re: [Numpy-discussion] My GSoC Proposal to Implement a Subset of NumPy for PyPy
by Dan Roberts 19 Apr '10

19 Apr '10

On Apr 18, 2010 6:46 PM, "Dan Roberts" <ademan555(a)gmail.com> wrote: I've been trying my best to take my time formulating my replies but I need to respond eventually. :-) This is embarassing, but I'm actually not sure where I talked about an interface specifically. I did rather nebulously talk about interfacing with C code and LAPACK and the interfaces there would be provided by the respective code and consumed by micronumpy, at least as I currently see it. I haven't consulted maciej about this yet, but I think working backwards from a complete C NumPy depends on a great deal of "ifs", many of which I think aren't satisfied. I need to look into this, but I assume NumPy operates on array structures directly, rather than through an interface. If it's through an interface, there's a real possibility that approach is possible, though it would require me to write some adaptors, I think it would be ok, and a low enough time investment. Like I said I'm currently speaking from ignorance so I need to look into it and get back to you. Cheers, Dan P.S. I agree about the sparse matrices, I've bugged fijal a small bit about that. P.P.S. Forgot to CC the mailing list... assumed this mail client would do it for me.. lol > > On Apr 17, 2010 12:25 AM, "Stéfan van der Walt" <stefan(a)sun.ac.za> wrote: > > Hi Dan > > > On 17 April 2010 06:50, Dan Roberts <ademan555(a)gmail.com> wrote: > Hi everybody, my name is Dan... > > Thanks for the introduction, and welcome to NumPy! > > > > I hadn't prepared for review by the NumPy mentors, but this can make my > proposal stronger t... > > This proposal builds a bridge between two projects, so even if it > technically falls under the... > > > > Why should we bother reimplimenting anything? PyPy, for those who are > unfamiliar, has the ... > > Your code has a fairly specialised application and it's worth > discussing exactly where it wou...

1 0