应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】用fortran求解大型线性方程组时出现的错误【已解决】
121/2返回列表12下一页
查看: 2387  |  回复: 11
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
求放心！
1楼2010-06-24 22:45:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csyky2007

木虫 (著名写手)
呵呵,估计很少人涉足此类问题吧,那我想先请教下个具体问题吧:为解线性方程组分配空间的L2...此类的函数是怎么使用呢?希望有用过的热心人能给予指教,不胜感激!!!
赞一下 回复此楼
求放心！
2楼2010-06-26 17:58:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxczxc0417

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
fortran
很久以前的故事了
赞一下(1人) 回复此楼
3楼2010-06-27 03:40:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csyky2007

木虫 (著名写手)
引用回帖:
Originally posted by zxczxc0417 at 2010-06-27 03:40:49:
fortran
很久以前的故事了

那求解一个大型线性方程组，一般推荐用什么方法或语言呢？
赞一下 回复此楼
求放心！
4楼2010-06-27 10:55:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxczxc0417

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by csyky2007 at 2010-06-27 10:55:29:


那求解一个大型线性方程组，一般推荐用什么方法或语言呢？

matlab
回复此楼
5楼2010-06-27 11:30:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csyky2007

木虫 (著名写手)
引用回帖:
Originally posted by zxczxc0417 at 2010-06-27 11:30:47:

matlab

我最开始就是使用matlab的，正因为求解不了100万阶的我才转向fortran

[ Last edited by csyky2007 on 2010-6-27 at 13:17 ]
赞一下 回复此楼
求放心！
6楼2010-06-27 13:13:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzxin111

铜虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+1):谢谢参与应助！ 2010-07-02 22:55:19
求解稀疏对称方程组，一般采用ICCG法
赞一下(2人) 回复此楼
7楼2010-07-02 19:53:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hxz

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
64位环境么?
回复此楼
8楼2010-07-03 12:46:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csyky2007

木虫 (著名写手)
引用回帖:
Originally posted by zzxin111 at 2010-07-02 19:53:57:
求解稀疏对称方程组，一般采用ICCG法

这是迭代方法吧，谢谢你的建议。
我现在找到了一个HSL库把问题解决了，呵呵
它有专门解决线性方程组的fortran程序
http://www.hsl.rl.ac.uk/catalog.html
赞一下 回复此楼
求放心！
9楼2010-07-05 22:26:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csyky2007

木虫 (著名写手)
引用回帖:
Originally posted by hxz at 2010-07-03 12:46:34:
64位环境么?

是32位，我问题已经解决了，谢谢大家的参与，谢谢！
赞一下 回复此楼
求放心！
10楼2010-07-05 22:27:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 csyky2007 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭