Re: [問題] 請教矩陣LU分解法

看板Programming作者時間18年前 (2006/11/27 02:01), 編輯推噓0(000)
留言0則, 0人參與, 最新討論串1/1
==> 在 bmw520i.bbs@ptt.cc (燃燒我的青春) 的文章中提到: > 最近上課上到關於矩陣解法的LU分解, > 我想這是滿普遍存在於一般數值分析教材, > 關於整個過程的演算程序我也知道, > 因為作業有個問題是求解二維Laplace方程式。 > 之前曾做過另一個解法是使用"三對角矩陣"解法, > 把方程式經過排列以"逐線"方式求解, > 最近要以LU分解去做該問題,在想LU分解若是以逐線方式求解去做並不難, > 但如此作法好像就跟"三對角矩陣"解法的精神似乎差不多!? > 可是若是"全域"方式去解,即一口氣解一個平面上所有變數, > 好像比較符合LU分解所謂一次求解的精神?! > 但是若是以這樣方式作,好像變數排列問題, > 在程式撰寫上是不是比較麻煩呢? > 嘗試以一個小維度矩陣去手寫看其排列, > 每條方程式求解變數並非完全有一規則, > 不知道說LU分解法在一般問題應用時,也是屬於以"逐線"方式求解嗎? > ps.聽老師曾講過是說,在大型平行電腦計算,以LU分解法為佳, > 若是在一般PC上,LU分解法會儲存龐大陣列資料, > 所以我會想說LU分解法在應用時,是"逐線"或是"全區域"去解。 我想在你的 matrix size 很小的時候,用 LU 分解反而花了更多時間 但是當 matrix size 大到一定程度 (ex. 100 x 100) LU 的方法一定比較快吧(純理論推測) 我沒有實際寫程式實驗這兩者的差別 也許你可以把 matrix size 加大,比較兩種 algorithm 的差別 不過根據之前的經驗,當 size 超過 15 x 15,程式跑起來 就開始懷疑自己是不是寫了一個 dead loop 應該是之前的寫法有點問題+電腦太爛吧 -- * Origin: ★ 交通大學資訊科學系 BBS ★ <bbs.cis.nctu.edu.tw: 140.113.23.3>
文章代碼(AID): #15QTPa00 (Programming)
文章代碼(AID): #15QTPa00 (Programming)