Re: [Numpy-discussion] f2py: ram usage

2016-04-11 Thread George Nurser
fortran array indices are reversed, but this is the most natural way in any case. --George Nurser On 10 April 2016 at 11:53, Sebastian Berg wrote: > On So, 2016-04-10 at 12:04 +0200, Vasco Gervasi wrote: > > Hi all, > > I am trying to write some code to do calculation onto an arra

[Numpy-discussion] Building numpy with OpenBLAS using bento

2012-11-20 Thread George Nurser
t doesn't seem to be appearing on it as far as I can see.) Best regards, George Nurser ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] Scipy dot

2012-11-10 Thread George Nurser
Note also that OpenBlas claims performance as good as MKL with Sandy Bridge processors. https://github.com/xianyi/OpenBLAS/wiki/faq#wiki-sandybridge_perf George Nurser. On 10 November 2012 00:38, Dag Sverre Seljebotn wrote: > On 11/09/2012 11:57 PM, Matthieu Brucher wrote: > > Hi, &g

Re: [Numpy-discussion] Scipy dot

2012-11-09 Thread George Nurser
Hi, It's really good to see this work being done. It would be great if this could somehow be put also into the np.einsum function, which currently doesn't even use blas, and is consequently substantially slower than current np.dot (timings on http://mail.scipy.org/pipermail/numpy-discussion/2012-O

Re: [Numpy-discussion] einsum slow vs (tensor)dot

2012-10-26 Thread George Nurser
On 25 October 2012 22:54, David Warde-Farley wrote: > On Wed, Oct 24, 2012 at 7:18 AM, George Nurser wrote: > > Hi, > > > > I was just looking at the einsum function. > > To me, it's a really elegant and clear way of doing array operations, > which > > i

[Numpy-discussion] einsum slow vs (tensor)dot

2012-10-24 Thread George Nurser
ation seems ~ 4-6x slower than dot or tensordot for decent size arrays. I suspect it is because the implementation does not use blas/lapack calls. cheers, George Nurser. E.g. (in ipython on Mac OS X 10.6, python 2.7.3, numpy 1.6.2 from macports) a = np.arange(60.).reshape(1500,400) b = np.a

Re: [Numpy-discussion] f2py with allocatable arrays

2012-07-03 Thread George Nurser
argument arrays clumsy, I know. George Nurser On 3 July 2012 02:17, Casey W. Stark wrote: > Hi numpy. > > Does anyone know if f2py supports allocatable arrays, allocated inside > fortran subroutines? The old f2py docs seem to indicate that the allocatable > array must be created

Re: [Numpy-discussion] problems with multiple outputs with numpy.nditer

2011-08-10 Thread George Nurser
Works fine with the [...]s. Thanks very much. --George On 10 August 2011 17:15, Mark Wiebe wrote: > On Wed, Aug 10, 2011 at 3:45 AM, George Nurser wrote: >> >> Hi, >> I'm running numpy 1.6.1rc2 + python 2.7.1 64-bit from python.org on OSX >> 10.6.8. >> &

[Numpy-discussion] problems with multiple outputs with numpy.nditer

2011-08-10 Thread George Nurser
o. I'm not sure whether I've just called it incorrectly, or whether perhaps it's only supposed to work with one output array. --George Nurser. ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] ANN: Numpy 1.6.0 beta 2

2011-04-05 Thread George Nurser
SUSE 11.3, python 2.7, gcc 4.3.4, gfortran from gcc 4.6.0, I get two failures on commit 1439a8ddcb2eda20fa102aa44e846783f29c0af3 (head of 1.6.x maintenance branch). --George. == FAIL: Test basic arithmetic function errors --

[Numpy-discussion] Einstein summation convention

2011-01-27 Thread George Nurser
undefined, since a.shape[0] != f.shape[0] E(f,'ij',f,'ji') = f_{ij}*f_{ji} = f*f.T Broadcasting would be explicit E(f,'ij',d,'j') =a_{ij}*d_j = a*d E(f,'ij',d,'i') =a_{ij}*d_i = a*d[:,None] If a definite order of precedence of sums were est

[Numpy-discussion] Fwd: [Matplotlib-users] Vectorization

2010-07-30 Thread George Nurser
-- Forwarded message -- From: George Nurser Date: 30 July 2010 22:37 Subject: Re: [Matplotlib-users] Vectorization To: Nicolas Bigaouette , Discussion of Numerical Python > I want to do the same for the calculation of the kinetic energy: > /2m. There is a laplacian

Re: [Numpy-discussion] Calling routines from a Fortran library using python

2010-02-18 Thread George Nurser
wrote: > If header files are provided, the work done by f2py is almost done. > But you don't know the real Fortran interface, so you still have to > use ctypes over f2py. > > Matthieu > > 2010/2/18 George Nurser : >> Hi Nils, >> I've not tried it, but you

Re: [Numpy-discussion] Calling routines from a Fortran library using python

2010-02-18 Thread George Nurser
Hi Nils, I've not tried it, but you might be able to interface with f2py your own fortran subroutine that calls the library. Then issue the f2py command with extra arguments -l -L. See section 5 of http://cens.ioc.ee/projects/f2py2e/usersguide/index.html#command-f2py --George. On 18 February 20

Re: [Numpy-discussion] f2py-users list not working

2009-11-04 Thread George Nurser
2009/11/4 Robin : > On Wed, Nov 4, 2009 at 2:38 PM, George Nurser wrote: >> Fortran can accept preprocessor directives, but f2py cannot. >> You first need to preprocess a .F (or .F90) file to create a .f (or >> .f90) file which you then pass to f2py >> The way I prepro

Re: [Numpy-discussion] f2py-users list not working

2009-11-04 Thread George Nurser
gfortran -E -DINTSIZE=8 file.F -o outdir/file.f The outdir is necessary in a case-insensitive file system (like default mac OSX) to prevent the .f files overwriting the .F file. Alternatively, it may be possible to use some other suffix than .f, but I've not tried that. Then f2py file.f G

Re: [Numpy-discussion] strange performance on mac 2.5/2.6 32/64 bit

2009-11-04 Thread George Nurser
n 90 introduced the INTERFACE block, which allows you to use different variable types as arguments to what appears externally to be the same routine. It then feeds the arguments to the appropriate version of the routine. I don't think f2py supports this, but i

Re: [Numpy-discussion] numpy build/installation problems ?

2009-10-19 Thread George Nurser
I had the same 4 errors in genfromtext yesterday when I upgraded numpy r 7539. mac os x python 2.5.2. --George. 2009/10/19 Pierre GM : > > On Oct 19, 2009, at 10:40 AM, josef.p...@gmail.com wrote: > >> I wanted to finally upgrade my numpy, so I can build scipy trunk >> again, but I get test fai

Re: [Numpy-discussion] MFDatasets and NetCDF4

2009-09-25 Thread George Nurser
2009/9/25 David Huard : > Hi George, > > On Fri, Sep 25, 2009 at 6:55 AM, George Nurser > wrote: >> >> Hi, >> I hope this is the right place to ask this. >> I've found the MFDataset works well in reading NetCDF3 files, but it >> appears that it doesn

[Numpy-discussion] MFDatasets and NetCDF4

2009-09-25 Thread George Nurser
ataset for NetCDF4 files sometime? It would be very useful. --George Nurser. ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

[Numpy-discussion] numpy/scipy/matplotlib + 10.6 + Apple python 2.6.1

2009-09-07 Thread George Nurser
There are some interesting instructions on how to make this work at http://blog.hyperjeff.net/?p=160. However I'm not sure that the recommendation to rename the Apple-supplied version of numpy is consistent with previous advice I've seen on this mailing list. --Geo

Re: [Numpy-discussion] svn numpy not building on osx 10.5.6, python.org python 2.5.2

2009-06-07 Thread George Nurser
Thanks for the quick fix. 2009/6/7 David Cournapeau : > George Nurser wrote: >> Sorry, I should have said that I'd always deleted the build directories. >> I now have a better idea about what the problem is. >> >> python setup.py config_fc --fcompiler=gnu95

Re: [Numpy-discussion] svn numpy not building on osx 10.5.6, python.org python 2.5.2

2009-06-07 Thread George Nurser
architecture binaries? But I'm puzzled as to why specifying the fortran compiler should make any difference -- I understood it isn't used to compile numpy. --George. 2009/6/7 David Cournapeau : > George Nurser wrote: >> running config_fc >> unifing config_fc, config,

[Numpy-discussion] svn numpy not building on osx 10.5.6, python.org python 2.5.2

2009-06-06 Thread George Nurser
Hi, the current svn version 7039 isn't compiling for me. Clean checkout, old numpy directories removed from site-packages.. Same command did work for svn r 6329 [george-nursers-macbook-pro-15:~/src/numpy] agn% python setup.py config_fc --fcompiler=gnu95 build_clib --fcompiler=gnu95 build_ext --fco

Re: [Numpy-discussion] numpy and the ACML

2009-01-25 Thread George Nurser
ries is deep in distutils: Configuration in numpy/distutils/misc_util.py George. 2009/1/24 Gideon Simpson : > That's not working for me. Any thoughts on how to troubleshoot it? > -gideon > > On Jan 24, 2009, at 6:18 PM, George Nurser wrote: > >> I did manage to get it working.

Re: [Numpy-discussion] numpy and the ACML

2009-01-24 Thread George Nurser
I did manage to get it working. I remember that both libcblas.a (or a link to it) and libacml.so had to be in the same directory. Also I had to comment out lines 399-400 of setup.py: # if ('NO_ATLAS_INFO',1) in blas_info.get('define_macros',[]): # return None # dotblas

Re: [Numpy-discussion] URGENT: Re: 1.1.0rc1, Mac Installer: please test it

2008-05-21 Thread George Nurser
2008/5/21 Pierre GM <[EMAIL PROTECTED]>: > Mmh, wait a minute: > * There shouldn't be any mstats.py nor morestats.py in numpy.ma any longer: I > moved the packages to scipy.stats along their respective unittests. Right. I hadn't deleted the previous /Library/Frameworks/Python.framework/Versions/2.

Re: [Numpy-discussion] URGENT: Re: 1.1.0rc1, Mac Installer: please test it

2008-05-21 Thread George Nurser
Hmm. I also get some problems with test(all=True) 2 failures (though they look spurious to me) + 18 errors. Intel MBP, 10.5.2, macPython 2.5.2, apple gcc 4.0.1 George Nurser. >>> numpy.test(all=True) Numpy is installed in /Library/Frameworks/Python.framework/Versions/2.5/lib/pytho

Re: [Numpy-discussion] 1.1.0rc1 OSX Installer - please test

2008-05-21 Thread George Nurser
for me. Intel, MBP, 10.5.2. >>> import numpy >>> numpy.__version__ '1.1.0rc1' >>> numpy.test(10) ran 1005 tests in 2.290s OK George Nurser. ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] numpy setup.py too restrictive, prevents use of fblas with cblas

2008-04-17 Thread George Nurser
as library himself, either. I have filed a ticket on this. Yes, it is the cblas interface from http://www.netlib.org/blas/blast-forum/cblas.tgz. Not having to build the interface would indeed make it a whole lot easier. George Nurser. ___ Numpy-di

[Numpy-discussion] numpy setup.py too restrictive, prevents use of fblas with cblas

2008-04-16 Thread George Nurser
blas be removed for v1.1? Should I file a bug here? Regards, George Nurser. ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] array allocation using tuples gives views of same array

2007-11-15 Thread George Nurser
On 15/11/2007, Timothy Hochberg <[EMAIL PROTECTED]> wrote: > > > On Nov 15, 2007 9:11 AM, Hans Meine <[EMAIL PROTECTED]> wrote: > > Am Donnerstag, 15. November 2007 16:29:12 schrieb Warren Focke: > > > > > On Thu, 15 Nov 2007, George Nurser wrote: > &

[Numpy-discussion] array allocation using tuples gives views of same array

2007-11-15 Thread George Nurser
I tried the (as I thought) nice compact form In [60]: a,b = (zeros((2,)),)*2 But... In [61]: b[0] = 2 In [62]: a Out[62]: array([ 2., 0.]) a and b are the _same_ array But In [68]: a,b = (zeros((2,)),zeros((2,))) In [69]: b[0] = 2 In [70]: a Out[70]: array([ 0., 0.]) is OK. a & b are

Re: [Numpy-discussion] interrupted svn updates

2007-05-11 Thread George Nurser
svn.scipy.org) I tried a fresh checkout in a new directory, so the problems can't be here. Regards, George Nurser. ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] interrupted svn updates

2007-05-11 Thread George Nurser
I should clear? Many thanks for loooking at this. George Nurser. ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion

[Numpy-discussion] interrupted svn updates

2007-05-11 Thread George Nurser
tion was closed by server. (http://svn.scipy.org) 2nd try produced 4 more & then died with above error 3rd & 4th just dies with above error. Any ideas -- is thw problem here or there? Many thanks, George Nurser. ___ Numpy-discussion mailing lis

Re: [Numpy-discussion] New release of pycdf package ported to NumPy

2007-02-21 Thread George Nurser
Hi Andre, I've downloaded bpycdf and it works very nicely with numpy; thanks very much for all your effort. One small problem; I'm probably being stupid, but I cannot see how to set a _Fillvalue as Float32. regards, George Nurser. On 12/02/07, Andre Gosselin <[EMAIL PROTECTED

Re: [Numpy-discussion] building NumPy with Intel CC & MKL (solved!)

2007-01-25 Thread George Nurser
ing python, and for Perhaps compiling python itself with icc might give a useful speedup. Apparently somebody managed this for python 2.3 in 2003: http://mail.python.org/pipermail/c++-sig/2003-October/005824.html --George Nurser. ___ Numpy-discussion mailing