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 wrote: > On Thu, Aug 11, 2011 at 4:43 PM, Jose Borreguero > wrote: > > a = random.randn(3,3) > > b = a.reshape(1,3,3).repeat(50,axis=0) > > scipy.linalg.block_diag( *b ) > > > > slightly simpler,

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

2011-08-11 Thread Jose Borreguero
Dear numpy users, I have a 3x3 matrix which I want to repeat 50 times along a diagonal, thus creating a 150x150 block diagonal matrix. I know of a method usin scipy.linalg.block_diag, but I don't know if this is the best one: a = random.randn(3,3) b = a.reshape(1,3,3).repeat(50,axis=0) scipy.lina

Re: [Numpy-discussion] ImportError: libatlas.so: cannot open shared object file: No such file or directory

2011-03-17 Thread Jose Borreguero
or, I found some comments on mixing fortran compilers. I haven't advance much because go from error to error. Maybe I'm mixing things up by adding directories to LD_LIBRARY_PATH :( Jose On Thu, Mar 17, 2011 at 1:35 PM, Ilan Schnell wrote: > It looks like atlas wasn't linked ri

[Numpy-discussion] ImportError: libatlas.so: cannot open shared object file: No such file or directory

2011-03-17 Thread Jose Borreguero
Dear Numpy/Scipy users, I just installed numpy but I have an error when importing >>> import numpy from numpy.linalg import lapack_lite ImportError: libatlas.so: cannot open shared object file: No such file or directory I have ATLAS libraries under /usr/local/atlas/lib libatlas.a libcb

Re: [Numpy-discussion] Build ERROR ConfigParser.MissingSectionHeaderError: File contains no section headers

2011-03-17 Thread Jose Borreguero
Thanks Josef, that worked! I was confused because I was thinking of site.cfg as some sort of bash script :) Jose On Thu, Mar 17, 2011 at 12:40 AM, wrote: > On Thu, Mar 17, 2011 at 12:23 AM, Jose Borreguero > wrote: > > Dear Numpy/SciPy users, > > > > I have

[Numpy-discussion] Build ERROR ConfigParser.MissingSectionHeaderError: File contains no section headers

2011-03-16 Thread Jose Borreguero
Dear Numpy/SciPy users, I have a build error with Numpy: $ /usr/local/bin/python2.7 setup.py build File "/usr/local/lib/python2.7/ConfigParser.py", line 504, in _read raise MissingSectionHeaderError(fpname, lineno, line) ConfigParser.MissingSectionHeaderError: File contains no section

Re: [Numpy-discussion] adding two pseudoinverse matrixes

2010-09-10 Thread Jose Borreguero
in a molecular dynamics simulation, A = . Same goes for B. -Jose On Fri, Sep 10, 2010 at 5:02 PM, Charles R Harris wrote: > > > On Fri, Sep 10, 2010 at 2:39 PM, Jose Borreguero wrote: > >> Dear Numpy users, >> >> I have to solve for Z in the following equation Z

[Numpy-discussion] adding two pseudoinverse matrixes

2010-09-10 Thread Jose Borreguero
brute force approach linalg.pinv( linalg.pinv(A) - lingal.pinv(B) ) gives me a matrix with all entries equal to 'infinity'. -Jose Borreguero ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] how to multiply the rows of a matrix by a different number?

2009-03-04 Thread Jose Borreguero
Sweet! I found that *M*b.reshape(1,1)* will also do the trick. Any guess which method is faster? On Tue, Mar 3, 2009 at 9:11 PM, wrote: > On Tue, Mar 3, 2009 at 8:53 PM, Jose Borreguero > wrote: > > I guess there has to be an easy way for this. I have: > > M.shape=(10

[Numpy-discussion] how to multiply the rows of a matrix by a different number?

2009-03-03 Thread Jose Borreguero
I guess there has to be an easy way for this. I have: M.shape=(1,3) N.shape=(1,) I want to do this: for i in range(1): M[i]*=N[i] without the explicit loop -Jose ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.s

Re: [Numpy-discussion] tensordot bug when summing over same axis indexes?

2009-02-24 Thread Jose Borreguero
haha, so it was a stupid error...my stupid error. [?] I incorrectly understood ([0,0],[1,1]) as index 0 of *a* summed with index 0 of *b*, and analogously for [1,1]. Gthanks, Josef On Tue, Feb 24, 2009 at 3:20 PM, wrote: > On Tue, Feb 24, 2009 at 2:55 PM, Jose Borreguero > wrote: &

[Numpy-discussion] tensordot bug when summing over same axis indexes?

2009-02-24 Thread Jose Borreguero
The following example: from numpy import * a=arange(12).reshape(2,3,2) b=arange(24).reshape(2,3,2,2) c=tensordot( a,b,axes=([0,0],[1,1]) ) defaults: c=tensordot( a,b,axes=([0,0],[1,1]) ) File "/usr/lib/python2.4/site-packages/numpy/core/numeric.py", line 359, in tensordot raise ValueError, "shape

Re: [Numpy-discussion] creating a Numeric array from a numpy array LEAKS memory

2008-10-24 Thread Jose Borreguero
=Numeric.array(na,Numeric.Float)* -Jose On Fri, Oct 24, 2008 at 2:16 PM, Jose Borreguero <[EMAIL PROTECTED]>wrote: > numpy 1.1.0 (from /usr/lib/python2.4/site-packages/numpy/version.py) > Numeric 24.2 (from > /usr/lib/python2.4/site-packages/Numeric/numeric_version.py) > > I also tried

Re: [Numpy-discussion] creating a Numeric array from a numpy array LEAKS memory

2008-10-24 Thread Jose Borreguero
have memory leaks if I use something like: *mylist=[0.0]*BIGNUMBER* *na=Numeric.array( mylist, Numeric.Float)* -Jose On Fri, Oct 24, 2008 at 1:54 PM, Travis E. Oliphant <[EMAIL PROTECTED]>wrote: > Jose Borreguero wrote: > > Dear numpy users, > > > > I need to pass a

[Numpy-discussion] creating a Numeric array from a numpy array LEAKS memory

2008-10-24 Thread Jose Borreguero
Dear numpy users, I need to pass a Numeric array to some oldie code from a numpy array. I decided to go like this: for i in range(BIGNUMER): my_numpy_array=grabArray(i) na=Numeric.array( my_numpy_array, Numeric.Float) oldie_code(na) The constructor line: na=Numeric.array( my_numpy_ar