 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2514 | 回復(fù): 11 | |||
csyky2007木蟲 (著名寫手)
|
[交流]
【求助】用fortran求解大型線性方程組時出現(xiàn)的錯誤【已解決】 已有6人參與
|
|
我用fortran進(jìn)行有限元分析，其中要涉及求解F=KU這個線性方程組。其中K為100萬階稀疏對稱正定矩陣，所以我用IMSL函數(shù)庫中的LSLXD進(jìn)行求解。當(dāng)程序進(jìn)行到LSLXD函數(shù)時就出現(xiàn)圖中所示的錯誤。我用這個程序?qū)嶒灹艘粋�7000階的系數(shù)矩陣時就沒有問題，所以說程序應(yīng)該沒有錯。 按照圖中的提示，其中提到workspace不足，我試著將stack按本站中的帖子調(diào)到了最大0xffffffff，基本沒有作用。 后我又指定itwksp=2000000000，可系統(tǒng)提示說：the required storage cannot be allocated. the workspace is based on n=1000000(矩陣的階數(shù)) and nz(下三角矩陣中非零元素的個數(shù)) 圖中的提示說是L2FXD出錯，這個函數(shù)在IMSL幫助里有說明，就是利用矩陣的階數(shù)和非零元素的個數(shù)精確指定worksapce。這個函數(shù)沒有提供例子，我還真不會使用，其中有幾個參數(shù)我實在搞不懂怎么定義。 最后值得一提是，當(dāng)我用7000階矩陣計算時，最后的命令提示我應(yīng)該設(shè)置itwksp=13636849（具體數(shù)值），而這次給的是itwksp=********，好無奈的。 見這個論壇上牛人蠻多的，我真希望你們能幫我解決這個困擾我一個多月的問題！  [ Last edited by csyky2007 on 2010-7-5 at 22:32 ] |
回復(fù)此樓» 猜你喜歡
0856調(diào)劑，是學(xué)校就去
已經(jīng)有9人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有9人回復(fù)
-
265求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
已經(jīng)有7人回復(fù)
-
材料學(xué)碩297已過四六級求調(diào)劑推薦
已經(jīng)有6人回復(fù)
-
295復(fù)試調(diào)劑
已經(jīng)有7人回復(fù)
-
本人考085602 化學(xué)工程 專碩
已經(jīng)有20人回復(fù)
-
一志愿中南化學(xué)337求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有23人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助：非線性方程組的求解（郁悶中）
已經(jīng)有26人回復(fù)
- 請教一個非線性schrodinger方程的解法
已經(jīng)有30人回復(fù)
- 求助：含多個未知系數(shù)的線性/矩陣方程組的解法
已經(jīng)有4人回復(fù)
- 用matlab解非線性方程組
已經(jīng)有4人回復(fù)
- fortran的array bounds exceeded錯誤
已經(jīng)有5人回復(fù)
- 用fortran解關(guān)于x=tanx方程的程序
已經(jīng)有5人回復(fù)
- 關(guān)于應(yīng)用fortran編寫高斯消去法程序來求解方程組的問題
已經(jīng)有5人回復(fù)
- 用fortran程序遺傳算法解非線性方程組
已經(jīng)有7人回復(fù)
- 第5章_解線性方程組的直接方法.ppt
已經(jīng)有13人回復(fù)
- 求Matlab高手解決線性方程組的迭代求解問題
已經(jīng)有17人回復(fù)
- matlab求解非線性方程組，錯誤提示怎么解決
已經(jīng)有5人回復(fù)
- 求高人幫忙看一下我這個Fortran程序的錯誤出在哪里？
已經(jīng)有12人回復(fù)
- 求matlab求解一個方程組的問題
已經(jīng)有6人回復(fù)
- 【求助】fortran求解矩陣
已經(jīng)有7人回復(fù)
- 【求助】用fortran怎么描述下面這樣的情況，我是新手，請大家?guī)兔?/a>
已經(jīng)有8人回復(fù)
- 【求助】解方程組的錯誤
已經(jīng)有7人回復(fù)
csyky2007
木蟲 (著名寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1476.9
- 散金: 3427
- 紅花: 42
- 帖子: 1050
- 在線: 622.3小時
- 蟲號: 906584
- 注冊: 2009-11-18
- 性別: GG
- 專業(yè): 固體力學(xué)
2樓2010-06-26 17:58:36
zxczxc0417
木蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1989.1
- 散金: 552
- 帖子: 925
- 在線: 544.3小時
- 蟲號: 326712
- 注冊: 2007-03-18
- 專業(yè): 催化化學(xué)
3樓2010-06-27 03:40:49
csyky2007
木蟲 (著名寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1476.9
- 散金: 3427
- 紅花: 42
- 帖子: 1050
- 在線: 622.3小時
- 蟲號: 906584
- 注冊: 2009-11-18
- 性別: GG
- 專業(yè): 固體力學(xué)
4樓2010-06-27 10:55:29
zxczxc0417
木蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1989.1
- 散金: 552
- 帖子: 925
- 在線: 544.3小時
- 蟲號: 326712
- 注冊: 2007-03-18
- 專業(yè): 催化化學(xué)
5樓2010-06-27 11:30:47
csyky2007
木蟲 (著名寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1476.9
- 散金: 3427
- 紅花: 42
- 帖子: 1050
- 在線: 622.3小時
- 蟲號: 906584
- 注冊: 2009-11-18
- 性別: GG
- 專業(yè): 固體力學(xué)
6樓2010-06-27 13:13:05
zzxin111
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2248.6
- 帖子: 47
- 在線: 15.4小時
- 蟲號: 571956
- 注冊: 2008-06-11
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
7樓2010-07-02 19:53:57
hxz
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2105
- 帖子: 370
- 在線: 16.7小時
- 蟲號: 176464
- 注冊: 2006-02-05
- 性別: GG
- 專業(yè): 永動機
8樓2010-07-03 12:46:34
csyky2007
木蟲 (著名寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1476.9
- 散金: 3427
- 紅花: 42
- 帖子: 1050
- 在線: 622.3小時
- 蟲號: 906584
- 注冊: 2009-11-18
- 性別: GG
- 專業(yè): 固體力學(xué)
|
這是迭代方法吧，謝謝你的建議。 我現(xiàn)在找到了一個HSL庫把問題解決了，呵呵 它有專門解決線性方程組的fortran程序 http://www.hsl.rl.ac.uk/catalog.html |
9樓2010-07-05 22:26:31
csyky2007
木蟲 (著名寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1476.9
- 散金: 3427
- 紅花: 42
- 帖子: 1050
- 在線: 622.3小時
- 蟲號: 906584
- 注冊: 2009-11-18
- 性別: GG
- 專業(yè): 固體力學(xué)
10樓2010-07-05 22:27:15
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 csyky2007 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +8 | S.H1 2026-03-18 | 8/400 |
|
|
|
[考研] 290求調(diào)劑 +4 | ^O^乜 2026-03-19 | 4/200 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 8/400 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +6 | 一顆meteor 2026-03-13 | 6/300 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|