| 查看: 2854 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
csyky2007木虫 (著名写手)
|
[交流]
【求助】用fortran求解大型线性方程组时出现的错误【已解决】 已有6人参与
|
||
|
我用fortran进行有限元分析，其中要涉及求解F=KU这个线性方程组。其中K为100万阶稀疏对称正定矩阵，所以我用IMSL函数库中的LSLXD进行求解。当程序进行到LSLXD函数时就出现图中所示的错误。我用这个程序实验了一个7000阶的系数矩阵时就没有问题，所以说程序应该没有错。 按照图中的提示，其中提到workspace不足，我试着将stack按本站中的帖子调到了最大0xffffffff，基本没有作用。 后我又指定itwksp=2000000000，可系统提示说：the required storage cannot be allocated. the workspace is based on n=1000000(矩阵的阶数) and nz(下三角矩阵中非零元素的个数) 图中的提示说是L2FXD出错，这个函数在IMSL帮助里有说明，就是利用矩阵的阶数和非零元素的个数精确指定worksapce。这个函数没有提供例子，我还真不会使用，其中有几个参数我实在搞不懂怎么定义。 最后值得一提是，当我用7000阶矩阵计算时，最后的命令提示我应该设置itwksp=13636849（具体数值），而这次给的是itwksp=********，好无奈的。 见这个论坛上牛人蛮多的，我真希望你们能帮我解决这个困扰我一个多月的问题！  [ Last edited by csyky2007 on 2010-7-5 at 22:32 ] |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
78506174
至尊木虫 (知名作家)
- 应助: 70 (初中生)
- 金币: 26353.9
- 散金: 50
- 红花: 10
- 帖子: 5712
- 在线: 830.9小时
- 虫号: 1331748
- 注册: 2011-06-26
- 专业: 金属结构材料
这个网址打不开了，我遇到了和你一样的问题 
12楼2015-10-05 22:50:12
csyky2007
木虫 (著名写手)
- 应助: 8 (幼儿园)
- 金币: 1476.9
- 散金: 3427
- 红花: 42
- 帖子: 1050
- 在线: 622.5小时
- 虫号: 906584
- 注册: 2009-11-18
- 性别: GG
- 专业: 固体力学
2楼2010-06-26 17:58:36
zxczxc0417
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1989.1
- 散金: 552
- 帖子: 925
- 在线: 544.3小时
- 虫号: 326712
- 注册: 2007-03-18
- 专业: 催化化学
3楼2010-06-27 03:40:49
csyky2007
木虫 (著名写手)
- 应助: 8 (幼儿园)
- 金币: 1476.9
- 散金: 3427
- 红花: 42
- 帖子: 1050
- 在线: 622.5小时
- 虫号: 906584
- 注册: 2009-11-18
- 性别: GG
- 专业: 固体力学
4楼2010-06-27 10:55:29
