Python：如何用Jacobian有效解决DAE？

问题描述

我正在尝试使用Assimulo软件包来求解一组微分代数方程（DAE）。我正在尝试使用类似于shown here的算法。但是，似乎没有传递稀疏矩阵的选项。我的Jacobian矩阵非常大，大约为3000 x3000。您知道是否有一种方法可以更有效地解决我的DAE？

解决方法

根据我使用稀疏 ODE 系统（更准确地说是使用半离散 PDE 系统）的经验，使用迭代线性求解器可以大大提高数值效率。据我所知，Assimulo 不允许提供雅可比稀疏模式，但更改线性求解器是解决此问题的另一种方法。

你会做这样的事情：

model = Explicit_Problem(ode_function,y0=y_init,t0=t_init)
simulator = CVode(model)
sim.linear_solver = 'SPGMR'

我不确定这是否也适用于 DAE 系统，但我认为值得一试。