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
and the
rationale for that decision was performance. What follows is that masked
array with the default nomask attribute behaves a regular array (hence
the nan), having a placeholder for mask to be set later, if needed.
That tripped me recently, I had Cython code which relied on shapes of
d
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
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
ssion/2012-October/064259.html ).
I really do think that the Einstein summation convention is a really clear
way of doing array operations, and its use should be encouraged.
--George
On 9 November 2012 06:18, Nicolas SCHEFFER wrote:
> Fred,
>
> Thanks for the advice.
> The code wil
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
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
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
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.
>>
&
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
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
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
I couldn't find an answer to my newbie question, so I'm posting it here.
I have:
a=numpy.array([[1,2],[3,4]])
b=numpy.array([[5,6],[7,8]])
Via broadcasting, I know that
a*[[5],[7]]=numpy.array([[5,10],[21,28]])
Being a recent convert from MATLAB, I expected the same result from a*b[:,0],
assum
-- 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
I'm suggesting writing a *new* Fortran interface, coupled with f2py.
The original library just needs to be linked to the new .so generated
by f2py. I am hoping (perhaps optimistically) that can be done in the
Fortran compilation...
--George.
On 18 February 2010 10:56, Matthieu Brucher
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 Feb
So with bincount I can exchange a loop over P for a loop over M? I
guess for me that is still really helpful.
Thanks!
- George
On Mon, Dec 21, 2009 at 6:35 AM, Pauli Virtanen wrote:
> Mon, 21 Dec 2009 09:35:08 +, Neil wrote:
> [clip]
>> I'm also interested to see if there a
vectorized and efficient way to do this without making
the obvious python for loop I included above? For what I am doing, P
is usually quite large. I am most interested in a clever use of
numpy/scipy commands and Python and not Cython.
Thanks in advance for any suggestions.
- George
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
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
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
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
>>
ere any drawbacks to this approach?
George
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion
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
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
Sturla Molden molden.no> writes:
>
> Erik Tollerud skrev:
> >> NumPy arrays on the GPU memory is an easy task. But then I would have to
> >> write the computation in OpenCL's dialect of C99?
> > This is true to some extent, but also probably difficult to do given
> > the fact that paralellizabl
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
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
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,
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 buil
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.
library_dirs = /noc/users/agn/ext/AMD64/acml/ifort64/lib
include_dirs = /noc/users/agn/ext/AMD64/acml/ifort64/include
Both libcblas.a (or a link to it) and libacml.so are in
/noc/users/agn/ext/AMD64/acml/ifort64/lib
HTH. George.
2009/1/24 Gideon Simpson :
> I've tried building CBLAS, wh
using sizeof on different platforms highlights
the problem.
Thanks.
George.
___
Numpy-discussion mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion
, sob, sob
sob). The problem was in the file MPL_isnan.h on line 26.
Consider the issue closed.
Thanks.
George.
___
Numpy-discussion mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion
--
Ran 1592 tests in 10.704s
OK (KNOWNFAIL=3, SKIP=1)
Thanks.
George.
___
Numpy-discussion
. I think the AMD ACML is
compiled using the Intel FORTRAN compiler, but will that effect it??
Anyway, I'll put an effort into constructing an example, but it will have to be
when I am back at the office from my vacation.
Cheers. Thanks.
George.
___
David Cournapeau gmail.com> writes:
>
> On Tue, Dec 9, 2008 at 12:50 AM, George Goussard emss.co.za> wrote:
> > Hello.
> >
> >
> >
> > I have been battling with the following error for the past week. The output
> > from the terminal is:
>
--
Ran 1592 tests in 10.704s
OK (KNOWNFAIL=3, SKIP=1)
Thanks.
George.
Message: 3
Date: Tue, 9 Dec 2008 02:43:25 +0900
From: "David Cournapeau
Initially MPL plots a graph but when you try to interact with the widget(for
example resize) then the output is displayed and the MPL figure is not updated.
Everything works with Windows 32-bit. Linux 32-bit and 64-bit are working
correctly.
Any ideas would be helpful.
Thanks.
George.
Python.framework/Versions/2.5/lib/python2.5/site-packages/numpy
After doing that, I get just 12 errors; no failures.
The errors are related to
* ma/tests/test_mrecords.py
* ma/tests/test_old_ma.py
* ma/tests/test_subclassing.py
-George.
Numpy version 1.1.0rc1
Python version 2.5.2 (r252:60911, Fe
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
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
MKL and the GOTO BLAS and allow them,
> too. It might also be worth including the appropriate subset of the
> cblas code provided in the tarball such that we can use any
> accelerated FORTRAN BLAS without the standard cblas interface. Then
> George wouldn't have the build the cbl
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
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:
> &
cate many arrays?
Regards, George.
___
Numpy-discussion mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion
On Sep 5, 12:29 pm, Francesc Altet <[EMAIL PROTECTED]> wrote:
> A Wednesday 05 September 2007, George Sakkis escrigué:
>
>
>
> > I was surprised to see that an in-place modification of a 2-d array
> > turns out to be slower from the respective non-mutating operatio
4701226007, 0.19477587235249172]
1 [0.073787590322233698, 1.9234369172618306]
So the 2-d in-place modification time grows linearly with the array
size but the 1-d operations are much more efficient, despite
allocating new arrays while doing so. What gives ?
George
___
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
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
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
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
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
53 matches
Mail list logo