| 5 | 5/1 | 返回列表 |
| 查看: 2516 | 回復(fù): 50 | |||
| 本帖產(chǎn)生 1 個(gè) 數(shù)學(xué)EPI ,點(diǎn)擊這里進(jìn)行查看 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
xjw0413銅蟲(chóng) (初入文壇)
|
[交流]
【討論】預(yù)調(diào)件共軛梯度法(PCG) 已有4人參與
|
||
|
有限元計(jì)算經(jīng)常碰到大型稀疏矩陣,由于此類(lèi)線性方程組通常條件數(shù)是比較大的,方程組的性態(tài)不好,所以最好用迭代方法求解,比方說(shuō)是預(yù)調(diào)件共軛梯度法,但此方法在選擇預(yù)調(diào)件矩陣時(shí)似乎沒(méi)有一個(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不就解出來(lái)了嗎,又何必再來(lái)PCG迭代呢? 2. 采用PCG方法的前提應(yīng)該是系數(shù)矩陣對(duì)稱(chēng)、正定吧,因?yàn)槠湓硎且粋(gè)相當(dāng)于勢(shì)函數(shù)的東西取極小值。那對(duì)于非正定的系數(shù)矩陣能求解嗎,我構(gòu)造了幾個(gè)非正定的,有的似乎是能夠收斂到正確結(jié)果的。 希望各位蟲(chóng)用解答和討論。 |

用戶注銷(xiāo) (知名作家)
鐵桿木蟲(chóng) (正式寫(xiě)手)
銅蟲(chóng) (初入文壇)

鐵桿木蟲(chóng) (正式寫(xiě)手)
|
在特定精度下終止程序?qū)嶋H上就是在一定的子空間內(nèi)取得了解的近似。我的意思是你得到看似正確的解可能是因?yàn)檫@個(gè)子空間還沒(méi)有包含任何系數(shù)矩陣的負(fù)曲率方向,這種情況下能得到滿足精度要求的一個(gè)近似解似乎也很自然。一旦偵測(cè)到一個(gè)負(fù)的方向,迭代序列很可能因此發(fā)散。具體的表現(xiàn)是什么,還是需要比較嚴(yán)謹(jǐn)?shù)耐茖?dǎo)。總之這種情況下用CG得到的結(jié)果是沒(méi)有說(shuō)服力的,這個(gè)不是改變終止準(zhǔn)則就能修正的。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-16 | 3/150 |
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
[考研] 290求調(diào)劑 +5 | 孔志浩 2026-03-12 | 10/500 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-15 | 6/300 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 材料工程327求調(diào)劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +5 | ZY-05 2026-03-13 | 6/300 |
|
|
[考研] 341求調(diào)劑 +4 | 番茄頭--- 2026-03-10 | 4/200 |
|
|
[考研] 考研調(diào)劑 +4 | 芬達(dá)46 2026-03-12 | 4/200 |
|
|
[論文投稿]
投稿問(wèn)題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|
|
[考博] 2026年博士申請(qǐng) +3 | QwQwQW10 2026-03-11 | 3/150 |
|
|
[考研] 調(diào)劑 +5 | 呵唔哦豁 2026-03-10 | 5/250 |
|