Mailman 3 January 2012 - SciPy-User

matplotlib: Simple legend code no longer works after upgrade to Ubuntu 11.10
by C Barrington-Leigh March 4, 2012

March 4, 2012

Oops; I just posted this to comp.lang.python, but I wonder whether matplotlib questions are supposed to go to scipy-user? Here it is: """ Before I upgraded to 2.7.2+ / 4 OCt 2011, the following code added a comment line to an axis legend using matplotlib / pylab. Now, the same code makes the legend appear "off-screen", ie way outside the axes limits. Can anyone help? And/or is there a new way to add a title and footer to the legend? Thanks! """ from pylab import * plot([0,0],[1,1],label='Ubuntu 11.10') lh=legend(fancybox=True,shadow=False) lh.get_frame().set_alpha(0.5) from matplotlib.offsetbox import TextArea, VPacker fontsize=lh.get_texts()[0].get_fontsize() legendcomment=TextArea('extra comments here', textprops=dict(size=fontsize)) show() # Looks fine here lh._legend_box = VPacker(pad=5, sep=0, children=[lh._legend_box,legendcomment], align="left") lh._legend_box.set_figure(gcf()) draw()

3 2

masked recarray, recarray with one field of type "ndarray"
by Gustavo Goretkin Feb. 19, 2012

Feb. 19, 2012

Does a recarray support masking? Can I have a recarray where one of the fields is an M-by-N ndarray (not recarray) of some dtype? ex: a = np.recarray(shape=(10),formats=['i4','f8','3-by-3 ndarray of dtype=float64'])

2 3

Accumulation sum using indirect indexes
by Alexander Kalinin Feb. 7, 2012

Feb. 7, 2012

Hello! I use SciPy in computer graphics applications. My task is to calculate vertex normals by averaging faces normals. In other words I want to accumulate vectors with the same ids. For example, ids = numpy.array([0, 1, 1, 2]) n = numpy.array([ [0.1, 0.1, 0.1], [0.1, 0.1, 0.1], [0.1, 0.1, 0.1], [0.1, 0.1 0.1] ]) I need result: nv = ([ [0.1, 0.1, 0.1], [0.2, 0.2, 0.2], [0.1, 0.1, 0.1]]) The most simple code: nv[ids] += n does not work, I know about this. For 1D arrays I use numpy.bincount(...) function. But this function does not work for 2D arrays. So, my question. What is the best way calculate accumulation sum for 2D arrays using indirect indexes? Sincerely, Alexander

7 16

scipy interpolate issues
by Travis Vaught Jan. 31, 2012

Jan. 31, 2012

Greetings, I'm experimenting with getting a smooth surface from a set of 3D points. With the attached CME (requires mayavi for the visualization) I'm getting unexpected results. The generated surface does not pass through the "known" points. When I transpose the results, it looks correct. Note: I also show the points connected by a Delaunay triangulation, whose surface is shown as contour lines. Am I using scipy.interpolate.Rbf inappropriately? TIA, Travis With transposing: without transposing:

2 1

[scipy-user] How to remove a value from an np array?
by Fabien Lafont Jan. 31, 2012

Jan. 31, 2012

I have that [3,2,4,8,7,8,9] and I want [3,2,8,7,8] how can I do? I've tried remove() but it works only on lists not on np.array.

7 11

distribution confidence interval for any percentile
by Christian K. Jan. 31, 2012

Jan. 31, 2012

Hi, dist.lognorm.interval provides the confidence interval for the median of a lognorm distribution. How can I get the confidence interval for any percentile, not only the median? Thanks in advance, Christian

1 0

[ANN] Summer School "Advanced Scientific Programming in Python" in Kiel, Germany
by Tiziano Zito Jan. 31, 2012

Jan. 31, 2012

Advanced Scientific Programming in Python ========================================= a Summer School by the G-Node and the Institute of Experimental and Applied Physics, Christian-Albrechts-Universität zu Kiel Scientists spend more and more time writing, maintaining, and debugging software. While techniques for doing this efficiently have evolved, only few scientists actually use them. As a result, instead of doing their research, they spend far too much time writing deficient code and reinventing the wheel. In this course we will present a selection of advanced programming techniques, incorporating theoretical lectures and practical exercises tailored to the needs of a programming scientist. New skills will be tested in a real programming project: we will team up to develop an entertaining scientific computer game. We use the Python programming language for the entire course. Python works as a simple programming language for beginners, but more importantly, it also works great in scientific simulations and data analysis. We show how clean language design, ease of extensibility, and the great wealth of open source libraries for scientific computing and data visualization are driving Python to become a standard tool for the programming scientist. This school is targeted at Master or PhD students and Post-docs from all areas of science. Competence in Python or in another language such as Java, C/C++, MATLAB, or Mathematica is absolutely required. Basic knowledge of Python is assumed. Participants without any prior experience with Python should work through the proposed introductory materials before the course. Date and Location ================= September 2—7, 2012. Kiel, Germany. Preliminary Program =================== Day 0 (Sun Sept 2) — Best Programming Practices - Best Practices, Development Methodologies and the Zen of Python - Version control with git - Object-oriented programming & design patterns Day 1 (Mon Sept 3) — Software Carpentry - Test-driven development, unit testing & quality assurance - Debugging, profiling and benchmarking techniques - Best practices in data visualization - Programming in teams Day 2 (Tue Sept 4) — Scientific Tools for Python - Advanced NumPy - The Quest for Speed (intro): Interfacing to C with Cython - Advanced Python I: idioms, useful built-in data structures, generators Day 3 (Wed Sept 5) — The Quest for Speed - Writing parallel applications in Python - Programming project Day 4 (Thu Sept 6) — Efficient Memory Management - When parallelization does not help: the starving CPUs problem - Advanced Python II: decorators and context managers - Programming project Day 5 (Fri Sept 7) — Practical Software Development - Programming project - The Pelita Tournament Every evening we will have the tutors' consultation hour: Tutors will answer your questions and give suggestions for your own projects. Applications ============ You can apply on-line at http://python.g-node.org Applications must be submitted before 23:59 UTC, May 1, 2012. Notifications of acceptance will be sent by June 1, 2012. No fee is charged but participants should take care of travel, living, and accommodation expenses. Candidates will be selected on the basis of their profile. Places are limited: acceptance rate last time was around 20%. Prerequisites: You are supposed to know the basics of Python to participate in the lectures. You are encouraged to go through the introductory material available on the website. Faculty ======= - Francesc Alted, Continuum Analytics Inc., USA - Pietro Berkes, Enthought Inc., UK - Valentin Haenel, Blue Brain Project, École Polytechnique Fédérale de Lausanne, Switzerland - Zbigniew Jędrzejewski-Szmek, Faculty of Physics, University of Warsaw, Poland - Eilif Muller, Blue Brain Project, École Polytechnique Fédérale de Lausanne, Switzerland - Emanuele Olivetti, NeuroInformatics Laboratory, Fondazione Bruno Kessler and University of Trento, Italy - Rike-Benjamin Schuppner, Technologit GbR, Germany - Bartosz Teleńczuk, Unité de Neurosciences Information et Complexité, Centre National de la Recherche Scientifique, France - Stéfan van der Walt, Helen Wills Neuroscience Institute, University of California Berkeley, USA - Bastian Venthur, Berlin Institute of Technology and Bernstein Focus Neurotechnology, Germany - Niko Wilbert, TNG Technology Consulting GmbH, Germany - Tiziano Zito, Institute for Theoretical Biology, Humboldt-Universität zu Berlin, Germany Organized by Christian T. Steigies and Christian Drews of the Institute of Experimental and Applied Physics, Christian-Albrechts-Universität zu Kiel , and by Zbigniew Jędrzejewski-Szmek and Tiziano Zito for the German Neuroinformatics Node of the INCF. Website: http://python.g-node.org Contact: python-info(a)g-node.org

1 0

[scipy_user] How to keep only the x first terms of an 1D array?
by Fabien Lafont Jan. 30, 2012

Jan. 30, 2012

Hello, Do somebody knows how to keep only the x first terms of an array? Fabien

2 1

Assimulo
by Christian Andersson Jan. 30, 2012

Jan. 30, 2012

Hello All, I would like to make you aware of that we have released a new version of Assimulo, which is a Python package for solving both ordinary differential equations and differential algebraic equations. Current interfaces include access to the Sundials solvers CVode and IDA with both sensitivity capabilities and event handling. Additionally, access to Radau5, Dopri5 and Rodas [by Hairer] are available, where Radau5 can also be used to solve differential algebraic equations. There are also plans to implement event handling for solvers that are missing the feature and include interfaces for more solvers. Assimulo is used as a teaching tool at Lund University and included as the default simulation environment in the platform JModelica.org where it has been used to solve industrial related problems with hundreds of states. Also, it can be used together with PyFMI, http://www.pyfmi.org, to simulate "Functional Mock-up Units". Assimulo can be found for download on http://www.assimulo.org and the documentation can be found on http://www.jmodelica.org/assimulo. Any suggestions and comments are very welcome! Best /Christian

1 0

SciPy 0.10: 8 failures on OS X Lion
by Samuel John Jan. 26, 2012

Jan. 26, 2012

Hi SciPy devs, I just want to report the following failures. Does anybody on OSX get the same/similar? Are they "known" ? Should I try scipy 0.11 because they are fixed there, perhaps? Six out of eight are related to arpack. I have installed numpy, scipy with gfortran provided by homebrew and I set the CC, CXX to use the non-llvm compilers of the 4.2 series. > Python 2.7.2 (default, Jan 13 2012, 15:00:18) > [GCC 4.2.1 (Based on Apple Inc. build 5658) (LLVM build 2336.1.00)] on darwin > >>> import scipy > >>> scipy.test() > Running unit tests for scipy > NumPy version 2.0.0.dev-55472ca > NumPy is installed in /usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy > SciPy version 0.10.0 > SciPy is installed in /usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy > Python version 2.7.2 (default, Jan 13 2012, 15:00:18) [GCC 4.2.1 (Based on Apple Inc. build 5658) (LLVM build 2336.1.00)] > nose version 1.1.2 > . > . > . > ====================================================================== > FAIL: test_arpack.test_symmetric_modes(True, <std-symmetric>, 'f', 2, 'LM', None, 0.5, <function aslinearoperator at 0x103a58f50>, None, 'normal') > ---------------------------------------------------------------------- > Traceback (most recent call last): > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest > self.test(*self.arg) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 235, in eval_evec > assert_allclose(LHS, RHS, rtol=rtol, atol=atol, err_msg=err) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 1213, in assert_allclose > verbose=verbose, header=header) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 677, in assert_array_compare > raise AssertionError(msg) > AssertionError: > Not equal to tolerance rtol=0.00178814, atol=0.000357628 > error for eigsh:standard, typ=f, which=LM, sigma=0.5, mattype=aslinearoperator, OPpart=None, mode=normal > (mismatch 100.0%) > x: array([[ 0.23815642, 0.1763755 ], > [-0.10785346, -0.32103487], > [ 0.12468303, -0.11230416],... > y: array([[ 0.23815642, 0.24814051], > [-0.10785347, -0.15634772], > [ 0.12468302, 0.05671416],... > > ====================================================================== > FAIL: test_arpack.test_symmetric_modes(True, <std-symmetric>, 'f', 2, 'LM', None, 0.5, <function aslinearoperator at 0x103a58f50>, None, 'cayley') > ---------------------------------------------------------------------- > Traceback (most recent call last): > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest > self.test(*self.arg) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 235, in eval_evec > assert_allclose(LHS, RHS, rtol=rtol, atol=atol, err_msg=err) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 1213, in assert_allclose > verbose=verbose, header=header) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 677, in assert_array_compare > raise AssertionError(msg) > AssertionError: > Not equal to tolerance rtol=0.00178814, atol=0.000357628 > error for eigsh:standard, typ=f, which=LM, sigma=0.5, mattype=aslinearoperator, OPpart=None, mode=cayley > (mismatch 100.0%) > x: array([[ 0.23815693, -0.33630507], > [-0.10785286, 0.02168 ], > [ 0.12468344, -0.11036437],... > y: array([[ 0.23815643, -0.2405392 ], > [-0.10785349, 0.14390968], > [ 0.12468311, -0.04574991],... > > ====================================================================== > FAIL: test_arpack.test_symmetric_modes(True, <std-symmetric>, 'f', 2, 'LA', None, 0.5, <function aslinearoperator at 0x103a58f50>, None, 'normal') > ---------------------------------------------------------------------- > Traceback (most recent call last): > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest > self.test(*self.arg) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 235, in eval_evec > assert_allclose(LHS, RHS, rtol=rtol, atol=atol, err_msg=err) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 1213, in assert_allclose > verbose=verbose, header=header) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 677, in assert_array_compare > raise AssertionError(msg) > AssertionError: > Not equal to tolerance rtol=0.00178814, atol=0.000357628 > error for eigsh:standard, typ=f, which=LA, sigma=0.5, mattype=aslinearoperator, OPpart=None, mode=normal > (mismatch 100.0%) > x: array([[ 28.80129188, -0.6379945 ], > [ 34.79312355, 0.27066791], > [-270.23255444, 0.4851834 ],... > y: array([[ 3.93467650e+03, -6.37994494e-01], > [ 3.90913859e+03, 2.70667916e-01], > [ -3.62176382e+04, 4.85183382e-01],... > > ====================================================================== > FAIL: test_arpack.test_symmetric_modes(True, <std-symmetric>, 'f', 2, 'SA', None, 0.5, <function aslinearoperator at 0x103a58f50>, None, 'normal') > ---------------------------------------------------------------------- > Traceback (most recent call last): > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest > self.test(*self.arg) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 235, in eval_evec > assert_allclose(LHS, RHS, rtol=rtol, atol=atol, err_msg=err) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 1213, in assert_allclose > verbose=verbose, header=header) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 677, in assert_array_compare > raise AssertionError(msg) > AssertionError: > Not equal to tolerance rtol=0.00178814, atol=0.000357628 > error for eigsh:standard, typ=f, which=SA, sigma=0.5, mattype=aslinearoperator, OPpart=None, mode=normal > (mismatch 100.0%) > x: array([[ 0.26260981, 0.23815559], > [-0.09760907, -0.10785484], > [ 0.06149647, 0.12468203],... > y: array([[ 0.23744165, 0.2381564 ], > [-0.13633069, -0.10785359], > [ 0.03132561, 0.12468301],... > > ====================================================================== > FAIL: test_arpack.test_symmetric_modes(True, <std-symmetric>, 'f', 2, 'SA', None, 0.5, <function aslinearoperator at 0x103a58f50>, None, 'cayley') > ---------------------------------------------------------------------- > Traceback (most recent call last): > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest > self.test(*self.arg) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 235, in eval_evec > assert_allclose(LHS, RHS, rtol=rtol, atol=atol, err_msg=err) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 1213, in assert_allclose > verbose=verbose, header=header) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 677, in assert_array_compare > raise AssertionError(msg) > AssertionError: > Not equal to tolerance rtol=0.00178814, atol=0.000357628 > error for eigsh:standard, typ=f, which=SA, sigma=0.5, mattype=aslinearoperator, OPpart=None, mode=cayley > (mismatch 100.0%) > x: array([[ 0.29524244, -0.2381569 ], > [-0.08169955, 0.10785299], > [ 0.06645597, -0.12468332],... > y: array([[ 0.24180251, -0.23815646], > [-0.14191195, 0.10785349], > [ 0.03568392, -0.12468307],... > > ====================================================================== > FAIL: test_arpack.test_symmetric_modes(True, <gen-symmetric>, 'f', 2, 'SM', None, 0.5, <function aslinearoperator at 0x103a58f50>, None, 'buckling') > ---------------------------------------------------------------------- > Traceback (most recent call last): > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest > self.test(*self.arg) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 235, in eval_evec > assert_allclose(LHS, RHS, rtol=rtol, atol=atol, err_msg=err) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 1213, in assert_allclose > verbose=verbose, header=header) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 677, in assert_array_compare > raise AssertionError(msg) > AssertionError: > Not equal to tolerance rtol=0.00178814, atol=0.000357628 > error for eigsh:general, typ=f, which=SM, sigma=0.5, mattype=aslinearoperator, OPpart=None, mode=buckling > (mismatch 100.0%) > x: array([[-0.10940548, 0.01676016], > [-0.07154097, 0.4628113 ], > [ 0.06895222, 0.49206394],... > y: array([[-0.10940547, 0.05459438], > [-0.07154103, 0.31407543], > [ 0.06895217, 0.37578294],... > > ====================================================================== > FAIL: test_arpack.test_symmetric_modes(True, <gen-symmetric>, 'f', 2, 'SA', None, 0.5, <function aslinearoperator at 0x103a58f50>, None, 'cayley') > ---------------------------------------------------------------------- > Traceback (most recent call last): > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest > self.test(*self.arg) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 235, in eval_evec > assert_allclose(LHS, RHS, rtol=rtol, atol=atol, err_msg=err) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 1213, in assert_allclose > verbose=verbose, header=header) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 677, in assert_array_compare > raise AssertionError(msg) > AssertionError: > Not equal to tolerance rtol=0.00178814, atol=0.000357628 > error for eigsh:general, typ=f, which=SA, sigma=0.5, mattype=aslinearoperator, OPpart=None, mode=cayley > (mismatch 100.0%) > x: array([[-0.4404992 , -0.01935683], > [-0.25650678, -0.11053132], > [-0.36893024, -0.13223556],... > y: array([[-0.44017013, -0.0193569 ], > [-0.25525379, -0.11053158], > [-0.36818443, -0.13223571],... > > ====================================================================== > FAIL: test_contingency.test_expected_freq > ---------------------------------------------------------------------- > Traceback (most recent call last): > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest > self.test(*self.arg) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/stats/tests/test_contingency.py", line 45, in test_expected_freq > assert_array_equal(e, np.ones_like(observed)) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 753, in assert_array_equal > verbose=verbose, header='Arrays are not equal') > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 677, in assert_array_compare > raise AssertionError(msg) > AssertionError: > Arrays are not equal > > (mismatch 100.0%) > x: array([[[ 0., 0.], > [ 0., 0.]], > ... > y: array([[[1, 1], > [1, 1]], > ... > > ====================================================================== > FAIL: test_distributions.test_frozen_fit_ticket_1536 > ---------------------------------------------------------------------- > Traceback (most recent call last): > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest > self.test(*self.arg) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/stats/tests/test_distributions.py", line 747, in test_frozen_fit_ticket_1536 > assert_almost_equal(params, true, decimal=2) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 451, in assert_almost_equal > return assert_array_almost_equal(actual, desired, decimal, err_msg) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 846, in assert_array_almost_equal > header=('Arrays are not almost equal to %d decimals' % decimal)) > File "/usr/local/Cellar/python/2.7.2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/testing/utils.py", line 677, in assert_array_compare > raise AssertionError(msg) > AssertionError: > Arrays are not almost equal to 2 decimals > > (mismatch 66.6666666667%) > x: array([ 1.416, 0. , 0.061]) > y: array([ 0.25, 0. , 0.5 ]) > > ---------------------------------------------------------------------- > Ran 5095 tests in 54.587s > > FAILED (KNOWNFAIL=12, SKIP=36, failures=9) > <nose.result.TextTestResult run=5095 errors=0 failures=9> > And some warning that look strange (at least to me): > site-packages/scipy/spatial/__init__.py:26: RuntimeWarning: numpy.ndarray size changed, may indicate binary incompatibility > site-packages/scipy/spatial/__init__.py:26: RuntimeWarning: numpy.ufunc size changed, may indicate binary incompatibility cheers, Samuel

2 2