程序的一部分, 目標(biāo)是求解一個(gè)類似泊松方程的場(chǎng)方程。原作者不是我, 程序?qū)懙煤苤苯,但很緩慢。?jīng)過(guò)我修改后, 速度雖然提上去了,但貌似兩個(gè)程序的結(jié)果不是很對(duì),按說(shuō)我都使用的相同的算法, 區(qū)別只在寫法上。
方程:
周期性邊界條件。
原始代碼: 注意這里的P1只是給了一個(gè)形式, 為的是說(shuō)明他的結(jié)構(gòu), 計(jì)算中的P1是有一個(gè)分布的,不是零。
修改后的代碼: 原始代碼迭代次數(shù)為498次,修改后要多迭代100次,而且結(jié)果還對(duì)不上。注意程序里的C1,C2,C3和方程里的不是一回事。
問(wèn)題:
1、第二段程序什么地方不正確,導(dǎo)致結(jié)果不同。
2、這樣的方程是不是有什么方法使用庫(kù)函數(shù)求解。
3、這種迭代格式是否正確。
4、還沒(méi)想到…… |