程序的一部分, 目标是求解一个类似泊松方程的场方程。原作者不是我, 程序写得很直接,但很缓慢。经过我修改后, 速度虽然提上去了,但貌似两个程序的结果不是很对,按说我都使用的相同的算法, 区别只在写法上。
方程:
周期性边界条件。
原始代码: 注意这里的P1只是给了一个形式, 为的是说明他的结构, 计算中的P1是有一个分布的,不是零。
修改后的代码: 原始代码迭代次数为498次,修改后要多迭代100次,而且结果还对不上。注意程序里的C1,C2,C3和方程里的不是一回事。
问题:
1、第二段程序什么地方不正确,导致结果不同。
2、这样的方程是不是有什么方法使用库函数求解。
3、这种迭代格式是否正确。
4、还没想到…… |