Mailman 3 December 2015 - NumPy-Discussion

f2py, numpy.distutils and multiple Fortran source files
by David Verelst Dec. 5, 2015

Dec. 5, 2015

Hi, For the wafo [1] package we are trying to include the extension compilation process in setup.py [2] by using setuptools and numpy.distutils [3]. Some of the extensions have one Fortran interface source file, but it depends on several other Fortran sources (modules). The manual compilation process would go as follows: gfortran -fPIC -c source_01.f gfortran -fPIC -c source_02.f f2py -m module_name -c source_01.o source_02.o source_interface.f Can this workflow be incorporated into … [View More]

3 3

NumPy Pyston support
by Boxiang Sun Dec. 5, 2015

Dec. 5, 2015

Hi all, I am an open source contributor of Pyston, and working on NumPy Pyston support(Pyston: https://github.com/dropbox/pyston) I disabled some tests which cause segfaults. This is current situation(in my local branch): ---------------------------------- Ran 3240 tests in 80.811s FAILED (KNOWNFAIL=3, SKIP=6, errors=334, failures=74) ---------------------------------- Indeed, the time consuming was longer than CPython. I think this is because of the additional exceptions / frame … [View More]

1 0

future of f2py and Fortran90+
by Eric Firing Dec. 4, 2015

Dec. 4, 2015

F2py is a great tool, but my impression is that it is being left behind by the evolution of Fortran from F90 onward. This is unfortunate; it would be nice to be able to easily wrap new Fortran libraries. I'm curious: has anyone been looking into what it would take to enable f2py to handle modern Fortran in general? And into prospects for getting such an effort funded? Eric

5 6

Recognizing a cycle in a vector
by Manolo Martínez Dec. 3, 2015

Dec. 3, 2015

Dear all, Suppose that I have a vector with the numerical solution of a differential equation -- more concretely, I am working with evolutionary game theory models, and the solutions are frequencies of types in a population that follows the replicator dynamics; but this is probably irrelevant. Sometimes these solutions are cyclical, yet I sample at points which do not correspond with the period of the cycle, so that np.allclose() cannot be directly applied. Is there any way to check for … [View More]

5 11

ANN: SciExp^2 (1.1.2)
by Lluís Vilanova Dec. 2, 2015

Dec. 2, 2015

Hi all! I'm really happy to make the first public announcement of SciExp^2 (actually, it's release 1.1.2). Home page: https://projects.gso.ac.upc.edu/projects/sciexp2 Description ----------- SciExp² (aka SciExp square or simply SciExp2) stands for Scientific Experiment Exploration, which contains a comprehensive framework for easing the workflow of creating, executing and evaluating experiments. The driving idea behind SciExp² is the need for quick and effortless design-space exploration. … [View More]

1 0

Indexing structured masked arrays with multidimensional fields; what with fill_value?
by Gerrit Holl Dec. 1, 2015

Dec. 1, 2015

Hello, usually, a masked array's .fill_value attribute has ndim=0 and the same dtype as the data attribute: In [27]: ar = array((0, [[0.0, 0.0, 0.0], [0.0, 0.0, 0.0]], 0.0), dtype="int, (2,3)float, float") In [28]: arm = ma.masked_array(ar) In [29]: arm.fill_value.ndim Out[29]: 0 In [31]: arm.fill_value.dtype Out[31]: dtype([('f0', '<i8'), ('f1', '<f8', (2, 3)), ('f2', '<f8')]) What would be the optimal approach for .fill_value if I address the member "f1" in this case? The … [View More]

1 0

ANN: SfePy 2015.4
by Robert Cimrman Dec. 1, 2015

Dec. 1, 2015

I am pleased to announce release 2015.4 of SfePy. Description ----------- SfePy (simple finite elements in Python) is a software for solving systems of coupled partial differential equations by the finite element method or by the isogeometric analysis (preliminary support). It is distributed under the new BSD license. Home page: http://sfepy.org Mailing list: http://groups.google.com/group/sfepy-devel Git (source) repository, issue tracker, wiki: http://github.com/sfepy Highlights of this … [View More]

1 0

Inconsistent/unexpected indexing semantics
by Lluís Vilanova Nov. 30, 2015

Nov. 30, 2015

Hi, TL;DR: There's a pending pull request deprecating some behaviour I find unexpected. Does anyone object? Some time ago I noticed that numpy yields unexpected results in some very specific cases. An array can be used to index multiple elements of a single dimension: >>> a = np.arange(8).reshape((2,2,2)) >>> a[ np.array([[0], [0]]) ] array([[[[0, 1], [2, 3]]], [[[0, 1], [2, 3]]]]) Nonetheless, if a list is used … [View More]

3 2