24小时热门版块排行榜    

查看: 2288  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhang86

木虫 (小有名气)

[求助] Fortran中高效计算矩阵的秩

最近使用Fortran处理大型矩阵(大概3000*3000),希望能判断矩阵是否奇异。
原本想用行列式判断,即IMSL中的det函数,但无奈的发现计算结果不靠谱,因为矩阵太大,虽然已经有特征值很少(1.0e-17),但计算出来的行列式还是很大的,无法由此判断奇异与否。
后来改用svd函数计算矩阵的秩,这样还是很靠谱的。不过新的问题又来了,因为矩阵太大,计算的速度较慢。
请版上的各位大牛们帮忙出出主意,看看有什么办法可以更为高效的求解矩阵的秩。
~~十分感谢~~
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhang86

木虫 (小有名气)

引用回帖:
2楼: Originally posted by chembetsey at 2012-09-16 09:48:19
IMSL中有rank函数,应该比你自己写的函数算的快,你可以试试。要想算得更快,可以用并行的rank。

谢谢。
但如何给定IMSL中的rank函数的误差限,也就是当特征值小于多少时可以认为是0。
还有,请问如何并行rank以提高效率?
我需要在依次迭代的过程中评价矩阵是否奇异,这样是不是没法用并行计算?
3楼2012-09-16 10:01:39
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

chembetsey

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
zhang86: 金币+2, 有帮助 2012-09-16 10:02:37
xzhdty: 谢谢 2012-09-16 21:14:26
IMSL中有rank函数,应该比你自己写的函数算的快,你可以试试。要想算得更快,可以用并行的rank。
2楼2012-09-16 09:48:19
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

chembetsey

木虫 (小有名气)

你自己查IMSL的手册吧,具体的我也不是很清楚。
并行只是rank函数并行,和你自己的程序没关系的。
4楼2012-09-16 10:44:41
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-24 3/150 2026-03-24 19:17 by pswait
[考研] 材料学硕,求调剂 6+3 糖葫芦888ll 2026-03-22 7/350 2026-03-24 17:11 by hello七七
[考研] 086003食品工程求调剂 +4 淼淼111 2026-03-24 4/200 2026-03-24 11:53 by 544594351
[考研] 【双一流院校新能源、环境材料,材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 7/350 2026-03-24 11:23 by 种大树
[考研] 材料与化工328分调剂 +4 。,。,。,。i 2026-03-23 4/200 2026-03-24 11:03 by 544594351
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 316求调剂 +7 梁茜雯 2026-03-19 7/350 2026-03-23 16:21 by lingjue
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
信息提示
请填处理意见