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

[Numpy-discussion] Setuptools leftover junk

[Fernando P. <fpe...@gm...>]

Hi all,

I recently noticed one of my in-house projects started leaving aroun
.egg-info directories after I ran its setup.py, even though I don't
use setuptools for anything at all.  For now I just added an extra
clean rule to my makefile and forgot about it, but it kind of annoyed
me.

Today I looked at the temp directory where I've been making my
numpy/scipy installs from SVN, and here's what I saw:

longs[site-packages]> d
/home/fperez/tmp/local/lib/python2.4/site-packages
total 228
drwxr-xr-x  2 fperez   4096 2006-06-21 22:16 dateutil/
drwxr-xr-x  7 fperez   4096 2006-06-28 02:50 matplotlib/
drwxr-xr-x 13 fperez   4096 2006-06-28 02:38 numpy/
drwxr-xr-x  2 fperez   4096 2006-06-21 21:28 numpy-0.9.9.2660-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-22 21:29 numpy-0.9.9.2665-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-24 11:33 numpy-0.9.9.2674-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-24 15:08 numpy-0.9.9.2675-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-25 12:40 numpy-0.9.9.2677-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-26 23:32 numpy-0.9.9.2691-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-28 02:38 numpy-0.9.9.2696-py2.4.egg-info/
-rw-r--r--  1 fperez     31 2006-03-18 20:11 pylab.py
-rw-r--r--  1 fperez    178 2006-06-24 13:29 pylab.pyc
drwxr-xr-x 20 fperez   4096 2006-06-28 11:20 scipy/
drwxr-xr-x  2 fperez   4096 2006-06-21 21:36 scipy-0.5.0.1990-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-22 21:36 scipy-0.5.0.1998-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-24 15:15 scipy-0.5.0.1999-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-25 12:46 scipy-0.5.0.2000-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-26 23:37 scipy-0.5.0.2004-py2.4.egg-info/
drwxr-xr-x  2 fperez   4096 2006-06-28 02:48 scipy-0.5.0.2012-py2.4.egg-info/


Is it really necessary to have all that setuptools junk left around,
for those of us who aren't asking for it explicitly?  My personal
opinions on setuptools aside, I think it's just a sane practice not to
create this kind of extra baggage unless explicitly requested.

I scoured my home directory for any .file which might be triggering
this inadvertedly, but I can't seem to find any, so I'm going to guess
this is somehow being caused by numpy's own setup.  If it's my own
mistake, I'll be happy to be shown how to coexist peacefully with
setuptools.

Since this also affects user code (I think via f2py or something
internal to numpy, since all I'm calling is f2py in my code), I really
think it would be nice to clean it.

Opinions?

f

Re: [Numpy-discussion] indexing bug in numpy r2694

[Pau G. <pau...@gm...>]

i don't know why 'where' is returning matrices.
if you use:

>>> idx = where(y.A > 0.5)[0]

everything will work fine (I guess)

pau


On 6/28/06, Keith Goodman <kwg...@gm...> wrote:
> >> x = asmatrix(rand(3,2))
>
> >> y = asmatrix(rand(3,1))
>
> >> y
>
> matrix([[ 0.77952062],
>        [ 0.97110465],
>        [ 0.77450218]])
>
> >> idx = where(y > 0.5)[0]
>
> >> idx
> matrix([[0, 1, 2]])
>
> >> x[idx,:]
>
> matrix([[ 0.24837887,  0.52988253],
>        [ 0.28661085,  0.43053076],
>        [ 0.05360893,  0.22668509]])
>
> So far everything works as it should. Now the problem:
>
> >> y[idx,:]
> ---------------------------------------------------------------------------
> exceptions.ValueError                                Traceback (most
> recent call last)
>
> /usr/local/lib/python2.4/site-packages/numpy/core/defmatrix.py in
> __getitem__(self, index)
>     120
>     121     def __getitem__(self, index):
> --> 122         out = N.ndarray.__getitem__(self, index)
>     123         # Need to swap if slice is on first index
>     124         retscal = False
>
> /usr/local/lib/python2.4/site-packages/numpy/core/defmatrix.py in
> __array_finalize__(self, obj)
>     116             self.shape = (1,1)
>     117         elif ndim == 1:
> --> 118             self.shape = (1,self.shape[0])
>     119         return
>     120
>
> ValueError: total size of new array must be unchanged
>
>
> And, on a related note, shouldn't this be a column vector?
>
> >> x[idx,0]
> matrix([[ 0.24837887,  0.28661085,  0.05360893]])
>
> 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] indexing bug in numpy r2694

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

>> x = asmatrix(rand(3,2))

>> y = asmatrix(rand(3,1))

>> y

matrix([[ 0.77952062],
       [ 0.97110465],
       [ 0.77450218]])

>> idx = where(y > 0.5)[0]

>> idx
matrix([[0, 1, 2]])

>> x[idx,:]

matrix([[ 0.24837887,  0.52988253],
       [ 0.28661085,  0.43053076],
       [ 0.05360893,  0.22668509]])

So far everything works as it should. Now the problem:

>> y[idx,:]
---------------------------------------------------------------------------
exceptions.ValueError                                Traceback (most
recent call last)

/usr/local/lib/python2.4/site-packages/numpy/core/defmatrix.py in
__getitem__(self, index)
    120
    121     def __getitem__(self, index):
--> 122         out = N.ndarray.__getitem__(self, index)
    123         # Need to swap if slice is on first index
    124         retscal = False

/usr/local/lib/python2.4/site-packages/numpy/core/defmatrix.py in
__array_finalize__(self, obj)
    116             self.shape = (1,1)
    117         elif ndim == 1:
--> 118             self.shape = (1,self.shape[0])
    119         return
    120

ValueError: total size of new array must be unchanged


And, on a related note, shouldn't this be a column vector?

>> x[idx,0]
matrix([[ 0.24837887,  0.28661085,  0.05360893]])

Re: [Numpy-discussion] fread codes versus numpy types

[Fernando P. <fpe...@gm...>]

On 6/28/06, Robert Kern <rob...@gm...> wrote:

> The Capitalized versions are actually old typecodes for backwards compatibility
> with Numeric. In recent development versions of numpy, they are no longer
> exposed except through the numpy.oldnumeric compatibility module. A decision was
> made for numpy to use the actual width of a type in its name instead of the
> width of its component parts (when it has parts).
>
> Code in scipy which still requires actual string typecodes is a bug. Please
> report such cases on the Trac:
>
>    http://projects.scipy.org/scipy/scipy

Well, an easy way to make all those poke their ugly heads in a hurry
would be to remove line 32 in scipy's init:

longs[Lib]> grep -n oldnum *py
__init__.py:31:import numpy.oldnumeric as _num
__init__.py:32:from numpy.oldnumeric import *


If we really want to push for the new api, I think it's fair to change
those two lines by simply

from numpy import oldnumeric

so that scipy also exposes oldnumeric, and let all deprecated names be
hidden there.

I just tried this change:

Index: __init__.py
===================================================================
--- __init__.py (revision 2012)
+++ __init__.py (working copy)
@@ -29,9 +29,8 @@

 # Import numpy symbols to scipy name space
 import numpy.oldnumeric as _num
-from numpy.oldnumeric import *
-del lib
-del linalg
+from numpy import oldnumeric
+
 __all__ += _num.__all__
 __doc__ += """
 Contents


and scipy's test suite still passes (modulo the test_cobyla thingie
Nils is currently fixing, which is not related to this).

Should I apply this patch, so we push the cleaned-up API even a bit harder?

Cheers,

f

Re: [Numpy-discussion] what happened to numarray type names ?

[Christopher B. <Chr...@no...>]

Erin Sheldon wrote:
> OK, I have changed all the examples that used dtype=Float or
> dtype=Int to float and int.

They are also available as:

numpy.float_
numpy.int_

-Chris


-- 
Christopher Barker, Ph.D.
Oceanographer
                                     		
NOAA/OR&R/HAZMAT         (206) 526-6959   voice
7600 Sand Point Way NE   (206) 526-6329   fax
Seattle, WA  98115       (206) 526-6317   main reception

Chr...@no...

Re: [Numpy-discussion] fread codes versus numpy types

[Robert K. <rob...@gm...>]

Glen W. Mabey wrote:
> Hello,
> 
> I see the following character codes defined in scipy (presumably) for
> use with scipy.io.fread() :
> 
> In [20]:scipy.Complex
> Out[20]:'D'
> 
> In [21]:scipy.Complex0
> Out[21]:'D'
> 
> In [22]:scipy.Complex128
> Out[22]:'G'
> 
> In [23]:scipy.Complex16
> Out[23]:'F'
> 
> In [24]:scipy.Complex32
> Out[24]:'F'
> 
> In [25]:scipy.Complex64
> Out[25]:'D'
> 
> In [26]:scipy.Complex8
> Out[26]:'F'
> 
> Then I see the following scalar types also defined:
> 
> In [27]:scipy.complex64
> Out[27]:<type 'complex64scalar'>
> 
> In [28]:scipy.complex128
> Out[28]:<type 'complex128scalar'>
> 
> In [29]:scipy.complex256
> Out[29]:<type 'complex256scalar'>
> 
> which correspond to types that exist within the numpy module.  These
> names seem to conflict in that (unless I misunderstand what's going on)
> scipy.complex64 actually occupies 64 bits of data (a 32-bit float for
> each of {real, imag}) whereas scipy.Complex64 looks like it occupies 128
> bits of data (a 64-bit double for each of {real, imag}).
> 
> Is there something I'm missing, or is this a naming inconsistency?

The Capitalized versions are actually old typecodes for backwards compatibility 
with Numeric. In recent development versions of numpy, they are no longer 
exposed except through the numpy.oldnumeric compatibility module. A decision was 
made for numpy to use the actual width of a type in its name instead of the 
width of its component parts (when it has parts).

Code in scipy which still requires actual string typecodes is a bug. Please 
report such cases on the Trac:

   http://projects.scipy.org/scipy/scipy

-- 
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless enigma
  that is made terrible by our own mad attempt to interpret it as though it had
  an underlying truth."
   -- Umberto Eco

[Numpy-discussion] matlab -> python translation

[Stefan v. d. W. <st...@su...>]

Hi all,

I recently saw discussions on the list regarding Matlab/Octave to
Python translation.  I brought this under John Eaton's attention (he
is the original author of Octave) -- below is his response.

Regards
St=E9fan

----- Forwarded message from "John W. Eaton" <jw...@be...> ----=
-

From: "John W. Eaton" <jw...@be...>

On 21-Jun-2006, Stefan van der Walt wrote:

| I'd like to bring this thread under your attention, in case you want
| to comment:
|=20
| http://aspn.activestate.com/ASPN/Mail/Message/numpy-discussion/3174978

Would you please pass along the following comments?


Translating the syntax might not be too hard, but to have a really
effective tool, you have to get all the details of the Matlab/Octave
function calls the same as well.  So would you do that by linking to
Octave's run-time libraries as well?  That could pobably be made to
work, but it would probably drag in a lot more code that some poeple
would expect when they just want to translate and run a relatively
small number of lines of Matlab code.

Another semantic detail that would likely cause trouble is the
(apparent) pass-by-value semantics of Matlab.  How would you reconcile
this with the mutable types of Python?

Finally, I would encourage anyone who wants to work on a Matlab/Octave
to Python translator using Octave's parser and run-time libraries to
work on this in a way that can be integrated with Octave.  Please
consider discuss your ideas about this project on the
mai...@oc... mailing list.


Thanks,

jwe

----- End forwarded message -----

[Numpy-discussion] fread codes versus numpy types

[Glen W. M. <Gle...@sw...>]

Hello,

I see the following character codes defined in scipy (presumably) for
use with scipy.io.fread() :


In [20]:scipy.Complex
Out[20]:'D'

In [21]:scipy.Complex0
Out[21]:'D'

In [22]:scipy.Complex128
Out[22]:'G'

In [23]:scipy.Complex16
Out[23]:'F'

In [24]:scipy.Complex32
Out[24]:'F'

In [25]:scipy.Complex64
Out[25]:'D'

In [26]:scipy.Complex8
Out[26]:'F'


Then I see the following scalar types also defined:


In [27]:scipy.complex64
Out[27]:<type 'complex64scalar'>

In [28]:scipy.complex128
Out[28]:<type 'complex128scalar'>

In [29]:scipy.complex256
Out[29]:<type 'complex256scalar'>


which correspond to types that exist within the numpy module.  These
names seem to conflict in that (unless I misunderstand what's going on)
scipy.complex64 actually occupies 64 bits of data (a 32-bit float for
each of {real, imag}) whereas scipy.Complex64 looks like it occupies 128
bits of data (a 64-bit double for each of {real, imag}).

Is there something I'm missing, or is this a naming inconsistency?

Glen

Re: [Numpy-discussion] int64 wierdness

[Vincent S. <sc...@sa...>]

Andrew Straw wrote:
> An SVN checkout from a week or two ago looks OK on my amd64 machine:
>
> astraw@hdmg:~$ python
> Python 2.4.3 (#2, Apr 27 2006, 14:43:32)
> [GCC 4.0.3 (Ubuntu 4.0.3-1ubuntu5)] on linux2
> Type "help", "copyright", "credits" or "license" for more information.
>  >>> import numpy
>  >>> numpy.__version__
> '0.9.9.2631'
>  >>> numpy.int64(9)**2
> 81
>  >>>
>   
Confirmed to be fixed on my gentoo amd64 machine, numpy svn of couple of
days ago:
>>> numpy.int64(9)**2
81
>>> numpy.__version__
'0.9.9.2665'

Cheers,
Vincent.
>
> EI wrote:
>
>   
>> numpy.__version__ says 0.9.8.
>>
>> Python 2.4.2, GCC 4.1, OpenSuSE 10.1 (x86_64).
>>
>> Thanks Travis,
>> Eugene
>>
>> On 6/27/06, *Travis Oliphant* < oli...@ie... 
>> <mailto:oli...@ie...>> wrote:
>>
>>     EI wrote:
>>     > Hi,
>>     >
>>     > I'm running python 2.4 on a 64bit linux and get strange results:
>>     > (int(9))**2 is equal to 81, as it should, but
>>     > (int64(9))**2 is equal to 0
>>
>>     Thanks for the bug-report.  Please provide the version of NumPy
>>     you are
>>     using so we can track it down, or suggest an upgrade.
>>
>>     -Travis
>>
>>
>> ------------------------------------------------------------------------
>>
>> 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
>>  
>>
>>     
>
>
> 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] Numpy Benchmarking

[Filip W. <fi...@ft...>]

Jens wrote:

> Dennis V. Perepelitsa wrote:

>>Hi, all.
>>
>>I've run some benchmarks comparing the performance of scipy, numpy, 
>>Numeric and numarray vs. MATLAB. There's also the beginnings of a 
>>benchmark framework included. The results are online at:
>>
>>     http://web.mit.edu/jonas/www/bench/
>>  
>>

> It's a little hard to see the curves for small matrix size, N.  How 
> about plotting the time divided by the theoretical number of operations
> - which would be N^2 or N^3.

Or use some logarithmic scale (one or both axis) where applicable.

fw

Re: [Numpy-discussion] Numpy Benchmarking

[<je...@fy...>]

Dennis V. Perepelitsa wrote:

>Hi, all.
>
>I've run some benchmarks comparing the performance of scipy, numpy,=20
>Numeric and numarray vs. MATLAB. There's also the beginnings of a=20
>benchmark framework included. The results are online at:
>
>     http://web.mit.edu/jonas/www/bench/
> =20
>

It's a little hard to see the curves for small matrix size, N.  How=20
about plotting the time divided by the theoretical number of operations=20
- which would be N^2 or N^3.

Jens J=F8rgen

Re: [Numpy-discussion] Numpy Benchmarking

[Arnd B. <arn...@we...>]

Hi,

On Wed, 28 Jun 2006, Jon Wright wrote:

>
> >>This strikes me as a little bit odd. Why not just provide the best-performing
> >>function to both SciPy and NumPy? Would NumPy be more difficult to install
> >>if the SciPy algorithm for inv() was incorporated?
> >>
> >>
> Having spent a few days recently trying out various different
> eigenvector routines in Lapack I would have greatly appreciated having a
> choice of which one to use

which routine are you trying to use?

> from without having to create my own
> wrappers, compiling atlas and lapack under windows (ouch). I noted that
> Numeric (24.2) seemed to be converting Float32 to double meaning my
> problem no longer fits in memory, which was the motivation for the work.
> Poking around in the svn of numpy.linalg appears to find the same lapack
> routine as Numeric (dsyevd). Perhaps I miss something in the code logic?

if you can convince the code to  get ssyevd instead of dsyevd
it might do what you want>


> The divide and conquer (*evd) uses more memory than the (*ev), as well
> as a factor of 2 for float/double, hence my problem, and the reason why
> "best performing" is a hard choice. I thought matlab has a look at the
> matrix dimensions and problem before deciding what to do (eg: the \
> operator).

Hmm, this is a hard choice, which might better left
in the hands of the knowledgeable user.
(e.g., aren't the divide and conquer routines substantially faster?)

Best, Arnd

Re: [Numpy-discussion] Numpy Benchmarking

[Jon W. <wr...@es...>]

>>This strikes me as a little bit odd. Why not just provide the best-performing 
>>function to both SciPy and NumPy? Would NumPy be more difficult to install
>>if the SciPy algorithm for inv() was incorporated?
>>    
>>
Having spent a few days recently trying out various different 
eigenvector routines in Lapack I would have greatly appreciated having a 
choice of which one to use from without having to create my own 
wrappers, compiling atlas and lapack under windows (ouch). I noted that 
Numeric (24.2) seemed to be converting Float32 to double meaning my 
problem no longer fits in memory, which was the motivation for the work. 
Poking around in the svn of numpy.linalg appears to find the same lapack 
routine as Numeric (dsyevd). Perhaps I miss something in the code logic? 
The divide and conquer (*evd) uses more memory than the (*ev), as well 
as a factor of 2 for float/double, hence my problem, and the reason why 
"best performing" is a hard choice. I thought matlab has a look at the 
matrix dimensions and problem before deciding what to do (eg: the \ 
operator).

Jon

Re: [Numpy-discussion] Numpy Benchmarking

[Robert K. <rob...@gm...>]

jo...@st... wrote:
> Hi,
> 
>  [TO]: NumPy uses Numeric's old wrapper to lapack algorithms. 
>  [TO]: 
>  [TO]: SciPy uses it's own f2py-generated wrapper (it doesn't rely on the
>  [TO]: NumPy wrapper).
>  [TO]: 
>  [TO]: The numpy.dual library exists so you can use the SciPy calls if the 
>  [TO]: person has SciPy installed or the NumPy ones otherwise.  It exists 
>  [TO]: precisely for the purpose of seamlessly taking advantage of 
>  [TO]: algorithms/interfaces that exist in NumPy but  are improved in SciPy.
> 
> This strikes me as a little bit odd. Why not just provide the best-performing 
> function to both SciPy and NumPy? Would NumPy be more difficult to install
> if the SciPy algorithm for inv() was incorporated?

That's certainly the case for the FFT algorithms. Scipy wraps more (and more 
complicated) FFT libraries that are faster than FFTPACK.

Most of the linalg functionality should probably be wrapping the same routines 
if an optimized LAPACK is available. However, changing the routine used in numpy 
in the absence of an optimized LAPACK would require reconstructing the f2c'ed 
lapack_lite library that we include with the numpy source. That hasn't been 
touched in so long that I would hesitate to do so. If you are willing to do the 
work and the testing to ensure that it still works everywhere, we'd probably 
accept the change.

-- 
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless enigma
  that is made terrible by our own mad attempt to interpret it as though it had
  an underlying truth."
   -- Umberto Eco

Re: [Numpy-discussion] Numpy Benchmarking

[<jo...@st...>]

Hi,

 [TO]: NumPy uses Numeric's old wrapper to lapack algorithms. 
 [TO]: 
 [TO]: SciPy uses it's own f2py-generated wrapper (it doesn't rely on the
 [TO]: NumPy wrapper).
 [TO]: 
 [TO]: The numpy.dual library exists so you can use the SciPy calls if the 
 [TO]: person has SciPy installed or the NumPy ones otherwise.  It exists 
 [TO]: precisely for the purpose of seamlessly taking advantage of 
 [TO]: algorithms/interfaces that exist in NumPy but  are improved in SciPy.

This strikes me as a little bit odd. Why not just provide the best-performing 
function to both SciPy and NumPy? Would NumPy be more difficult to install
if the SciPy algorithm for inv() was incorporated?

Joris




Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm

Re: [Numpy-discussion] (no subject)

[<ka...@ya...>]

:―― INFORMATION  ―――――――――――――――――――――――――:
 不正・悪質なサイトを一切排除しておりますので
 安心してご利用ください。   http://love-match.bz/pc/?03
:――――――――――――――――――――――――――――――――――:

*・゜゜・*:.。. .。.:*・゜゜・*:.。..。:*・゜゜・*:.。..。:**・゜゜・*


━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ 不倫・ワリキリ専門の無料出会いサイト『Love☆Match』
　-----------------------------------------------------------------
　　登録料・利用料　・・・・・・・・・【無料】
　　メールの送受信　・・・・・・・・・【無料】
　　ユーザーの検索　・・・・・・・・・【無料】
　　掲示板の閲覧・書込み　・・・・・・【無料】
　　画像交換・アップロード　・・・・・【無料】
　　アドレス交換・電話番号交換　・・・【無料】
　-----------------------------------------------------------------
　　どれだけ使っても全て無料！　http://love-match.bz/pc/?03

━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ いつでも女性ユーザーがいっぱい。その理由は？

　-----------------------------------------------------------------
　ＰＣ＆モバイルに対応！いつでもどこでも気軽に楽しめる！
　-----------------------------------------------------------------
　仕事中は携帯電話から、プライベートは自宅からのんびりと。
　気になる相手といつでも繋がっているから、新密度も急速にＵＰ。
　http://love-match.bz/pc/?03

　-----------------------------------------------------------------
　ＰＣから簡単プロフィール作成。ネット初心者でもラクラク参加ＯＫ
　-----------------------------------------------------------------
　面倒な登録は一切不要。パソコンから簡単なプロフィールを作成して
　初心者の方や女性でもすぐに参加できます。
　http://love-match.bz/pc/?03
　-----------------------------------------------------------------
　自由恋愛対応！直電・直メ交換支援ツール
　-----------------------------------------------------------------
　基本的にメールアドレスや電話番号は非公開ですが
　仲良くなった人だけにメールアドレスや電話番号を教える事ができます。
　http://love-match.bz/pc/?03

　-----------------------------------------------------------------
　写真アップロードに対応！好みの相手を素早くＣＨＥＣＫ！
　-----------------------------------------------------------------
　待ち合わせ場所にイメージとまったく違う人が来たら…。
　ピュアックスなら会う前に写真交換ができるから、そんな不安も解消。
　http://love-match.bz/pc/?03

　-----------------------------------------------------------------
　スレッド掲示板で秘密のパートナー検索も効率ＵＰ！
　-----------------------------------------------------------------
　メインの掲示板のほかにスレッド型の掲示板を設置。
　メル友から秘密のパートナーまで目的別のユーザーが集う掲示板です。
　http://love-match.bz/pc/?03


━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ 毎日500人近くのユーザーが続々参加中！！

□-----------------------------------------------------------------

　リエ（２１歳/会社員）
　いつも１人でエッチなことを考えてます。
　メールだといろいろ話せるんだけど、実際に会うとあまりしゃべれなく
　なっちゃうので、盛り上げてくれるような楽しい男の人いないかな？
　引っ込み思案のせいか、男性経験はあまり無いんです。
　優しく＆楽しくリードしてくれる男性からのメール待ってます。
　[写真有り] http://love-match.bz/pc/?03

□-----------------------------------------------------------------

　真菜（２４歳/フリーター）
　彼氏が浮気してて超アタマきたっ！まなだって遊びたい盛りだし、ずっと
　ガマンしてたのにさ！かっこいい人見つけて思いっきりふってやるつもりで
　登録してみた（笑）
　[写真有り] http://love-match.bz/pc/?03

□-----------------------------------------------------------------

　みさ（３４歳/専業主婦）
　殆ど家に帰ってこない仕事人間のだんなさまと二人きりの毎日で、ほんと
　寂しい思いをしています。年下の男の子がいれば仲良くなりたいな。
　年下の人とは付き合ったことがないので興味津々です(^^)
　[写真無し] http://love-match.bz/pc/?03

□-----------------------------------------------------------------

　恭子（２８歳/会社員）
　彼氏とはいつも同じようなセックスばかりでかなり冷め気味です。
　誰かあたしと熱いセックスを楽しみませんか？めんどくさい事は
　言いません。ただ、いつもと違うドキドキするような事がしたい
　だけなんです。
　[写真無し] http://love-match.bz/pc/?03

□-----------------------------------------------------------------

　ななこ（２８歳/主婦）
　半年前にだんなと別れて今は×1です。
　夜のお仕事なので、昼間まったりと過ごしませんか？
　心身ともに疲れ気味で、今、激しく癒されたいです。
　[写真有り] http://love-match.bz/pc/?03

□-----------------------------------------------------------------

　祥子（３１歳/クリエイター）
　平日は18時くらいまでは大体仕事してるので、その後に食事したり
　楽しく飲んだりできるパートナー希望です。年上でも年下でも
　かまいませんので気軽にメールを送って頂けると嬉しいです。
　[写真有り] http://love-match.bz/pc/?03

□-----------------------------------------------------------------

　ゅヵ`（２０歳/学生）
　まずゎ会ってみないとはじまらなぃょね？！
　横浜近辺の人で、いろんな意味でオトナな人は
　プロフ付きでめぇる送って☆
　[写真有り] http://love-match.bz/pc/?03

□-----------------------------------------------------------------


   出会いサイトのサクラに騙されないように↓
　┏━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┓
　　【裏】無料の出会い情報
　　-------------------------------------------------------------
　　お金と時間を持て余している人妻の間で、噂になってるあのサイト
　　[登録・利用料全て無料]　http://love-match.bz/pc/?03
　　-------------------------------------------------------------
　　彼女達が求めているのはこんな男性です。
　　?年上女性にリードしてもらいたい、経験少なめの男性
　　?体力・テクニックに自信が有る男性
　　男性会員が不足しています。我こそは、と思う方は今すぐ参加！
　　[登録・利用料全て無料]　http://love-match.bz/pc/?03
　┗━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┛

広東省茂名市人民大街3-6-4-533
友誼網絡公司
139-3668-7892

Re: [Numpy-discussion] (no subject)

[<em...@ya...>]

:―― INFORMATION  ―――――――――――――――――――――――――:
 不正・悪質なサイトを一切排除しておりますので
 安心してご利用ください。   http://love-match.bz/pc/?010
:――――――――――――――――――――――――――――――――――:

*・゜゜・*:.。. .。.:*・゜゜・*:.。..。:*・゜゜・*:.。..。:**・゜゜・*

お金と時間を持て余している人妻の間で、噂になってるあのサイト
　　[登録・利用料全て無料]　http://love-match.bz/pc/?010
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ 不倫・ワリキリ専門の無料出会いサイト『Love☆Match』
　-----------------------------------------------------------------
　　登録料・利用料　・・・・・・・・・【無料】
　　メールの送受信　・・・・・・・・・【無料】
　　ユーザーの検索　・・・・・・・・・【無料】
　　掲示板の閲覧・書込み　・・・・・・【無料】
　　画像交換・アップロード　・・・・・【無料】
　　アドレス交換・電話番号交換　・・・【無料】
　-----------------------------------------------------------------
　　どれだけ使っても全て無料！　http://love-match.bz/pc/?010

━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ いつでも女性ユーザーがいっぱい。その理由は？

　-----------------------------------------------------------------
　ＰＣ＆モバイルに対応！いつでもどこでも気軽に楽しめる！
　-----------------------------------------------------------------
　仕事中は携帯電話から、プライベートは自宅からのんびりと。
　気になる相手といつでも繋がっているから、新密度も急速にＵＰ。
　http://love-match.bz/pc/?010

　-----------------------------------------------------------------
　ＰＣから簡単プロフィール作成。ネット初心者でもラクラク参加ＯＫ
　-----------------------------------------------------------------
　面倒な登録は一切不要。パソコンから簡単なプロフィールを作成して
　初心者の方や女性でもすぐに参加できます。
　http://love-match.bz/pc/?010
　-----------------------------------------------------------------
　自由恋愛対応！直電・直メ交換支援ツール
　-----------------------------------------------------------------
　基本的にメールアドレスや電話番号は非公開ですが
　仲良くなった人だけにメールアドレスや電話番号を教える事ができます。
　http://love-match.bz/pc/?010

　-----------------------------------------------------------------
　写真アップロードに対応！好みの相手を素早くＣＨＥＣＫ！
　-----------------------------------------------------------------
　待ち合わせ場所にイメージとまったく違う人が来たら…。
　ピュアックスなら会う前に写真交換ができるから、そんな不安も解消。
　http://love-match.bz/pc/?010

　-----------------------------------------------------------------
　スレッド掲示板で秘密のパートナー検索も効率ＵＰ！
　-----------------------------------------------------------------
　メインの掲示板のほかにスレッド型の掲示板を設置。
　メル友から秘密のパートナーまで目的別のユーザーが集う掲示板です。
　http://love-match.bz/pc/?010


━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ 毎日500人近くのユーザーが続々参加中！！

□-----------------------------------------------------------------

　リエ（２１歳/会社員）
　いつも１人でエッチなことを考えてます。
　メールだといろいろ話せるんだけど、実際に会うとあまりしゃべれなく
　なっちゃうので、盛り上げてくれるような楽しい男の人いないかな？
　引っ込み思案のせいか、男性経験はあまり無いんです。
　優しく＆楽しくリードしてくれる男性からのメール待ってます。
　[写真有り] http://love-match.bz/pc/?010

□-----------------------------------------------------------------

　真菜（２４歳/フリーター）
　彼氏が浮気してて超アタマきたっ！まなだって遊びたい盛りだし、ずっと
　ガマンしてたのにさ！かっこいい人見つけて思いっきりふってやるつもりで
　登録してみた（笑）
　[写真有り] http://love-match.bz/pc/?010

□-----------------------------------------------------------------

　みさ（３４歳/専業主婦）
　殆ど家に帰ってこない仕事人間のだんなさまと二人きりの毎日で、ほんと
　寂しい思いをしています。年下の男の子がいれば仲良くなりたいな。
　年下の人とは付き合ったことがないので興味津々です(^^)
　[写真無し] http://love-match.bz/pc/?010

□-----------------------------------------------------------------

　恭子（２８歳/会社員）
　彼氏とはいつも同じようなセックスばかりでかなり冷め気味です。
　誰かあたしと熱いセックスを楽しみませんか？めんどくさい事は
　言いません。ただ、いつもと違うドキドキするような事がしたい
　だけなんです。
　[写真無し] http://love-match.bz/pc/?010

□-----------------------------------------------------------------

　ななこ（２８歳/主婦）
　半年前にだんなと別れて今は×1です。
　夜のお仕事なので、昼間まったりと過ごしませんか？
　心身ともに疲れ気味で、今、激しく癒されたいです。
　[写真有り] http://love-match.bz/pc/?010

□-----------------------------------------------------------------

　祥子（３１歳/クリエイター）
　平日は18時くらいまでは大体仕事してるので、その後に食事したり
　楽しく飲んだりできるパートナー希望です。年上でも年下でも
　かまいませんので気軽にメールを送って頂けると嬉しいです。
　[写真有り] http://love-match.bz/pc/?010

□-----------------------------------------------------------------

　ゅヵ`（２０歳/学生）
　まずゎ会ってみないとはじまらなぃょね？！
　横浜近辺の人で、いろんな意味でオトナな人は
　プロフ付きでめぇる送って☆
　[写真有り] http://love-match.bz/pc/?010

□-----------------------------------------------------------------


   出会いサイトのサクラに騙されないように↓
　┏━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┓
　　【裏】無料の出会い情報
　　-------------------------------------------------------------
　　お金と時間を持て余している人妻の間で、噂になってるあのサイト
　　[登録・利用料全て無料]　http://love-match.bz/pc/?010
　　-------------------------------------------------------------
　　彼女達が求めているのはこんな男性です。
　　?年上女性にリードしてもらいたい、経験少なめの男性
　　?体力・テクニックに自信が有る男性
　　男性会員が不足しています。我こそは、と思う方は今すぐ参加！
　　[登録・利用料全て無料]　http://love-match.bz/pc/?010
　┗━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┛

広東省茂名市人民大街3-6-4-533
友誼網絡公司
139-3668-7892

Re: [Numpy-discussion] (no subject)

[<mik...@ya...>]

:―― INFORMATION  ―――――――――――――――――――――――――:
 不正・悪質なサイトを一切排除しておりますので
 安心してご利用ください。   http://love-match.bz/pc/?07
:――――――――――――――――――――――――――――――――――:

*・゜゜・*:.。. .。.:*・゜゜・*:.。..。:*・゜゜・*:.。..。:**・゜゜・*


━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ 不倫・ワリキリ専門の無料出会いサイト『Love☆Match』
　-----------------------------------------------------------------
　　登録料・利用料　・・・・・・・・・【無料】
　　メールの送受信　・・・・・・・・・【無料】
　　ユーザーの検索　・・・・・・・・・【無料】
　　掲示板の閲覧・書込み　・・・・・・【無料】
　　画像交換・アップロード　・・・・・【無料】
　　アドレス交換・電話番号交換　・・・【無料】
　-----------------------------------------------------------------
　　どれだけ使っても全て無料！　http://love-match.bz/pc/?07

━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ いつでも女性ユーザーがいっぱい。その理由は？

　-----------------------------------------------------------------
　ＰＣ＆モバイルに対応！いつでもどこでも気軽に楽しめる！
　-----------------------------------------------------------------
　仕事中は携帯電話から、プライベートは自宅からのんびりと。
　気になる相手といつでも繋がっているから、新密度も急速にＵＰ。
　http://love-match.bz/pc/?07

　-----------------------------------------------------------------
　ＰＣから簡単プロフィール作成。ネット初心者でもラクラク参加ＯＫ
　-----------------------------------------------------------------
　面倒な登録は一切不要。パソコンから簡単なプロフィールを作成して
　初心者の方や女性でもすぐに参加できます。
　http://love-match.bz/pc/?07
　-----------------------------------------------------------------
　自由恋愛対応！直電・直メ交換支援ツール
　-----------------------------------------------------------------
　基本的にメールアドレスや電話番号は非公開ですが
　仲良くなった人だけにメールアドレスや電話番号を教える事ができます。
　http://love-match.bz/pc/?07

　-----------------------------------------------------------------
　写真アップロードに対応！好みの相手を素早くＣＨＥＣＫ！
　-----------------------------------------------------------------
　待ち合わせ場所にイメージとまったく違う人が来たら…。
　ピュアックスなら会う前に写真交換ができるから、そんな不安も解消。
　http://love-match.bz/pc/?07

　-----------------------------------------------------------------
　スレッド掲示板で秘密のパートナー検索も効率ＵＰ！
　-----------------------------------------------------------------
　メインの掲示板のほかにスレッド型の掲示板を設置。
　メル友から秘密のパートナーまで目的別のユーザーが集う掲示板です。
　http://love-match.bz/pc/?07


━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
　□■ 毎日500人近くのユーザーが続々参加中！！

□-----------------------------------------------------------------

　リエ（２１歳/会社員）
　いつも１人でエッチなことを考えてます。
　メールだといろいろ話せるんだけど、実際に会うとあまりしゃべれなく
　なっちゃうので、盛り上げてくれるような楽しい男の人いないかな？
　引っ込み思案のせいか、男性経験はあまり無いんです。
　優しく＆楽しくリードしてくれる男性からのメール待ってます。
　[写真有り] http://love-match.bz/pc/?07

□-----------------------------------------------------------------

　真菜（２４歳/フリーター）
　彼氏が浮気してて超アタマきたっ！まなだって遊びたい盛りだし、ずっと
　ガマンしてたのにさ！かっこいい人見つけて思いっきりふってやるつもりで
　登録してみた（笑）
　[写真有り] http://love-match.bz/pc/?07

□-----------------------------------------------------------------

　みさ（３４歳/専業主婦）
　殆ど家に帰ってこない仕事人間のだんなさまと二人きりの毎日で、ほんと
　寂しい思いをしています。年下の男の子がいれば仲良くなりたいな。
　年下の人とは付き合ったことがないので興味津々です(^^)
　[写真無し] http://love-match.bz/pc/?07

□-----------------------------------------------------------------

　恭子（２８歳/会社員）
　彼氏とはいつも同じようなセックスばかりでかなり冷め気味です。
　誰かあたしと熱いセックスを楽しみませんか？めんどくさい事は
　言いません。ただ、いつもと違うドキドキするような事がしたい
　だけなんです。
　[写真無し] http://love-match.bz/pc/?07

□-----------------------------------------------------------------

　ななこ（２８歳/主婦）
　半年前にだんなと別れて今は×1です。
　夜のお仕事なので、昼間まったりと過ごしませんか？
　心身ともに疲れ気味で、今、激しく癒されたいです。
　[写真有り] http://love-match.bz/pc/?07

□-----------------------------------------------------------------

　祥子（３１歳/クリエイター）
　平日は18時くらいまでは大体仕事してるので、その後に食事したり
　楽しく飲んだりできるパートナー希望です。年上でも年下でも
　かまいませんので気軽にメールを送って頂けると嬉しいです。
　[写真有り] http://love-match.bz/pc/?07

□-----------------------------------------------------------------

　ゅヵ`（２０歳/学生）
　まずゎ会ってみないとはじまらなぃょね？！
　横浜近辺の人で、いろんな意味でオトナな人は
　プロフ付きでめぇる送って☆
　[写真有り] http://love-match.bz/pc/?07

□-----------------------------------------------------------------


   出会いサイトのサクラに騙されないように↓
　┏━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┓
　　【裏】無料の出会い情報
　　-------------------------------------------------------------
　　お金と時間を持て余している人妻の間で、噂になってるあのサイト
　　[登録・利用料全て無料]　http://love-match.bz/pc/?07
　　-------------------------------------------------------------
　　彼女達が求めているのはこんな男性です。
　　?年上女性にリードしてもらいたい、経験少なめの男性
　　?体力・テクニックに自信が有る男性
　　男性会員が不足しています。我こそは、と思う方は今すぐ参加！
　　[登録・利用料全て無料]　http://love-match.bz/pc/?07
　┗━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┛

広東省茂名市人民大街3-6-4-533
友誼網絡公司
139-3668-7892

Re: [Numpy-discussion] Should cholesky return upper or lower triangular matrix?

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

On 6/27/06, Robert Kern <rob...@gm...> wrote:
> Keith Goodman wrote:
> > Isn't the Cholesky decomposition by convention an upper triangular
> > matrix? I noticed, by porting Octave code, that linalg.cholesky
> > returns the lower triangular matrix.
> >
> > References:
> >
> > http://mathworld.wolfram.com/CholeskyDecomposition.html
> > http://www.mathworks.com/access/helpdesk/help/techdoc/ref/chol.html
>
> Lower:
> http://en.wikipedia.org/wiki/Cholesky_decomposition
> http://www.math-linux.com/spip.php?article43
> http://planetmath.org/?op=getobj&from=objects&id=1287
> http://rkb.home.cern.ch/rkb/AN16pp/node33.html#SECTION000330000000000000000
> http://www.riskglossary.com/link/cholesky_factorization.htm
> http://www.library.cornell.edu/nr/bookcpdf/c2-9.pdf
>
> If anything, the convention appears to be lower-triangular.

If you give me a second, I'll show you that the wikipedia supports my claim.

OK. Lower it is. It will save me a transpose when I calculate joint
random variables.

Re: [Numpy-discussion] Should cholesky return upper or lower triangular matrix?

[Robert K. <rob...@gm...>]

Keith Goodman wrote:
> Isn't the Cholesky decomposition by convention an upper triangular
> matrix? I noticed, by porting Octave code, that linalg.cholesky
> returns the lower triangular matrix.
> 
> References:
> 
> http://mathworld.wolfram.com/CholeskyDecomposition.html
> http://www.mathworks.com/access/helpdesk/help/techdoc/ref/chol.html

Lower:
http://en.wikipedia.org/wiki/Cholesky_decomposition
http://www.math-linux.com/spip.php?article43
http://planetmath.org/?op=getobj&from=objects&id=1287
http://rkb.home.cern.ch/rkb/AN16pp/node33.html#SECTION000330000000000000000
http://www.riskglossary.com/link/cholesky_factorization.htm
http://www.library.cornell.edu/nr/bookcpdf/c2-9.pdf

If anything, the convention appears to be lower-triangular.

-- 
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless enigma
  that is made terrible by our own mad attempt to interpret it as though it had
  an underlying truth."
   -- Umberto Eco

[Numpy-discussion] Should cholesky return upper or lower triangular matrix?

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

Isn't the Cholesky decomposition by convention an upper triangular
matrix? I noticed, by porting Octave code, that linalg.cholesky
returns the lower triangular matrix.

References:

http://mathworld.wolfram.com/CholeskyDecomposition.html
http://www.mathworks.com/access/helpdesk/help/techdoc/ref/chol.html

Re: [Numpy-discussion] Numpy Benchmarking

[Robert K. <rob...@gm...>]

Keith Goodman wrote:
> Scipy computes the inverse of a matrix faster than numpy (except if
> the dimensions of x are small). But scipy is slower than numpy for
> eigh (I only checked for symmetric positive definite matrices):

Looks like scipy uses *SYEV and numpy uses the better *SYEVD (the D stands for 
divide-and-conquer) routine. Both should probably be using the RRR versions 
(*SYEVR) if I'm reading the advice in the LUG correctly.

-- 
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless enigma
  that is made terrible by our own mad attempt to interpret it as though it had
  an underlying truth."
   -- Umberto Eco

Re: [Numpy-discussion] Numpy Benchmarking

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

On 6/27/06, Keith Goodman <kwg...@gm...> wrote:
> On 6/27/06, Travis Oliphant <oli...@ie...> wrote:
>
> > The numpy.dual library exists so you can use the SciPy calls if the
> > person has SciPy installed or the NumPy ones otherwise.  It exists
> > precisely for the purpose of seamlessly taking advantage of
> > algorithms/interfaces that exist in NumPy but  are improved in SciPy.
>
> That sounds very interesting. It would make a great addition to the
> scipy performance page:
>
> http://scipy.org/PerformanceTips
>
> So if I need any of the following functions I should import them from
> scipy or from numpy.dual? And all of them are faster?
>
> fft
> ifft
> fftn
> ifftn
> fft2
> ifft2
> norm
> inv
> svd
> solve
> det
> eig
> eigvals
> eigh
> eigvalsh
> lstsq
> pinv
> cholesky
>
> http://svn.scipy.org/svn/numpy/trunk/numpy/dual.py
>

Scipy computes the inverse of a matrix faster than numpy (except if
the dimensions of x are small). But scipy is slower than numpy for
eigh (I only checked for symmetric positive definite matrices):

from numpy import asmatrix, randn
from numpy.linalg import eigh as Neigh
from scipy.linalg import eigh as Seigh
import time

def test(N):

  x = asmatrix(randn(N,2*N))
  x = x * x.T

  t0 = time.time()
  eigval, eigvec = Neigh(x)
  t1 = time.time()

  t2 = time.time()
  eigval, eigvec = Seigh(x)
  t3 = time.time()

  print 'NumPy:', t1-t0, 'seconds'
  print 'SciPy:', t3-t2, 'seconds'

>> dual.test(10)
NumPy: 0.000217914581299 seconds
SciPy: 0.000226020812988 seconds

>> dual.test(100)
NumPy: 0.0123109817505 seconds
SciPy: 0.0321230888367 seconds

>> dual.test(200)
NumPy: 0.0793058872223 seconds
SciPy: 0.082535982132 seconds

>> dual.test(500)
NumPy: 0.59161400795 seconds
SciPy: 1.41600894928 seconds

Re: [Numpy-discussion] Numpy Benchmarking

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

On 6/27/06, Travis Oliphant <oli...@ie...> wrote:

> The numpy.dual library exists so you can use the SciPy calls if the
> person has SciPy installed or the NumPy ones otherwise.  It exists
> precisely for the purpose of seamlessly taking advantage of
> algorithms/interfaces that exist in NumPy but  are improved in SciPy.

That sounds very interesting. It would make a great addition to the
scipy performance page:

http://scipy.org/PerformanceTips

So if I need any of the following functions I should import them from
scipy or from numpy.dual? And all of them are faster?

fft
ifft
fftn
ifftn
fft2
ifft2
norm
inv
svd
solve
det
eig
eigvals
eigh
eigvalsh
lstsq
pinv
cholesky

http://svn.scipy.org/svn/numpy/trunk/numpy/dual.py

Re: [Numpy-discussion] Numpy Benchmarking

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

Keith Goodman wrote:
> On 6/27/06, Dennis V. Perepelitsa <dv...@mi...> wrote:
>
>   
>> I've run some benchmarks comparing the performance of scipy, numpy,
>> Numeric and numarray vs. MATLAB.
>>     
>
> I enjoyed looking at the results.
>
> The most interesting result, for me, was that inverting a matrix is
> much faster in scipy than numpy. How can that be? I would have guessed
> that numpy handled the inversion for scipy since numpy is the core.
>
> The two calls were scipy.linalg.inv(m) and numpy.linalg.inv(m).
>   
NumPy uses Numeric's old wrapper to lapack algorithms. 

SciPy uses it's own f2py-generated wrapper (it doesn't rely on the NumPy 
wrapper).

The numpy.dual library exists so you can use the SciPy calls if the 
person has SciPy installed or the NumPy ones otherwise.  It exists 
precisely for the purpose of seamlessly taking advantage of 
algorithms/interfaces that exist in NumPy but  are improved in SciPy.


-Travis