PyDSTool temporary web pages

The Front page and Main project page will temporarily redirect here while the original wiki server is being restored. Many of the wiki pages are not yet available here.

Project Overview

  1. Project overview
  2. Motivation
  3. Features implemented
  4. User interface
  5. Why Python?
  6. Design philosophy
  7. Coming attractions
  8. Historical context
  9. Comparison to other packages
    1. General packages
      1. Matlab / Simulink
      2. XPP
      3. AUTO and MATCONT
      4. Custom code in other languages
    2. Application-specific packages
      1. NEURON and GENESIS

1. Project overview

With PyDSTool we aim to provide a powerful suite of computational tools for the development, simulation, and analysis of dynamical systems that are used for the modeling of physical processes in many scientific disciplines. We place emphasis on the support of data analysis as part of the process of data-driven modeling. Our focus is on models involving ordinary differential equations (ODEs), differential-algebraic equations (DAEs), and discrete maps.

The software package is written primarily in the high-level interpreted language [WWW] Python, and is freely available under the BSD open-source license from the [WWW] SourceForge community website. For certain computationally intensive tasks, code written in C and Fortran can be invoked.

In pursuing our aims we have developed a hierarchy of object-oriented programming classes available to the modeler. These classes provide a means for “data abstraction”, that is, to provide the user with intuitive “objects” that encapsulate common data structures and methods for manipulating that data through a well-defined interface. These classes naturally make the software highly modular, and they range from defining hybrid model simulators, to individual data points, to symbolic mathematical expressions, to bifurcation analyzers. Some of these classes are built over classes exported by SciPy. We believe that many of these classes have wider scientific utility.

The PyDSTool software is "research code" in a Beta stage of development, and should not be treated as a complete or comprehensive dynamical systems package, with the associated expectation that its design and implementation have thoroughly stabilized and have been well tested. We have added features as and when we have had a use for them in our own research, and have omitted many important features that we would love to add if time permits us, or if our research so demands.

You might like to submit feature requests, or you may also like to contribute to the code yourself. We are also interested to hear your opinions about the possibility of adding some of our classes to SciPy (perhaps in modified form). Please contact us at the [WWW] SourceForge open discussion forum or via email.

2. Motivation

A [WWW] dynamical system is defined by a set of rules or transformations for determining how points in a multidimensional space move in time. Time may be either discrete or continuous. The traces of the points as they move in time are called trajectories. The goal of Dynamical Systems Theory is to provide a comprehensive description of the geometric structures arising from these trajectories. In addition to elucidating the dynamics associated with an individual dynamical ssystem, bifurcation theory may be used to describe how the dynamics of a system varies with changes in parameter values.

Interactive numerical and graphical exploration are important tools in dynamical systems research for several reasons:

Consequently, there is a critical need for computational environments that provide effective tools for exploring dynamical systems with minimial effort on the part of the user. The PyDSTool project attempts to provide such an environment for use by both casual and advanced practitioners of Dynamical Systems Theory.

3. Features implemented

Our package incorporates:

Most of these features are simply not available in the native SciPy ODE integrators.

4. User interface

Users of PyDSTool need to be familiar with working in interactive, command-line environments such as UNIX and Matlab, including the writing of simple command scripts. There is presently no graphical interface for PyDSTool. Our emphasis is on the interactivity of a command-line and the rapid prototyping possibilities of script-based computing.

In building a core library of Python classes, supporting many fundamental concepts in dynamical systems modeling, we provide more than just a glue with which to interface multiple tools. Our classes involve storing and maintaining a "context" that carries a lot of useful mathematical baggage. Through interaction with our Python environment at the script level, users can build complex models in a structured way, and have access to mathematically intuitive information about the models, using the intrinsic context of all the Python objects at the heart of their computations.

Our UI model is for users to interactively "query" objects for basic information (known in Python-speak as introspection), and also to be able to treat them as unitary objects of computation for use with tools and utilities such as optimizers, parameter estimators, and so forth.

We believe it is crucial for users to be able to combine the application of tools in a nested or interleaved fashion, in order to make the most flexible and dynamic manipulations of a model. Such rich combinations are practically impossible in disjointed software environments, and we believe our community is eager to be able to smoothly set-up and maintain such situations for their own modeling projects. It is a challenge to cleanly and efficiently interface different legacy algorithms with the core Python code in order to maximize the use and re-use of the context associated with the core objects.

Users are provided with an interface for the specification of both simple and complex dynamical systems models, using minimal programming syntax, and a range of options in converting these abstract specifications into instantiated numerical solvers for a specific system. Within the same interactive session, users have immediate access to analysis tools for continuation, parameter estimation, optimization, and so on. These tools are each tailored for use with the core PyDSTool structures to ensure the user has to write as little additional computer code as possible. Extensive documentation for the project has been provided online on this wiki.

A key aspect in the design of PyDSTool is the provision of adequate diagnostic information and querying utilities for data structures and computations. Users can expect helpful information regarding the status of their model development and computations beyond the guidance of the online documentation, through in-built querying commands and detailed error messages. The object-oriented nature of the software also provides inherent protection of the users’ conceptualization of data-flow and control in their PyDSTool scripts. See the UserDocumentation for more details.

5. Why Python?

Python has several benefits as a platform for our package. It is an open-source language, enthusiastically developed by a large community of people drawn from both academic and commercial backgrounds. The powerful object-oriented and dynamic typing features of the language greatly aid in the development of flexible and intuitive data structures and user interface (UI) elements. Being a high level language its interactive mode of UI, via a prompt, follows the same principle as that of Matlab, and shares the same advantages for quick prototyping, querying and scripting of complex object manipulations and computations.

For more about what Python has to offer in scientific computing see the PythonResources page or a page at the [WWW] python wiki.

6. Design philosophy

In our design we have emphasized modularized data structures and interface design that facilitates data-driven approaches to the modeling of physical processes, and we have built upon standard numerical, scientific and graphics libraries for Python (for instance, SciPy and MatPlotLib). These, in turn, make use of well-established and efficient legacy codes for numerical integration of ODEs, and for dealing with linear algebra, optimization, and root solving (for instance, the LAPACK and MINPACK Fortran libraries). These legacy codes are typically interfaced using SWIG. The low-level languages of these codes provide the computational speed that Python itself lacks, in the places for which computation is most intensive.

All of the code involved in the PyDSTool project is open source, and we have aimed to create as few dependencies on external software packages as possible. In particular, the package can be used with Microsoft Windows, Mac OS X and Linux machines.

On top of the third-party libraries we have added several new tools and capabilities. We have enhanced legacy numerical integration code for ordinary differential equations to perform various additional tasks of use in hybrid systems modeling, implemented at the C-code level for maximal efficiency. This includes supporting discrete event detection during dynamical evolution. Adding arbitrary user-specified event detection to a model permits ODEs and maps to be used in combination as "hybrid" dynamical systems. Also, the inclusion of data-based time series inputs to a dynamical system’s evolution equations is a feature that aids data-driven modeling.

Utilities have been added that allow the movement of data and model specifications both in and out of PyDSTool, for sharing in other software environments. As well as basic importing and exporting of numerical data as text files, this also includes more systems-level interfacing. For instance, a user can export a dynamical model’s definition to a Matlab environment in which Automatic Differentiation is available for parameter sensitivity calculations (via the package ADMC++). Also, PyDSTool can be interfaced with the systems biology modeling package SloppyCell, through which PyDSTool inherits access to an interface with the Systems Biology Markup Language (SBML) for model specification, and the LaTeX mathematical markup language. Further interfaces to packages are in active development, such as to the original DsTool and to other simulation tools such as XPP or NEURON.

7. Coming attractions

For the most up to date information about these, see the page ComingAttractions.

Support for delayed systems is currently limited to dirac delta function pulses, as part of a hybrid system (see HybridSystems). DDE integrators written by Hairer and Wanner follow the same interfacing pattern as our existing integrators and will be added soon.

Automatic differentiation and Taylor-series integration capabilities will be added once we can get a working interface to an open-source AD package such as ADOL-C. In the meantime, we provide a model export tool to interface PyDSTool models to work with the ADMC++ environment in Matlab, a package that was also developed at Cornell (by Eric Phipps). See [WWW] this paper about ADMC++ and its application to parameter estimation for periodic orbits.

Although adequately supported in the underlying numerics of Python and SciPy, PyDSTool does not presently support dynamical systems having phase spaces ranging over the complex numbers. The application of PyDSTool to date has been to certain types of model in the physical sciences which have not required expressions of complex dynamics. It is hoped that at some point this can be redressed with explicitly complex dynamical systems. In the meantime, complex dynamics must be represented using two-dimensional real number phase spaces, and complex arithmetic implemented using auxiliary user-defined functions on 2D reals.

We have successfully tested an update of our automatic generation and compilation of C-code, although it does not yet avoid Python's distutils package by using Scons or autoconf / automake. This updated improves the efficiency of C-based vector field support and will make them more easily configured in new installations of Python/SciPy/PyDSTool.

We look forward to adding a range of other standard tools that are available in other packages for dynamical systems, such as phase-plane analysis, Poincare maps, Lyapunov exponent calculations, support for symbolic dynamics and complex dynamics, BVP solvers, averaging tools, invariant manifold computations, and so on.

8. Historical context

Approximately ten years ago the original [WWW] DsTool package was written here at Cornell, and has provided many years of functionality to the applied mathematics community. We have decided to write a new package from the ground up, as our concept of what is needed in our community has evolved. Our new project is different from existing approaches because it attempts to combine a large array of tools under one, centralized umbrella. We wish to support both straightforward simulation tools, such as numerical integrators of differential equations, and sophisticated analysis tools such as continuators and parameter estimation / sensitivity tools.

Individually, several examples of the software tools to implement these algorithms exist in the world of scientific computing. However, there has not yet been a one-stop solution to providing mathematically sophisticated users with a combined suite of such tools that are highly integrated and interactively operated. For instance, Matlab is geared much less towards applied mathematicians in this regard as it is to physicists and engineers. It is still typical for applied scientists to need to use a variety of separate programs to achieve their dynamical systems modeling goals, but increasingly sophisticated questions are being asked of the models, which are ever-harder to answer when the tools are not integrated. Example domains that will benefit from PyDSTool are computational neuroscience, biomechanics and robotics, and advanced engineering projects.

9. Comparison to other packages

In the broadest terms, we believe PyDSTool combines or extends a range of features from packages such as Matlab, XPP, AUTO, and NEURON into an integrated and open environment.

9.1. General packages

9.1.1. Matlab / Simulink

[WWW] Matlab and Simulink.

See also [WWW] this comparison of the NumPy library to Matlab.

9.1.2. XPP

[WWW] XPP (X-PhasePlane)

9.1.3. AUTO and MATCONT

[WWW] AUTO, [WWW] Matcont.

9.1.4. Custom code in other languages

9.2. Application-specific packages


[WWW] Neuron at Duke, [WWW] Neuron at Yale,[WWW] Genesis.

PyDSTool source code is hosted by SourceForge:

Get PyDSTool at Fast, secure and Free Open Source software downloads