使用solve_ivp求解ODE时并行计算
我有一台带有64GB RAM的i7 6540HQ计算机。我正在尝试解决一个PDE，该PDE转换为具有100个离散点的ODE，如下代码所示，总计1002个方程。当我检查任务管…
我如何避免在numpy中具有nan值
我有一个非线性的PDE系统，并且我试图使用有限差异方法来解决它。但是我一直遇到以下错误，并且图形不正确： Runtime Warning：在Double_scalars中遇…
在Gekko上使用scipy.integrate.trapz时出错
我正在尝试使用Gekko库解决PDE，该PDE涉及通过在一定时间与空间集成速度来计算体积流量。 Qgap==scipy.integrate.trapz(2*np.pi*xpos*np.array(u1)) …
如何在Java或PDE中打印实例的内容？
class bouncingBall(): def __init__(self, bounce, window, position): self.bounce = bounce self.window = window self.position = position ball …
如何在广义加权剩余方法（Chebyshev Galerkin Modal形式）中施加边界条件
我正在努力将广义加权剩余方法与RK4相结合。 GWRM零件将PDE分解到光谱空间，在该空间中，未知数是Chebyshev系数A_K。但是，我很难看到如何在这种情况…
解决1D Navier Stokes具有可压缩质量保护（液压阻尼器）的问题
我想在圆柱施加的管（笛卡尔山脉）上求解一个1D Navier方程。 沿y方向的流动为p1，右腔室最初都在8e5 pa处。圆柱壁是固定的，宽度'L'的活塞正在用A移…
用Gekko Python求解部分微分方程的问题
我在尝试求解下面附加的部分微分方程时得到了收敛的解决方案。 在我的代码中，我想通过集成2 pi r*v（r）*DR来计算量流量。我已经使用scipy.integrate…
如何在 py-pde (python) 中输入一般非线性 pde？
我正在尝试使用 py-pde 来模拟 2d 图 z=u(x,y) 的平均曲率流。 平均曲率流方程有一个函数 u 的 Hessian 项 请参阅此处的方程（D_i 表示相对于 x 和 y …
在 FiPy 中求解总体平衡 PDE 时遇到问题
我是使用FIPY以及通常使用有限体积方法的初学者。因此，如果我的问题很幼稚，请原谅我 ，我试图解决形式的人口平衡方程 pde 边界条件如下 边界条件 x…
将MATLAB代码转换用于使用Kronecker产品生成2D Laplacian矩阵的代码
以下MATLAB代码使用Kronecker产品生成2D Laplacian矩阵 方法。 function A=A(N) % Assemble the system matrix A e = ones(N,1); D = spdiags([e -2*e…
所需插件“org.eclipse.pde.junit.runtime”找不到
我正在使用 Eclipse Indigo，需要使用 Junit PDE 来测试我的插件功能。我安装了 Junit 和这两个插件： http://market.eclipsesource.com/yoxos/node/j…
AJDT 无头 PDE 构建期间未设置 ajdt.pdebuild.scripts
我有一个在 eclipse 3.5 中开发的 eclipse rcp 应用程序。我能够通过以下目标条目在 ant 中成功执行 pde headless build（从 eclipse 外部的命令 shel…
