Re: [Numpy-discussion] numpy dot returns [nan nan nan]

2013-08-24 Thread Tom Bennett
Hi Warren, Yes you are absolutely right. I had some values close to log(x), where x is almost 0. That caused the problem. Thanks, Tom On Sat, Aug 24, 2013 at 12:39 PM, Warren Weckesser < warren.weckes...@gmail.com> wrote: > On 8/24/13, Warren Weckesser wrote: > > On 8/24/

[Numpy-discussion] numpy dot returns [nan nan nan]

2013-08-24 Thread Tom Bennett
Hi All, I have two arrays, A and B.A is 3 x 100,000 and B is 100,000. If I do np.dot(A,B), I get [nan, nan, nan]. However, np.any(np.isnan(A))==False and np.any(no.isnan(B))==False. And also np.seterr(all='print') does not print anything. I am not wondering what is going on and how to avoid. In

Re: [Numpy-discussion] How to Keep An Array Two Dimensional

2012-11-25 Thread Tom Bennett
< warren.weckes...@gmail.com> wrote: > > > On Sun, Nov 25, 2012 at 8:24 PM, Tom Bennett > wrote: > >> Hi, >> >> I am trying to extract n columns from an 2D array and then operate on the >> extracted columns. Below is the code: >> >> A is an MxN 2D a

[Numpy-discussion] How to Keep An Array Two Dimensional

2012-11-25 Thread Tom Bennett
Hi, I am trying to extract n columns from an 2D array and then operate on the extracted columns. Below is the code: A is an MxN 2D array. u = A[:,:n] #extract the first n columns from A B = np.dot(u, u.T) #take outer product. This code works when n>1. However, when n=1, u becomes an 1D array i