python - 显式方案 RuntimeWarning:长标量中遇到溢出
问题描述
今天是个好日子。
由于代码太大,我附上repl的链接:
https://replit.com/@PabloJunior1/Explicit-scheme-ENG
我正在尝试为二维方程实现显式差分方案。我处理了所有公式并且程序运行良好并且在 256 次迭代后收敛,如果我使用时间步长 τ = 0.25 * h ^ 2。但是如果我使用步长 τ = h ^ 2,那么在大约第 370 次迭代之后这个错误出去,这个过程进入一个永恒的循环。不幸的是,这是需要使用的值。.
我会注意到,在当前参数 N1 = N2 = 10 和误差精度 eps = 10 ^ (- 5) 的情况下,步长 τ = h ^ 2 程序应该在大约 75-85 步内收敛。我用隐式方案和交替方向隐式方案检查了它。如有必要,我可以附上原始公式。
此外,我已经尝试将 float64 类型分配给所有数组,正如在类似问题中所建议的那样。但是我的代码仍然不起作用。我不知道如何解决它。
解决方案
推荐阅读
- c# - 如何在 Oracle 身份验证中使用 MS Access 连接字符串
- reactjs - 除非我在 React 中刷新页面,否则状态不会改变
- html - 在自定义滑块 Jquery 上隐藏/显示下一个和上一个按钮
- python - 如何在函数中检测变量赋值?
- typescript - 如果仅定义了字符串,则 Nest.js 验证字符串数组
- reactjs - 为什么 Material-UI 日期与占位符重叠值并抱怨受控输入?
- r - 如何使用 errbar 从数据集中绘制多列
- python - 如何使用 isinstance() 根据对象的类型将列表分成两个列表?
- python - 仍然可以安装 Python 2 以与 Spyder 一起使用吗?
- android - 在不调用 recompose 的情况下,如何更新 jetpack 中的特定 ui 组件?