From david.grant at telus.net  Wed Apr  6 22:22:49 2005
From: david.grant at telus.net (David Grant)
Date: Wed, 06 Apr 2005 19:22:49 -0700
Subject: [SciPy-dev] kron function in scipy.linalg.basic
Message-ID: <42549979.1020802@telus.net>

Re-posting to scipy-dev as it might be a bug.

I'm using the kron() function in scipy.linalg.basic but I get an error

I do this:


>>>>>> from scipy.linalg.basic import kron
>>>>>> from Numeric import *
>>>>>> a=ones((2,2))
>>>>>> b=ones((2,2))
>>>>>> kron(a,b)
>>>      
>>>
Traceback (most recent call last):
  File "<stdin>", line 1, in ?
  File "/usr/lib/python2.3/site-packages/scipy/linalg/basic.py", line
506, in kron
    o = outerproduct(a,b)
NameError: global name 'outerproduct' is not defined

I cannot understand why outerproduct would not be imported in basic.py.... 

Doing a "from Numeric import outerproduct, concatenate" seems to fix the problem, but I haven't tested with any real numbers or anything. concatenate is missing as well.

Dave


-- 
David J. Grant
http://www.davidandnasha.ca

-------------- next part --------------
A non-text attachment was scrubbed...
Name: david.grant.vcf
Type: text/x-vcard
Size: 180 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20050406/f033e226/attachment.vcf>

From oliphant at ee.byu.edu  Thu Apr  7 02:42:36 2005
From: oliphant at ee.byu.edu (Travis Oliphant)
Date: Thu, 07 Apr 2005 00:42:36 -0600
Subject: [SciPy-dev] kron function in scipy.linalg.basic
In-Reply-To: <42549979.1020802@telus.net>
References: <42549979.1020802@telus.net>
Message-ID: <4254D65C.70307@ee.byu.edu>

David Grant wrote:

>Re-posting to scipy-dev as it might be a bug.
>
>I'm using the kron() function in scipy.linalg.basic but I get an error
>  
>
Bugs like this are the most common, because they don't show up when you 
lazily test code you write from an interactive shell like emacs provides 
(which has all the names already imported).

Thanks for finding them.  The fix is easy (just edit the file and import 
the name that is missing).

-Travis



From david.grant at telus.net  Thu Apr  7 02:53:33 2005
From: david.grant at telus.net (David Grant)
Date: Wed, 06 Apr 2005 23:53:33 -0700
Subject: [SciPy-dev] kron function in scipy.linalg.basic
In-Reply-To: <4254D65C.70307@ee.byu.edu>
References: <42549979.1020802@telus.net> <4254D65C.70307@ee.byu.edu>
Message-ID: <4254D8ED.90602@telus.net>

Travis Oliphant wrote:

> David Grant wrote:
>
>> Re-posting to scipy-dev as it might be a bug.
>>
>> I'm using the kron() function in scipy.linalg.basic but I get an error
>>  
>>
> Bugs like this are the most common, because they don't show up when
> you lazily test code you write from an interactive shell like emacs
> provides (which has all the names already imported).
>
> Thanks for finding them.  The fix is easy (just edit the file and
> import the name that is missing).
>
>
Thanks Travis. This is the kind of thing that happens to me when using
IDLE as well, so I know what you mean.

I did so some tests and it looks like the function works as intended,
just by importing those two missing Numeric functions.

Will you do this in CVS? I don't mind doing it, but of course, I'm not a
scipy developer.

Dave

-- 
David J. Grant
http://www.davidandnasha.ca

-------------- next part --------------
A non-text attachment was scrubbed...
Name: david.grant.vcf
Type: text/x-vcard
Size: 180 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20050406/6d24ea0b/attachment.vcf>

From david.grant at telus.net  Thu Apr  7 04:59:30 2005
From: david.grant at telus.net (David Grant)
Date: Thu, 07 Apr 2005 01:59:30 -0700
Subject: [SciPy-dev] kron function in scipy.linalg.basic
In-Reply-To: <4254D65C.70307@ee.byu.edu>
References: <42549979.1020802@telus.net> <4254D65C.70307@ee.byu.edu>
Message-ID: <4254F672.7020709@telus.net>

Travis Oliphant wrote:

> David Grant wrote:
>
>> Re-posting to scipy-dev as it might be a bug.
>>
>> I'm using the kron() function in scipy.linalg.basic but I get an error
>>  
>>
> Bugs like this are the most common, because they don't show up when
> you lazily test code you write from an interactive shell like emacs
> provides (which has all the names already imported).
>
> Thanks for finding them.  The fix is easy (just edit the file and
> import the name that is missing).
>
>
Well it turned out it needed one more addition, the reshape function, so
here's the diff:


-(/usr/lib/python2.3/site-packages/scipy/linalg:$)-> diff -u
basic.py.orig basic.py

--- basic.py.orig       2005-04-07 00:52:23.000000000 -0700
+++ basic.py    2005-04-07 00:52:31.000000000 -0700
@@ -18,6 +18,8 @@
 from scipy_base import asarray_chkfinite
 import calc_lwork
 
+#added by DG (April 6, 2005)
+from Numeric import outerproduct, concatenate, reshape
 
 class LinAlgError(Exception):
     pass


-- 
David J. Grant
http://www.davidandnasha.ca

-------------- next part --------------
A non-text attachment was scrubbed...
Name: david.grant.vcf
Type: text/x-vcard
Size: 180 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20050407/9057e9af/attachment.vcf>

From david.grant at telus.net  Fri Apr  8 19:01:40 2005
From: david.grant at telus.net (David Grant)
Date: Fri, 08 Apr 2005 16:01:40 -0700
Subject: [SciPy-dev] matrixmultiply()
Message-ID: <42570D54.2090509@telus.net>

Can anyone tell me why the * operation on matrices is for element-wise
multiplcation and not matrix multiplication for Numeric/numpy arrays?

Or at least tell me a nice way to do A*B*C*D*E*F, where the * denotes
matrix multiplication. I once wrote a function where I passed those to a
function as arrary of arrays, so mult([A,B,C,D,E,F]) and the function
multiplied them all together. But I can't remember where I put this
function, and I was wondering if there was an easier way.

-- 
David J. Grant
http://www.davidandnasha.ca

-------------- next part --------------
A non-text attachment was scrubbed...
Name: david.grant.vcf
Type: text/x-vcard
Size: 180 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20050408/33618360/attachment.vcf>

From aisaac at american.edu  Fri Apr  8 19:11:53 2005
From: aisaac at american.edu (Alan G Isaac)
Date: Fri, 8 Apr 2005 19:11:53 -0400 (Eastern Daylight Time)
Subject: [SciPy-dev] matrixmultiply()
In-Reply-To: <42570D54.2090509@telus.net>
References: <42570D54.2090509@telus.net>
Message-ID: <Mahogany-0.66.0-116-20050408-191201.00@american.edu>

On Fri, 08 Apr 2005, David Grant apparently wrote: 
> Can anyone tell me why the * operation on matrices is for element-wise 
> multiplcation and not matrix multiplication for Numeric/numpy arrays? 

Because they're arrays.

> Or at least tell me a nice way to do A*B*C*D*E*F, where the * denotes 
> matrix multiplication.

Try matrices.
>>> from scipy import *
>>> help(mat)

hth,
Alan Isaac






From joe at enthought.com  Sat Apr  9 20:03:32 2005
From: joe at enthought.com (Joe Cooper)
Date: Sat, 09 Apr 2005 19:03:32 -0500
Subject: [SciPy-dev] New Enthon Test Release
Message-ID: <42586D54.3090109@enthought.com>

Hi all,

I posted a new test release of Enthon yesterday, and I think this one is 
approaching useful.  I'm running out of things I can test, in my limited 
knowledge of all of the components that it contains.

So, I'd like to ask anyone who has a vested interest in an easy to 
install bundle of Python stuff for Windows to download this release and 
give it a try.  A fresh install is probably the best way to test it, as 
I haven't figured out how to make Wise remove the extraneous menu items 
on upgrade (the uninstaller seems to leave a few sticking around 
too...will probably get that fixed before the final release).

Grab it here:

http://download.enthought.com/enthought_python-2.3.5-1069.exe

Submit bug reports to the roundup instance, or send me an email, if you 
find problems.

Thanks!



From flory at fdu.edu  Sun Apr 10 02:05:37 2005
From: flory at fdu.edu (David Flory)
Date: Sun, 10 Apr 2005 02:05:37 -0400
Subject: [SciPy-dev] New Enthon Test Release
In-Reply-To: <42586D54.3090109@enthought.com>
Message-ID: <!~!UENERkVCMDkAAQACAAAAAAAAAAAAAAAAABgAAAAAAAAA48bIjzpX0xGAzADAT0HvxiIvKwAQAAAAYV9QKFkbEkiRz2lpjO1d0QEAAAAA@fdu.edu>

 
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Is Enthon planning to do a Python 2.4.x release?  If so, when?

- -------------------------------------------
David Flory, Ph.D.
Professor of Physics 
 

> -----Original Message-----
> From: scipy-dev-bounces at scipy.net 
> [mailto:scipy-dev-bounces at scipy.net] On Behalf Of Joe Cooper
> Sent: Saturday, April 09, 2005 8:04 PM
> To: Scipy-dev at scipy.net
> Subject: [SciPy-dev] New Enthon Test Release
> 
> Hi all,
> 
> I posted a new test release of Enthon yesterday, and I think 
> this one is approaching useful.  I'm running out of things I 
> can test, in my limited knowledge of all of the components 
> that it contains.
> 
> So, I'd like to ask anyone who has a vested interest in an 
> easy to install bundle of Python stuff for Windows to 
> download this release and give it a try.  A fresh install is 
> probably the best way to test it, as I haven't figured out 
> how to make Wise remove the extraneous menu items on upgrade 
> (the uninstaller seems to leave a few sticking around 
> too...will probably get that fixed before the final release).
> 
> Grab it here:
> 
> http://download.enthought.com/enthought_python-2.3.5-1069.exe
> 
> Submit bug reports to the roundup instance, or send me an 
> email, if you find problems.
> 
> Thanks!
> 
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-dev
> 

-----BEGIN PGP SIGNATURE-----
Version: PGP 8.1

iQA/AwUBQljCMVe2/rcN3lp8EQJDIgCfVyw3U1dmlkpeLzn5RYuFUzb+nfUAn2+A
A0+eeLljx3O87l+ZvGtdSI5Q
=XWX/
-----END PGP SIGNATURE-----



From prabhu_r at users.sf.net  Tue Apr 12 14:58:07 2005
From: prabhu_r at users.sf.net (Prabhu Ramachandran)
Date: Wed, 13 Apr 2005 00:28:07 +0530
Subject: [SciPy-dev] New Enthon Test Release
In-Reply-To: <42586D54.3090109@enthought.com>
References: <42586D54.3090109@enthought.com>
Message-ID: <16988.6719.719911.214344@monster.linux.in>

Hi Joe,

>>>>> "Joe" == Joe Cooper <joe at enthought.com> writes:

    Joe> Hi all, I posted a new test release of Enthon yesterday, and
    Joe> I think this one is approaching useful.  I'm running out of
    Joe> things I can test, in my limited knowledge of all of the
    Joe> components that it contains.
[...]
    Joe> http://download.enthought.com/enthought_python-2.3.5-1069.exe

I gave it a whirl.  Install is smooth!  I have a few VTK gripes
though.

 1. On March 1st I sent you a mail on how I needed to reorganize the
    site-packages/vtk_python directory along with a directory listing.
    It does not look like this has been done for 1069.  Basically
    there are the following problems:

    a. You have two copies of the vtk*Python.dll files, one in
       vtk_python and one inside vtk_python/vtk.  You really should
       put all the dll's into the vtk_python dir and leave the
       vtk_python/vtk directory with pure Python modules.

    b. You don't need to build the Tcl libraries.  You do need to
       point VTK to the Tcl/Tk libraries in order to build the Tk
       widgets used by Tkinter but you don't need to build the Tcl
       bindings to get them built.

    I also posted some of these comments on Jan 28th.   See here
    (specifically look at point 6):

http://www.scipy.org/mailinglists/mailman?fn=scipy-dev/2005-January/002800.html


 2. It looks like you have not enabled GL2PS support for VTK.  I'd
    appreciate this very much.  To get this option do the following:
    when you run CMake, turn on the "show advanced values" and set
    VTK_USE_GL2PS on.  You can turn off the "show advanced values"
    once you do this.  Then re-configure and rebuild.  This allows
    users to save VTK scenes to vector PS/EPS/PDF files which is very
    handy.  I'd appreciate it very much if this feature could be
    incorporated in the next intermediate release of Enthon.  Then can
    I cut a 1.5 MayaVi release with a binary that supports gl2ps.


Thanks!


cheers,
prabhu



From schofield at ftw.at  Tue Apr 19 12:48:19 2005
From: schofield at ftw.at (Ed Schofield)
Date: Tue, 19 Apr 2005 18:48:19 +0200
Subject: [SciPy-dev] Logical operations, Numeric.sum() and overflow
Message-ID: <42653653.90608@ftw.at>

Hi all,

Importing 'scipy' changes the output of the following code:

>>>import random
>>>import Numeric, RandomArray
>>>(m,n) = (1000,3)
>>>x = RandomArray.random((m,n))
>>>y = x < 0.5
>>>assert sum(y) == Numeric.sum(y)

from nothing to an AssertionError.

This is random behaviour: the error occurs about 90% of the time with
this value of m on my PC (NumPy 23.1, SciPy 0.3.2, Python 2.4.1, Linux
2.6.11) , but reducing m to 500 makes it occur only about 20% of the time.

There appear to be two causes:
(1) importing scipy changes the behaviour of "y = x < 0.5" to return an
array of typecode 'b' rather than 'l'.
(2) Numeric.sum() is prone to overflow errors, returning an array of
type 'b' rather than increasing precision:

>>> a = Numeric.array([[253,254,255],[1,1,1]],'b')
>>> Numeric.sum(a)
array([254, 255,   0],'b')

Here's my two cents.  On point (1), unit tests are needed to ensure a
simple 'import scipy' can't change the behaviour of unrelated Numpy
code.  (How is this even possible?)

Point (2) seems to indicate a design flaw with Numeric.  How do Octave
and Matlab deal with this?  Whatever they do, it "just works", whereas
Numeric feels "broken" in this respect; this overflow propagates through
other operations (Numeric.average() in my case), and finding such bugs
can take hours.  Any suggestions / ideas?


-- Ed




From charles.harris at sdl.usu.edu  Tue Apr 19 13:03:00 2005
From: charles.harris at sdl.usu.edu (Charles R Harris)
Date: Tue, 19 Apr 2005 11:03:00 -0600
Subject: [SciPy-dev] Logical operations, Numeric.sum() and overflow
In-Reply-To: <42653653.90608@ftw.at>
References: <42653653.90608@ftw.at>
Message-ID: <1113930180.10073.3.camel@E011704>

On Tue, 2005-04-19 at 18:48 +0200, Ed Schofield wrote:
> Hi all,
> 
> Importing 'scipy' changes the output of the following code:
> 
> >>>import random
> >>>import Numeric, RandomArray
> >>>(m,n) = (1000,3)
> >>>x = RandomArray.random((m,n))
> >>>y = x < 0.5
> >>>assert sum(y) == Numeric.sum(y)
> 
> from nothing to an AssertionError.
> 
> This is random behaviour: the error occurs about 90% of the time with
> this value of m on my PC (NumPy 23.1, SciPy 0.3.2, Python 2.4.1, Linux
> 2.6.11) , but reducing m to 500 makes it occur only about 20% of the time.
> 
> There appear to be two causes:
> (1) importing scipy changes the behaviour of "y = x < 0.5" to return an
> array of typecode 'b' rather than 'l'.
> (2) Numeric.sum() is prone to overflow errors, returning an array of
> type 'b' rather than increasing precision:
> 
> >>> a = Numeric.array([[253,254,255],[1,1,1]],'b')
> >>> Numeric.sum(a)
> array([254, 255,   0],'b')
> 
> Here's my two cents.  On point (1), unit tests are needed to ensure a
> simple 'import scipy' can't change the behaviour of unrelated Numpy
> code.  (How is this even possible?)
> 
> Point (2) seems to indicate a design flaw with Numeric.  How do Octave
> and Matlab deal with this?  Whatever they do, it "just works", whereas
> Numeric feels "broken" in this respect; this overflow propagates through
> other operations (Numeric.average() in my case), and finding such bugs
> can take hours.  Any suggestions / ideas?
> 
> 
> -- Ed

Numarray uses the byte boolean types also and has given me similar
problems. I haven't felt the need to complain, but it is annoying. I
suppose the byte type saves a bit of space, but frankly I think it would
be better to stick to plain old integers as c does.

chuck



From joe at enthought.com  Tue Apr 19 15:49:31 2005
From: joe at enthought.com (Joe Cooper)
Date: Tue, 19 Apr 2005 14:49:31 -0500
Subject: [SciPy-dev] SciPy.org brief outages planned
Message-ID: <426560CB.2040100@enthought.com>

Hi all,

I'm preparing to migrate SciPy.org to a new server, and so I will be 
taking the SciPy.org community website offline for a few minutes a 
couple of times over the next few days.  Given the current state of the 
old system (Plone eats way more than her fair share of memory, leaving 
the box slow and easily confused at time), it's likely no one will 
notice these outages, but I thought I'd give everyone a heads up just in 
case it takes longer than I expect.

Anyway, no other services should be adversely impacted at this time, 
though there will be an hour or two of downtime when the final switch 
occurs, as I'll need to shutdown everything to do the final rsync of 
mail, web, and CVS data.  I do not have a firm schedule for when the big 
switch will occur, but it is likely to be this week.

Thanks!



From eric at enthought.com  Wed Apr 20 01:47:27 2005
From: eric at enthought.com (eric jones)
Date: Wed, 20 Apr 2005 00:47:27 -0500
Subject: [SciPy-dev] job openings at Enthought
Message-ID: <4265ECEF.6050004@enthought.com>

Hey group,

We have a number of scientific/python related jobs open.  If you have 
any interest, please see:

    http://www.enthought.com/careers.htm

thanks,
eric



From venkatraghavan at gmail.com  Wed Apr 20 05:34:01 2005
From: venkatraghavan at gmail.com (Venkat Raghavan V.C.)
Date: Wed, 20 Apr 2005 04:34:01 -0500
Subject: [SciPy-dev] Control Systems Module
Message-ID: <f4956e0205042002342909a121@mail.gmail.com>

Hi all,

I am new to scipy and python. I have some experience in programming in
C, and compiling stuff. I do have good control systems knowledge,
including the math!

I would like to write a Control Systems module for scipy, similiar to
the one in Matlab. I have done quite a bit of programming in Matlab,
though I don't think that's going to help me out in this regard,
except maybe regarding design.

The problem is, how do I get started off? I have about 3 months of
time during summer which I intend to have finished this off. I would
be grateful if someone gives me a general overview on how to go about
it.

I have found this netlib library,
http://www.netlib.org/ieeecss/cascade/ after minimal searching.
I also know of the scipy api docs. I just wanted to build a roadmap
before I get started off!

Regards,
Venkat



From rkern at ucsd.edu  Wed Apr 20 06:52:54 2005
From: rkern at ucsd.edu (Robert Kern)
Date: Wed, 20 Apr 2005 03:52:54 -0700
Subject: [SciPy-dev] Control Systems Module
In-Reply-To: <f4956e0205042002342909a121@mail.gmail.com>
References: <f4956e0205042002342909a121@mail.gmail.com>
Message-ID: <42663486.8000308@ucsd.edu>

Venkat Raghavan V.C. wrote:
> Hi all,
> 
> I am new to scipy and python. I have some experience in programming in
> C, and compiling stuff. I do have good control systems knowledge,
> including the math!
> 
> I would like to write a Control Systems module for scipy, similiar to
> the one in Matlab. I have done quite a bit of programming in Matlab,
> though I don't think that's going to help me out in this regard,
> except maybe regarding design.

Fantastic!

> The problem is, how do I get started off? I have about 3 months of
> time during summer which I intend to have finished this off. I would
> be grateful if someone gives me a general overview on how to go about
> it.

Not knowing anything specific about control theory codes, I can really 
only give some general advice. And in fact, most of the advice I can 
give is to use f2py wherever possible! I think your code would probably 
go into its own scipy.control subpackage. Other than that, look around 
scipy for analagues to mimic.

> I have found this netlib library,
> http://www.netlib.org/ieeecss/cascade/ after minimal searching.
> I also know of the scipy api docs. I just wanted to build a roadmap
> before I get started off!

-- 
Robert Kern
rkern at ucsd.edu

"In the fields of hell where the grass grows high
  Are the graves of dreams allowed to die."
   -- Richard Harter



From nwagner at mecha.uni-stuttgart.de  Wed Apr 20 11:02:47 2005
From: nwagner at mecha.uni-stuttgart.de (Nils Wagner)
Date: Wed, 20 Apr 2005 17:02:47 +0200
Subject: [SciPy-dev] An addition to basic.py
Message-ID: <42666F17.1080405@mecha.uni-stuttgart.de>

Hi all,

it would be nice if the following function could be added to basic.py in
/Lib/linalg

from scipy import *

def vec(a):
    """vec function of the matrix a is the vector formed by stacking
       the columns of a on top of one another
       Reference: Alan J. Laub, Matrix Analysis for Scientists & Engineers
       SIAM 2005
    """
    if not a.iscontiguous():
        a = reshape(a, a.shape)

    return concatenate(transpose(a))
#
# A short test
#
a = rand(3,3)
b = rand(3,3)
c = rand(3,3)
d = dot(a,dot(b,c))
x1 = vec(d)
x2 = dot(linalg.kron(transpose(c),a),vec(b))
print 'x1-x2 should be zero',linalg.norm(x1-x2)

Nils



From aisaac at american.edu  Wed Apr 20 13:21:40 2005
From: aisaac at american.edu (Alan G Isaac)
Date: Wed, 20 Apr 2005 13:21:40 -0400 (Eastern Daylight Time)
Subject: [SciPy-dev] An addition to basic.py
In-Reply-To: <42666F17.1080405@mecha.uni-stuttgart.de>
References: <42666F17.1080405@mecha.uni-stuttgart.de>
Message-ID: <Mahogany-0.66.0-944-20050420-132145.00@american.edu>

1. I prefer the solution you can find in pyGAUSS:
def vec(x): return transpose([ravel(transpose(x))])
def vecr(x): return transpose([ravel(x)])
In particular, I prefer to get back a 2-D array.
2. I suppose someone familiar with tensors should
address the proper generalization beyond 2-D inputs.

fwiw,
Alan Isaac






From nwagner at mecha.uni-stuttgart.de  Wed Apr 20 13:53:05 2005
From: nwagner at mecha.uni-stuttgart.de (Nils Wagner)
Date: Wed, 20 Apr 2005 19:53:05 +0200
Subject: [SciPy-dev] An addition to basic.py
In-Reply-To: <Mahogany-0.66.0-944-20050420-132145.00@american.edu>
References: <42666F17.1080405@mecha.uni-stuttgart.de>
	<Mahogany-0.66.0-944-20050420-132145.00@american.edu>
Message-ID: <web-19396980@uni-stuttgart.de>

On Wed, 20 Apr 2005 13:21:40 -0400 (Eastern Daylight Time)
  Alan G Isaac <aisaac at american.edu> wrote:
> 1. I prefer the solution you can find in pyGAUSS:
> def vec(x): return transpose([ravel(transpose(x))])
> def vecr(x): return transpose([ravel(x)])
> In particular, I prefer to get back a 2-D array.
> 2. I suppose someone familiar with tensors should
> address the proper generalization beyond 2-D inputs.
> 
> fwiw,
> Alan Isaac
> 
> 
Just a suggestion. And how about the reverse ?

a = rand(3,3)
b = vec(a)
a = transpose(reshape(a,(3,3))) # not the best way

Nils

> 
> 
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-dev

  



From nwagner at mecha.uni-stuttgart.de  Thu Apr 21 07:53:57 2005
From: nwagner at mecha.uni-stuttgart.de (Nils Wagner)
Date: Thu, 21 Apr 2005 13:53:57 +0200
Subject: [SciPy-dev] An addition to basic.py
In-Reply-To: <Mahogany-0.66.0-944-20050420-132145.00@american.edu>
References: <42666F17.1080405@mecha.uni-stuttgart.de>
	<Mahogany-0.66.0-944-20050420-132145.00@american.edu>
Message-ID: <42679455.4040000@mecha.uni-stuttgart.de>

Alan G Isaac wrote:

>1. I prefer the solution you can find in pyGAUSS:
>def vec(x): return transpose([ravel(transpose(x))])
>def vecr(x): return transpose([ravel(x)])
>In particular, I prefer to get back a 2-D array.
>2. I suppose someone familiar with tensors should
>address the proper generalization beyond 2-D inputs.
>
>fwiw,
>Alan Isaac
>
>
>
>
>_______________________________________________
>Scipy-dev mailing list
>Scipy-dev at scipy.net
>http://www.scipy.net/mailman/listinfo/scipy-dev
>  
>


Please can you add these functions to basic.py.

Thanks in advance
                                           Nils

 




From aisaac at american.edu  Thu Apr 21 10:12:14 2005
From: aisaac at american.edu (Alan Isaac)
Date: Thu, 21 Apr 2005 10:12:14 -0400 (Eastern Standard Time)
Subject: [SciPy-dev] An addition to basic.py
In-Reply-To: <42679455.4040000@mecha.uni-stuttgart.de>
References: <42666F17.1080405@mecha.uni-stuttgart.de>
	<Mahogany-0.66.0-944-20050420-132145.00@american.edu><42679455.4040000@mecha.uni-stuttgart.de>
Message-ID: <Mahogany-0.66.0-2468-20050421-124950.00@CASC64043305.american.edu>

> Alan G Isaac wrote: 
>> 1. I prefer the solution you can find in pyGAUSS:
>> def vec(x): return transpose([ravel(transpose(x))])
>> def vecr(x): return transpose([ravel(x)])
>> In particular, I prefer to get back a 2-D array.
>> 2. I suppose someone familiar with tensors should
>> address the proper generalization beyond 2-D inputs.


On Thu, 21 Apr 2005, Nils Wagner wrote:
> Please can you add these functions to basic.py.


I apologize if I gave the impression that I am a SciPy developer.
But pyGAUSS is under the MIT license, so this (trivial) code is
unencumbered if you can persuade a developer to include it.

I agree that vec and vecr are basic functionality.
But I do not know how they should generalize beyond 2-D
inputs.

Cheers,
Alan Isaac





From oliphant at ee.byu.edu  Fri Apr 22 07:03:06 2005
From: oliphant at ee.byu.edu (Travis Oliphant)
Date: Fri, 22 Apr 2005 05:03:06 -0600
Subject: [SciPy-dev] Re: [SciPy-user] Vector mean calculation?
In-Reply-To: <89f3084c0a6f6705488057af7d9b987e@cwazy.co.uk>
References: <4268475A.6020803@haas.berkeley.edu>
	<1114130503.426848471de1d@webmail.colorado.edu>
	<89f3084c0a6f6705488057af7d9b987e@cwazy.co.uk>
Message-ID: <4268D9EA.60309@ee.byu.edu>

Lance Boyle wrote:

>
> On Apr 21, 2005, at 5:41 PM, Fernando.Perez at colorado.edu wrote:
>
>> Quoting Thomas Davidoff <davidoff at haas.berkeley.edu>:
>>
>>> I have a vector with min equal to 412.  Why am I getting a negative 
>>> mean?
>>
>>
>> Without more details, it's quite difficult to be sure.  Most likely, 
>> because
>> it's very long and of integer type (ints are 32 bit objects, so you are
>> effectively doing mod 2^31 arithmetic with negative wraparound).
>>
> Why does Python allow modular arithmetic when the user doesn't request 
> modular arithmetic?


To get speed, scipy must assume that if you are using an integer array 
then you actually are requesting modular arithmetic (to the size of the 
integers in question).  The other option would be to silently and 
automatically upcast when you get too large. 

I don't think that is the right way to go.  It would be a good idea, 
however, to allow the user to specify a larger type within which to 
accomplish a reduction.   Look for this in scipy.base (the new Numeric).

I'm thinking of something like

add.reduce(a, Int32) or add.reduce(a, Float64) 

which gives a return type for the array.



-Travis




From venkatraghavan at gmail.com  Sat Apr 23 20:12:58 2005
From: venkatraghavan at gmail.com (Venkat Raghavan)
Date: Sat, 23 Apr 2005 19:12:58 -0500
Subject: [SciPy-dev] Control Systems Module, Issues
In-Reply-To: <20050420175346.710023EBD5@www.scipy.com>
References: <20050420175346.710023EBD5@www.scipy.com>
Message-ID: <e08f11e188fea1ca7e70febb6a31b2e7@gmail.com>

Hi

Is it okay to implement the functions, name-wise, as it is done in 
Matlab. For example, say the series function in Matlab is implemented 
as : series(first param, second param). Is it okay to use the same 
function names?

I think it's okay, but as I have never done this before, I am not sure. 
Also, this has the added advantage of easy portability of code across 
from Matlab, and lesser learning curve (if any).

I've got my books out, and hope to have some results soon. The whole 
control systems module seems 'very doable' right now :)

Regards,
Venkat 



From prasansamtani at gmail.com  Mon Apr 25 21:03:56 2005
From: prasansamtani at gmail.com (Prasan Samtani)
Date: Mon, 25 Apr 2005 18:03:56 -0700
Subject: [SciPy-dev] Fuzzy logic library for SciPy
Message-ID: <5da497b705042518035e185dc3@mail.gmail.com>

Greetings scipy-devers,

I am currently beginning work on a Fuzzy Logic toolkit in Python
(modeled after the MatLab fuzzy logic toolbox
http://www.mathworks.com/products/fuzzylogic/). Before I get ahead of
myself, I was wondering if someone has already been working on
something like this, I want to make sure I'm not reinventing the
wheel.

Thanks,
Prasan



From schofield at ftw.at  Tue Apr 26 13:54:39 2005
From: schofield at ftw.at (Ed Schofield)
Date: Tue, 26 Apr 2005 19:54:39 +0200
Subject: [SciPy-dev] Bug fix in optimize/lbfgsb.py
Message-ID: <426E805F.9060300@ftw.at>

Hi all,

This patch fixes a bug in SciPy 0.3.2 when passing user-defined
arguments to the LBFGSB minimization routine.

File: optimize/lbfgsb.py

179c179
<             f[0], g = func_and_grad(x, *args)
---
>             f[0], g = func_and_grad(x)

This bug is outstanding from version 0.3 (see
http://www.scipy.net/pipermail/scipy-user/2004-June/002935.html).

Best wishes,
Ed



From cookedm at physics.mcmaster.ca  Tue Apr 26 14:05:34 2005
From: cookedm at physics.mcmaster.ca (David M. Cooke)
Date: Tue, 26 Apr 2005 14:05:34 -0400
Subject: [SciPy-dev] Bug fix in optimize/lbfgsb.py
In-Reply-To: <426E805F.9060300@ftw.at> (Ed Schofield's message of "Tue, 26
	Apr 2005 19:54:39 +0200")
References: <426E805F.9060300@ftw.at>
Message-ID: <qnk64y9za4x.fsf@arbutus.physics.mcmaster.ca>

Ed Schofield <schofield at ftw.at> writes:

> Hi all,
>
> This patch fixes a bug in SciPy 0.3.2 when passing user-defined
> arguments to the LBFGSB minimization routine.
>
> File: optimize/lbfgsb.py
>
> 179c179
> <             f[0], g = func_and_grad(x, *args)
> ---
>>             f[0], g = func_and_grad(x)
>
> This bug is outstanding from version 0.3 (see
> http://www.scipy.net/pipermail/scipy-user/2004-June/002935.html).

That's been fixed in CVS since the end of November.

-- 
|>|\/|<
/--------------------------------------------------------------------------\
|David M. Cooke                      http://arbutus.physics.mcmaster.ca/dmc/
|cookedm at physics.mcmaster.ca



From t.zito at biologie.hu-berlin.de  Wed Apr 27 06:41:59 2005
From: t.zito at biologie.hu-berlin.de (Tiziano Zito)
Date: Wed, 27 Apr 2005 12:41:59 +0200
Subject: [SciPy-dev] EuroPython conference 2005
Message-ID: <20050427104159.GA23626@itb.biologie.hu-berlin.de>

Is someone of the SciPy developers coming to the EuroPython
conference this year? I think that Numeric3 and SciPy deserve
a talk in the Science track.

The announce on python-announce list follows:

EuroPython is held -for the second time- in G?teborg, Sweden, during
the week of June 27-29, 2005.  Registration is open and now is your
chance to submit a proposal (or more).

Several topics are handled in parallel tracks at EuroPython2005, these
include: Business, Education, Lightning talks, Python Frameworks,
Python Language, Refereed papers, Science, Social skills/General
topics, Tutorials/Neopythes, Zope and Zope Related, Zope Lightning
talks. If you would like to tell the community something about your
work, interest or experience, then then submit button is your way to
go.

The Refereed Track is intended for prestigious peer-reviewed talks for
presenting technical and scientific papers. To our knowledge,
EuroPython is the only Python conference, with peer-reviewed
talks. Everybody can submit a paper for the Refereed Track, as long as
it follows the Refereed paper instructions. For detailed information
concerning this subject, please have a look at the Tracks & Talks
webpage of EuroPython. To give you yet another opportunity to submit a
refereed paper, we have extended the deadline to coincide with the
deadline for regular talks. Lat day for submission is 1 May 2005!

If you don't want to participate as a speaker, you might be interested
in hearing the talks. Registration for EuroPython2005 is now
open. Several registration alternatives are possible. Pre-Registration
will be available until 17th of June, but Early Bird registration will
end on May the 15th. Everybody is welcome: from novice to advanced.

Our keynote speakers will be Guido van Rossum and Steve Pemberton.
You may not have heard of Steve, since his focus is in an area that is
not directly connected to Python. However, we expect him to give you
an interesting perspective on what is going on in the industry.
Find out more at http://homepages.cwi.nl/~steven/.
If you don't know who Guido van Rossum is, you really need to come to
Europython to find out.

Feel free to spread the word about EuroPython2005 to your friends, family,
colleagues, ...

For all conference details, go to:
http://www.europython.org

See you at G?teborg.

EuroPython Team.

About EuroPython: Europython is an annual conference for the Python
and Zope communities. It circulates between different sites in
Europe. Having started in Charleroi, Belgium, it is now in G?teborg,
Sweden and will move to CERN in Swizerland next year. Europython is a
community conference run by volunteers.