| 查看: 2293 | 回复: 4 | ||
zhang86木虫 (小有名气)
|
[求助]
Fortran中高效计算矩阵的秩
|
|
最近使用Fortran处理大型矩阵（大概3000*3000），希望能判断矩阵是否奇异。 原本想用行列式判断，即IMSL中的det函数，但无奈的发现计算结果不靠谱，因为矩阵太大，虽然已经有特征值很少（1.0e-17），但计算出来的行列式还是很大的，无法由此判断奇异与否。 后来改用svd函数计算矩阵的秩，这样还是很靠谱的。不过新的问题又来了，因为矩阵太大，计算的速度较慢。 请版上的各位大牛们帮忙出出主意，看看有什么办法可以更为高效的求解矩阵的秩。 ~~十分感谢~~ |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
fortran能求带字母矩阵的本征值吗？（矩阵元素带有字母，能得到形式解吗）
已经有5人回复
- fortran程序怎样调用lapack库进行矩阵对角化？
已经有7人回复
- 【求助】Fortran语言赋值问题？
已经有3人回复
- 【求助】fortran求解矩阵
已经有7人回复
- 【求助】fortran：如何求解 复数矩阵的秩【已完结】
已经有5人回复
- 【分享】fortran入门视频
已经有20人回复
- 【求助】fortran计算时出现NaN
已经有11人回复
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
2楼2012-09-16 09:48:19
zhang86
木虫 (小有名气)
- 应助: 37 (小学生)
- 金币: 3152.9
- 帖子: 218
- 在线: 117.3小时
- 虫号: 1354456
- 注册: 2011-07-24
- 性别: GG
- 专业: 固体力学
3楼2012-09-16 10:01:39
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
4楼2012-09-16 10:44:41
 |
5楼2012-09-16 11:16:01
