Re: [Numpy-discussion] Scipy 0.6.0 to 0.7.0, sparse matrix change

On Mon, Jun 15, 2009 at 6:47 AM, Fadhley Salim wrote: > > I'm trying to track down a numerical discrepancy in our proejct. We > noticed that a certain set of results are different having upgraded from > scipy 0.6.0 to 0.7.0. > > The following item from the Scipy change-log is our current number-on

Re: [Numpy-discussion] Interleaved Arrays and

I'm not sure it's worth having a function to replace a one-liner (column_stack followed by reshape). But if you're going to implement this with slice assignment, you should take advantage of the flexibility this method allows and offer the possibility of interleaving "raggedly", that is, where the

Re: [Numpy-discussion] Interleaved Arrays and

On 2009-06-16 16:05 , Robert wrote: > Neil Martinsen-Burrell wrote: >> On 06/16/2009 02:18 PM, Robert wrote: >>>>>> n = 10 >>>>>> xx = np.ones(n) >>>>>> yy = np.arange(n) >>>>>> aa = np.column_stack((xx,yy)) >>>>>> bb = np.column_stack((xx+1,yy)) >>>>>> aa >>> ar

Re: [Numpy-discussion] Interleaved Arrays and

Neil Martinsen-Burrell wrote: > On 06/16/2009 02:18 PM, Robert wrote: >> >>> n = 10 >> >>> xx = np.ones(n) >> >>> yy = np.arange(n) >> >>> aa = np.column_stack((xx,yy)) >> >>> bb = np.column_stack((xx+1,yy)) >> >>> aa >> array([[ 1., 0.], >> [ 1., 1.], >> [ 1.,

Re: [Numpy-discussion] Interleaved Arrays and

On 06/16/2009 02:18 PM, Robert wrote: > >>> n = 10 > >>> xx = np.ones(n) > >>> yy = np.arange(n) > >>> aa = np.column_stack((xx,yy)) > >>> bb = np.column_stack((xx+1,yy)) > >>> aa > array([[ 1., 0.], > [ 1., 1.], > [ 1., 2.], > [ 1., 3.], >

Re: [Numpy-discussion] ticket #1096

On Tue, Jun 16, 2009 at 11:47 AM, Pauli Virtanen wrote: > Anyway, fixed in r7059. > Thanks! ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] Interleaved Arrays and

Ian Mallett wrote: > > n = #blah > testlist = [] > for x in xrange(n): > for y in xrange(n): > testlist.append([x,y]) > testlist.append([x+1,y]) > > If "testlist" is an array (i.e., I could go: "array(testlist)"), it > works nicely. However, my Python method is certainly im

Re: [Numpy-discussion] ticket #1096

On 2009-06-16, Brian Lewis wrote: > http://projects.scipy.org/numpy/ticket/1096 > > Is the fix to this to check if (line 95 of > trunk/numpy/core/src/umathmodule.c.src > ) > >const @type@ tmp = x - y; > > is

[Numpy-discussion] ticket #1096

http://projects.scipy.org/numpy/ticket/1096 Is the fix to this to check if (line 95 of trunk/numpy/core/src/umathmodule.c.src ) const @type@ tmp = x - y; is -inf or not. And if it is, just to return -inf.

Re: [Numpy-discussion] state of G3F2PY?

I would be interested in testing your GSOC project and will do what I can in the mean time. I do develop on windows, but the library lives on linux, macos, and windows, so we can test on anyg--it also binds with ifort, gfortran, etc. so seems rather robust. Cheers, William On Tue, Jun 16, 2009 a

Re: [Numpy-discussion] state of G3F2PY?

On Tue, Jun 16, 2009 at 8:21 AM, william ratcliff wrote: > Hi!  I'm looking at trying to bind a rather large (>150K lines of code) > crystallography library to python and would like to know what the state of > F2py is.  Are allocatable arrays supported?  Derived types?  Modules, > Pointers, etc.? 

Re: [Numpy-discussion] Array resize question

On 2009-06-16 09:42 , Cristi Constantin wrote: > Good day. > I have this array: > > a = array([[u'0', u'0', u'0', u'0', u'0', u' '], > [u'1', u'1', u'1', u'1', u'1', u' '], > [u'2', u'2', u'2', u'2', u'2', u' '], > [u'3', u'3', u'3', u'3', u'3', u' '], > [u'4', u'4', u'4', u'4', u'4', u'']], > dtyp

[Numpy-discussion] Array resize question

Good day. I have this array: a = array([[u'0', u'0', u'0', u'0', u'0', u' '],    [u'1', u'1', u'1', u'1', u'1', u' '],    [u'2', u'2', u'2', u'2', u'2', u' '],    [u'3', u'3', u'3', u'3', u'3', u' '],    [u'4', u'4', u'4', u'4', u'4', u'']],   dtype='http://docs.scipy.org/doc/

[Numpy-discussion] state of G3F2PY?

Hi! I'm looking at trying to bind a rather large (>150K lines of code) crystallography library to python and would like to know what the state of F2py is. Are allocatable arrays supported? Derived types? Modules, Pointers, etc.? Is there a list somewhere? Has anyone else looked into wrapping

Re: [Numpy-discussion] More on doc-ing new functions

Mon, 15 Jun 2009 09:52:05 -0700, David Goldsmith kirjoitti: [clip] > Is there a protocol for making sure these things get done? (Just don't > want to reinvent the wheel.) I don't think so. The current way is that people forget to do it, and then someone fixes it afterwards :) I'm not sure how