Mailman 3 February 2006 - SciPy-User

TypeError: unsupported operand type(s) for +: 'slice' and 'int'
by Nils Wagner Feb. 21, 2006

Feb. 21, 2006

Hi all, The operation (see below for details) works fine with dense matrices but not with sparse matrices. shape(V) (317, 60) <type 'instance'> shape(R0), type(R0) (317, 10) <type 'instance'> Traceback (most recent call last): File "ggsl_sparse.py", line 116, in ? V,H,P = g_arnoldi(A,B,R0,n,s,k) File "ggsl_sparse.py", line 68, in g_arnoldi V[:,:s] = R0/beta File "/usr/lib64/python2.4/site-packages/scipy/sparse/sparse.py", line 1191, in __setitem__ self.indptr = resize1d(self.indptr, row+2) TypeError: unsupported operand type(s) for +: 'slice' and 'int' How can I fix this problem ? Nils

1 0

Basic linear algebra for sparse matrices
by Nils Wagner Feb. 21, 2006

Feb. 21, 2006

Hi Ed, *Also, SparseM* provides some basic R functionality for linear algebra with sparse matrices. http://cran.r-project.org/src/contrib/Descriptions/SparseM.html Nils

1 0

Matrix exponential for sparse matrices
by Nils Wagner Feb. 21, 2006

Feb. 21, 2006

Hi Ed, Here you can find some details on matrix exponentials for both, dense and sparse matrices. http://www.maths.uq.edu.au/expokit/guide.html#Usage Nils

1 0

Right way to import scipy
by Bill Baxter Feb. 21, 2006

Feb. 21, 2006

I posted this on the numpy list, but it seems more a scipy issue. If you do this: import numpy as N import scipy as S Then S.linalg looks to be identical to N.linalg. In particular it doesn't have things like lu, qr, gmres, etc. If you subsequently do 'from scipy.linalg import lu', then suddenly S.linalglooks different from N.linalg, and it has not just lu, but all the other scipy linalg goodies too. Is this a bug? I would like to add info about getting at these functions to the Scipy for Matlab Users wiki page, but I'd like to know first if A) this is the expected behavior, or just a bug and B) if these import statements are really the best way to get at scipy.linalg.lu et al. import numpy as N import scipy as S from scipy.linalg import lu (This is with Numpy 0.9.5 and Scipy 0.4.6. With numpy 0.9.4/scipy 0.4.4 the above imports make scipy.linalg.lu et al accessible, but they throw exceptions trying to access ndarray.dtypechar when you try to use them.) Thanks, Bill Baxter

3 4

mtrand no longer in scipy.stats
by David Huard Feb. 21, 2006

Feb. 21, 2006

Hi, I just updated from svn to the new release and I can't import mtrand anymore from scipy.stats. Does someone know where it went ? I simply want to set the seed before generating random values using the rvs method. Thanks, David

2 2

ANN: SciPy 0.4.6 released
by Ed Schofield Feb. 20, 2006

Feb. 20, 2006

=========================== SciPy 0.4.6 Scientific tools for Python =========================== I'm pleased to announce the release of SciPy 0.4.6. This release adds support for the latest NumPy (version 0.9.5), features a new maximum entropy package, improved spline evaluation, and various bug fixes for signaltools, sparse matrices, weave, conjugate gradients, and more. It is available for download from http://www.scipy.org/Download as a source tarball for Linux/Solaris/OS X/BSD/Windows (64-bit and 32-bit) and as an executable installer for Win32. More information on SciPy is available at http://www.scipy.org/ =========================== SciPy is an Open Source library of scientific tools for Python. It contains a variety of high-level science and engineering modules, including modules for statistics, optimization, integration, linear algebra, Fourier transforms, signal and image processing, genetic algorithms, ODE solvers, special functions, and more.

2 1

problem with test check_algebraic_log_weight in 0.4.6?
by Erick Tryzelaar Feb. 20, 2006

Feb. 20, 2006

I just tried scipy.test on the just-released 0.4.6, and it ends up exiting my interactive python with this message: >>> scipy.test(verbosity=5) ... [snip] ... check_algebraic_log_weight (scipy.integrate.quadpack.test_quadpack.test_quad) Adjust D1MACH by uncommenting data statements appropriate for your machine. STOP 779 Any idea what is happening? I'm running on mac os x, if that helps. Thanks, -e

1 0

Fwd: io.read_array and columns of strings
by Eric Nodwell Feb. 20, 2006

Feb. 20, 2006

Hi, I notice that io.read_array (scipy version 0.4.4, numpy version 0.9.4) does not seem to support reading a column of strings. For example: >>> A = io.read_array ('r1683.log.layers', atype=[(string,10), Float]) Traceback (most recent call last): File "<stdin>", line 1, in ? File "/sw/lib/python2.4/site-packages/scipy/io/array_import.py", line 358, in read_array if not atype[k] in typecodes: TypeError: 'in <string>' requires string as left operand Perhaps this is intended behaviour, although standard numpy arrays do support string types, for example A = array (('one', 'two'), (string,10)) Of course, it's easy to work around this by reading the data with a loop instead of read_array. So this is "wish-list" kind of item for convenience and consistency. Any thoughts? thanks, Eric

1 0

scipy cvs : bug in fftconvolve ?
by Xavier Gnata Feb. 20, 2006

Feb. 20, 2006

Hi, (If it is a repost, please just ignore). Current scipy CVS complains about a "NameError: global name 'fftn' is not defined" calling scipy.signal.fftconvolve(I,K) where I and K are 2d numpy arrays. Looking at signaltools.py, you can see that we only import the following fft related stuff: from scipy.fftpack import fft, ifft, ifftshift, fft2, ifft2 Looks like fftn is just missing. Anyway, the code is buggy because fftn is not defined anywhere. BTW, is fftconvole the best way to perform a 2D fft-based convolution with scipy? (If not, naming conventions are not very intuitive ;)) Xavier -- ############################################ Xavier Gnata CRAL - Observatoire de Lyon 9, avenue Charles André 69561 Saint Genis Laval cedex Phone: +33 4 78 86 85 58 Fax: +33 4 78 86 83 86 E-mail: gnata(a)obs.univ-lyon1.fr ############################################

1 0

ANN: Python Enthought Edition Version 0.9.2 Released
by Travis N. Vaught Feb. 18, 2006

Feb. 18, 2006

Enthought is pleased to announce the release of Python Enthought Edition Version 0.9.2 (http://code.enthought.com/enthon/) -- a python distribution for Windows. This is a kitchen-sink-included Python distribution including the following packages/tools out of the box: Numeric 24.2 SciPy 0.3.3 IPython 0.6.15 Enthought Tool Suite 1.0.2 wxPython 2.6.1.0 PIL 1.1.4 mingw 20030504-1 f2py 2.45.241_1926 MayaVi 1.5 Scientific Python 2.4.5 VTK 4.4 and many more... 0.9.2 Release Notes Summary --------------------------- Version 0.9.2 of Python Enthought Edition is the first to include the Enthought Tool Suite Package (http://code.enthought.com/ets/). Other changes include upgrading to Numeric 24.2, including MayaVi 1.5 (rather than 1.3) and removing a standalone PyCrust package in favor of the one included with wxPython. Also, elementtree and celementtree have been added to the distribution. Notably, this release is still based on Python 2.3.5 and still includes SciPy 0.3.3. You'll also notice that we have changed the version numbering to a major.minor.point format (from a build number format). see full release notes at: http://code.enthought.com/release/changelog-enthon0.9.2.shtml Best, Travis N. Vaught

3 4