| 查看: 695 | 回复: 1 | ||
[求助]
奇异值分解法求广义逆代码求助 已有1人参与
|
|
以下是从《FORTRAN常用算法程序集》里面找到的求广义逆的奇异值分解法子程序代码，但是用例题测试之后发现计算有误，想请教下大神是不是子程序有误？ SUBROUTINE BGINV(M,N,A,AA,L,EPS,U,V,KA,S,E,WORK) DIMENSION A(M,N),U(M,M),V(N,N),AA(N,M) DIMENSION S(KA),E(KA),WORK(KA) DOUBLE PRECISION A,U,V,AA,S,E,WORK CALL BMUAV(A,M,N,U,V,L,EPS,KA,S,E,WORK) IF(L.EQ.0) THEN K=1 10 IF(A(K,K).NE.0.0) THEN K=K+1 IF(K.LE.MIN(M,N)) GO TO 10 END IF K=K-1 IF(K.NE.0) THEN DO 40 I=1,N DO 40 J=1,M AA(I,J)=0.0 DO 30 II=1,K 30 AA(I,J)=A(I,J)+V(II,I)*U(J,II)/A(II,II) 40 CONTINUE END IF END IF RETURN END |
回复此楼» 猜你喜欢
博士去二本高校当辅导员，如何调整心态？
已经有8人回复
-
27细胞生物学申博
已经有5人回复
-
面上有专家说实验设备不是我们单位的
已经有5人回复
-
江苏省自然基金 什么时候出结果
已经有6人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有11人回复
-
闲聊
已经有3人回复
-
如何查询2026年获批项目名称？
已经有3人回复
-
两块石头
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有15人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
2楼2015-08-19 10:02:30