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

2016-06-29 Thread Andrew Jaffe
On 28/06/2016 18:50, Ralf Gommers wrote: On Tue, Jun 28, 2016 at 5:50 PM, Chris Barker mailto:chris.bar...@noaa.gov>> wrote: > This doesn't really matter too much imho, we have to support Accelerate > either way. do we? -- so if we go OpenBlas, and someone want to do a sim

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

2015-10-15 Thread Andrew Jaffe
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, @@?

2014-03-19 Thread Andrew Jaffe
On 16/03/2014 01:31, josef.p...@gmail.com wrote: > > > > On Sat, Mar 15, 2014 at 8:47 PM, Warren Weckesser > mailto:warren.weckes...@gmail.com>> wrote: > > > On Sat, Mar 15, 2014 at 8:38 PM, > wrote: > > I think I wouldn't use anything like @@ often eno

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

2013-06-13 Thread Andrew Jaffe
On 11/06/2013 22:11, Chris Barker - NOAA Federal wrote: > On Tue, Jun 11, 2013 at 1:28 PM, Ralf Gommers 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 able to have it

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

2013-04-02 Thread Andrew Jaffe
>> >> 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

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

2011-12-21 Thread Andrew Jaffe
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]

2011-05-16 Thread Andrew Jaffe
On 16/05/2011 18:45, Ralf Gommers wrote: > > On Mon, May 16, 2011 at 12:41 PM, Andrew Jaffe <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. > >

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

2011-05-16 Thread Andrew Jaffe
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 ar

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

2011-02-08 Thread Andrew Jaffe
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 +

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

2010-09-22 Thread Andrew Jaffe
g version is not yet available"! Andrew > > David H. > > On Tue, Sep 21, 2010 at 8:25 AM, Andrew Jaffe <mailto:a.h.ja...@gmail.com>> wrote: > > Hi all, > > I've got an ascii file with a relatively complicated structure, > originally written by

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

2010-09-21 Thread Andrew Jaffe
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 i2,1x,f5

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

2010-04-08 Thread Andrew Jaffe
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

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

2010-04-08 Thread Andrew Jaffe
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 few

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

2009-09-11 Thread Andrew Jaffe
On 11/09/2009 08:33, Robert Kern wrote: > On Fri, Sep 11, 2009 at 10:24, Andrew Jaffe 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. >> ther

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

2009-09-11 Thread Andrew Jaffe
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, surve

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

2009-07-23 Thread Andrew Jaffe
Hi, Charles R Harris wrote: > > On Thu, Jul 23, 2009 at 7:14 AM, per freem > 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 multivariate

Re: [Numpy-discussion] PEP: named axis

2009-02-11 Thread Andrew Jaffe
Pauli Virtanen wrote: > Wed, 11 Feb 2009 22:21:30 +0000, 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 al

Re: [Numpy-discussion] PEP: named axis

2009-02-11 Thread Andrew Jaffe
Robert Kern wrote: > On Fri, Feb 6, 2009 at 03:22, Stéfan van der Walt wrote: >> Hi Robert >> >> 2009/2/6 Robert Kern : This could be implemented but would require adding information to the NumPy array. >>> More than that, though. Every function and method that takes an axis >>> or reduc

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

2007-07-02 Thread Andrew Jaffe
I wrote: > 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)? But, sorry, of course

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

2007-07-02 Thread Andrew Jaffe
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 n

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

2007-06-04 Thread Andrew Jaffe
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 i

Re: [Numpy-discussion] test_multiarray.test_clip fails on Solaris 8 system

2007-04-02 Thread Andrew Jaffe
> The following test fails on a Solaris 8 system: > == > FAIL: check_basic (numpy.core.tests.test_multiarray.test_clip) > -- > Traceback (most