 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2519 | 回復(fù): 50 | |||
| 本帖產(chǎn)生 1 個(gè) 數(shù)學(xué)EPI ，點(diǎn)擊這里進(jìn)行查看 | |||
xjw0413銅蟲 (初入文壇)
|
[交流]
【討論】預(yù)調(diào)件共軛梯度法（PCG） 已有4人參與
|
||
|
有限元計(jì)算經(jīng)常碰到大型稀疏矩陣，由于此類線性方程組通常條件數(shù)是比較大的，方程組的性態(tài)不好，所以最好用迭代方法求解，比方說是預(yù)調(diào)件共軛梯度法，但此方法在選擇預(yù)調(diào)件矩陣時(shí)似乎沒有一個(gè)同一的標(biāo)準(zhǔn)，大多推薦的是采用incomplete LU decomposition做為預(yù)調(diào)件矩陣。incomplete LU decomposition的計(jì)算方法似乎又有很多種。 1. incomplete LU decomposition 的計(jì)算時(shí)間應(yīng)該比 LU decomposition要快速的多吧，不然直接用LU decomposition不就解出來了嗎，又何必再來PCG迭代呢？ 2. 采用PCG方法的前提應(yīng)該是系數(shù)矩陣對稱、正定吧，因?yàn)槠湓硎且粋�(gè)相當(dāng)于勢函數(shù)的東西取極小值。那對于非正定的系數(shù)矩陣能求解嗎，我構(gòu)造了幾個(gè)非正定的，有的似乎是能夠收斂到正確結(jié)果的。 希望各位蟲用解答和討論。 |
回復(fù)此樓» 猜你喜歡
268求調(diào)劑
已經(jīng)有9人回復(fù)
-
301求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)336分求調(diào)劑
已經(jīng)有5人回復(fù)
-
301求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有9人回復(fù)
-
304求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有5人回復(fù)
-
344求調(diào)劑
已經(jīng)有4人回復(fù)
-
326求調(diào)劑
已經(jīng)有6人回復(fù)
-
296求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 關(guān)于渤海大學(xué)應(yīng)用化學(xué)研究生調(diào)劑
已經(jīng)有13人回復(fù)
- 不知道選擇調(diào)劑還是再考一年
已經(jīng)有47人回復(fù)
- 面上基金項(xiàng)目預(yù)算里"勞務(wù)費(fèi)"自動(dòng)生成的 "博士后"能刪掉嗎?
已經(jīng)有5人回復(fù)
- [2012-04-15]第12屆多尺度、多功能梯度材料國際會(huì)議（12th International Symposium o
已經(jīng)有56人回復(fù)
- 如果用高級氧化方法處理生活污水，都需要測什么指標(biāo)呢？
已經(jīng)有28人回復(fù)
- 溶脹平衡法測定交聯(lián)聚合物的交聯(lián)密度
已經(jīng)有433人回復(fù)
- 關(guān)于 共軛梯度法 求助～～～好心人幫忙 有金幣的哦
已經(jīng)有11人回復(fù)
- 請教：Anglient 1200 數(shù)據(jù)分析方法
已經(jīng)有12人回復(fù)
- 【求助】求四氫呋喃的除水干燥方法
已經(jīng)有20人回復(fù)
- 求助N,N-二甲基乙酰胺縮二甲醇的合成方法，金幣可加，還有沒有毒性小一點(diǎn)的方法。
已經(jīng)有8人回復(fù)
- 【求助】求共軛梯度法中的Armijo線性搜索程序
已經(jīng)有5人回復(fù)
匿名
用戶注銷 (知名作家)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3228.2
- 散金: 24
- 紅花: 1
- 帖子: 6711
- 在線: 384.6小時(shí)
- 蟲號: 0
- 注冊: 2009-11-17
- 專業(yè): 大氣環(huán)境與全球氣候變化
27樓2011-11-14 10:10:48
匿名
用戶注銷 (知名作家)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3228.2
- 散金: 24
- 紅花: 1
- 帖子: 6711
- 在線: 384.6小時(shí)
- 蟲號: 0
- 注冊: 2009-11-17
- 專業(yè): 大氣環(huán)境與全球氣候變化
49樓2012-09-02 14:29:27
saladin983
鐵桿木蟲 (正式寫手)
- 數(shù)學(xué)EPI: 4
- 應(yīng)助: 6 (幼兒園)
- 貴賓: 0.2
- 金幣: 9197.9
- 紅花: 3
- 帖子: 687
- 在線: 115.7小時(shí)
- 蟲號: 448168
- 注冊: 2007-11-01
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
2樓2010-09-06 22:58:01
xjw0413
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 224.1
- 帖子: 41
- 在線: 11.5小時(shí)
- 蟲號: 1073607
- 注冊: 2010-08-12
- 性別: GG
- 專業(yè): 固體力學(xué)
3樓2010-09-07 14:04:33
saladin983
鐵桿木蟲 (正式寫手)
- 數(shù)學(xué)EPI: 4
- 應(yīng)助: 6 (幼兒園)
- 貴賓: 0.2
- 金幣: 9197.9
- 紅花: 3
- 帖子: 687
- 在線: 115.7小時(shí)
- 蟲號: 448168
- 注冊: 2007-11-01
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
lovibond(金幣+1): 鼓勵(lì)耐心交流 2011-06-18 10:35:05
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
lovibond(金幣+1): 鼓勵(lì)耐心交流 2011-06-18 10:35:05
|
在特定精度下終止程序?qū)嶋H上就是在一定的子空間內(nèi)取得了解的近似。我的意思是你得到看似正確的解可能是因?yàn)檫@個(gè)子空間還沒有包含任何系數(shù)矩陣的負(fù)曲率方向，這種情況下能得到滿足精度要求的一個(gè)近似解似乎也很自然。一旦偵測到一個(gè)負(fù)的方向，迭代序列很可能因此發(fā)散。具體的表現(xiàn)是什么，還是需要比較嚴(yán)謹(jǐn)?shù)耐茖?dǎo)。總之這種情況下用CG得到的結(jié)果是沒有說服力的，這個(gè)不是改變終止準(zhǔn)則就能修正的。 |
4樓2010-09-07 16:06:00
匿名
用戶注銷 (知名作家)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3228.2
- 散金: 24
- 紅花: 1
- 帖子: 6711
- 在線: 384.6小時(shí)
- 蟲號: 0
- 注冊: 2009-11-17
- 專業(yè): 大氣環(huán)境與全球氣候變化
5樓2011-06-14 10:37:26
saladin983
鐵桿木蟲 (正式寫手)
- 數(shù)學(xué)EPI: 4
- 應(yīng)助: 6 (幼兒園)
- 貴賓: 0.2
- 金幣: 9197.9
- 紅花: 3
- 帖子: 687
- 在線: 115.7小時(shí)
- 蟲號: 448168
- 注冊: 2007-11-01
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
6樓2011-06-15 03:32:09
匿名
用戶注銷 (知名作家)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3228.2
- 散金: 24
- 紅花: 1
- 帖子: 6711
- 在線: 384.6小時(shí)
- 蟲號: 0
- 注冊: 2009-11-17
- 專業(yè): 大氣環(huán)境與全球氣候變化
7樓2011-06-15 08:57:22
saladin983
鐵桿木蟲 (正式寫手)
- 數(shù)學(xué)EPI: 4
- 應(yīng)助: 6 (幼兒園)
- 貴賓: 0.2
- 金幣: 9197.9
- 紅花: 3
- 帖子: 687
- 在線: 115.7小時(shí)
- 蟲號: 448168
- 注冊: 2007-11-01
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖
nono2009(金幣+2): 鼓勵(lì)應(yīng)助 2011-07-01 07:10:32
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖
nono2009(金幣+2): 鼓勵(lì)應(yīng)助 2011-07-01 07:10:32
|
通常意義下的PCG，preconditioner只要求對稱正定，這里的CC(i)只要都是正數(shù)，顯然P可以用作preconditioner。但是需要知道的是，對于不同的線性方程組，preconditioner只能依照系數(shù)矩陣來做具體的選擇。這里為什么這么用，針對這個(gè)代碼要解決的問題，你可以看看P的應(yīng)用有沒有達(dá)到這樣的效果： (1) 使系數(shù)矩陣的特征值聚集； (2) 改善條件數(shù)。 因?yàn)镻是個(gè)簡單的對角矩陣，除非系數(shù)矩陣也近似于相關(guān)的對角矩陣，否則，我更偏向于(2)。 另外，你可以讓P=I（單位矩陣），這樣PCG就變成標(biāo)準(zhǔn)CG了，可以對比一下，這個(gè)P有沒有改善你的計(jì)算結(jié)果。 |
8樓2011-06-18 00:29:42
匿名
用戶注銷 (知名作家)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3228.2
- 散金: 24
- 紅花: 1
- 帖子: 6711
- 在線: 384.6小時(shí)
- 蟲號: 0
- 注冊: 2009-11-17
- 專業(yè): 大氣環(huán)境與全球氣候變化
9樓2011-06-18 08:58:08
saladin983
鐵桿木蟲 (正式寫手)
- 數(shù)學(xué)EPI: 4
- 應(yīng)助: 6 (幼兒園)
- 貴賓: 0.2
- 金幣: 9197.9
- 紅花: 3
- 帖子: 687
- 在線: 115.7小時(shí)
- 蟲號: 448168
- 注冊: 2007-11-01
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖
nono2009(金幣+1, 數(shù)學(xué)EPI+1): 鼓勵(lì)應(yīng)助 2011-07-01 07:11:01
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖
nono2009(金幣+1, 數(shù)學(xué)EPI+1): 鼓勵(lì)應(yīng)助 2011-07-01 07:11:01
|
注意一點(diǎn)，通常選擇preconditioner的時(shí)候一個(gè)常規(guī)的做法就是選擇一個(gè)系數(shù)矩陣的近似，這樣的話，條件預(yù)優(yōu)后得到的矩陣就會(huì)接近于單位矩陣。 假定使用LU分解，然后用L和U的乘積作為P，那么實(shí)際上就是就是用系數(shù)矩陣自身作為P，亦即條件預(yù)優(yōu)后的矩陣嚴(yán)格地等于單位矩陣。這是最理想的情況，但是，事實(shí)上這么做沒有任何意義，因?yàn)橐呀?jīng)有LU分解，那么可以直接用這個(gè)分解求解原方程組，只需要簡單的向前向后代入就足夠，無需使用CG這樣的迭代方法。 這里使用ILU（不完全的LU分解）是因?yàn)長U分解代價(jià)太大，尤其是對于大規(guī)模線性方程組而言。ILU分解得到的L和U作乘積之后可以充當(dāng)系數(shù)矩陣的近似，按照之前所說的，也就是P的一種理想選擇。 你可以看看ILU分解的具體做法，Yousef Saad的Iterative methods for sparse linear systems一書里面有專門的章節(jié)介紹，有空可以看看。這是本很好的書，在作者的個(gè)人主頁可以直接下載（http://www-users.cs.umn.edu/~saad/PS/all_pdf.zip）。 |
10樓2011-06-18 15:21:12
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xjw0413 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【0856】化學(xué)工程（085602）313 分，本科學(xué)科評估A類院�；瘜W(xué)工程與工藝，誠求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 303求調(diào)劑 +3 | 睿08 2026-03-17 | 3/150 |
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
|
[考研] 材料工程327求調(diào)劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
|
[考研] 材料371求調(diào)劑 +9 | 鱷魚? 2026-03-11 | 11/550 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 快樂的香蕉 2026-03-11 | 4/200 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] 求調(diào)劑（材料與化工327） +4 | 愛吃香菜啦 2026-03-11 | 4/200 |
|
|
|
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|