Mailman 3 April 2010 - SciPy-Dev

SLSQP Constrained Optimizer Status
by Rob Falck Nov. 24, 2013

Nov. 24, 2013

I'm currently implementing the Sequential Least Squares Quadratic Programming (SLSQP) optimizer, by Dieter Kraft, for use in Scipy. The Fortran code being wrapped with F2PY is here: http://www.netlib.org/toms/733 (its use within Scipy has been cleared) SLSQP provides for bounds on the independent variables, as well as equality and inequality constraint functions, which is a capability that doesn't exist in scipy.optimize. Currently the code works, although the constraint normals are being generated by approximation. I'm working on a way to pass these in. I think the most elegant way will be a single function that returns the matrix of constraint normals. For a demonstration of what the code can do, here is an optimization of f(x,y) = 2xy + 2x - x**2 - 2y**2 Example 14.2 in Chapra & Canale gives the maximum as x=2.0, y=1.0. The unbounded optimization tests find this solution. As expected, its faster when derivatives are provided rather than approximated. Unbounded optimization. Derivatives approximated. Elapsed time: 1.45792961121 ms Results [[1.9999999515712266, 0.99999996181577444], -1.9999999999999984, 4, 0, 'Optimization terminated successfully.'] Unbounded optimization. Derivatives provided. Elapsed time: 1.03211402893 ms Results [[1.9999999515712266, 0.99999996181577444], -1.9999999999999984, 4, 0, 'Optimization terminated successfully.'] The following example uses an equality constraint to find the optimal when x=y. Bound optimization. Derivatives approximated. Elapsed time: 1.384973526 ms Results [[0.99999996845920858, 0.99999996845920858], -0.99999999999999889, 4, 0, 'Optimization terminated successfully.'] I've tried to conform to the syntax used by the other optimizers in scipy.optimize. The function definition and doc string are below. If anyone is interested in testing it out, let me know. def fmin_slsqp( func, x0 , eqcons=[], f_eqcons=None, ieqcons=[], f_ieqcons=None, bounds = [], fprime = None, fprime_cons=None,args = (), iter = 100, acc = 1.0E-6, iprint = 1, full_output = 0, epsilon = _epsilon ): """ Minimize a function using Sequential Least SQuares Programming Description: Python interface function for the SLSQP Optimization subroutine originally implemented by Dieter Kraft. Inputs: func - Objective function (in the form func(x, *args)) x0 - Initial guess for the independent variable(s). eqcons - A list of functions of length n such that eqcons[j](x0,*args) == 0.0 in a successfully optimized problem f_eqcons - A function of the form f_eqcons(x, *args) that returns an array in which each element must equal 0.0 in a successfully optimized problem. If f_eqcons is specified, eqcons is ignored. ieqcons - A list of functions of length n such that ieqcons[j](x0,*args) >= 0.0 in a successfully optimized problem f_ieqcons - A function of the form f_ieqcons(x0, *args) that returns an array in which each element must be greater or equal to 0.0 in a successfully optimized problem. If f_ieqcons is specified, ieqcons is ignored. bounds - A list of tuples specifying the lower and upper bound for each independent variable [ (xl0, xu0), (xl1, xu1), ...] fprime - A function that evaluates the partial derivatives of func fprime_cons - A function of the form f(x, *args) that returns the m by n array of constraint normals. If not provided, the normals will be approximated. Equality constraint normals precede inequality constraint normals. args - A sequence of additional arguments passed to func and fprime iter - The maximum number of iterations (int) acc - Requested accuracy (float) iprint - The verbosity of fmin_slsqp. iprint <= 0 : Silent operation iprint == 1 : Print summary upon completion (default) iprint >= 2 : Print status of each iterate and summary full_output - If 0, return only the minimizer of func (default). Otherwise, output final objective function and summary information. epsilon - The step size for finite-difference derivative estimates. Outputs: ( x, { fx, gnorm, its, imode, smode }) x - The final minimizer of func. fx - The final value of the objective function. its - The number of iterations. imode - The exit mode from the optimizer, as an integer. smode - A string describing the exit mode from the optimizer. Exit modes are defined as follows: -1 : Gradient evaluation required (g & a) 0 : Optimization terminated successfully. 1 : Function evaluation required (f & c) 2 : Number of equality constraints larger than number of independent variables 3 : More than 3*n iterations in LSQ subproblem 4 : Inequality constraints incompatible 5 : Singular matrix E in LSQ subproblem 6 : Singular matrix C in LSQ subproblem 7 : Rank-deficient equality constraint subproblem HFTI 8 : Positive directional derivative for linesearch 9 : Iteration limit exceeded -- - Rob Falck

4 16

np.savetxt: apply patch in enhancement ticket 1079 to add headers?
by Skipper Seabold June 8, 2010

June 8, 2010

Hi all, I am assuming that this is ok to request via the list... Could we discuss or could someone apply the patch in enhancement ticket 1079? http://projects.scipy.org/numpy/ticket/1079 I needed this functionality recently, and this is a quick and easy fix that may have been overlooked. There is also another enhancement request about this here: http://projects.scipy.org/numpy/ticket/1236 The only thing that I can think of that might need to be added is a test to see that the header length is the same as the number of columns, but really that might just be up to the user to supply the right headers. It might also be nice to have a header = True, that uses the field names for a structured array, but I can live without that. Cheers, Skipper

3 11

preparing for scipy 0.8 release
by Ralf Gommers May 25, 2010

May 25, 2010

Hi all, Now that scipy 0.7.2 and numpy 1.4.1 are out the door, it's time to look at what needs to be done for scipy 0.8. Some questions - What needs to be done and what do you all want included that's not in trunk yet in terms of features / new code / code cleanups? - Right now scipy trunk needs to be built against numpy trunk. Is it feasible to build against numpy 1.3 still, to get one backwards-compatible release? If not, 0.8 should be released in parallel with numpy 2.0. - Time line: do we want a release before SciPy2010 (end of June)? If so, there's only a month left before a freeze. If not, we could just provide snapshot binaries in June. Here's some things that I personally think should be done: - clean up tests so nothing is printed when they are run (deprecation warnings etc) - clean up build warnings as far as possible - create a little bit of order in Trac Cheers, Ralf

11 26

namespace recommendation in scipy tutorial and docs for np.emath
by josef.pktd＠gmail.com May 15, 2010

May 15, 2010

I'm going through some scipy tickets and I found http://projects.scipy.org/scipy/ticket/882 I don't know much about the old usage of mixing (mixed up) namespaces. But I think it would be good if someone could check the namespace convention/recommendation in the basic scipy tutorial http://docs.scipy.org/doc/scipy/reference/tutorial/basic.html When I tried to figure out the difference between scipy.* and numpy.*, I also found that the functions in numpy.emath (numpy.scipy) have good docstrings, but they are not included in the docs (at least I don't find them), because http://docs.scipy.org/numpy/docs/numpy-docs/reference/routines.emath.rst/#n… is still a stub page without function links. http://docs.scipy.org/doc/numpy/reference/routines.emath.html Josef

2 2

io -> npyio move: doc wiki update
by Ralf Gommers May 13, 2010

May 13, 2010

Hi, Here's a reminder that after the rename of numpy.io, the pages and history of all its objects in the doc wiki need to be carried over. Cheers, Ralf

2 1

Update on the "Needs editing" numpy docstrings
by David Goldsmith May 4, 2010

May 4, 2010

Hi, all! I've been working on the remaining "Needs editing" status numpy docstrings, of which there are presently 74. I've compiled two lists: those I think might be classifiable as "Unimportant," and those I think might require somewhat lengthy docstrings (the latter list is meant to try to recruit helpers). ;-) Here are those lists: List of potentially "Unimportant"s timeinteger (autogenerated subclass of signedinteger) timedelta64 (autogenerated subclass of timeinteger) polytemplate.rel_import Everything in distutils that still has status "Needs editing" linalg (docstring - which consists of tables of objects - auto-generated?) lib.shape_base.get_array_wrap f2py.diagnose and its "Needs editing" status attributes numpy.lib._iotools.StringConverter.iterupgrade dtype.base and dtype.metadata List of potentially "long" ones All "Needs editing" status doc.* core.umath, core.shape_base, core.function_base numpy-docs/* If you have opinions on these (or if you look at the whole list of remaining "Needs editing" docstrings and feel that I've left some off the "Unimportant" list) please do chime in. Thanks, DG -- Mathematician: noun, someone who disavows certainty when their uncertainty set is non-empty, even if that set has measure zero.

2 1

Fortran cleanup advice?
by Ralf Gommers May 4, 2010

May 4, 2010

Hi, I'm trying to make a start cleaning up the build warnings. One of the reasons it took me a while to get going with building scipy was the huge amount of harmless(?) warnings that hide the actually significant ones. Most warnings are from Fortran code, I count six types (examples at end of mail): 1. Nonconforming tab character. This one's straightforward, change tabs to spaces. 2. Unused variable. 3. Obsolete ASSIGN statement. 4. Label defined but not used. 5. Obsolete assigned GOTO statement 6. ' uninitialized. For 2-6, is there a good way to fix these without hiding actual bugs? For example 2 can be fixed by assigning to variable at top of subroutine with a comment saying what's going on there, but I'm not sure it is a good idea. Ralf Warning: Nonconforming tab character in column 1 of line 510 build/scons/scipy/integrate/dop/dop853.f:364.42: & SOLOUT,IOUT,IDID,NMAX,UROUND,METH,NSTIFF,SAFE,BETA,FAC1,FAC2, 1 Warning: Unused variable meth declared at (1) build/scons/scipy/optimize/slsqp/slsqp_optmz.f:1933.72: 10 assign 30 to next 1 Warning: Obsolete: ASSIGN statement at (1) build/scons/scipy/sparse/linalg/eigen/arpack/ARPACK/SRC/znaup2.f:322.5: 10 continue 1 Warning: Label 10 at (1) defined but not used build/scons/scipy/special/cdflib/dinvr.f:346.18: GO TO i99999 1 Warning: Obsolete: Assigned GOTO statement at (1) build/scons/scipy/special/specfun/specfun.f:6598: warning: ' may be used uninitialized in this function 'REALPART_EXPR <cf>build/scons/scipy/special/specfun/specfun.f:6598: warning: ' may be used uninitialized in this function

4 5

Failure to build superlumodule using scons, rev 6364
by Christoph Gohlke April 30, 2010

April 30, 2010

Dear Scipy developers, scipy svn trunk rev 6364, specifically the superlumodule extension, fails to build for me on Windows/Python2.6.5/numpy1.4.1/numscons-0.12.0dev/msvc9/ifort11.1/mkl using the setupscons script: python setupscons.py scons -b --fcompiler=ifort --compiler=msvc bdist_wininst --user-access-control=auto scons: *** [build\scons\scipy\sparse\linalg\dsolve\_csuperlumodule.obj] Source `build\scons\scipy\sparse\linalg\dsolve\_csuperlumodule.c' not found, needed by target `build\scons\scipy\sparse\linalg\dsolve\_csuperlumodule.obj'. scons: building terminated because of errors. Thank you, Christoph

1 0

scipy.test bails with segmentation error
by hannes April 30, 2010

April 30, 2010

Dear List, I installed numpy and scipy from svn source as of yesterday (30 of april), however get an segmentation fault in the scipy.test().. Running unit tests for scipy NumPy version 1.4.1 NumPy is installed in /usr/local/lib/python2.6/dist-packages/numpy SciPy version 0.8.0.dev6362 SciPy is installed in /usr/local/lib/python2.6/dist-packages/scipy Python version 2.6.2 (release26-maint, Apr 19 2009, 01:58:18) [GCC 4.3.3] nose version 0.10.4 ...... ..saxpy:n=3 ....scopy:n=4 ..scopy:n=3 .............sscal:n=4 ....sswap:n=4 ..sswap:n=3 .....zaxpy:n=4 ..zaxpy:n=3 ....zcopy:n=4 ..zcopy:n=3 .............zscal:n=4 ....zswap:n=4 ..zswap:n=3 .......Result may be inaccurate, approximate err = 6.98199990784e-09 ...Result may be inaccurate, approximate err = 1.25504040833e-10 ..................................................................................................................... ...................................................................................................................... ......................................................................................................................... ........................................................................................................................ ..............ESSSSSSSSSSS.......Segmentation fault Anyone can shed some light ? python -c 'import os,sys;print os.name,sys.platform' posix linux2 uname -a Linux hirubu 2.6.28-18-generic #60-Ubuntu SMP Fri Mar 12 04:26:47 UTC 2010 x86_64 GNU/Linux hannes@hirubu:~$ gcc -v Using built-in specs. Target: x86_64-linux-gnu Configured with: ../src/configure -v --with-pkgversion='Ubuntu 4.3.3-5ubuntu4' --with-bugurl=file:///usr/share/doc/gcc-4.3/README.Bugs --enable-languages=c,c++,fortran,objc,obj-c++ --prefix=/usr --enable-shared --with-system-zlib --libexecdir=/usr/lib --without-included-gettext --enable-threads=posix --enable-nls --with-gxx-include-dir=/usr/include/c++/4.3 --program-suffix=-4.3 --enable-clocale=gnu --enable-libstdcxx-debug --enable-objc-gc --enable-mpfr --with-tune=generic --enable-checking=release --build=x86_64-linux-gnu --host=x86_64-linux-gnu --target=x86_64-linux-gnu Thread model: posix gcc version 4.3.3 (Ubuntu 4.3.3-5ubuntu4) hannes@hirubu:~$ g77 --version bash: g77: command not found hannes@hirubu:~$ python -c 'import sys;print sys.version' 2.6.2 (release26-maint, Apr 19 2009, 01:58:18) [GCC 4.3.3] for point 5 cd scipy/Lib/linalg python setup_atlas_version.py build_ext --inplace --force python -c 'import atlas_version' i cant give any information as the scipy dir does not contain any directory lib for python /usr/lib/python2.6/site-packages/numpy/distutils/system_info.py lapack_info: libraries lapack not found in /usr/local/lib FOUND: libraries = ['lapack'] library_dirs = ['/usr/lib'] language = f77 lapack_opt_info: lapack_mkl_info: mkl_info: libraries mkl,vml,guide not found in /usr/local/lib libraries mkl,vml,guide not found in /usr/lib NOT AVAILABLE NOT AVAILABLE atlas_threads_info: Setting PTATLAS=ATLAS libraries ptf77blas,ptcblas,atlas not found in /usr/local/lib libraries lapack_atlas not found in /usr/local/lib libraries ptf77blas,ptcblas,atlas not found in /usr/lib/atlas libraries lapack_atlas not found in /usr/lib/atlas libraries ptf77blas,ptcblas,atlas not found in /usr/lib __main__.atlas_threads_info NOT AVAILABLE atlas_info: libraries f77blas,cblas,atlas not found in /usr/local/lib libraries lapack_atlas not found in /usr/local/lib libraries f77blas,cblas,atlas not found in /usr/lib/atlas libraries lapack_atlas not found in /usr/lib/atlas __main__.atlas_info FOUND: libraries = ['lapack', 'f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib'] language = f77 include_dirs = ['/usr/include'] /usr/local/lib/python2.6/dist-packages/numpy/distutils/command/config.py:394: DeprecationWarning: +++++++++++++++++++++++++++++++++++++++++++++++++ Usage of get_output is deprecated: please do not use it anymore, and avoid configuration checks involving running executable on the target machine. +++++++++++++++++++++++++++++++++++++++++++++++++ DeprecationWarning) customize GnuFCompiler Could not locate executable g77 Found executable /usr/bin/f77 gnu: no Fortran 90 compiler found gnu: no Fortran 90 compiler found customize IntelFCompiler Could not locate executable ifort Could not locate executable ifc customize LaheyFCompiler Could not locate executable lf95 customize PGroupFCompiler Could not locate executable pgf90 Could not locate executable pgf77 customize AbsoftFCompiler Could not locate executable f90 absoft: no Fortran 90 compiler found absoft: no Fortran 90 compiler found absoft: no Fortran 90 compiler found absoft: no Fortran 90 compiler found absoft: no Fortran 90 compiler found absoft: no Fortran 90 compiler found customize NAGFCompiler Found executable /usr/bin/f95 customize VastFCompiler customize GnuFCompiler gnu: no Fortran 90 compiler found gnu: no Fortran 90 compiler found customize CompaqFCompiler Could not locate executable fort customize IntelItaniumFCompiler Could not locate executable efort Could not locate executable efc customize IntelEM64TFCompiler customize Gnu95FCompiler Found executable /usr/bin/gfortran customize Gnu95FCompiler customize Gnu95FCompiler using config FOUND: libraries = ['lapack', 'f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib'] language = f77 define_macros = [('ATLAS_INFO', '"\\"?.?.?\\""')] include_dirs = ['/usr/include'] wx_info: Found executable /usr/bin/wx-config FOUND: libraries = ['wx_baseu_xml-2.8', 'wx_baseu_net-2.8', 'wx_baseu-2.8'] extra_link_args = ['-pthread', '-Wl,-Bsymbolic-functions'] define_macros = [('WX_INFO', '"\\"2.8.9\\""'), ('WX_VERSION_2_8_9', None), ('WX_RELEASE_2_8', None), ('_FILE_OFFSET_BITS', '64'), ('_LARGE_FILES', None), ('wxUSE_GUI', '0')] extra_compile_args = ['-pthread'] include_dirs = ['/usr/lib/wx/include/base-unicode-release-2.8', '/usr/include/wx-2.8'] lapack_atlas_info: libraries lapack_atlas,f77blas,cblas,atlas not found in /usr/local/lib libraries lapack_atlas not found in /usr/local/lib libraries lapack_atlas,f77blas,cblas,atlas not found in /usr/lib/atlas libraries lapack_atlas not found in /usr/lib/atlas __main__.lapack_atlas_info FOUND: libraries = ['lapack', 'lapack_atlas', 'f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib'] language = f77 include_dirs = ['/usr/include'] umfpack_info: libraries umfpack not found in /usr/local/lib libraries umfpack not found in /usr/lib NOT AVAILABLE _pkg_config_info: Found executable /usr/bin/pkg-config NOT AVAILABLE lapack_atlas_threads_info: Setting PTATLAS=ATLAS libraries lapack_atlas,ptf77blas,ptcblas,atlas not found in /usr/local/lib libraries lapack_atlas not found in /usr/local/lib libraries lapack_atlas,ptf77blas,ptcblas,atlas not found in /usr/lib/atlas libraries lapack_atlas not found in /usr/lib/atlas libraries lapack_atlas,ptf77blas,ptcblas,atlas not found in /usr/lib __main__.lapack_atlas_threads_info NOT AVAILABLE x11_info: libraries X11 not found in /usr/local/lib FOUND: libraries = ['X11'] library_dirs = ['/usr/lib'] include_dirs = ['/usr/include'] blas_info: libraries blas not found in /usr/local/lib FOUND: libraries = ['blas'] library_dirs = ['/usr/lib'] language = f77 fftw_info: libraries fftw3 not found in /usr/local/lib libraries fftw3 not found in /usr/lib fftw3 not found libraries rfftw,fftw not found in /usr/local/lib libraries rfftw,fftw not found in /usr/lib fftw2 not found NOT AVAILABLE f2py_info: FOUND: sources = ['/usr/local/lib/python2.6/dist-packages/numpy/f2py/src/fortranobject.c'] include_dirs = ['/usr/local/lib/python2.6/dist-packages/numpy/f2py/src'] gdk_pixbuf_xlib_2_info: NOT AVAILABLE dfftw_threads_info: libraries drfftw_threads,dfftw_threads not found in /usr/local/lib libraries drfftw_threads,dfftw_threads not found in /usr/lib dfftw threads not found NOT AVAILABLE atlas_blas_info: libraries f77blas,cblas,atlas not found in /usr/local/lib libraries f77blas,cblas,atlas not found in /usr/lib/atlas FOUND: libraries = ['f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib'] language = c include_dirs = ['/usr/include'] fftw3_info: libraries fftw3 not found in /usr/local/lib libraries fftw3 not found in /usr/lib fftw3 not found NOT AVAILABLE blas_opt_info: blas_mkl_info: libraries mkl,vml,guide not found in /usr/local/lib libraries mkl,vml,guide not found in /usr/lib NOT AVAILABLE atlas_blas_threads_info: Setting PTATLAS=ATLAS libraries ptf77blas,ptcblas,atlas not found in /usr/local/lib libraries ptf77blas,ptcblas,atlas not found in /usr/lib/atlas libraries ptf77blas,ptcblas,atlas not found in /usr/lib NOT AVAILABLE customize GnuFCompiler gnu: no Fortran 90 compiler found gnu: no Fortran 90 compiler found customize IntelFCompiler customize LaheyFCompiler customize PGroupFCompiler customize AbsoftFCompiler absoft: no Fortran 90 compiler found absoft: no Fortran 90 compiler found absoft: no Fortran 90 compiler found absoft: no Fortran 90 compiler found absoft: no Fortran 90 compiler found absoft: no Fortran 90 compiler found customize NAGFCompiler customize VastFCompiler customize GnuFCompiler gnu: no Fortran 90 compiler found gnu: no Fortran 90 compiler found customize CompaqFCompiler customize IntelItaniumFCompiler customize IntelEM64TFCompiler customize Gnu95FCompiler customize Gnu95FCompiler customize Gnu95FCompiler using config FOUND: libraries = ['f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib'] language = c define_macros = [('ATLAS_INFO', '"\\"?.?.?\\""')] include_dirs = ['/usr/include'] sfftw_info: libraries srfftw,sfftw not found in /usr/local/lib libraries srfftw,sfftw not found in /usr/lib sfftw not found NOT AVAILABLE xft_info: FOUND: libraries = ['Xft'] define_macros = [('XFT_INFO', '"\\"2.1.13\\""'), ('XFT_VERSION_2_1_13', None)] include_dirs = ['/usr/include/freetype2'] fft_opt_info: fftw2_info: libraries rfftw,fftw not found in /usr/local/lib libraries rfftw,fftw not found in /usr/lib fftw2 not found NOT AVAILABLE dfftw_info: libraries drfftw,dfftw not found in /usr/local/lib libraries drfftw,dfftw not found in /usr/lib dfftw not found NOT AVAILABLE djbfft_info: NOT AVAILABLE NOT AVAILABLE gdk_x11_2_info: NOT AVAILABLE agg2_info: NOT AVAILABLE numarray_info: NOT AVAILABLE blas_src_info: NOT AVAILABLE fftw_threads_info: libraries rfftw_threads,fftw_threads not found in /usr/local/lib libraries rfftw_threads,fftw_threads not found in /usr/lib fftw threads not found NOT AVAILABLE _numpy_info: FOUND: define_macros = [('NUMERIC_VERSION', '"\\"24.2\\""'), ('NUMERIC', None)] gdk_info: NOT AVAILABLE gtkp_x11_2_info: NOT AVAILABLE sfftw_threads_info: libraries srfftw_threads,sfftw_threads not found in /usr/local/lib libraries srfftw_threads,sfftw_threads not found in /usr/lib sfftw threads not found NOT AVAILABLE boost_python_info: NOT AVAILABLE freetype2_info: FOUND: libraries = ['freetype'] define_macros = [('FREETYPE2_INFO', '"\\"9.20.3\\""'), ('FREETYPE2_VERSION_9_20_3', None)] include_dirs = ['/usr/include/freetype2'] gdk_2_info: NOT AVAILABLE lapack_src_info: NOT AVAILABLE gtkp_2_info: NOT AVAILABLE gdk_pixbuf_2_info: NOT AVAILABLE amd_info: libraries amd not found in /usr/local/lib libraries amd not found in /usr/lib NOT AVAILABLE Numeric_info: FOUND: define_macros = [('NUMERIC_VERSION', '"\\"24.2\\""'), ('NUMERIC', None)] numerix_info: numpy_info: FOUND: define_macros = [('NUMPY_VERSION', '"\\"1.4.1\\""'), ('NUMPY', None)] FOUND: define_macros = [('NUMPY_VERSION', '"\\"1.4.1\\""'), ('NUMPY', None)]

3 6

numpy 1.4.1 f2py on Snow Leopard.
by Charles سمير Doutriaux April 30, 2010

April 30, 2010

Hello, I jsut built numpy 1.4.1 on my mac (64bit only no universal) Now f2py fails building any fortran code... my python: 2.6.5 meryem /svn/cdat/trunk >gcc --version i686-apple-darwin10-gcc-4.2.1 (GCC) 4.2.1 (Apple Inc. build 5646) (dot 1) Copyright (C) 2007 Free Software Foundation, Inc. This is free software; see the source for copying conditions. There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. meryem /svn/cdat/trunk >gfortran --version GNU Fortran (GCC) 4.5.0 20100107 (experimental) Copyright (C) 2009 Free Software Foundation, Inc. GNU Fortran comes with NO WARRANTY, to the extent permitted by law. You may redistribute copies of GNU Fortran under the terms of the GNU General Public License. For more information about these matters, see the file named COPYING I'm also attaching bzipped logs of my numpy and scipy (as an example that would be familiar to you guys) installs. scipy obviously fails because the f2py doesn't work. Thx, C.

2 2