[Numpy-discussion] GPU blas with numpy?
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
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?
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
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
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
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