 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (425)
- >導(dǎo)師招生 (42)
- >考研 (38)
- >基金申請(qǐng) (8)
- >考博 (7)
- >蟲(chóng)友互識(shí) (5)
- >碩博家園 (4)
- >招聘信息布告欄 (2)
- >博后之家 (2)
- >專(zhuān)業(yè)外語(yǔ) (1)
- >材料綜合 (1)
- >有機(jī)交流 (1)
- >找工作 (1)
- >論文投稿 (1)
- >休閑灌水 (1)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2516 | 回復(fù): 11 | |||
csyky2007木蟲(chóng) (著名寫(xiě)手)
|
[交流]
【求助】用fortran求解大型線性方程組時(shí)出現(xiàn)的錯(cuò)誤【已解決】 已有6人參與
|
|
我用fortran進(jìn)行有限元分析，其中要涉及求解F=KU這個(gè)線性方程組。其中K為100萬(wàn)階稀疏對(duì)稱(chēng)正定矩陣，所以我用IMSL函數(shù)庫(kù)中的LSLXD進(jìn)行求解。當(dāng)程序進(jìn)行到LSLXD函數(shù)時(shí)就出現(xiàn)圖中所示的錯(cuò)誤。我用這個(gè)程序?qū)嶒?yàn)了一個(gè)7000階的系數(shù)矩陣時(shí)就沒(méi)有問(wèn)題，所以說(shuō)程序應(yīng)該沒(méi)有錯(cuò)。 按照?qǐng)D中的提示，其中提到workspace不足，我試著將stack按本站中的帖子調(diào)到了最大0xffffffff，基本沒(méi)有作用。 后我又指定itwksp=2000000000，可系統(tǒng)提示說(shuō)：the required storage cannot be allocated. the workspace is based on n=1000000(矩陣的階數(shù)) and nz(下三角矩陣中非零元素的個(gè)數(shù)) 圖中的提示說(shuō)是L2FXD出錯(cuò)，這個(gè)函數(shù)在IMSL幫助里有說(shuō)明，就是利用矩陣的階數(shù)和非零元素的個(gè)數(shù)精確指定worksapce。這個(gè)函數(shù)沒(méi)有提供例子，我還真不會(huì)使用，其中有幾個(gè)參數(shù)我實(shí)在搞不懂怎么定義。 最后值得一提是，當(dāng)我用7000階矩陣計(jì)算時(shí)，最后的命令提示我應(yīng)該設(shè)置itwksp=13636849（具體數(shù)值），而這次給的是itwksp=********，好無(wú)奈的。 見(jiàn)這個(gè)論壇上牛人蠻多的，我真希望你們能幫我解決這個(gè)困擾我一個(gè)多月的問(wèn)題！  [ Last edited by csyky2007 on 2010-7-5 at 22:32 ] |
回復(fù)此樓» 猜你喜歡
【考研調(diào)劑】化學(xué)專(zhuān)業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿華南師大 070300（化學(xué)）304分求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。
已經(jīng)有3人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境專(zhuān)碩，總分308求調(diào)劑
已經(jīng)有3人回復(fù)
-
354求調(diào)劑
已經(jīng)有8人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿西南交大，求調(diào)劑
已經(jīng)有5人回復(fù)
-
299求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料專(zhuān)業(yè)求調(diào)劑
已經(jīng)有6人回復(fù)
-
311求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助：非線性方程組的求解（郁悶中）
已經(jīng)有26人回復(fù)
- 請(qǐng)教一個(gè)非線性schrodinger方程的解法
已經(jīng)有30人回復(fù)
- 求助：含多個(gè)未知系數(shù)的線性/矩陣方程組的解法
已經(jīng)有4人回復(fù)
- 用matlab解非線性方程組
已經(jīng)有4人回復(fù)
- fortran的array bounds exceeded錯(cuò)誤
已經(jīng)有5人回復(fù)
- 用fortran解關(guān)于x=tanx方程的程序
已經(jīng)有5人回復(fù)
- 關(guān)于應(yīng)用fortran編寫(xiě)高斯消去法程序來(lái)求解方程組的問(wèn)題
已經(jīng)有5人回復(fù)
- 用fortran程序遺傳算法解非線性方程組
已經(jīng)有7人回復(fù)
- 第5章_解線性方程組的直接方法.ppt
已經(jīng)有13人回復(fù)
- 求Matlab高手解決線性方程組的迭代求解問(wèn)題
已經(jīng)有17人回復(fù)
- matlab求解非線性方程組，錯(cuò)誤提示怎么解決
已經(jīng)有5人回復(fù)
- 求高人幫忙看一下我這個(gè)Fortran程序的錯(cuò)誤出在哪里？
已經(jīng)有12人回復(fù)
- 求matlab求解一個(gè)方程組的問(wèn)題
已經(jīng)有6人回復(fù)
- 【求助】fortran求解矩陣
已經(jīng)有7人回復(fù)
- 【求助】用fortran怎么描述下面這樣的情況，我是新手，請(qǐng)大家?guī)兔?/a>
已經(jīng)有8人回復(fù)
- 【求助】解方程組的錯(cuò)誤
已經(jīng)有7人回復(fù)
csyky2007
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1476.9
- 散金: 3427
- 紅花: 42
- 帖子: 1050
- 在線: 622.3小時(shí)
- 蟲(chóng)號(hào): 906584
- 注冊(cè): 2009-11-18
- 性別: GG
- 專(zhuān)業(yè): 固體力學(xué)
2樓2010-06-26 17:58:36
zxczxc0417
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1989.1
- 散金: 552
- 帖子: 925
- 在線: 544.3小時(shí)
- 蟲(chóng)號(hào): 326712
- 注冊(cè): 2007-03-18
- 專(zhuān)業(yè): 催化化學(xué)
3樓2010-06-27 03:40:49
csyky2007
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1476.9
- 散金: 3427
- 紅花: 42
- 帖子: 1050
- 在線: 622.3小時(shí)
- 蟲(chóng)號(hào): 906584
- 注冊(cè): 2009-11-18
- 性別: GG
- 專(zhuān)業(yè): 固體力學(xué)
4樓2010-06-27 10:55:29
zxczxc0417
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1989.1
- 散金: 552
- 帖子: 925
- 在線: 544.3小時(shí)
- 蟲(chóng)號(hào): 326712
- 注冊(cè): 2007-03-18
- 專(zhuān)業(yè): 催化化學(xué)
5樓2010-06-27 11:30:47
csyky2007
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1476.9
- 散金: 3427
- 紅花: 42
- 帖子: 1050
- 在線: 622.3小時(shí)
- 蟲(chóng)號(hào): 906584
- 注冊(cè): 2009-11-18
- 性別: GG
- 專(zhuān)業(yè): 固體力學(xué)
6樓2010-06-27 13:13:05
zzxin111
銅蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2248.6
- 帖子: 47
- 在線: 15.4小時(shí)
- 蟲(chóng)號(hào): 571956
- 注冊(cè): 2008-06-11
- 性別: GG
- 專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
7樓2010-07-02 19:53:57
hxz
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2105
- 帖子: 370
- 在線: 16.7小時(shí)
- 蟲(chóng)號(hào): 176464
- 注冊(cè): 2006-02-05
- 性別: GG
- 專(zhuān)業(yè): 永動(dòng)機(jī)
8樓2010-07-03 12:46:34
csyky2007
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1476.9
- 散金: 3427
- 紅花: 42
- 帖子: 1050
- 在線: 622.3小時(shí)
- 蟲(chóng)號(hào): 906584
- 注冊(cè): 2009-11-18
- 性別: GG
- 專(zhuān)業(yè): 固體力學(xué)
|
這是迭代方法吧，謝謝你的建議。 我現(xiàn)在找到了一個(gè)HSL庫(kù)把問(wèn)題解決了，呵呵 它有專(zhuān)門(mén)解決線性方程組的fortran程序 http://www.hsl.rl.ac.uk/catalog.html |
9樓2010-07-05 22:26:31
csyky2007
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1476.9
- 散金: 3427
- 紅花: 42
- 帖子: 1050
- 在線: 622.3小時(shí)
- 蟲(chóng)號(hào): 906584
- 注冊(cè): 2009-11-18
- 性別: GG
- 專(zhuān)業(yè): 固體力學(xué)
10樓2010-07-05 22:27:15
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 csyky2007 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +8 | 申子申申 2026-03-19 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專(zhuān)業(yè)調(diào)劑） +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶(hù)573181 2026-03-20 | 5/250 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4 | 守候夕陽(yáng)CF 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚(yú) 2026-03-18 | 6/300 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽(tīng)風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專(zhuān)業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|