24小时热门版块排行榜

>论坛更新日志 (3034)
>虫友互识 (606)
>休闲灌水 (105)
>基金申请 (91)
>博后之家 (85)
>导师招生 (82)
>论文投稿 (78)
>文献求助 (68)
>教师之家 (44)
>硕博家园 (39)
>考博 (23)
>药学 (19)
>招聘信息布告栏 (13)
>医学 (11)
>公派出国 (10)
>考研 (9)

返回列表

kimileegdut

捐助贵宾 (小有名气)

应助: 1 (幼儿园)
帖子: 56
在线: 24小时
虫号: 3526305
注册: 2014-11-08

[求助] 奇异值分解法求广义逆代码求助已有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

回复此楼