FLASH Code傲娇记
近日因为要计算相对论性激波的缘故，需要进行数值模拟，于是就面临一个选择：是自力更生编写代码，还是直接拿来现有程序？考虑一开始要求解的问题还算简单，使用现有大型模拟程序颇有牛刀杀鸡的感觉，而且这些大型程序的学习使用也非易事，有必要专门为此抽出两三个月的时间，本人最初是计划走前一路线的，相关代码主体也很快就编写完毕了。但是数值方法多少会带来误差，就本人所用的算法而言，人工耗散的懒往往是不可以偷的。不过相对论性流体基本方程组在分离变量之后造型奇特且复杂，自己一下子就因为耗散的添加方法而卡了壳。这时候合作导师又提出了新的课题，必须通过复杂的数值模拟工作才能解决，于是心一横，还是拿来主义好。 (more…)
