Hello numpy fans,
a patch and pull request for ticket 1573 has been posted, giving numpy.savetxt
the possibility to save complex arrays to text files. Formatting of the output
is supported roughly along the same lines as for real numbers.
https://github.com/numpy/numpy/pull/172
Share, test,
Hi,
I don't know about documentation, but I always use syntax like
zeros(10, dtype=numpy.float64)
where you have dtypes like
numpy.int8
numpy.uint32
numpy.complex128 # == two numpy.float64, one for real, one for imag.
etc. This is usually less confusing to my mind. One caveat: floats above 64
On 30. nov. 2011, at 12:09, Giovanni Plantageneto wrote:
Dear all,
I have a simple question. I would like to have all the parameters of a
model written in a configuration file (text), and I would like to have
all the parameters in the file automatically defined inside a program.
I find
On 6. des. 2011, at 17:32, Roger Binns wrote:
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
On 06/12/11 01:58, Pauli Virtanen wrote:
I think this cannot be helped --- it does not make sense to explain
basic Numpy concepts in every docstring, especially `axis` and `shape`
are very
On 15. jan. 2012, at 01:21, josef.p...@gmail.com wrote:
On Sat, Jan 14, 2012 at 5:25 PM, Benjamin Root ben.r...@ou.edu wrote:
On Sat, Jan 14, 2012 at 4:16 PM, Benjamin Root ben.r...@ou.edu wrote:
On Sat, Jan 14, 2012 at 4:12 PM, Charles R Harris
charlesr.har...@gmail.com wrote:
This
If by store you mean store on disk, I recommend h5py datasets and
attributes. Reportedly pytables is also good but I don't have any
first hand experience there. Both python modules use the hdf5 library,
written in C/C++/Fortran.
Paul
On Wed, Jan 25, 2012 at 7:47 PM, Val Kalatsky
On 27. jan. 2012, at 14:52, Chao YUE wrote:
Dear all,
suppose I have a ndarray a:
In [66]: a
Out[66]: array([0, 1, 2, 3, 4])
how can use it as 5X1 array without doing a=a.reshape(5,1)?
Several ways, this is one, although not much simpler.
In [6]: a
Out[6]: array([0, 1, 2, 3, 4])
In
An example I really like is LibreOffice's get involved page.
http://www.libreoffice.org/get-involved/
Producing something similar for NumPy will take some work, but I believe it's
needed.
Speaking as someone who has contributed to numpy in a microscopic fashion, I
agree completely. I
In the language wars, I have one question. Why is Fortran not being considered?
Fortran already implements many of the features that we want in NumPy:
- slicing and similar operations, at least some of the fancy indexing kind
- element-wise array operations and function calls
- array
On 20. feb. 2012, at 16:29, Sturla Molden wrote:
Den 20.02.2012 08:35, skrev Paul Anton Letnes:
In the language wars, I have one question. Why is Fortran not being
considered? Fortran already implements many of the features that we want in
NumPy:
Yes ... but it does not make Fortran
As others on this list, I've also been confused a bit by the prolific numpy
interfaces to reading text. Would it be an idea to create some sort of object
oriented solution for this purpose?
reader = np.FileReader('my_file.txt')
reader.loadtxt() # for backwards compat.; np.loadtxt could
Hi,
this probably does not help with your problem. However, I would recommend
changing your fortran code to:
subroutine print_bit_array(bits)
use iso_fortran_env
integer(kind=int8),intent(in),dimension(:)::bits
print*,'bits = ',bits
end subroutine print_bit_array
In that way you could
Ah, come to think of it, I think that f2py only supports literal kind values.
Maybe that's your problem.
Paul
On 17. apr. 2012, at 07:58, Sameer Grover wrote:
On Tuesday 17 April 2012 11:02 AM, John Mitchell wrote:
Hi,
I am using f2py to pass a numpy array of type numpy.int8 to fortran.
you to mean that
'int8' is not a literal kind value while 1 and 8 are examples of literal kind
values.
Thanks,
John
On Tue, Apr 17, 2012 at 10:12 AM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote:
Ah, come to think of it, I think that f2py only supports literal kind values
On 21. apr. 2012, at 00:16, Drew Frank wrote:
On Fri, Apr 20, 2012 at 11:45 AM, Chris Barker chris.bar...@noaa.gov wrote:
On Fri, Apr 20, 2012 at 11:39 AM, Dag Sverre Seljebotn
d.s.seljeb...@astro.uio.no wrote:
Oh, right. I was thinking small as in fits in L2 cache, not small as
in a few
On 3. mai 2012, at 19:33, Moroney, Catherine M (388D) wrote:
A quick recap of the problem: a 128x512 array of 7-element vectors
(element), and a 5000-vector
training dataset (targets). For each vector in element, I want to find the
best-match in targets,
defined as minimizing the
On Fri, May 4, 2012 at 12:49 AM, Keith Goodman kwgood...@gmail.com wrote:
On Thu, May 3, 2012 at 3:12 PM, Moroney, Catherine M (388D)
catherine.m.moro...@jpl.nasa.gov wrote:
Here is the python code:
def single(element, targets):
if (isinstance(element, tuple)):
xelement =
Hi,
I'm getting a couple of errors when testing. System:
Arch Linux (updated today)
Python 3.2.3
gcc 4.7.0
(Anything else?)
I think that this error:
AssertionError: selectedrealkind(19): expected -1 but got 16
is due to the fact that newer versions of gfortran actually supports
precision this
All tests for 1.6.2rc1 pass on
Mac OS X 10.7.3
python 2.7.2
gcc 4.2 (Apple)
Great!
Paul
On 6. mai 2012, at 00:12, Charles R Harris wrote:
On Sat, May 5, 2012 at 2:56 PM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote:
Hi,
I'm getting a couple of errors when testing. System
On Sat, May 12, 2012 at 9:50 PM, Ralf Gommers
ralf.gomm...@googlemail.com wrote:
On Sun, May 6, 2012 at 12:12 AM, Charles R Harris
charlesr.har...@gmail.com wrote:
On Sat, May 5, 2012 at 2:56 PM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote:
Hi,
I'm getting a couple of errors
On Mon, May 14, 2012 at 9:47 PM, Ralf Gommers
ralf.gomm...@googlemail.com wrote:
On Sun, May 13, 2012 at 1:14 PM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote:
On Sat, May 12, 2012 at 9:50 PM, Ralf Gommers
ralf.gomm...@googlemail.com wrote:
On Sun, May 6, 2012 at 12:12 AM
On 4. juni 2012, at 16:27, Thouis (Ray) Jones wrote:
On Fri, Jun 1, 2012 at 6:56 PM, Chris Withers ch...@simplistix.co.uk wrote:
On 01/06/2012 16:39, Benjamin Root wrote:
import numpy
numpy.zeros(10)[-123]
Traceback (most recent call last):
File stdin, line
On 7. juni 2012, at 10:30, Thouis (Ray) Jones wrote:
I've opened a PR at https://github.com/numpy/numpy/pull/296 for discussion.
A typical result
np.zeros((3,3))[[1,2,3]]
Traceback (most recent call last):
File stdin, line 1, in module
IndexError: index 3 is out of bounds for axis 0:
Hello,
I don't know exactly what went wrong. I'd start out my debugging by
1) which python # see whether you're running apple's python in /usr/bin/python,
or the one you tried to install
2) which easy_install # did you run Apple-python's easy_install, or the one you
(tried to) installed?
3) If
Hi,
are you sure that you want g77 and not gfortran? If you want gfortran, you
should pass the
--fcompiler=gnu95
flag to setup.py.
Which redhat version are you building on? (I don't know red hat well enough to
comment, but perhaps someone else do...)
Paul
On 6. juli 2012, at 03:00,
However, I got the following error message:
error: Command /usr/bin/g77 -g -Wall -g -Wall -shared
build/temp.linux-x86_64-2.6/build/src.linux-x86_64-2.6/scipy/integrate/vodemodule.o
build/temp.linux-x86_64-2.6/build/src.linux-x86_64-2.6/fortranobject.o
-L/home/username/lib/ -L/usr/lib64
On 10. aug. 2012, at 09:54, Mark Bakker wrote:
I am giving this a second try. Can anybody help me out?
I think there is a problem with assigning a 1D complex array of length one
to a position in another complex array.
Example:
a = ones(1,'D')
b = ones(1,'D')
a[0] = b
For convenience and clarity, this is the diff in question:
https://github.com/Dynetrekk/numpy-1/commit/5bde67531a2005ef80a2690a75c65cebf97c9e00
And this is my numpy fork:
https://github.com/Dynetrekk/numpy-1/
Paul
On Sun, Sep 30, 2012 at 4:14 PM, Paul Anton Letnes
paul.anton.let...@gmail.com
On 3. okt. 2012, at 17:48, Wes McKinney wrote:
On Monday, October 1, 2012, Chris Barker wrote:
Paul,
Nice to see someone working on these issues, but:
I'm not sure the problem you are trying to solve -- accumulating in a
list is pretty efficient anyway -- not a whole lot overhead.
On 1. okt. 2012, at 21:07, Chris Barker wrote:
Paul,
Nice to see someone working on these issues, but:
I'm not sure the problem you are trying to solve -- accumulating in a
list is pretty efficient anyway -- not a whole lot overhead.
Oh, there's significant overhead, since we're not
On 3. okt. 2012, at 18:22, Chris Barker wrote:
On Wed, Oct 3, 2012 at 9:05 AM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote:
I'm not sure the problem you are trying to solve -- accumulating in a
list is pretty efficient anyway -- not a whole lot overhead.
Oh, there's significant
Hi Pierre,
first I'd like to congratulate you on a product which seems to finally solve
all the problems I have at work regarding python scripting. A portable, fully
featured python distribution is indeed very useful.
One question though: Is there a way to run winpython's IPython under cygwin,
On 29. okt. 2012, at 11:29, Radek Machulka wrote:
Hi,
is there a way how to save more arrays into single npy (npz if possible) file
in loop? Something like:
fw = open('foo.bar', 'wb')
while foo:
arr = np.array(bar)
np.savez_compressed(fw, arr)
fw.close()
Or some
On 4. jan. 2013, at 21:42, m...@eml.cc wrote:
Hiall,
I am trying to embed numerical code in a mexFunction,
as called by MATLAB, written as a Cython function.
NumPy core functions and BLAS work fine, but calls to LAPACK
function such as SVD seem to be made against to MATLAB's linked
On 17.01.2013 04:43, Patrick Marsh wrote:
Thanks, everyone for chiming in. Now that I know this behavior
exists, I can explicitly prevent it in my code. However, it would be
nice if a warning or something was generated to alert users about the
inconsistency between var += ... and var = var
On Tue, Nov 23, 2010 at 3:37 AM, Sebastian Walter
sebastian.wal...@gmail.com wrote:
On Tue, Nov 23, 2010 at 11:17 AM, Gael Varoquaux
gael.varoqu...@normalesup.org wrote:
On Tue, Nov 23, 2010 at 11:13:23AM +0100, Sebastian Walter wrote:
I'm not familiar with dichotomy optimization.
On 3. des. 2010, at 16.24, Fabian Pedregosa wrote:
Hi all.
Macports installs gfortran as part of the gcc package, but names it
gfortran-mp-$version, without providing a symbolic link to a default
gcfortran executable, and thus numpy.distutils is unable to find the
right executable.
The
Mabe I am wrong somehow, but in my experience the easiest install of scipy
is 'port install py26-scipy'. For new users, I do not see why one would
recommend to build manually from source? Macports can do it for you,
automagically...
Paul
4. des.. 2010 15.04 Ralf Gommers
On 12. jan. 2011, at 16.40, dstaley wrote:
Zachary Pincus-2 wrote:
Is it possible to use a variable in an array name? I am looping
through a
bunch of calculations, and need to have each array as a separate
entity.
I'm pretty new to python and numpy, so forgive my ignorance. I'm
Hi!
This little snippet of code tricked me (in a more convoluted form). The *=
operator does not change the datatype of the left hand side array. Is this
intentional? It did fool me and throw my results quite a bit off. I always
assumed that 'a *= b' means exactly the same as 'a = a * b' but
Heisann!
On 16. mars 2011, at 14.30, Dag Sverre Seljebotn wrote:
On 03/16/2011 02:24 PM, Paul Anton Letnes wrote:
Hi!
This little snippet of code tricked me (in a more convoluted form). The *=
operator does not change the datatype of the left hand side array. Is this
intentional? It did
On 16. mars 2011, at 15.49, Chris Barker wrote:
On 3/16/11 6:34 AM, Charles R Harris wrote:
On Wed, Mar 16, 2011 at 7:24 AM, Paul Anton Letnes
Yes, it is intentional. Numpy is more C than Python in this case,
I don't know that C has anything to do with it -- the *= operators were
added
On 16. mars 2011, at 15.57, Dag Sverre Seljebotn wrote:
On 03/16/2011 02:35 PM, Paul Anton Letnes wrote:
Heisann!
Hei der,
On 16. mars 2011, at 14.30, Dag Sverre Seljebotn wrote:
On 03/16/2011 02:24 PM, Paul Anton Letnes wrote:
Hi!
This little snippet of code tricked me (in a more
This comes up for discussion on a fairly regular basis. I tend towards the
more warnings side myself, but you aren't going to get the current behavior
changed unless you can convince a large bunch of people that it is the right
thing to do, which won't be easy. For one thing, a lot of
On 20. mars 2011, at 16.08, Ben Smith wrote:
So, in addition to my computer science work, I'm a PhD student in econ. Right
now, the class is using GAUSS for almost everything. This sort of pisses me
off because it means people are building libraries of code that become
valueless when
I'm no expert, but I just pulled off the scipy+numpy+GotoBLAS2 installation.
From what I gather, the Makefile for libgoto2 downloads and compiles the
generic lapack from netlib. It also wraps lapack into libgoto2.so/.a. I
believe the idea is as long as the BLAS implementation is fast(TM), the
One way is to use python dicts. Pseudocode:
arrays = {}
integer = 0
for file in allfiles:
integer +=1
data = file.read()
arrays['a' + '%d' % integer] = data
I would consider using the filename as the key instead of a1, a2, etc. That way
you have that information readily
Hi!
I have had a look at the list of numpy.loadtxt tickets. I have never
contributed to numpy before, so I may be doing stupid things - don't be afraid
to let me know!
My opinions are my own, and in detail, they are:
1752:
I attach a possible patch. FWIW, I agree with the request. The
Hi Derek!
On 26. mars 2011, at 15.48, Derek Homeier wrote:
Hi again,
On 26 Mar 2011, at 15:20, Derek Homeier wrote:
1562:
I attach a possible patch. This could also be the default
behavior to my mind, since the function caller can simply call
numpy.squeeze if needed. Changing default
On 26. mars 2011, at 21.44, Derek Homeier wrote:
Hi Paul,
having had a look at the other tickets you dug up,
My opinions are my own, and in detail, they are:
1752:
I attach a possible patch. FWIW, I agree with the request. The
patch is written to be compatible with the fix in
Hi.
When looking at the loadtxt/savetxt tickets, I noticed that the 're' module is
imported in an odd place. I therefore suggest that this import is moved to the
top of the file, in order to gather these as much as possible. I find the code
easier to read then. After all, there is no 'try /
On 4. apr. 2011, at 15.34, Charles R Harris wrote:
On Sun, Apr 3, 2011 at 4:35 AM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote:
Hi.
When looking at the loadtxt/savetxt tickets, I noticed that the 're' module
is imported in an odd place. I therefore suggest that this import
On 4. apr. 2011, at 16.42, Charles R Harris wrote:
On Mon, Apr 4, 2011 at 8:29 AM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote:
On 4. apr. 2011, at 15.34, Charles R Harris wrote:
On Sun, Apr 3, 2011 at 4:35 AM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote
On 5. apr. 2011, at 15.39, Alan G Isaac wrote:
On 4/5/2011 9:26 AM, François Steinmetz wrote:
It does not change the dtype, 'int' is just interpreted as 'int64' :
So the meaning of 'int' is system specific?
import numpy as np; a=np.eye(2,dtype='int'); a.dtype
dtype('int32')
On 25. apr. 2011, at 19.57, Pauli Virtanen wrote:
On Mon, 25 Apr 2011 10:16:13 -0700, Rob Beezer wrote:
[clip]
Many more details and complete transcripts are at:
http://trac.sagemath.org/sage_trac/ticket/11248
Any thoughts or advice to help us understand this would be greatly
appreciated.
On Mac OS X 10.6.7 with macports python 2.7 I get as follows:
OK (KNOWNFAIL=3, SKIP=1)
Python version:
Python 2.7.1 (r271:86832, Jan 12 2011, 10:44:27)
[GCC 4.2.1 (Apple Inc. build 5664)] on darwin
% gfortran --version
GNU Fortran (GCC) 4.5.3
Well done to everyone working on numpy 1.6!
Paul.
On 4. mai 2011, at 17.34, Derek Homeier wrote:
Hi Paul,
I've got back to your suggestion re. the ndmin flag for loadtxt from a few
weeks ago...
On 27.03.2011, at 12:09PM, Paul Anton Letnes wrote:
1562:
I attach a possible patch. This could also be the default
behavior to my mind
On 4. mai 2011, at 20.33, Benjamin Root wrote:
On Wed, May 4, 2011 at 7:54 PM, Derek Homeier
de...@astro.physik.uni-goettingen.de wrote:
On 05.05.2011, at 2:40AM, Paul Anton Letnes wrote:
But: Isn't the numpy.atleast_2d and numpy.atleast_1d functions written for
this? Shouldn't we
On 5. mai 2011, at 08.49, Benjamin Root wrote:
On Wed, May 4, 2011 at 11:08 PM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote:
On 4. mai 2011, at 20.33, Benjamin Root wrote:
On Wed, May 4, 2011 at 7:54 PM, Derek Homeier
de...@astro.physik.uni-goettingen.de wrote
Hi,
it seems that I have found a bug in numpy.ndarray. numpy 1.5.1, python 2.7.1
from macports on mac os x 10.6.7. I got the same error on Fedora 14 with numpy
1.4.1 and python 2.7. Appending a [0] to the last line solves the problem.
% python testcrash.py
On 8. mai 2011, at 17.32, Warren Weckesser wrote:
On Sun, May 8, 2011 at 7:23 PM, Charles R Harris charlesr.har...@gmail.com
wrote:
On Sun, May 8, 2011 at 3:15 PM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote:
Hi,
it seems that I have found a bug in numpy.ndarray. numpy
Dear numpy developers, users, and fans:
I wish for a numpy that will magically detect all my errors. Barring the
invention of a numpy version which knows what I need, not what I ask of it, I
have the following, simpler wish: if one takes myarray.real or myarray.imag of
a float (or possibly
Hi,
I had the same problem. I think this might work:
FFLAGS='-fopenmp' f2py -c (etc)
The thing is that f2py doesn't let you pass the -fopenmp flag at the right time
to the compiler, so you have to use some sort of environment variable trick. By
the way, as far as I know, this is the case also
Hi!
In my quest for a bug-free numpy I have, I think, fixed ticket 1793.
https://github.com/numpy/numpy/pull/123
Enjoy - feedback welcome, of course.
Cheers,
Paul.
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
(A pull request has been submitted on github, but I'm posting here so people
can discuss the user interface issues.)
As of now, the fmt= kwarg kan be (for complex dtype):
a) a single specifier, fmt='%.4e', resulting in numbers formatted like '
(%s+%sj)' % (fmt, fmt)
b) a full string specifying
Looks like you have done some great work! I've been using f2py in the past, but
I always liked the idea of cython - gradually wrapping more and more code as
the need arises. I read somewhere that fortran wrapping with cython was coming
- dare I ask what the status on this is? Is it a goal for
On 10. aug. 2011, at 21.03, Gael Varoquaux wrote:
On Wed, Aug 10, 2011 at 04:01:37PM -0400, Anne Archibald wrote:
A 1 Gb text file is a miserable object anyway, so it might be desirable
to convert to (say) HDF5 and then throw away the text file.
+1
G
+1 and a very warm recommendation
On 13. aug. 2011, at 20.42, Charles R Harris wrote:
2011/8/11 Dmitrey tm...@ukr.net
bug in KUBUNTU 11.04, latest numpy git snapshot build with Python3
import numpy
Traceback (most recent call last):
File stdin, line 1, in module
File
On 19. aug. 2011, at 19.57, Ian wrote:
Right. I'm new to NumPy so I figured I'd check if there was some nifty way of
preserving the shape without storing it in the database that I hadn't
discovered yet. No worries, I'll store the shape alongside the array. Thanks
for the reply.
I love
Hi!
On 21. aug. 2011, at 00.18, Chris Withers wrote:
Hi All,
I've got a tree of nested dicts that at their leaves end in numpy arrays
of identical sizes.
What's the easiest way to persist these to disk so that I can pick up
with them where I left off?
Probably giving them names like
On 25. aug. 2011, at 23.49, David Warde-Farley wrote:
On 2011-08-25, at 2:42 PM, Chris.Barker wrote:
On 8/24/11 9:22 AM, Anthony Scopatz wrote:
You can use Python pickling, if you do *not* have a requirement for:
I can't recall why, but it seem pickling of numpy arrays has been
I recently got into trouble with these calculations (although I used scipy). I
actually got segfaults and bus errors. The solution for me was to not link
against ATLAS, but rather link against Apple's blas/lapack libraries. That got
everything working again. I would suggest trying to install
?
On Mon, Aug 29, 2011 at 9:31 AM, Paul Anton Letnes
paul.anton.let...@gmail.com wrote:
I recently got into trouble with these calculations (although I used
scipy). I actually got segfaults and bus errors. The solution for me was
to not link against ATLAS, but rather link against Apple's blas
You can use the BLAS and LAPACK environment variables.
export BLAS=/path/to/libatlas.so
export LAPACK=/path/to/libatlas.so
python setup.py build
I've recently had problems with ATLAS solving equation systems
incorrectly for certain inputs with no adequate explanation.
Re-running the same
74 matches
Mail list logo