[交流] 【求助】用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 ]

回复此楼