numpy-discussion — Discussion list for all users of Numerical Python

Re: [Numpy-discussion] log can't handle big ints

[Travis O. <oli...@ie...>]

David Grant wrote:
> I am using numpy-0.9.8 and it seems that numpy's log2 function can't
> handle large integers?
>
> In [19]: a=11111111111111111111111111111111111111111111111111111111111111111111111111111111111111111
>
> In [20]: math.log(a,2)
> Out[20]: 292.48167544353294
>
> In [21]: numpy.log2(a)
>   

Ufuncs on objects (like the long object) work by looking for the 
corresponding method.  It's not found for long objects.

Convert the long object to a float first.   I'm not sure of any other 
way to "fix" it.  I suppose if no method is found an attempt to convert 
them to floats could be performed under the covers on all object array 
inputs.

-Travis

Re: [Numpy-discussion] Installation and Uninstallation

[Bryce H. <bhe...@en...>]

python setup.py install --prefix=your_path

You shouldn't have to clean up the previous install, if it got to the 
point where it was copy files, the first one would have failed. Next 
time you run setup.py with the --prefix option, it will pick up where 
the previous install left off.

Bryce

Elijah Gregory wrote:
> Dear NumPy Users,
>
>    I am attempting to install numpy-0.9.8 as a user on unix system. 
> When I install numpy by typing "python setup.py install" as per the 
> (only) instructions in the README.txt file everything proceeds 
> smoothly until some point where the script attempts to write a file to 
> the root-level /usr/lib64.  How can I configure the setup.py script to 
> use my user-level directories which I do have access to?  Also, given 
> that the install exited with an error, how do I clean up the aborted 
> installation?  Thank you for your help,
>
> regards,
>
>    Elijah Gregory
> ------------------------------------------------------------------------
>
> -------------------------------------------------------------------------
> Using Tomcat but need to do more? Need to support web services, security?
> Get stuff done quickly with pre-integrated technology to make your job easier
> Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
> http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
> ------------------------------------------------------------------------
>
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion
>   

[Numpy-discussion] Installation and Uninstallation

[Elijah G. <eli...@gm...>]

Dear NumPy Users,

   I am attempting to install numpy-0.9.8 as a user on unix system. When I
install numpy by typing "python setup.py install" as per the (only)
instructions in the README.txt file everything proceeds smoothly until some
point where the script attempts to write a file to the root-level
/usr/lib64.  How can I configure the setup.py script to use my user-level
directories which I do have access to?  Also, given that the install exited
with an error, how do I clean up the aborted installation?  Thank you for
your help,

regards,

   Elijah Gregory

[Numpy-discussion] log can't handle big ints

[David G. <dav...@gm...>]

I am using numpy-0.9.8 and it seems that numpy's log2 function can't
handle large integers?

In [19]: a=11111111111111111111111111111111111111111111111111111111111111111111111111111111111111111

In [20]: math.log(a,2)
Out[20]: 292.48167544353294

In [21]: numpy.log2(a)
---------------------------------------------------------------------------
exceptions.AttributeError                            Traceback (most
recent call last)

/home/david/<ipython console>

/usr/lib/python2.4/site-packages/numpy/lib/ufunclike.py in log2(x, y)
     52     x = asarray(x)
     53     if y is None:
---> 54         y = umath.log(x)
     55     else:
     56         umath.log(x, y)

AttributeError: 'long' object has no attribute 'log'

Does anyone else get this in numpy? if not, what version are you using?

-- 
David Grant
http://www.davidgrant.ca

Re: [Numpy-discussion] some work on arpack

[Nils W. <nw...@ia...>]

On Wed, 16 Aug 2006 08:08:16 -0700
  "David Grant" <dav...@gm...> wrote:
> On 8/16/06, Nils Wagner <nw...@ia...>=20
>wrote:
>>
>> Keith Goodman wrote:
>> > On 8/15/06, David Grant <dav...@gm...> wrote:
>> >
>> >
>> >> My idea is (if I have time) to write an eigs-like=20
>>function in python
>> >> that will only perform a subset of what Matlab's eigs=20
>>does for. It
>> >> will, for example, compute a certain number of=20
>>eigenvalues and
>> >> eigenvectors for a real, sparse, symmetric matrix=20
>>(the case I'm
>> >> interested in)
>> >>
>> >
>> >
>> AFAIK, pysparse (in the sandbox) includes a module that=20
>>implements a
>> Jacobi-Davidson
>> eigenvalue solver for the symmetric, generalised matrix=20
>>eigenvalue
>> problem (JDSYM).
>> Did someone test pysparse ?
>>
>> I did try pysparse a few years ago (I think right before=20
>>sparse stuff came
> into scipy). I think there is probably an old post=20
>asking the list about
> sparse stuff and I think Travis had just written it and=20
>told me about it...
> can't remember. Can JDSYM just return the k lowest=20
>eigenvalues/eigenvectors?
>=20
> --=20
> David Grant
> http://www.davidgrant.ca

Yes. See
  http://people.web.psi.ch/geus/pyfemax/pysparse_examples.html

for details.

Nils

Re: [Numpy-discussion] some work on arpack

[David G. <dav...@gm...>]

On 8/16/06, Nils Wagner <nw...@ia...> wrote:
>
> Keith Goodman wrote:
> > On 8/15/06, David Grant <dav...@gm...> wrote:
> >
> >
> >> My idea is (if I have time) to write an eigs-like function in python
> >> that will only perform a subset of what Matlab's eigs does for. It
> >> will, for example, compute a certain number of eigenvalues and
> >> eigenvectors for a real, sparse, symmetric matrix (the case I'm
> >> interested in)
> >>
> >
> >
> AFAIK, pysparse (in the sandbox) includes a module that implements a
> Jacobi-Davidson
> eigenvalue solver for the symmetric, generalised matrix eigenvalue
> problem (JDSYM).
> Did someone test pysparse ?
>
> I did try pysparse a few years ago (I think right before sparse stuff came
into scipy). I think there is probably an old post asking the list about
sparse stuff and I think Travis had just written it and told me about it...
can't remember. Can JDSYM just return the k lowest eigenvalues/eigenvectors?

-- 
David Grant
http://www.davidgrant.ca

Re: [Numpy-discussion] some work on arpack

[David G. <dav...@gm...>]

On 8/16/06, Albert Strasheim <fu...@gm...> wrote:
>
> Hello all
>
> > -----Original Message-----
> > From: num...@li... [mailto:numpy-
> > dis...@li...] 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 <kwg...@gm...> wrote:
> >
> >       On 8/15/06, David Grant <dav...@gm...> wrote:
> >
> >       > My idea is (if I have time) to write an eigs-like function in
> > python
> >       > that will only perform a subset of what Matlab's eigs does for.
> It
> >       > will, for example, compute a certain number of eigenvalues and
> >       > eigenvectors for a real, sparse, symmetric matrix (the case I'm
> >       > interested in)
> >
> >       Will it also work for a real, dense, symmetric matrix? That's the
> > case
> >       I'm interested in. But even if it doesn't, your work is great news
> > for
> >       numpy.
> >
> > Real, dense, symmetric, well doesn't scipy already have something for
> > this? I'm honestly not sure on the arpack side of things, I thought
> arpack
> > was only useful (over other tools) for sparse matrices, I could be
> wrong.
>
> Maybe SciPy can also do this, but what makes ARPACK useful is that it can
> get you a few eigenvalues and eigenvectors of a massive matrix without
> having to have the whole thing in memory. Instead, you provide ARPACK with
> a
> function that does A*x on your matrix. ARPACK passes a few x's to your
> function and a few eigenvalues and eigenvectors fall out.


Cool, thanks for the info.

-- 
David Grant
http://www.davidgrant.ca

Re: [Numpy-discussion] array equivalent to string.split(sep)

[Alan G I. <ai...@am...>]

On Wed, 16 Aug 2006, David Huard apparently wrote:=20
> I have a time series that I want to split into contiguous=20
> groups differentiated by a condition.

Perhaps itertools.groupby()?

fwiw,
Alan Isaac

Re: [Numpy-discussion] some work on arpack

[Albert S. <fu...@gm...>]

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

Re: [Numpy-discussion] some work on arpack

[Albert S. <fu...@gm...>]

Hello all

> -----Original Message-----
> From: num...@li... [mailto:numpy-
> dis...@li...] 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 <kwg...@gm...> wrote:
> 
> 	On 8/15/06, David Grant <dav...@gm...> wrote:
> 
> 	> My idea is (if I have time) to write an eigs-like function in
> python
> 	> that will only perform a subset of what Matlab's eigs does for. It
> 	> will, for example, compute a certain number of eigenvalues and
> 	> eigenvectors for a real, sparse, symmetric matrix (the case I'm
> 	> interested in)
> 
> 	Will it also work for a real, dense, symmetric matrix? That's the
> case
> 	I'm interested in. But even if it doesn't, your work is great news
> for
> 	numpy.
> 
> Real, dense, symmetric, well doesn't scipy already have something for
> this? I'm honestly not sure on the arpack side of things, I thought arpack
> was only useful (over other tools) for sparse matrices, I could be wrong.

Maybe SciPy can also do this, but what makes ARPACK useful is that it can
get you a few eigenvalues and eigenvectors of a massive matrix without
having to have the whole thing in memory. Instead, you provide ARPACK with a
function that does A*x on your matrix. ARPACK passes a few x's to your
function and a few eigenvalues and eigenvectors fall out.

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
keep A and A' in memory I could quickly calculate A'A*x, which is what
ARPACK needs.

Cheers,

Albert

Re: [Numpy-discussion] some work on arpack

[David G. <dav...@gm...>]

On 8/16/06, Keith Goodman <kwg...@gm...> wrote:
>
> On 8/15/06, David Grant <dav...@gm...> wrote:
>
> > My idea is (if I have time) to write an eigs-like function in python
> > that will only perform a subset of what Matlab's eigs does for. It
> > will, for example, compute a certain number of eigenvalues and
> > eigenvectors for a real, sparse, symmetric matrix (the case I'm
> > interested in)
>
> Will it also work for a real, dense, symmetric matrix? That's the case
> I'm interested in. But even if it doesn't, your work is great news for
> numpy.
>

Real, dense, symmetric, well doesn't scipy already have something for this?
I'm honestly not sure on the arpack side of things, I thought arpack was
only useful (over other tools) for sparse matrices, I could be wrong.

-- 
David Grant
http://www.davidgrant.ca

Re: [Numpy-discussion] some work on arpack

[Nils W. <nw...@ia...>]

Nils Wagner wrote:
> Keith Goodman wrote:
>  =20
>> On 8/15/06, David Grant <dav...@gm...> wrote:
>>
>>  =20
>>    =20
>>> My idea is (if I have time) to write an eigs-like function in python
>>> that will only perform a subset of what Matlab's eigs does for. It
>>> will, for example, compute a certain number of eigenvalues and
>>> eigenvectors for a real, sparse, symmetric matrix (the case I'm
>>> interested in)
>>>    =20
>>>      =20
>>  =20
>>    =20
> AFAIK, pysparse (in the sandbox) includes a module that implements a
> Jacobi-Davidson
> eigenvalue solver for the symmetric, generalised matrix eigenvalue
> problem (JDSYM).
> Did someone test pysparse ?
>
> Nils
>
>  =20
>> Will it also work for a real, dense, symmetric matrix? That's the case
>> I'm interested in. But even if it doesn't, your work is great news for
>> numpy.
>>
>> ----------------------------------------------------------------------=
---
>> Using Tomcat but need to do more? Need to support web services, securi=
ty?
>> Get stuff done quickly with pre-integrated technology to make your job=
 easier
>> Download IBM WebSphere Application Server v.1.0.1 based on Apache Gero=
nimo
>> http://sel.as-us.falkag.net/sel?cmd=3Dlnk&kid=3D120709&bid=3D263057&da=
t=3D121642
>> _______________________________________________
>> Numpy-discussion mailing list
>> Num...@li...
>> https://lists.sourceforge.net/lists/listinfo/numpy-discussion
>>  =20
>>    =20
>
>
> =20
>
>
> -----------------------------------------------------------------------=
--
> Using Tomcat but need to do more? Need to support web services, securit=
y?
> Get stuff done quickly with pre-integrated technology to make your job =
easier
> Download IBM WebSphere Application Server v.1.0.1 based on Apache Geron=
imo
> http://sel.as-us.falkag.net/sel?cmd=3Dlnk&kid=3D120709&bid=3D263057&dat=
=3D121642
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion
>  =20
Ok it's not ready...

gcc: Lib/sandbox/pysparse/src/spmatrixmodule.c
In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:17:
Lib/sandbox/pysparse/src/ll_mat.c: In function =91LLMat_matvec_transp=92:
Lib/sandbox/pysparse/src/ll_mat.c:760: error: =91CONTIGUOUS=92 undeclared
(first use in this function)
Lib/sandbox/pysparse/src/ll_mat.c:760: error: (Each undeclared
identifier is reported only once
Lib/sandbox/pysparse/src/ll_mat.c:760: error: for each function it
appears in.)
Lib/sandbox/pysparse/src/ll_mat.c: In function =91LLMat_matvec=92:
Lib/sandbox/pysparse/src/ll_mat.c:797: error: =91CONTIGUOUS=92 undeclared
(first use in this function)
In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:18:
Lib/sandbox/pysparse/src/csr_mat.c: In function =91CSRMat_matvec_transp=92=
:
Lib/sandbox/pysparse/src/csr_mat.c:119: error: =91CONTIGUOUS=92 undeclare=
d
(first use in this function)
Lib/sandbox/pysparse/src/csr_mat.c: In function =91CSRMat_matvec=92:
Lib/sandbox/pysparse/src/csr_mat.c:146: error: =91CONTIGUOUS=92 undeclare=
d
(first use in this function)
In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:19:
Lib/sandbox/pysparse/src/sss_mat.c: In function =91SSSMat_matvec=92:
Lib/sandbox/pysparse/src/sss_mat.c:83: error: =91CONTIGUOUS=92 undeclared
(first use in this function)
In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:17:
Lib/sandbox/pysparse/src/ll_mat.c: In function =91LLMat_matvec_transp=92:
Lib/sandbox/pysparse/src/ll_mat.c:760: error: =91CONTIGUOUS=92 undeclared
(first use in this function)
Lib/sandbox/pysparse/src/ll_mat.c:760: error: (Each undeclared
identifier is reported only once
Lib/sandbox/pysparse/src/ll_mat.c:760: error: for each function it
appears in.)
Lib/sandbox/pysparse/src/ll_mat.c: In function =91LLMat_matvec=92:
Lib/sandbox/pysparse/src/ll_mat.c:797: error: =91CONTIGUOUS=92 undeclared
(first use in this function)
In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:18:
Lib/sandbox/pysparse/src/csr_mat.c: In function =91CSRMat_matvec_transp=92=
:
Lib/sandbox/pysparse/src/csr_mat.c:119: error: =91CONTIGUOUS=92 undeclare=
d
(first use in this function)
Lib/sandbox/pysparse/src/csr_mat.c: In function =91CSRMat_matvec=92:
Lib/sandbox/pysparse/src/csr_mat.c:146: error: =91CONTIGUOUS=92 undeclare=
d
(first use in this function)
In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:19:
Lib/sandbox/pysparse/src/sss_mat.c: In function =91SSSMat_matvec=92:
Lib/sandbox/pysparse/src/sss_mat.c:83: error: =91CONTIGUOUS=92 undeclared
(first use in this function)
error: Command "gcc -pthread -fno-strict-aliasing -DNDEBUG -O2
-fmessage-length=3D0 -Wall -D_FORTIFY_SOURCE=3D2 -g -fPIC
-ILib/sandbox/pysparse/include/
-I/usr/lib64/python2.4/site-packages/numpy/core/include
-I/usr/include/python2.4 -c Lib/sandbox/pysparse/src/spmatrixmodule.c -o
build/temp.linux-x86_64-2.4/Lib/sandbox/pysparse/src/spmatrixmodule.o"
failed with exit status 1

Nils

[Numpy-discussion] array equivalent to string.split(sep)

[David H. <dav...@gm...>]

Hi,

I have a time series that I want to split into contiguous groups
differentiated by a condition. I didn't find a vectorized way to that, so I
ended up doing a for loop...

I know there are split functions that split arrays into equal lengths
subarrays, but is there a swell trick to return a sequence of arrays
separated by a condition ?

For instance, I would like to do something like:

>>> a = array([1,1,1,1,1,5,1,1,1,1,1,1,6,2,1,1])
>>> a.argsplit(a>1)
[[0,1,2,3,4], [6,7,8,9,10,11], [14,15]]

Thanks,

David

Re: [Numpy-discussion] some work on arpack

[Nils W. <nw...@ia...>]

Keith Goodman wrote:
> On 8/15/06, David Grant <dav...@gm...> wrote:
>
>   
>> My idea is (if I have time) to write an eigs-like function in python
>> that will only perform a subset of what Matlab's eigs does for. It
>> will, for example, compute a certain number of eigenvalues and
>> eigenvectors for a real, sparse, symmetric matrix (the case I'm
>> interested in)
>>     
>
>   
AFAIK, pysparse (in the sandbox) includes a module that implements a
Jacobi-Davidson
eigenvalue solver for the symmetric, generalised matrix eigenvalue
problem (JDSYM).
Did someone test pysparse ?

Nils

> Will it also work for a real, dense, symmetric matrix? That's the case
> I'm interested in. But even if it doesn't, your work is great news for
> numpy.
>
> -------------------------------------------------------------------------
> Using Tomcat but need to do more? Need to support web services, security?
> Get stuff done quickly with pre-integrated technology to make your job easier
> Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
> http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion
>   


 

Re: [Numpy-discussion] some work on arpack

[Keith G. <kwg...@gm...>]

On 8/15/06, David Grant <dav...@gm...> wrote:

> My idea is (if I have time) to write an eigs-like function in python
> that will only perform a subset of what Matlab's eigs does for. It
> will, for example, compute a certain number of eigenvalues and
> eigenvectors for a real, sparse, symmetric matrix (the case I'm
> interested in)

Will it also work for a real, dense, symmetric matrix? That's the case
I'm interested in. But even if it doesn't, your work is great news for
numpy.

Re: [Numpy-discussion] Numpy-discussion Digest, Vol 3, Issue 42

[Matthew T. <mm...@cs...>]

For starters, wow.  I'm sorry.  I didn't mean to spam my problem 5  
times.  My mail server decided to fritz out today and I thought it  
was Sourceforge rejecting my emails since they didn't originate the  
address I'm registered as.  My apologies.

> Do you need to specify --prefix if you've already got Python installed
> somewhere?
>
> Are you missing it.

I tried it again without setting it.  No more luck.

> This  seems to be a path issue.  Can you give us
>
> import sys
> print sys.path()

[
'',
'/home/m/mmt/local/lib/python2.4/site-packages',
'/home/m/mmt/local/lib/python2.4/site-packages/PIL',
'/home/m/mmt/local/lib/python2.4/site-packages/numpy',
'/grads/mmt/local/lib/python24.zip',
'/grads/mmt/local/lib/python2.4',
'/grads/mmt/local/lib/python2.4/plat-linux2',
'/grads/mmt/local/lib/python2.4/lib-tk',
'/grads/mmt/local/lib/python2.4/lib-dynload',
'/grads/mmt/local/lib/python2.4/site-packages',
'/grads/mmt/local/lib/python2.4/site-packages/PIL'
]

The top 3 are added to my python path by myself, the rest are  
included by default.  FYI: /grads/mmt and /home/m/mmt map to the same  
diretory.

Sorry again about the repeat emails.

Matt

Re: [Numpy-discussion] numpy 1.0b2 problems

[Travis O. <oli...@ie...>]

Matthew Trentacoste wrote:
> Hey.  I'm trying to get numpy up and running on SuSE 10.1 and not  
> having much luck.
>
> I've been working with 1.0b2 and can get it to install without any  
> errors, but can't do anything with it.  I run a local install of  
> python 2.4.3 just to keep out of whatever weirdness gets installed on  
> my machine by our sysadmins.  Pretty standard fare, untar the ball,  
> and './setup.py install --prefix=$HOME/local'
>
> It will complete that without issue, but when I try to run the test,  
> I get:
>
> Python 2.4.3 (#1, Aug 15 2006, 18:09:56)
> [GCC 4.1.0 (SUSE Linux)] on linux2
> Type "help", "copyright", "credits" or "license" for more information.
>  >>> import numpy
>  >>> numpy.test(1)
> Traceback (most recent call last):
>    File "<stdin>", line 1, in ?
>    File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
> __init__.py", line 77, in test
>      return NumpyTest().test(level, verbosity)
>    File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/testing/ 
> numpytest.py", line 285, in __init__
>      from numpy.distutils.misc_util import get_frame
>    File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
> distutils/__init__.py", line 5, in ?
>      import ccompiler
>    File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
> distutils/ccompiler.py", line 6, in ?
>      from distutils.ccompiler import *
>    File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
> distutils/__init__.py", line 5, in ?
>      import ccompiler
>    File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
> distutils/ccompiler.py", line 7, in ?
>      from distutils import ccompiler
> ImportError: cannot import name ccompiler
>
>   

This  seems to be a path issue.  Can you give us

import sys
print sys.path()

-Travis

Re: [Numpy-discussion] numpy 1.0b2 problems

[Travis O. <oli...@ie...>]

Matthew Trentacoste wrote:
> Hey.  I'm trying to get numpy up and running on SuSE 10.1 and not  
> having much luck.
>
> I've been working with 1.0b2 and can get it to install without any  
> errors, but can't do anything with it.  I run a local install of  
> python 2.4.3 just to keep out of whatever weirdness gets installed on  
> my machine by our sysadmins.  Pretty standard fare, untar the ball,  
> and './setup.py install --prefix=$HOME/local'
>   
Do you need to specify --prefix if you've already got Python installed 
somewhere?

Are you missing it.

[Numpy-discussion] Fwd: numpy 1.0b2 problems

[Matthew T. <mm...@cs...>]

Hey.  I'm trying to get numpy up and running on SuSE 10.1 and not  
having much luck.

I've been working with 1.0b2 and can get it to install without any  
errors, but can't do anything with it.  I run a local install of  
python 2.4.3 just to keep out of whatever weirdness gets installed on  
my machine by our sysadmins.  Pretty standard fare, untar the ball,  
and './setup.py install --prefix=$HOME/local'

It will complete that without issue, but when I try to run the test,  
I get:

Python 2.4.3 (#1, Aug 15 2006, 18:09:56)
[GCC 4.1.0 (SUSE Linux)] on linux2
Type "help", "copyright", "credits" or "license" for more information.
 >>> import numpy
 >>> numpy.test(1)
Traceback (most recent call last):
   File "<stdin>", line 1, in ?
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
__init__.py", line 77, in test
     return NumpyTest().test(level, verbosity)
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/testing/ 
numpytest.py", line 285, in __init__
     from numpy.distutils.misc_util import get_frame
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/__init__.py", line 5, in ?
     import ccompiler
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/ccompiler.py", line 6, in ?
     from distutils.ccompiler import *
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/__init__.py", line 5, in ?
     import ccompiler
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/ccompiler.py", line 7, in ?
     from distutils import ccompiler
ImportError: cannot import name ccompiler

Once installed, it messes up trying to install anything else, so I  
have to move it out of the way in the short term.

Any thoughts?

Thanks
Matt

[  matthew m trentacoste                 mm...@cs...                 ]
[                                                                      ]
[  graduate student                      lead software developer       ]
[  university of british columbia        brightside technologies       ]
[  http://www.cs.ubc.ca/~mmt             http://brightsidetech.com     ]
[  +1 (604) 827-3979                     +1 (604) 228-4624             ]

[Numpy-discussion] numpy 1.0b2 problems

[Matthew T. <mm...@cs...>]

Hey.  I'm trying to get numpy up and running on SuSE 10.1 and not  
having much luck.

I've been working with 1.0b2 and can get it to install without any  
errors, but can't do anything with it.  I run a local install of  
python 2.4.3 just to keep out of whatever weirdness gets installed on  
my machine by our sysadmins.  Pretty standard fare, untar the ball,  
and './setup.py install --prefix=$HOME/local'

It will complete that without issue, but when I try to run the test,  
I get:

Python 2.4.3 (#1, Aug 15 2006, 18:09:56)
[GCC 4.1.0 (SUSE Linux)] on linux2
Type "help", "copyright", "credits" or "license" for more information.
 >>> import numpy
 >>> numpy.test(1)
Traceback (most recent call last):
   File "<stdin>", line 1, in ?
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
__init__.py", line 77, in test
     return NumpyTest().test(level, verbosity)
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/testing/ 
numpytest.py", line 285, in __init__
     from numpy.distutils.misc_util import get_frame
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/__init__.py", line 5, in ?
     import ccompiler
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/ccompiler.py", line 6, in ?
     from distutils.ccompiler import *
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/__init__.py", line 5, in ?
     import ccompiler
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/ccompiler.py", line 7, in ?
     from distutils import ccompiler
ImportError: cannot import name ccompiler

Once installed, it messes up trying to install anything else, so I  
have to move it out of the way in the short term.

Any thoughts?

Thanks
Matt

[  matthew m trentacoste                 mm...@cs...                 ]
[                                                                      ]
[  graduate student                      lead software developer       ]
[  university of british columbia        brightside technologies       ]
[  http://www.cs.ubc.ca/~mmt             http://brightsidetech.com     ]
[  +1 (604) 827-3979                     +1 (604) 228-4624             ]

[Numpy-discussion] Numpy 1.0b2 install issues

[Matthew T. <mm...@cs...>]

Hey.  I'm trying to get numpy up and running on SuSE 10.1 and not  
having much luck.

I've been working with 1.0b2 and can get it to install without any  
errors, but can't do anything with it.  I run a local install of  
python 2.4.3 just to keep out of whatever weirdness gets installed on  
my machine by our sysadmins.  Pretty standard fare, untar the ball,  
and './setup.py install --prefix=$HOME/local'

It will complete that without issue, but when I try to run the test,  
I get:

Python 2.4.3 (#1, Aug 15 2006, 18:09:56)
[GCC 4.1.0 (SUSE Linux)] on linux2
Type "help", "copyright", "credits" or "license" for more information.
 >>> import numpy
 >>> numpy.test(1)
Traceback (most recent call last):
   File "<stdin>", line 1, in ?
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
__init__.py", line 77, in test
     return NumpyTest().test(level, verbosity)
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/testing/ 
numpytest.py", line 285, in __init__
     from numpy.distutils.misc_util import get_frame
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/__init__.py", line 5, in ?
     import ccompiler
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/ccompiler.py", line 6, in ?
     from distutils.ccompiler import *
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/__init__.py", line 5, in ?
     import ccompiler
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/ccompiler.py", line 7, in ?
     from distutils import ccompiler
ImportError: cannot import name ccompiler

This pretty much borks everything.  I have to remove it before I can  
try to install other packages and stuff.

Any thoughts?

Thanks
Matt

[Numpy-discussion] some work on arpack

[David G. <dav...@gm...>]

Building an arpack extension turned out to be surprisingly simple. For
example for dsaupd:

f2py -c dsaupd.f -m dsaupd -L/usr/lib/blas/atlas:/usr/lib/lapack/atlas
-llapack -lblas -larpack

It took me a long time to get the command down to something that
simple. Took me a while even to figure out I could just use the arpack
library on my computer rather than re-linking all of arpack! I was
able to import the dsaupd.so python module just fine and I was also
able to call it just fine. I'll have to tweak the pyf file in order to
get some proper output. But this gives me confidence that arpack is
easy to hook into which is what others have said in the past, but
without any experience with f2py I had no idea myself. f2py is
awesome, for anyone who doesn't know.

Matlab has interfaces for the arpack functions like dsaupd, dseupd,
dnaupd, znaupd, zneupd (the mex file documentation claims those are
the only ones, but they have more). Matlab has a C interface to these
functions in arpackc.mex* and the script eigs.m does the grunt work,
providing a very high-level interface as well as doing some linear
algebra (the same type of stuff that is done in arpack's examples
directory I gather) and various other things.

My idea is (if I have time) to write an eigs-like function in python
that will only perform a subset of what Matlab's eigs does for. It
will, for example, compute a certain number of eigenvalues and
eigenvectors for a real, sparse, symmetric matrix (the case I'm
interested in)... I hope that this subset-of-matlab's-eigs function
will not be too hard to write. Then more functionality can be added on
to eigs.py later... Does this make sense?

Has anyone else started work on arpack integration at all?

-- 
David Grant
http://www.davidgrant.ca

[Numpy-discussion] numpy 1.0b2 problems

[Matthew T. <mm...@cs...>]

Hey.  I'm trying to get numpy up and running on SuSE 10.1 and not  
having much luck.

I've been working with 1.0b2 and can get it to install without any  
errors, but can't do anything with it.  I run a local install of  
python 2.4.3 just to keep out of whatever weirdness gets installed on  
my machine by our sysadmins.  Pretty standard fare, untar the ball,  
and './setup.py install --prefix=$HOME/local'

It will complete that without issue, but when I try to run the test,  
I get:

Python 2.4.3 (#1, Aug 15 2006, 18:09:56)
[GCC 4.1.0 (SUSE Linux)] on linux2
Type "help", "copyright", "credits" or "license" for more information.
 >>> import numpy
 >>> numpy.test(1)
Traceback (most recent call last):
   File "<stdin>", line 1, in ?
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
__init__.py", line 77, in test
     return NumpyTest().test(level, verbosity)
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/testing/ 
numpytest.py", line 285, in __init__
     from numpy.distutils.misc_util import get_frame
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/__init__.py", line 5, in ?
     import ccompiler
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/ccompiler.py", line 6, in ?
     from distutils.ccompiler import *
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/__init__.py", line 5, in ?
     import ccompiler
   File "/home/m/mmt/local/lib/python2.4/site-packages/numpy/ 
distutils/ccompiler.py", line 7, in ?
     from distutils import ccompiler
ImportError: cannot import name ccompiler

Any thoughts?

Thanks
Matt

[  matthew m trentacoste                 mm...@cs...                 ]
[                                                                      ]
[  graduate student                      lead software developer       ]
[  university of british columbia        brightside technologies       ]
[  http://www.cs.ubc.ca/~mmt             http://brightsidetech.com     ]
[  +1 (604) 827-3979                     +1 (604) 228-4624             ]

Re: [Numpy-discussion] f2py --include_paths from command line

[Stephen W. <drs...@gm...>]

On 8/15/06, Ethan T Coon <et...@co...> wrote:
>
> Hi all,
>
> The following line:
>
> f2py -c -m _test --include_paths ./include test.f


Typing f2py alone seems to indicate the syntax should be

f2py -I./include [other args] test.f

I tried this and it seems to work here.

[Numpy-discussion] request for new array method: arr.abs()

[Sebastian H. <ha...@ms...>]

Hi!
numpy renamed the *function* abs to absolute.
Most functions like mean, min, max, average, ...
have an equivalent array *method*.

Why is absolute left out ?
I think it should be added .

Furthermore, looking at some line  of code that have multiple calls to 
absolute [ like f(absolute(a), absolute(b), absolute(c)) ]
I think "some people" might prefer less typing and less reading,
like f( a.abs(), b.abs(), c.abs() ).
One could even consider not requiring the "function call" parenthesis '()'  at 
all - but I don't know about further implications that might have.

Thanks,
Sebastian Haase

PS: is there any performace hit in using the built-in abs function ?