网上找了个求解矩阵的程序

shrine

用起来有点问题
        error #6633: The type of the actual argument differs from the type of the dummy argument.   [1.E-15]        

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	program main

02	implicit none

03	real(8)::a3(18,9),A3inv(9,18)

04	integer::log

05

a3=1.

06	call BGINV(18,9,A3,A3inv, Log, 1.e-15)

07

end

08

09	!!-----------------------------------------

10	!!         general inverse

11	!!-----------------------------------------

12	SUBROUTINE BGINV(M,N,A,AA,L,EPS)

13	double precision :: A(M,N),U(M,M),V(N,N),AA(N,M),EPS

14	double precision,ALLOCATABLE :: S(:), E(:), WORK(:)

15	KA=MAX(M,N)+1

16	ALLOCATE(S(KA))

17	ALLOCATE(E(KA))

18	ALLOCATE(WORK(KA))

19	CALL BMUAV(A,M,N,U,V,L,EPS,KA,S,E,WORK)

20	IF (L.EQ.0) THEN

21

K=1

22	10     IF (A(K,K).NE.0.0) THEN

23

K=K+1

24	IF (K.LE.MIN(M,N)) GOTO 10

25

END IF

26

K=K-1

27	IF (K.NE.0) THEN

28	DO 40 I=1,N

29	DO 40 J=1,M

30	AA(I,J)=0.0

31	DO 30 II=1,K

32	30        AA(I,J)=AA(I,J)+V(II,I)*U(J,II)/A(II,II)

33	40      CONTINUE

34

END IF

35

END IF

36	DEALLOCATE(S)

37	DEALLOCATE(E)

38	DEALLOCATE(WORK)

39

RETURN

40

END

41

42	!!-----------------------------------------

43	!!         singular decomposite

44	!!-----------------------------------------

45	SUBROUTINE BMUAV(A,M,N,U,V,L,EPS,KA,S,E,WORK)

46

!略去

47

end

kyra

1.E-15 改为 1.D-15