[SOSlib-discuss] Identifying algebraic parts of DAE system

[Randall B. <r.b...@au...>]

Hi

 

Andrew Miller has created code for decomposing a system of DAE's into
possibly many algebraic subsystems and one differential system, which can
then be solved in a decoupled manner.

 

This code is part of the open source CellML api, I'm not sure if it would be
useful to you, but it seems related to a comment on
http://www.tbi.univie.ac.at/~raim/odeSolver/doc/api/group__ida.html: "This
code is not working yet. It is not compiled with the current package. It is
included in the documentation merely to motivate people to help us implement
this functionality. The main problem is that ODE construction currently
can't decide which variables are to be defined via algebraic constraints.
Contact us, if you want to help!"

 

The code is available from
https://svn.physiomeproject.org/svn/physiome/CellML_DOM_API/

 

A visual diff showing the relevant code is visible at
http://svnviewer.bioeng.auckland.ac.nz/projects/physiome/changeset/2344,
though there have been minor changes to it since.

 

Regards,

 

Randall Britten

___________________________

Software Development Leader

Bioengineering Institute

Auckland University