使用 sfepy 求解简单的二维微分方程

Question

我正在努力学习 sfepy。 为此，我想求解微分方程

在三角域 (2D) 上。

阅读文档后有两件事我不明白。

1. 如何指定单个三角形网格？ 该代码似乎假设您已经拥有网格文件，但它没有提供指向它的链接，所以我不知道如何构造我的数据。

2. 我很难学习如何将等式 map 应用到 sfepy 的语法中。 即我不知道如何使用库指定问题，即使在遵循教程之后也是如此。

Answer 1

我有解决你的问题 1。

示例脚本linear_elastic_interactive.py确实提到了第 35 行使用的网格的路径，例如mesh = Mesh.from_file(data_dir + '/meshes/2d/rectangle_tri.mesh') 。 教程还提到“这个文件应该从顶级 SfePy 源目录运行，这样它才能正确找到网格文件”。 所以如果你浏览源代码，你可以在这里找到网格： 链接