Mailman 3 August 2003 - SciPy-User

Couple of Matrix() questions
by Jeffrey B. Layton 23 Sep '03

23 Sep '03

Good morning! I've got a couple of quick questions regarding the matrix() class (mat() class?). Anyway, I want to erase a column or a row of a two-dimensional matrix. I create the matrix using, >>> A=mat('[0,3,2];[5,1,2];[9,-6,7]') and I want to blank the second column using something like the following (complete with error message), >>> A[:,1]=[] Traceback (most recent call last): File "<stdin>", line 1, in ? File "/usr/lib/python2.2/site-packages/Numeric/Matrix.py", line 184, in __setitem__ self.array[index] = value How do I blank (or set to zero) the appropriate slice of the matrix? My second question is how to concatenate matrices create with mat(). I tried the following (complete with error messages), >>> A=mat('[0,1];[2,3]') >>> B=mat('[4,5];[6,7]') >>> C=mat('[8,9];[10,11]') >>> D=mat('[12,13];[14,15]') >>> E=mat('[A,B];[C,D]') Traceback (most recent call last): File "<stdin>", line 1, in ? File "/usr/lib/python2.2/site-packages/Numeric/Matrix.py", line 69, in __init__ data = _convert_from_string(data) File "/usr/lib/python2.2/site-packages/Numeric/Matrix.py", line 36, in _convert_from_string newrow.extend(map(_eval,temp)) File "/usr/lib/python2.2/site-packages/Numeric/Matrix.py", line 24, in _eval return eval(astr.translate(_table,_todelete)) File "<string>", line 0 ^ SyntaxError: unexpected EOF while parsing >>> E=[A,B; C,D] File "<stdin>", line 1 E=[A,B; C,D] I know this is a holdover from Matlab, but I would like to get is the following, E=[0, 1, 4 ,5; 2, 3, 6, 7; 8, 9, 12, 13; 10, 11, 14, 15] How do I do this in scipy? TIA! Jeff

6 13

user contributed documentation and code
by Andrew Straw 31 Aug '03

31 Aug '03

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I would like for us, as a community, to be able contribute documentation and sample code to the scipy project. I think a wiki would be a good, low-energy starting point, but unfortunately the link from the scipy "Documentation and Tutorials" page is broken. Could we get one of the Enthought folks to fix this? (If you're evaluating wikis, I recently found setting up a phpWiki ridiculously easy, even though I knew nothing about php. And yes, I know it's not in our favorite language, but I got it working within 30 minutes. :) There is also mention of user webpages and annotation of scipy.org pages if you become a scipy site member. However, I just tried to test-drive this functionality and also failed... It could be my browser -- OS X Safari is not the most compatible browser in the world, but that mainly seems to do with JavaScript support, which seems unlikely to be the problem. What do people think? Cheers! Andrew -----BEGIN PGP SIGNATURE----- iD8DBQE/UCH61xWcCSPVbpgRAthiAJ9eJ9M3lGTlP9MtnQMIT+cIFNRhDQCfS4CC 0eS5GQBJfTiLuMGF3FefErI= =uHan -----END PGP SIGNATURE-----

2 1

RE: [SciPy-user] tutorial problems
by William Griffin 31 Aug '03

31 Aug '03

I have had the same problem with tkplt in WinXP and Win2000. I can get chaco to work if I use PyCrust, and it works very well if I use it interactively, but I've not been able to get it to work reliably in module form. I can usually correct the problem with chaco but with tkplt, I've not been successful. The next time I work with either, I'll post the error messages. I figured that I was the only one not able to get it workly properly! bill -----Original Message----- From: Gary Pajer [mailto:pajer@iname.com] Sent: Sat 8/30/2003 4:38 PM To: scipy-user(a)scipy.net Cc: Subject: Re: [SciPy-user] tutorial problems It's a problem with something. It's been brought up here before, but it hasn't generated any response. Since things work fine using PyShell, perhaps it's a low priority issue. Another thing that doesn't work right is chaco.tkplt. It fails the same way on both of my systems: WinXP and WinNT I have to take the attitude that since I can't code C, I can't complain. Someday I'm sure all will be fixed. > HI, > I'm new to scipy and seem to have run into a problem with the plotting > tutorial available from the scipy site. When I follow the directions for > gplt, every thing seems to work ok, but as I get to the plt module, > python crashes. Using a wx based shell like pyshell works great every > time. When I'm using the dos shell or pythonwin it seem to crash at the > line plt.plot ((1,2,3)). > > import gui_thread > from scipy import plt > plt.plot((1,2,3)) > > Crash as soon as plt.plot((1,2,3)) is passed to the interpiter. Could > this be a problem with the gui_thread or something else? > > Can anyone help with this? > > Thanks > > _______________________________________________ > SciPy-user mailing list > SciPy-user(a)scipy.net > http://www.scipy.net/mailman/listinfo/scipy-user <http://www.scipy.net/mailman/listinfo/scipy-user> _______________________________________________ SciPy-user mailing list SciPy-user(a)scipy.net http://www.scipy.net/mailman/listinfo/scipy-user <http://www.scipy.net/mailman/listinfo/scipy-user>

1 0

tutorial problems
by SGD 30 Aug '03

30 Aug '03

HI, I'm new to scipy and seem to have run into a problem with the plotting tutorial available from the scipy site. When I follow the directions for gplt, every thing seems to work ok, but as I get to the plt module, python crashes. Using a wx based shell like pyshell works great every time. When I'm using the dos shell or pythonwin it seem to crash at the line plt.plot ((1,2,3)). import gui_thread from scipy import plt plt.plot((1,2,3)) Crash as soon as plt.plot((1,2,3)) is passed to the interpiter. Could this be a problem with the gui_thread or something else? Can anyone help with this? Thanks

2 1

ANN: PAIDA-3.0a1
by K.KISHIMOTO 29 Aug '03

29 Aug '03

Announcing PAIDA-3.0a1 (No GUI) http://paida.sourceforge.net/ PIDA is a Python implementation of AIDA (Abstract Interface for Data Analysis). I'm NOT a member of the AIDA project, but I have developed PAIDA for my use. The goals of the AIDA (http://aida.freehep.org/) project are to define abstract interfaces for common physics analysis objects, such as - histogram - ntuple - cloud - profile - fitter - IO etc. PAIDA uses SciPy in its core, and is planning to use Chaco or another technology for GUIs. Let me know of any bugs and suggestions. Have fun! K. KISHIMOTO korry(a)users.sourceforge.net

2 1

SciPy.org upgrades complete
by Joe Cooper 28 Aug '03

28 Aug '03

Hi all, The work on scipy.org is complete, and everything should be back up and running normally. If you spot problems with mail, Zope, or the roundup bug tracker, let me know. Thanks!

1 0

scipy.org mail server going down for upgrades
by Joe Cooper 28 Aug '03

28 Aug '03

Hi all, The scipy mail server will be taken down in just a few minutes for a system upgrade. Service should resume in under an hour, assuming no problems. Thanks!

1 0

xml.dom.minidom broken
by Robert Ferrell 28 Aug '03

28 Aug '03

The distribution of xml.dom.minidom in SciPy seems to be broken. I think this is not unique to SciPy, but rather common to all the Python distributions. The solution I'm using is to build PyXML myself. Has any one else encountered this problem? thanks, -robert

1 0

wxplt histogram + normal fit function (histfit)
by Sajec, Mike TQO 27 Aug '03

27 Aug '03

I wrote this function to create histograms with normal fits, and thought it may be useful to other scipy users. It requires matplotlib for the normal pdf function (http://matplotlib.sourceforge.net/). Here it is: Best Regards, Mike #*************************************************************************** **** # Module: tqplots.py # Author: Mike Sajec # Date: 8-13-03 # # Change Log: # 8-13-03 - added first function: histfit #*************************************************************************** **** from matplotlib import mlab from scipy import * from chaco import * from chaco import wxplt as plt import RandomArray as ra def histfit(data): """histfit(data) data = a 1-dimensional list or numpy array of numbers""" #data for testing... #data = ra.normal(0,3,100) mu = mlab.mean(data) sigma = mlab.std(data) mhist = mlab.hist(data) xhist = mhist[1] yhist = mhist[0] #put xhist & yhist into a list of tuple pairs histvalues=[] for ii in range(len(xhist)): histvalues.append((xhist[ii],yhist[ii])) histplot = plt.PlotValue(histvalues, type='bar') f1 = plt.figure() f1.canvas.add(histplot) plt.hold() # Plot the normal fit curve xnormfit = arange(mu-4*sigma,mu+4*sigma,(8*sigma/100)) ynormfit = mlab.normpdf(xnormfit,mu,sigma) # ynormfit = ynormfit/max(ynormfit)*max(yhist) #normalize #put xnormfit & ynormfit into a list of tuple pairs fitvalues=[] for ii in range(len(xnormfit)): fitvalues.append((xnormfit[ii],ynormfit[ii])) fitplot = plt.PlotValue(fitvalues,type='line') f1.canvas.add(fitplot) #Plot the limits #plt.plot((mu-4*sigma)*ones(2),[0,max(yhist)]) # plot lower spec limit #plt.plot((mu+4*sigma)*ones(2),[0,max(yhist)]) # plot upper spec limit return f1

1 0

Re: [SciPy-user] About SciPy '03 Workshop.
by Eric Jones 26 Aug '03

26 Aug '03

Hey Andrea, Sounds like a good idea. How about this: I will try and round up as many presentations as possible at the conference if someone else will handle putting them on the web at scipy.org. Any takers? eric Andrea Riciputi <ariciputi(a)pito.com> wrote .. > Hi there, > I've read Scipy '03 workshop announce and I'm very interested in the > topics that will be discussed. On the other side I live in Italy > (Europe) and I can't come to Pasadena to attend to the workshop. I > think many other Python/SciPy users have the same problem and I > wondering if it would be possible to put online something like a > workshop's proceeding, obviously after the workshop's end. > > Perhaps preparing a "real" proceding would turn in a time-comsuming > task, but, in my mind, just a web page to download all the presetations > (PDF form, please) would be a great help for users like me and it > wouldn't be to much time-consuming for the workshop's organization team. > > Any comment? > > Thanks, > Andrea. > > --- > Andrea Riciputi > > "Science is like sex: sometimes something useful comes out, > but that is not the reason we are doing it" -- (Richard Feynman) > > _______________________________________________ > SciPy-user mailing list > SciPy-user(a)scipy.net > http://www.scipy.net/mailman/listinfo/scipy-user

1 0