| 查看: 1958 | 回复: 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 ] |
回复此楼» 猜你喜欢
专科院校申请国自然
已经有22人回复
-
论文被硕导卡
已经有5人回复
-
第一篇论文投稿进程记录
已经有11人回复
-
博后网上办公系统维护
已经有22人回复
-
暑假孩子一个人在家可以干什么
已经有7人回复
-
offer选择
已经有3人回复
-
夏至，要不要硕博联谊
已经有8人回复
-
基金评审时，评审专家还回去看申请人代表作之外的文章吗？
已经有12人回复
-
CSC博士联培对将来就业有用吗
已经有5人回复
-
对象没有，还非常想读博，难以抉择
已经有52人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：非线性方程组的求解（郁闷中）
已经有26人回复
- 请教一个非线性schrodinger方程的解法
已经有30人回复
- 请大神求解复杂非线性方程组（有报酬）
已经有22人回复
- 求助：含多个未知系数的线性/矩阵方程组的解法
已经有4人回复
- 用matlab解非线性方程组
已经有4人回复
- 线性不等式组的解的存在性问题
已经有8人回复
- fortran的array bounds exceeded错误
已经有5人回复
- 用fortran解关于x=tanx方程的程序
已经有5人回复
- 关于应用fortran编写高斯消去法程序来求解方程组的问题
已经有5人回复
- 求解非线性方程组
已经有7人回复
- 用fortran程序遗传算法解非线性方程组
已经有7人回复
- 如何使用灵活求解非线性方程组
已经有22人回复
- 第5章_解线性方程组的直接方法.ppt
已经有13人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- fortran里的求根问题
已经有8人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 如何使用mathematica求解非线性微分方程组
已经有11人回复
- 求高人帮忙看一下我这个Fortran程序的错误出在哪里？
已经有12人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 【求助】fortran求解矩阵
已经有7人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】解方程组的错误
已经有7人回复
csyky2007
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1785.4
- 散金: 3163
- 红花: 41
- 帖子: 971
- 在线: 587.8小时
- 虫号: 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 (幼儿园)
- 金币: 1785.4
- 散金: 3163
- 红花: 41
- 帖子: 971
- 在线: 587.8小时
- 虫号: 906584
- 注册: 2009-11-18
- 性别: GG
- 专业: 固体力学
4楼2010-06-27 10:55:29
zxczxc0417
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1989.1
- 散金: 552
- 帖子: 925
- 在线: 544.3小时
- 虫号: 326712
- 注册: 2007-03-18
- 专业: 催化化学
5楼2010-06-27 11:30:47
csyky2007
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1785.4
- 散金: 3163
- 红花: 41
- 帖子: 971
- 在线: 587.8小时
- 虫号: 906584
- 注册: 2009-11-18
- 性别: GG
- 专业: 固体力学
6楼2010-06-27 13:13:05
zzxin111
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2248.6
- 帖子: 47
- 在线: 15.4小时
- 虫号: 571956
- 注册: 2008-06-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2010-07-02 19:53:57
8楼2010-07-03 12:46:34
csyky2007
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1785.4
- 散金: 3163
- 红花: 41
- 帖子: 971
- 在线: 587.8小时
- 虫号: 906584
- 注册: 2009-11-18
- 性别: GG
- 专业: 固体力学
|
这是迭代方法吧，谢谢你的建议。 我现在找到了一个HSL库把问题解决了，呵呵 它有专门解决线性方程组的fortran程序 http://www.hsl.rl.ac.uk/catalog.html |
9楼2010-07-05 22:26:31
csyky2007
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1785.4
- 散金: 3163
- 红花: 41
- 帖子: 971
- 在线: 587.8小时
- 虫号: 906584
- 注册: 2009-11-18
- 性别: GG
- 专业: 固体力学
10楼2010-07-05 22:27:15
5