今天对程序大修了一下,顺便把所有算例测试了一遍,突然发现二维浅水方程有些算例出现了明显的错误. 这次突然出现的错误让我有点措手不及,因为一直没有修改过浅水方程求解器,所以这些算例很久没有测试过了.硬着头皮从斜率限制器,数值通量,计算项等等所有计算函数里面寻找 bug,但是找了一下午仍然一无所获.看哪里都好像不太对,但是不管怎么改错误仍然会出现. 好在上个月我写算例时候当时是能够正确计算的,于是利用 git reset 命令开始返回过去提交的代码版本,开始一个个测试,最后终于在6月2号版本程序里得…