41得票1回答
收到"Signal 15"是什么意思?

可能是什么原因导致使用一个名为SUNDIALS/CVODE的库(一种数值ODE求解器)在Gentoo Linux集群上运行的C、MPI程序给我反复发出Signal 15 received.? 这个代码是由MPI、Sundials、Linux、C或其他谁发出的? 请注意,我对以下技术基本上是...

27得票5回答
使用自适应步长来处理scipy.integrate.ode

scipy.integrate.ode的(简要)文档称,有两种方法(dopri5和dop853)具有步长控制和密集输出。 通过查看示例和代码本身,我只能看到从积分器获取输出的非常简单的方法。 也就是说,您只需向前移动积分器一些固定的dt,获取该时间的函数值,然后重复。 我的问题具有相当可变的...

22得票6回答
传递参数给 solve_ivp(新的 SciPy ODE API)

要使用SciPy解决简单的ODE问题,我过去常常使用odeint函数,格式如下: scipy.integrate.odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rto...

18得票2回答
scipy.integrate.odeint和scipy.integrate.ode之间有什么区别?

我已经阅读了这些此处的文档,但它只是说明odeint是“基于ODEPACK中lsoda的简单接口的积分器”。 实际上有什么区别?在什么情况下使用一个比另一个更合适?

15得票2回答
有没有适用于普通微分方程(ODE)求解器的C++库?

更具体地说,我对8阶Dormand-Prince嵌入式方法感兴趣,它基于龙格-库塔法和刚性方程。 我使用Numerical Recipes 3,但我经常在编译其库时遇到麻烦。我想知道是否有其他替代方案。

13得票2回答
多个scipy.integrate.ode实例

我想在多个线程中使用scipy.integrate.ode(或scipy.integrate.odeint)实例,每个线程都使用一个CPU核心,以便同时解决多个IVP问题。然而文档上说:“这个积分器不可重入。您不能同时使用两个使用“vode”积分器的ode实例。 (同时,如果多次实例化ode...

11得票1回答
在Java中的微分方程

我正在尝试用Java创建一个简单的SIR流行病模型的仿真程序。 基本上,SIR由三个微分方程组成: S'(t) = - l(t) * S(t) I'(t) = l(t) * S(t) - g(t) * I(t) R'(t) = g(t) * I(t) 其中,S - 易感人群,I - 感染者...

10得票2回答
奇怪的SciPy ODE集成错误

我正在实现一个非常简单的易感-感染-恢复模型,该模型针对一个稳定人口进行建模,这是一个相当琐碎的任务。但是,我在使用PysCeS或SciPy时遇到了求解器错误,它们都使用lsoda作为其基础求解器。这仅发生在特定参数值下,我对此感到困惑。我使用的代码如下: import numpy as n...

9得票2回答
用零替换模型(ODE系统)中的负值

我目前正在使用 deSolve 解决一组普通微分方程组,想知道是否有任何方法可以防止微分变量值低于零。我看到了一些其他帖子关于在向量、数据框等中将负值设置为零,但由于这是生物模型(T细胞计数不可能为负),我需要阻止它发生,以便这些值不会影响结果,而不仅仅是在最终输出中替换负数。

9得票4回答
使用scipy.integrate.odeint求解一个带有可变常数的ODE系统?

我目前有一个带有时间变量的常数的ODE系统。例如: def fun(u, t, a, b, c): x = u[0] y = u[1] z = u[2] dx_dt = a * x + y * z dy_dt = b * (y-z) dz_dt...