 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2520 | 回復(fù): 50 | |||
| 本帖產(chǎn)生 1 個 數(shù)學(xué)EPI ，點(diǎn)擊這里進(jìn)行查看 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
xjw0413銅蟲 (初入文壇)
|
[交流]
【討論】預(yù)調(diào)件共軛梯度法（PCG） 已有4人參與
|
||
|
有限元計(jì)算經(jīng)常碰到大型稀疏矩陣，由于此類線性方程組通常條件數(shù)是比較大的，方程組的性態(tài)不好，所以最好用迭代方法求解，比方說是預(yù)調(diào)件共軛梯度法，但此方法在選擇預(yù)調(diào)件矩陣時似乎沒有一個同一的標(biāo)準(zhǔn)，大多推薦的是采用incomplete LU decomposition做為預(yù)調(diào)件矩陣。incomplete LU decomposition的計(jì)算方法似乎又有很多種。 1. incomplete LU decomposition 的計(jì)算時間應(yīng)該比 LU decomposition要快速的多吧，不然直接用LU decomposition不就解出來了嗎，又何必再來PCG迭代呢？ 2. 采用PCG方法的前提應(yīng)該是系數(shù)矩陣對稱、正定吧，因?yàn)槠湓硎且粋�相當(dāng)于勢函數(shù)的東西取極小值。那對于非正定的系數(shù)矩陣能求解嗎，我構(gòu)造了幾個非正定的，有的似乎是能夠收斂到正確結(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)"自動生成的 "博士后"能刪掉嗎?
已經(jīng)有5人回復(fù)
- [2012-04-15]第12屆多尺度、多功能梯度材料國際會議（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ù)
saladin983
鐵桿木蟲 (正式寫手)
- 數(shù)學(xué)EPI: 4
- 應(yīng)助: 6 (幼兒園)
- 貴賓: 0.2
- 金幣: 9197.9
- 紅花: 3
- 帖子: 687
- 在線: 115.7小時
- 蟲號: 448168
- 注冊: 2007-11-01
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
nono2009(金幣+2): 鼓勵應(yīng)助 2011-07-01 07:10:32
小木蟲(金幣+0.5):給個紅包，謝謝回帖
nono2009(金幣+2): 鼓勵應(yīng)助 2011-07-01 07:10:32
|
通常意義下的PCG，preconditioner只要求對稱正定，這里的CC(i)只要都是正數(shù)，顯然P可以用作preconditioner。但是需要知道的是，對于不同的線性方程組，preconditioner只能依照系數(shù)矩陣來做具體的選擇。這里為什么這么用，針對這個代碼要解決的問題，你可以看看P的應(yīng)用有沒有達(dá)到這樣的效果： (1) 使系數(shù)矩陣的特征值聚集； (2) 改善條件數(shù)。 因?yàn)镻是個簡單的對角矩陣，除非系數(shù)矩陣也近似于相關(guān)的對角矩陣，否則，我更偏向于(2)。 另外，你可以讓P=I（單位矩陣），這樣PCG就變成標(biāo)準(zhǔn)CG了，可以對比一下，這個P有沒有改善你的計(jì)算結(jié)果。 |
8樓2011-06-18 00:29:42
saladin983
鐵桿木蟲 (正式寫手)
- 數(shù)學(xué)EPI: 4
- 應(yīng)助: 6 (幼兒園)
- 貴賓: 0.2
- 金幣: 9197.9
- 紅花: 3
- 帖子: 687
- 在線: 115.7小時
- 蟲號: 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小時
- 蟲號: 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小時
- 蟲號: 448168
- 注冊: 2007-11-01
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
lovibond(金幣+1): 鼓勵耐心交流 2011-06-18 10:35:05
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
lovibond(金幣+1): 鼓勵耐心交流 2011-06-18 10:35:05
|
在特定精度下終止程序?qū)嶋H上就是在一定的子空間內(nèi)取得了解的近似。我的意思是你得到看似正確的解可能是因?yàn)檫@個子空間還沒有包含任何系數(shù)矩陣的負(fù)曲率方向，這種情況下能得到滿足精度要求的一個近似解似乎也很自然。一旦偵測到一個負(fù)的方向，迭代序列很可能因此發(fā)散。具體的表現(xiàn)是什么，還是需要比較嚴(yán)謹(jǐn)?shù)耐茖?dǎo)。總之這種情況下用CG得到的結(jié)果是沒有說服力的，這個不是改變終止準(zhǔn)則就能修正的。 |
4樓2010-09-07 16:06:00
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料與化工求調(diào)劑 +5 | 為學(xué)666 2026-03-16 | 5/250 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 332求調(diào)劑 +3 | zjy101327 2026-03-11 | 6/300 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
|
[考研] 求調(diào)劑 資源與環(huán)境 285 +3 | 未名考生 2026-03-10 | 3/150 |
|