[問題] 一維Ising model
開發平台(Platform): (Ex: VC++, GCC, Linux, ...)
Linux
額外使用到的函數庫(Library Used): (Ex: OpenGL, ...)
無
問題(Question)：
估計是61~68行，算trial Hamiltonian時發生錯誤
同樣計算 trail Hamiltonian如果以暴力法算，最後結果是正確的
當我想要避免用暴力法時，誤差變得非常大
一維Ising model大致上考慮兩個邊界問題就可以
H = -Σ S * S 此公式是計算Hamiltonian的公式
i i+1
計算 trail Hamiltonian時我考慮改變隨機選取的S 以及前後項
x
但是結果無法得到我所預期的結果
餵入的資料(Input)：
kT 可以是不大於20的正整數
預期的正確結果(Expected Output)：
輸出兩個值中，後者是理論預測值
前者是用蒙地卡羅方法計算的，應該要接近理論值
誤差大約在 2~3左右
錯誤結果(Wrong Output)：
算出的結果離理論值非常多
ex: -173.219 -197.375 when kT=5
程式碼(Code)：(請善用置底文網頁, 記得排版)
http://dpaste.com/731578/
補充說明(Supplement)：
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.167.79.118
C_and_CPP 近期熱門文章
PTT數位生活區 即時熱門文章
