Re: [Numpy-discussion] SVD does not converge on clean matrix

2011-08-12 Thread Charanpal Dhanjal
Thank Nadav for testing out the matrix. I wonder if you had a chance to check if the resulting decomposition contained NaN or Inf values? As far I understood, numpy.linalg.svd uses routines in LAPACK and ATLAS (if available) to compute the corresponding SVD. I did some complementary tests on

Re: [Numpy-discussion] how to create a block diagonal matrix by repeating the block?

2011-08-12 Thread Jose Borreguero
Thanks! Jose On Thu, Aug 11, 2011 at 8:15 PM, Fernando Perez fperez@gmail.comwrote: On Thu, Aug 11, 2011 at 4:43 PM, Jose Borreguero borregu...@gmail.com wrote: a = random.randn(3,3) b = a.reshape(1,3,3).repeat(50,axis=0) scipy.linalg.block_diag( *b ) slightly simpler, but

[Numpy-discussion] Statistical distributions on samples

2011-08-12 Thread Andrea Gavana
Hi All, I am working on something that appeared to be a no-brainer issue (at the beginning), by my complete ignorance in statistics is overwhelming and I got stuck. What I am trying to do can be summarized as follows Let's assume that I have to generate a sample of a 1,000 values for a

Re: [Numpy-discussion] SVD does not converge on clean matrix

2011-08-12 Thread Warren Weckesser
On Fri, Aug 12, 2011 at 4:03 AM, Charanpal Dhanjal dhan...@telecom-paristech.fr wrote: Thank Nadav for testing out the matrix. I wonder if you had a chance to check if the resulting decomposition contained NaN or Inf values? As far I understood, numpy.linalg.svd uses routines in LAPACK and

[Numpy-discussion] nditer confusion

2011-08-12 Thread Neal Becker
There'a a boatload of options for nditer. I need a simple explanation, maybe a few simple examples. Is there anything that might help? ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org

Re: [Numpy-discussion] Statistical distributions on samples

2011-08-12 Thread Christopher Jordan-Squire
Hi Andrea--An easy way to get something like this would be import numpy as np import scipy.stats as stats sigma = #some reasonable standard deviation for your application x = stats.norm.rvs(size=1000, loc=125, scale=sigma) x = x[x50] x = x[x200] That will give a roughly normal distribution to

Re: [Numpy-discussion] SVD does not converge on clean matrix

2011-08-12 Thread Nadav Horesh
I tested all the the result 3 matrices with alltrue(infinite(mat)) and got True answer for all of them. Nadav From: numpy-discussion-boun...@scipy.org [numpy-discussion-boun...@scipy.org] On Behalf Of Warren Weckesser [warren.weckes...@enthought.com] Sent:

Re: [Numpy-discussion] nditer confusion

2011-08-12 Thread Mark Wiebe
I'll write up some more introductory-style documentation, you're right that the examples I put in the reference page aren't a nice simple starting point. Will post back here for feedback when I have a draft for you to review. Cheers, Mark On Fri, Aug 12, 2011 at 7:30 AM, Neal Becker

[Numpy-discussion] Theano 0.4.1 released

2011-08-12 Thread Frédéric Bastien
=== Announcing Theano 0.4.1 === This is an important release, with lots of new features, bug fixes and some deprecation warning. The upgrade is recommended for everybody. For those using the bleeding edge version in the mercurial repository, we