Mailman 3 August 2004 - SciPy-User

FAIL: check_lpmn (scipy.special.basic.test_basic.test_lpmn)
by Nils Wagner Aug. 6, 2004

Aug. 6, 2004

scipy.test() with latest cvs ====================================================================== FAIL: check_lpmn (scipy.special.basic.test_basic.test_lpmn) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.3/site-packages/scipy/special/tests/test_basic.py", line 1625, in check_lpmn 1.50000]])),4) File "/usr/lib/python2.3/site-packages/scipy_test/testing.py", line 681, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 66.6666666667%): Array 1: [[ [ 1. 0. 0.]] [ [ 0. 0. 0.]]] Array 2: [[ [ 1. 0.5 -0.125]] [ [ 0. 1. 1.5 ]]] ---------------------------------------------------------------------- Ran 936 tests in 2.826s FAILED (failures=1)

1 0

vectorize question
by bulatov＠cs.orst.edu Aug. 6, 2004

Aug. 6, 2004

What's the best way to evaluate a function at every point of a space (ie, for plotting purposes)? For instance, here's how I'd do it in regular python values=zeros((5,5), Float) for i in range(5): for j in range(5): values[i,j]=some_custom_function(i,j) plot(values) Is there some SciPy way of doing it more concisely? I looked at vectorize, but I can only use it for similar task in 1 dimension. Yaroslav

2 1

RE: [SciPy-user] Re: plot_func
by Matt Kubilus Aug. 6, 2004

Aug. 6, 2004

Looking at the gplt modules, it appears that the necessary function to send commands directly is a hidden function. I ended up needing to use scipy.gplt with time series data that I was logging. There doesn't appear to be a way to do this (Gnuplot can, scipy.gplt makes this not possible) so anyways, I hacked into the pyPlot.py module and added a function that essentially just sends a command to the hidden function that actually calls the plotting routines. This allowed me to use the full functionality of gnuplot which can to way more than just what scipy.gplt does. If anyone is interested, I can post the actual code changes. --Matt Kubilus >From: Alan G Isaac <aisaac(a)american.edu> >Reply-To: aisaac(a)american.edu,SciPy Users List <scipy-user(a)scipy.net> >To: aisaac(a)american.edu, scipy-user(a)scipy.net >CC: scipy-dev(a)scipy.net >Subject: [SciPy-user] Re: plot_func >Date: Tue, 3 Aug 2004 12:29:40 -0400 (Eastern Daylight Time) > >On Sun, 1 Aug 2004 Alan G Isaac apparently wrote: > > I noticed that pyPlot.py has: > > def plot_func(self,func): > > self._init_plot() > > if(type(func) == StringType): > > self._send('plot '+ func) > > else: > > print 'error: requires a function as a string' > > but gplt.plot_func is not available. > > > > >OK, I've confirmed that if I make an instance of >scipy.gplt.pyPlot.Plot >then the plot_func method works fine for this instance. > >I don't know enough about Python to guess why gplt is >not giving us access to this. Can anyone help? > >Thank you, >Alan Isaac > > > >_______________________________________________ >SciPy-user mailing list >SciPy-user(a)scipy.net >http://www.scipy.net/mailman/listinfo/scipy-user _________________________________________________________________ Discover the best of the best at MSN Luxury Living. http://lexus.msn.com/

2 1

RE: [SciPy-user] pyPlot.py vs. Gnuplot.py
by Matt Kubilus Aug. 5, 2004

Aug. 5, 2004

Another new user commenting. I found that pyPlot was less reliable than gnuplot and had less features. For instance, with the wx based pyPlot, once I closed a window with a plot in my wx based application, I could no longer open the graph. The C object has already been destroyed was the error in a nut shell. Restarting the application was the only workaround. Has anyone else ran into, fixed a similar problem. Gnuplot is also able to do interesting things like 3d plots (not necessary for what I use it for) I use this in a manufacturing/data gathering setting. --Matt Kubilus >From: Alan G Isaac <aisaac(a)american.edu> >Reply-To: aisaac(a)american.edu,SciPy Users List <scipy-user(a)scipy.net> >To: gnuplot-py-users(a)lists.sourceforge.net,scipy-user ><scipy-user(a)scipy.net> >Subject: [SciPy-user] pyPlot.py vs. Gnuplot.py >Date: Tue, 3 Aug 2004 16:23:12 -0400 (Eastern Daylight Time) > >As a new SciPy user, I am interested in understanding >the relationship between the scipy.gplt module and the >Gnuplot.py project. > >What considerations influence which was chosen, and >which a new user should choose? > >This is especially important to me as I am planning >to use one of these in a teaching setting. > >Thank you, >Alan Isaac > > >_______________________________________________ >SciPy-user mailing list >SciPy-user(a)scipy.net >http://www.scipy.net/mailman/listinfo/scipy-user _________________________________________________________________ Discover the best of the best at MSN Luxury Living. http://lexus.msn.com/

2 1

Re: [SciPy-user] Kronecker product kron(a,b)
by Alan G Isaac Aug. 5, 2004

Aug. 5, 2004

On Wed Sep 24 16:52:28 2003 Nils Wagner scipy-user(a)scipy.net wrote > There is another useful function which might be interesting for scipy, > the so-called Kronecker product of two matrices. > def kron(a,b): > if not a.iscontiguous(): > a = reshape(a, a.shape) > if not b.iscontiguous(): > b = reshape(b, b.shape) > o = outerproduct(a,b) > o.shape = a.shape + b.shape > return concatenate(concatenate(o, axis=1), axis=1) > It would be great, if this becomes a standard function in scipy. I agree. It the meantime, thanks for posting this. Apparently the numarray folks are thinking of including a kroneckerproduct function. Cheers, Alan Isaac

1 0

Building Debug extrensions with scipy_distutils on WinXP
by Berthold Höllmann Aug. 5, 2004

Aug. 5, 2004

Hello, I'm trying to build a Windows Debug version of an inhouse extension module whose setup.py requires scipy_distutils. I call python_d.exe setup.py config_fc --fcompiler=compaqv build --debug but get ... creating build\temp.win32-2.3\lib\PyBrushMod\flib compile options: '-g -c' DF:f77: lib\PyBrushMod\flib\area3.f DF: error: Switch '-g' is ambiguous error: Command "DF /f77rtl /fixed /nologo /MD /WX /iface=(cref,nomixed_str_len_arg) /names:lowercase /assume:underscore /debug /Ox /fast /optimize:5 /unroll:0 /math_library:fast /threads -g -c /compil e_only lib\PyBrushMod\flib\area3.f /object:build\temp.win32-2.3\lib\PyBrushMod\flib\area3.o" failed with exit status 1 [31863 refs] The strange thing for me is, when I look into 'compaqfcompiler.py' I find the class 'CompaqVisualFCompiler' which defines def get_flags_debug(self): return ['/debug'] but for some strange reason '-g' is inserted. Kind regards Berthold Höllmann -- Germanischer Lloyd AG CAE Development Vorsetzen 35 20459 Hamburg Phone: +49(0)40 36149-7374 Fax: +49(0)40 36149-7320 e-mail: hoel(a)gl-group.com Internet: http://www.gl-group.com This e-mail contains confidential information for the exclusive attention of the intended addressee. Any access of third parties to this e-mail is unauthorised. Any use of this e-mail by unintended recipients such as copying, distribution, disclosure etc. is prohibited and may be unlawful. When addressed to our clients the content of this e-mail is subject to the General Terms and Conditions of GL's Group of Companies applicable at the date of this e-mail. GL's Group of Companies does not warrant and/or guarantee that this message at the moment of receipt is authentic, correct and its communication free of errors, interruption etc.

2 2

Constrained Optimizers (multivariate) in scipy
by Nils Wagner Aug. 5, 2004

Aug. 5, 2004

Dear all, I am looking for a small test example using one of the constrained optimizers currently available in scipy 1. fmin_l_bfgs_b 2. fmin_tnc 3. fmin_cobyla Any pointer would be appreciated. Thanks in advance. Nils

1 0

Bug in optimize.newton
by Martin Wiechert Aug. 4, 2004

Aug. 4, 2004

Hi, I've been bitten by a bug in the snippet from optimize.newton below (yesterday's CVS). I think, the actual bug is in line 317, which should read if abs(p-p1) < tol: instead of if abs(p-p0) < tol: This will typically only cause one superfluous iteration, which is no problem. But consider the following situation: alter_numeric () has been called. So, if one passes, say, a Float32 0d array as initial value, this will coerce all other scalars into Float32 0d arrays. So if you happen to get identical p's for two consecutive iterations, this will go unnoticed, because the ZeroDivisionError is not raised (rather the next p is nan), and line 317 doesn't compare the right p's. else: # Secant method p0 = x0 p1 = x0*(1+1e-4) q0 = apply(func,(p0,)+args) q1 = apply(func,(p1,)+args) for iter in range(maxiter): try: p = p1 - q1*(p1-p0)/(q1-q0) except ZeroDivisionError: if p1 != p0: print "Tolerance of %g reached" % (p1-p0) return (p1+p0)/2.0 if abs(p-p0) < tol: return p p0 = p1 q0 = q1 p1 = p q1 = apply(func,(p1,)+args) Regards, Martin Wiechert.

1 0

Re[2]: [SciPy-user] gplt
by Alan G Isaac Aug. 3, 2004

Aug. 3, 2004

On Tue, 03 Aug 2004, someone wrote: > Is there also a way to print the names of my > datalists instead of the strings "/tmp/(a)5745.4" as well as to > differentiate the colors of the lines that appear using your gplt.plot > function so I can tell which line corresponds to which data set without > the sole use of the markers? d=arange(-1.,1.,0.01) gplt.plot(d,d**2,"title 'test' with lines lt 3",d,d**2,"notitle with points lt 1") Color is controlled by the linetype, unfortunately. You'll have to see the gnuplot documentation. Cheers, Alan Isaac

1 0

plot_func
by Alan G Isaac Aug. 3, 2004

Aug. 3, 2004

I noticed that pyPlot.py has: def plot_func(self,func): self._init_plot() if(type(func) == StringType): self._send('plot '+ func) else: print 'error: requires a function as a string' but gplt.plot_func is not available. Why? I want to use this as part of gplt! How should I proceed? Thanks, Alan Isaac

1 1