Re: [Numpy-discussion] Accelerate or OpenBLAS for numpy / scipy wheels?

On 28/06/2016 18:50, Ralf Gommers wrote: On Tue, Jun 28, 2016 at 5:50 PM, Chris Barker > wrote: > This doesn't really matter too much imho, we have to support Accelerate > either way. do we? -- so if we go OpenBlas, and

[Numpy-discussion] how is toolchain selected for compiling (OS X with python.org build)?

This isn't, strictly speaking, a numpy question, but I suspect it's something that numpy devs and users have some insight into. I am trying to compile an extension that requires a fairly advanced c++ compiler. Using the built-in apple python, it defaults to the latest clang from apple, and it

Re: [Numpy-discussion] [RFC] should we argue for a matrix power operator, @@?

On 16/03/2014 01:31, josef.p...@gmail.com wrote: On Sat, Mar 15, 2014 at 8:47 PM, Warren Weckesser warren.weckes...@gmail.com mailto:warren.weckes...@gmail.com wrote: On Sat, Mar 15, 2014 at 8:38 PM, josef.p...@gmail.com mailto:josef.p...@gmail.com wrote: I think I

Re: [Numpy-discussion] Dropping support for, Accelerate/veclib?

On 11/06/2013 22:11, Chris Barker - NOAA Federal wrote: On Tue, Jun 11, 2013 at 1:28 PM, Ralf Gommers ralf.gomm...@gmail.com wrote: The binaries will still be built against python.org Python, so there shouldn't be an issue here. Same for building from source. My point was that it's nice to be

Re: [Numpy-discussion] Raveling, reshape order keyword unnecessarily confuses index and memory ordering

Proposal - * Deprecate the use of C and F meaning backwards and forwards index ordering for ravel, reshape * Prefer Z and N, being graphical representations of unraveling in 2 dimensions, axis1 first and axis0 first respectively (excellent naming idea by Paul Ivanov) What

Re: [Numpy-discussion] Getting non-normalized eigenvectors from generalized eigenvalue solution?

Just to be completely clear, there is no such thing as a non-normalized eigenvector. An eigenvector is only determined *up to a scalar normalization*, which is obvious from the eigenvalue equation: A v = l v where A is the matrix, l is the eigenvalue, and v is the eigenvector. Obviously v is

Re: [Numpy-discussion] ANN: NumPy 1.6.0 [Problems on OS X w/ Python 2.6 and XCode 4]

Dear all, I have OSX 10.6.7, XCode 4, and Python.org python 2.6.6 and 2.7.1, where 2.7 is 64-bit. With 2.7, easy_install successfully compiles and installs the package, both over the web and with an explicit download. With 2.6, there seems to be a problem with attempting to compile the PPC

Re: [Numpy-discussion] ANN: NumPy 1.6.0 [Problems on OS X w/ Python 2.6 and XCode 4]

On 16/05/2011 18:45, Ralf Gommers wrote: On Mon, May 16, 2011 at 12:41 PM, Andrew Jaffe a.h.ja...@gmail.com mailto:a.h.ja...@gmail.com wrote: Dear all, I have OSX 10.6.7, XCode 4, and Python.org python 2.6.6 and 2.7.1, where 2.7 is 64-bit. With 2.7, easy_install

Re: [Numpy-discussion] Computing the norm of an array of vectors

On 08/02/2011 16:44, Ben Gamari wrote: I have an array of (say, row) vectors, v = [ [ a1, a2, a3 ], [ b1, b2, b3 ], [ c1, c2, c3 ], ... ] What is the optimal way to compute the norm of each vector, norm(v)**2 = [ [ a1**2 + a2**2 + a3**2 ],

Re: [Numpy-discussion] read ascii file from complex fortran format() -- genfromtxt

David H. On Tue, Sep 21, 2010 at 8:25 AM, Andrew Jaffe a.h.ja...@gmail.com mailto:a.h.ja...@gmail.com wrote: Hi all, I've got an ascii file with a relatively complicated structure, originally written by fortran with the format: 135format(a12,1x,2(f10.5,1x),i3,1x

[Numpy-discussion] read ascii file from complex fortran format() -- genfromtxt

Hi all, I've got an ascii file with a relatively complicated structure, originally written by fortran with the format: 135format(a12,1x,2(f10.5,1x),i3,1x,4(f9.3,1x),4(i2,1x),3x, 1 16(f7.2,1x),i3,3x,f13.5,1x,f10.5,1x,f10.6,1x,i3,1x, 2 4(f10.6,1x), 2

[Numpy-discussion] long integers in loadtxt -- bad parsing?

Hi all, I am trying to read some 19-digit integers using loadtxt (or genfromtxt -- same problem). The numbers are smaller than the max int64 (and the max uint64 -- same problem with either one). Below, Out[184] shows that python has no problem with the conversion, but loadtxt gets the last

Re: [Numpy-discussion] long integers in loadtxt -- bad parsing?

Hi, I am trying to read some 19-digit integers using loadtxt (or genfromtxt -- same problem). The numbers are smaller than the max int64 (and the max uint64 -- same problem with either one). Below, Out[184] shows that python has no problem with the conversion, but loadtxt gets the last few

[Numpy-discussion] searchsorted for exact matches, not preserving order

Dear all, I've got two (integer) arrays, and I want to find the indices in the first one that have entries in the second. I.E. I want all idx s.t. there exists a j with a[idx]=b[j]. Here is my current implementation (with a = pixnums, b=surveypix) import numpy as np def matchPix(pixnums,

Re: [Numpy-discussion] searchsorted for exact matches, not preserving order

On 11/09/2009 08:33, Robert Kern wrote: On Fri, Sep 11, 2009 at 10:24, Andrew Jaffea.h.ja...@gmail.com wrote: Dear all, I've got two (integer) arrays, and I want to find the indices in the first one that have entries in the second. I.E. I want all idx s.t. there exists a j with a[idx]=b[j].

Re: [Numpy-discussion] pdf for multivariate normal function?

Hi, Charles R Harris wrote: On Thu, Jul 23, 2009 at 7:14 AM, per freem perfr...@gmail.com mailto:perfr...@gmail.com wrote: i'm trying to find the function for the pdf of a multivariate normal pdf. i know that the function multivariate_normal can be used to sample from the

Re: [Numpy-discussion] PEP: named axis

Pauli Virtanen wrote: Wed, 11 Feb 2009 22:21:30 +, Andrew Jaffe wrote: [clip] Maybe I misunderstand the proposal, but, actually, I think this is completely the wrong semantics for axis= anyway. axis= in numpy refers to what is also a dimension, not a column. I think the proposal

Re: [Numpy-discussion] how do I configure with gfortran

This is slightly off-topic, but probably of interest to anyone reading this thread: Is there any reason why we don't use --fcompiler=gfortran as an alias for --fcompiler=gfortran (and --fcompiler=g77 as an alias for --fcompiler=gnu, for that matter)? Those seem to me to be much more mnemonic

Re: [Numpy-discussion] flatten() without copy - is this possible?

dmitrey wrote: hi all. in the numpy for matlab users I read y = x.flatten(1) turn array into vector (note that this forces a copy) Is there any way to do the trick wthout copying? What are the problems here? Just other way of array elements indexing... One important question is