| 查看: 1952 | 回复: 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 ] |
回复此楼» 猜你喜欢
我标书代码变了
已经有65人回复
-
江南大学到瑞士招聘，称取消非升即走，改预聘+长聘
已经有7人回复
-
投稿求助
已经有7人回复
-
2024国社科通讯评审
已经有28人回复
-
药学硕士找不到工作，打算去做科研助理了
已经有17人回复
-
豫北虫友互识
已经有3人回复
-
最新消息：2024国自然下载文件名变了
已经有29人回复
-
试剂价格
已经有3人回复
-
34岁读博士晚吗
已经有32人回复
-
博后基金刷到的BUG,图片来的更直观
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：非线性方程组的求解（郁闷中）
已经有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