|
|
发表于 2003-10-13 16:19:10
|
显示全部楼层
Clawpack的问题
[这个贴子最后由supersonic在 2003/10/14 06:05pm 第 2 次编辑]
也许可以。clawpack是一个通用双曲系统求解器。关键在于方程和相应的黎曼求解方法,他们都是是可以自己构造的。你可以看看
http://www.math.tu-cottbus.de/~deiter/
http://www.math.tu-cottbus.de/~deiter/amroc/index.htm
这个是以clawpack为基础的,它的代码包里面对clawpack的源代码文件做了详尽的分类,其中专门包含一个equations目录,我想修改这个估计可以吧。另外source相关的函数也要改。
作者就是做爆炸模拟的,当然有化学反应。不过AMROC里面似乎没有实现多组分化学反应方程。你得自己写。
clawpack精度不低了,如果需要提高精度,在clawpack框架内修改比自己从头编写应该简单些。
不过clawpack程序结构不算很简单,特别是数据管理部分。由于f77没有动态内存分配功能,clawpack内建有一套存储管理机制,有些让人晕。
纯属建议,自己拿决定吧。 |
|