Hello all
These problems were fixed in SVN just after 1.0rc1 was released.
http://projects.scipy.org/scipy/numpy/changeset/3119
Regards,
Albert
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Achim Gaedke
> Sent: 17 September 2006 22:00
>
Hey Eric
> Hi again
>
> after some hours of debugging I finally (I think) found the problem:
>
> numpy.sum([[0,1,2],[2,3,4]])
> 24
>
> numpy.sum([[0,1,2],[2,3,4]],axis=0)
> array([2, 4, 6])
>
> numpy.sum([[0,1,2],[2,3,4]],axis=1)
> array([3, 9])
>
>
> Isn't the first line supposed to act as w
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Sebastian Haase
> Sent: 19 September 2006 19:24
> To: Discussion of Numerical Python
> Subject: [Numpy-discussion] arr.dtype.kind is 'i' for dtype=unit !?
>
> Hi,
> What are the possib
Hello all
I recently started looking at David Cournapeau's PyEM package, specifically
his implementation of the K-Means algorithm. He implemented part of this
algorithm with in pure Python version and also provided a Pyrex alternative
that is significantly faster (about 10 times with the data I te
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Travis Oliphant
> Sent: 03 October 2006 02:32
> To: Discussion of Numerical Python
> Subject: Re: [Numpy-discussion] Vectorizing code, for loops, and all that
>
> Travis Oliphant wrote
Hello all
> -Original Message-
> From: Travis Oliphant [mailto:[EMAIL PROTECTED]
> Sent: 03 October 2006 03:51
> To: Discussion of Numerical Python; Albert Strasheim
> Subject: Re: ***[Possible UCE]*** Re: [Numpy-discussion] Vectorizing
> code, for loops, and all that
&
Hello all
I recently started looking at David Cournapeau's PyEM package,
specifically his implementation of the K-Means algorithm. He
implemented part of this algorithm with in pure Python version and
also provided a Pyrex alternative that is significantly faster (about
10 times with the data I te
Hello all
> -Original Message-
> From: Travis Oliphant [mailto:[EMAIL PROTECTED]
> Sent: 03 October 2006 03:51
> To: Discussion of Numerical Python; Albert Strasheim
> Subject: Re: ***[Possible UCE]*** Re: [Numpy-discussion] Vectorizing code,
> for loops, and all that
&
Hello all
Some comments from a Windows user's perspective.
On Thu, 05 Oct 2006, Travis Oliphant wrote:
> John Hunter wrote:
>
> >>"Robert" == Robert Kern <[EMAIL PROTECTED]> writes:
> >
> >Robert> IMO, I'd rather see this and similar functions go into
> >Robert> scipy. New functions
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Karol Langner
> Sent: 07 October 2006 14:07
> To: NumPy List
> Subject: [Numpy-discussion] problem with lapack_lite / BLAS
>
> Dear list,
>
> I'm trying to compile ATLAS 3.6.0 and LA
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Tim Hochberg
> Sent: 12 October 2006 21:24
> To: numpy-discussion
> Subject: [Numpy-discussion] More SVN testing
>
>
> I just checked in a couple of changes to SVN. I was going to che
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Aric Hagberg
> Sent: Sunday, October 22, 2006 4:09 PM
> To: Discussion of Numerical Python
> Subject: Re: [Numpy-discussion] some work on arpack
>
> On Sat, Oct 21, 2006 at 02:05:42PM
Argh.
>
> In addition to the wrapper there is a Python interface (and some tests).
> I don't know if the interface is like "eigs" - I don't use Matlab.
http://www.mathworks.de/access/helpdesk/help/techdoc/ref/index.html?/access/
helpdesk/help/techdoc/ref/eigs.html
should give you some idea of
Hello all
I'm trying to sort an array with two fields, but I'm getting a result that
doesn't seem to make sense.
What I tried (first attempt): I have two 2-D arrays. I would like to sort
one based on the sort of the other. I managed to do this with argsort.
However, the fancy indexing required to
Hey Fernando
Maybe you can give the code a spin under Valgrind. It's going to be slow,
but if the crash is being caused by memory corruption that happens all the
time as the process is running, maybe Valgrind will show it.
You need some Valgrind suppressions for Python. It seems the 2.3 source tr
Hello all
I'm trying to generate random 32-bit integers. None of the following seem to
do the trick with NumPy 1.0.dev3383:
In [32]: N.random.randint(-2**31, 2**31-1)
ValueError: low >= high
In [43]: N.random.random_integers(-2**31, 2**31-1)
OverflowError: long int too large to convert to int
I
Hey Travis
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Travis Oliphant
> Sent: Tuesday, October 24, 2006 12:32 AM
> To: Discussion of Numerical Python
> Subject: [Numpy-discussion] Release of 1.0 coming
>
> The long awaited day is coming.
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Travis Oliphant
> Sent: Tuesday, October 24, 2006 12:04 AM
> To: Discussion of Numerical Python
> Subject: Re: [Numpy-discussion] Strange results when sorting array
> withfields
>
>
Hello all
I'm trying to get SciPy compiling with MSVC and the Intel Visual Fortran
compiler, version 9.1. I'm running into problems with the version
detection code failing.
The current version_cmd in intel.py is:
'version_cmd' : [fc_exe, "-FI -V -c %(fname)s.f -o %(fname)s.o" \
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Travis Oliphant
> Sent: Thursday, October 19, 2006 5:45 PM
> To: Discussion of Numerical Python
> Subject: Re: [Numpy-discussion] adding an attribute to an nd-array
>
> Stefan van der
Hello all
I managed to get NumPy to compile with the free for non-commercial use Linux
versions of the Intel C compiler, version 9.1 and Intel MKL 8.1.
Instructions at the bottom of this page:
http://www.scipy.org/Installing_SciPy/Linux
I ran the NumPy test suite, and it turned up a few test fai
Hey George
On Tue, 31 Oct 2006, George Sakkis wrote:
> Is there a more elegant and/or faster way to read some records from a
> file and then sort them by different fields ? What I have now is too
> specific and error-prone in general:
>
> import numpy as N
> records = N.fromfile(a_file, dtype=N.
Howdy
On Wed, 08 Nov 2006, Vincent Broman wrote:
> Building an rpm of numpy-1.0.1.dev3432-1 on fedora core 6 is failing for me.
> With either python-2.4.3 or 2.4.4 I try "python setup.py bdist_rpm"
> inside the source directory, and everything seems to go well except for
> many "File listed twice
Argh,
On Thu, 09 Nov 2006, Albert Strasheim wrote:
> %_unpackaged_files_terminate_build 1
Cut and paste error. Make that
%_unpackaged_files_terminate_build 0
Cheers,
Albert
-
Using Tomcat but need to do more? Need
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Travis Oliphant
> Sent: 08 June 2006 22:27
> To: numpy-discussion
> Subject: [Numpy-discussion] Array Protocol change for Python 2.6
>
> ...
>
> I would like to eliminate all the other
Hello all
For my Summer of Code project, I'm adding Support Vector Machine code to
SciPy. Underneath, I'm currently using libsvm. Thus far, I've been compiling
libsvm as a shared library (DLL on Windows) using SCons and doing the
wrapping with ctypes.
Now, I would like to integrate my code into t
Hey Glen
http://www.scipy.org/Cookbook/C_Extensions
covers most of the boilerplate you need to get started with extension
modules.
Regards,
Albert
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Glen W. Mabey
> Sent: 16 June 2006 18:24
> T
Hello all
On Fri, 16 Jun 2006, Travis Oliphant wrote:
> I just updated the array interface page to emphasize we now have version
> 3. NumPy still supports objects that expose (the C-side) of version 2
> of the array interface, though.
> Please voice concerns now if you have any.
In the d
Hey Sheldon
With NumPy you can use dtype's newbyteorder method to convert any dtype's
byte order to an order you specify:
In [1]: import numpy as N
In [2]: x = N.array([1],dtype='i4')
In [4]: xle = N.asarray(x, dtype=x.dtype.newbyteorder('<'))
In [5]: yle = N.asarray(y, dtype=y.dtype.newbyteorder
Hello all
Travis Oliphant wrote:
> Mathew Yeates wrote:
> > when I try and load a file with numpy.fromfile I keep getting a message
>
> > 7245092 items requested but only 3899 read. Its always the same number
> read.
> >
> >
> Which platform are you on? Could you show exactly how you are ca
Hey Chuck
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Charles R Harris
> Sent: 01 July 2006 19:57
> To: Robert Kern
> Cc: numpy-discussion@lists.sourceforge.net
> Subject: Re: [Numpy-discussion] Time for beta1 of NumPy 1.0
>
> All,
>
> T
Hello all
Travis Oliphant wrote:
> I've been playing a bit with ctypes and realized that with a little
> help, it could be made much easier to interface with NumPy arrays.
> Thus, I added a ctypes attribute to the NumPy array. If ctypes is
> installed, this attribute returns a "conversion" object
Hello all,
Travis Oliphant wrote:
> Hey Albert, I read the post you linked to on the ctypes mailing list.
> I hope I didn't step on any toes with what I did in NumPy. I was just
Certainly not. This is great stuff!
> working on a ctypes interface and realized that a lot of the cruft to
> conve
Hello all
Travis Oliphant wrote:
>
> Unfortunately, from the source code this is not true. It would be an
> improvement, but the source code shows that the from_param of each type
> does something special and only works with particular kinds of
> data-types --- basic Python types or ctypes types
Hello all
Travis Oliphant wrote:
> The ctypes-conversion object has attributes which return c_types aware
> objects so that the information can be passed directly to c-code (as an
> integer, the number of dimensions can already be passed using c-types).
>
> The information available and it's corr
Hello all
> > In this case, I think one wants strides[0] (the row stride) to return
> 40.
> >
>
> Why do you think that?
>
> All sliced arrays keep the same strides information as their
> "parents". This is the essence of a "view". The striding is exactly
> the same as before (the data hasn't
Hello all
On Mon, 03 Jul 2006, Travis Oliphant wrote:
> Fernando Perez wrote:
> > Hi all,
> >
> > but my lack of familiarity with all the details of new type creation got me
> > a
> > bit lost. I'm sure the information I need is all there, but right now I
> > don't
> > really see the forest
Hey Thomas
Thomas Heller wrote:
> Thomas Heller schrieb:
> > I've also played a little, and I think one important limitation in
> ctypes
> > is that items in the argtypes list have to be ctypes types.
>
> Thi swas misleading: I mean that this limitation should probably be
> removed, because it p
Hello all
On Tue, 04 Jul 2006, Thomas Heller wrote:
> Albert Strasheim schrieb:
> > Hey Thomas
> >
> > Thomas Heller wrote:
> >> Thomas Heller schrieb:
> >> > I've also played a little, and I think one important limitation in
> >> cty
Hello all
I noticed that arrays with dtypes with field names don't compare correctly
when using the equality operator.
Example:
In [17]: dt = N.dtype({'names' : ['f1', 'f2'],
'formats' : [N.intc, N.intc]})
In [18]: x = N.array([(0,0)], dtype=dt)
In [19]: y = N.array([(0,
Hello all
>
> Python 2.3b1+ (#2, Jun 10 2003, 20:53:51)
> [GCC 3.0.2 20010905 (Red Hat Linux 7.1 3.0.1-3)] on linux2
> Type "help", "copyright", "credits" or "license" for more information.
> >>> import Numeric
> >>> a=Numeric.arange(20)
> >>> b=a[::2]
> >>> c=a[1::2]
> >>> a
> array([ 0, 1
Hello all
> Please,
>
> try out the new C-API and let's get the bugs wrinkled out.
>
> Hopefully this will give us a more solid foundation for the future...
> I've already committed changes to matplotlib SVN that allow it to work
> with old and new NumPy.
What implications, if any, do these cha
Hello all
Various people wrote:
> > Im curious though: the several projects recently using ctypes
> > and numpy to wrap libraries (Pygame SDL, OpenGL, svm) must have come
> > across the issue of using a creating a numpy array from a ctypes
> > pointer. Ill have to look further.
> >
> It depends
Hey David
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of David Huard
> Sent: 17 July 2006 16:11
> To: numpy-discussion@lists.sourceforge.net
> Subject: Re: [Numpy-discussion] unique() should return a sorted array
>
> Hi,
>
> I attached a pa
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Keith Goodman
> Sent: 18 July 2006 15:55
> To: Thomas Heller
> Cc: numpy-discussion@lists.sourceforge.net
> Subject: Re: [Numpy-discussion] Cannot build numpy svn on Windows
>
> On 7/1
Hello all
In some situations, I have to work with very large matrices. My Windows
machine has 3 GB RAM, so I would expect to be able to use most of my
process's address space for my matrix.
Unfortunately, with matrices much larger than 700 or 800 MB, one starts
running into heap fragmentation pro
Hey Lars
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Lars Friedrich
> Sent: 25 July 2006 12:50
> To: numpy-discussion@lists.sourceforge.net
> Subject: [Numpy-discussion] ctypes, numpy-array
>
> Hello,
>
> I would like to work with some d
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Lars Friedrich
> Sent: 25 July 2006 13:55
> To: numpy-discussion@lists.sourceforge.net
> Subject: Re: [Numpy-discussion] ctypes, numpy-array
>
> > What's might be happening here is tha
Hey Mathew
The problem is that ATLAS doesn't provide all the LAPACK functions, only a
few that the ATLAS developers have optimized.
To get a complete LAPACK library, you need to build the Fortran LAPACK
library, and then put the ATLAS-optimized functions into this library.
Details here:
http://
Hello all
Another +1. When I build a Windows installer, I get:
numpy-1.0b2.dev2915.win32-py2.4.exe
This tells me everything I want to know.
Regards,
Albert
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Sasha
> Sent: 27 July 2006 21:26
>
Hello all
>
> Are you running on a 64 bit machine? I'll fix the test (which was
> written on a 32 bit machine), but I also noticed some strange
> behaviour:
>
> In [11]: N.intp('0xfff',16)
> Out[11]: 1152921504606846975
>
> In [12]: N.intp('0x',16)
> ---
Hello all
With NumPy 1.0 mere weeks away, I'm hoping we can improve the documentation
a bit before the final release. Some things we might want to think about:
1. Documentation Sprint
This page:
http://www.scipy.org/SciPy2006/CodingSprints
mentions a possible Documentation Sprint at SciPy 2006
Hello all
With the nice ctypes integration in NumPy, and with Python 2.5 which will
include ctypes around the corner, a remote possibility exists that within
the next year or two, I might not be the only person that wants to use NumPy
with ctypes.
This is probably going to mean that this someone
Hello all
Just a quick note on the ndpointer function that Travis recently added to
NumPy (thanks Travis!).
When wrapping functions with ctypes, one can specify the argument types of
the function. ctypes then checks that the parameters are valid before
invoking the C function.
This is described
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of David Grant
> Sent: 16 August 2006 17:11
> To: Discussion of Numerical Python
> Subject: Re: [Numpy-discussion] some work on arpack
>
>
>
> On 8/16/06, Keith Goodman <[EMAIL PROTECTE
Argh...
> I recently used MATLAB's eigs to do exactly this. I had a dense matrix A
> with dimensions m x n, where m >> n. I wanted the eigenvalues of A'A
> (which
> has dimensions m x m, which is too large to keep in memory). But I could
Make that AA'.
Cheers,
Albert
-
Hello all
>
> I decided to upgrade to 1.0b2 just to see what I get and now I get 7kB of
> "possibly lost" memory, coming from PyObject_Malloc (in
> /usr/lib/libpython2.4.so.1.0). This is a constant 7kB, however, and it
> isn't getting any larger if I increase the loop iterations. Looks good
> the
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of David Cournapeau
> Sent: 18 August 2006 06:55
> To: Discussion of Numerical Python
> Subject: [Numpy-discussion] ctypes: how does load_library work ?
>
>
> That works OK, but to avoid
Hello all
>
> Furthermore: the TRAC tickets do not seem to allow attachment of
> patches. Did I miss something there?
After submitting the initial report, you can attach files to the ticket.
Regards,
Albert
-
Using Tomca
Hello all
>
>
> >>> a = N.arange(4, dtype='>c8')
> >>> a.imag.max()
> 4.60060298822e-41
Confirmed on Windows 32-bit with 1.0b4.dev3050.
I created a ticket here:
http://projects.scipy.org/scipy/numpy/ticket/265
Regards,
Albert
---
A complete code snippet that reproduces the bug would be most helpful.
If there is a memory corruption problem, it might show up if we run the
problematic code under Valgrind.
Regards,
Albert
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of
Hello all
I just ran NumPy and SciPy through Valgrind, and everything looks clean on
that the NumPy side.
Some other things that could be fixed for RC1:
- GCC 4.1.1 warning in ufuncobject.c:
numpy/core/src/ufuncobject.c: In function âPyUFunc_RegisterLoopForTypeâ:
numpy/core/src/ufuncobject.c:32
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Sebastian Haase
> Sent: 15 September 2006 03:21
> To: numpy-discussion
> Subject: [Numpy-discussion] how to get info about internals of an
> arrayobject ?
>
> Hi,
> what I'm asking is if numpy ha
Hello all
In [1]: import numpy as N
In [3]: N.dtype({'names' : ['x', 'y'],
'formats' : [N.intc, N.float64]},
align=True)
Out[3]: dtype([('x', '
String Form:
Namespace: Interactive
File: c:\python24\lib\site-packages\numpy\__init__.py
Docstring:
Hello all
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:numpy-
> [EMAIL PROTECTED] On Behalf Of Martin Wiechert
> Sent: 15 September 2006 19:53
> To: Discussion of Numerical Python
> Subject: Re: [Numpy-discussion]PyArray_DescrConverter - alignment /
> trailingunused bytes
>
> On
Argh
> >
>
> > One final question. To me the repr of a dtype with gaps is a little bit
> > puzzling:
> >
> > >>> dtype ({'names': ['a', 'b', 'c'], 'formats': [' > 'offsets': [0, 16, 24]})
> > dtype([('a', '|S4'), ('', '|V12'), ('b', ' > ' >
> > There should be no gap between "b" and "c" but stil
66 matches
Mail list logo