24小时热门版块排行榜    

查看: 2393  |  回复: 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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

求放心!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxczxc0417

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流

csyky2007

木虫 (著名写手)

csyky2007

木虫 (著名写手)

zxczxc0417

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
普通表情 高级回复(可上传附件)
信息提示
请填处理意见