[Numpy-discussion] GPU blas with numpy?

2007-12-06 Thread Jon Wright
Hello,

Does anyone have any experience with numpy and using a GPU? eg:

http://developer.nvidia.com/object/cuda.html

...they seem to have a GPU based blas available.

Thanks for any info,

Jon

___
Numpy-discussion mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion


Re: [Numpy-discussion] multinomial question

2007-12-06 Thread Alan G Isaac
 Alan G Isaac wrote:
 I would think that 
 multinomial(1,prob,size=ntrials).sum(axis=0)
 would be equivalent to 
 multinomial(ntrials,prob)
 but the first gives a surprising result.  (See below.)
 Explanation? 


On Wed, 05 Dec 2007, Robert Kern apparently wrote:
 Pretty much anyone who derives their binomial distribution algorithm from 
 http://www.unc.edu/~gfish/fcmc.html is also wrong. 
 SVN now has a bound such that CDF(bound) is within 1e-16 (or so) of 1.0. 


Thanks!
Alan



___
Numpy-discussion mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion


Re: [Numpy-discussion] GPU blas with numpy?

2007-12-06 Thread Arnar Flatberg
Hi Jon

Have you looked at pystream?
http://code.google.com/p/pystream/

Arnar

On Dec 6, 2007 12:39 PM, Jon Wright [EMAIL PROTECTED] wrote:
 Hello,

 Does anyone have any experience with numpy and using a GPU? eg:

 http://developer.nvidia.com/object/cuda.html

 ...they seem to have a GPU based blas available.

 Thanks for any info,

 Jon

 ___
 Numpy-discussion mailing list
 Numpy-discussion@scipy.org
 http://projects.scipy.org/mailman/listinfo/numpy-discussion

___
Numpy-discussion mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion


[Numpy-discussion] Max on Complex Array

2007-12-06 Thread Tom Johnson
What is the meaning of numpy's maximum on a complex-valued array?
___
Numpy-discussion mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion


Re: [Numpy-discussion] Max on Complex Array

2007-12-06 Thread Robert Kern
Tom Johnson wrote:
 What is the meaning of numpy's maximum on a complex-valued array?

We impose a lexicographical ordering on the complex space. Complex numbers are
compared first by their real component and then, if the real components are
equal, by their imaginary component.


In [1]: from numpy import *

In [19]: c = random.randint(0, 5, 10) + random.random(10)*1j

In [20]: c
Out[20]:
array([ 0.+0.68448275j,  1.+0.97849291j,  3.+0.22114928j,  4.+0.65409519j,
3.+0.91550523j,  4.+0.50667105j,  1.+0.34576644j,  4.+0.97286048j,
1.+0.07268317j,  0.+0.52885086j])

In [21]: c.sort()

In [22]: c
Out[22]:
array([ 0.+0.52885086j,  0.+0.68448275j,  1.+0.07268317j,  1.+0.34576644j,
1.+0.97849291j,  3.+0.22114928j,  3.+0.91550523j,  4.+0.50667105j,
4.+0.65409519j,  4.+0.97286048j])

-- 
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 mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion


Re: [Numpy-discussion] Max on Complex Array

2007-12-06 Thread Tom Johnson
On Dec 6, 2007 3:23 PM, Robert Kern [EMAIL PROTECTED] wrote:
 Tom Johnson wrote:
  What is the meaning of numpy's maximum on a complex-valued array?

 We impose a lexicographical ordering on the complex space.

Good to know, but is this documented somewhere.  a.max? does not mention this.
___
Numpy-discussion mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion