fork(1) download

  PROGRAM HALLO3B
    ! IMPLICIT NONE
     PARAMETER(n=6)
     INTEGER::i,j
     DOUBLE PRECISION DA(n,n),A(n,n),X(n,n),Z(n,n),S(n),E,max
     DO i=1,n
      DO j=1,n
        READ(*,*) A(i,j)
      END DO
     END DO
     CALL HALLO3(A,X)
     CALL HALLO3(X,Z)
     DO i=1,n
      DO j=1,n
       DA(i,j)=Z(i,j)-A(i,j)
      END DO
     END DO 
     DO i=1,n
      S(i)=0
      DO j=1,n
       S(i)=S(i)+DABS(DA(i,j))
      END DO
     END DO
     max=S(1)
     DO i=2,n
     IF (S(i)>max) THEN
      max=S(i)
     END IF
     E=max
     END
 
 
 
 
 
 
   SUBROUTINE HALLO3(A,X)
   PARAMETER(n=6)
   INTEGER::i,j,k,IOPT,IREG,ND,MON(n,n)
    DOUBLE PRECISION A(n,n),AI(n,n+1),B(n),X(n,n),D
    ND=n
   DO i=1,n
    DO j=1,n
     IF (i==j) THEN
      MON(i,j)=1
     ELSE
      MON(i,j)=0
     END IF
      AI(i,j)=A(i,j)
    END DO
   END DO
   IOPT=1
   CALL DGAUSS(ND,n,A,B,IOPT,D,IREG)
   IOPT=2
   DO k=1,n
    DO i=1,n
     AI(i,n+1)=MON(i,k)
    END DO
    CALL DGAUSS(ND,n,A,B,IOPT,D,IREG)
    DO i=1,n
     X(i,k)=B(i)
    END DO
   END DO
  END
 
 
 
 
  SUBROUTINE DGAUSS(ND,N,A,B,IOPT,D,IREG)
   DOUBLE PRECISION A,B,D,S,EE
      DIMENSION A(ND,N),B(ND),L(ND)
      EE=1.D-12
      IREG=1
      NM=N-1
      IF(IOPT.EQ.2) GO TO 60
 
      D=1.D0
      DO 10 K=1,NM
      KP=K+1
      M=K
         DO 20 I=KP,N
         IF(DABS(A(I,K)).GT.DABS(A(M,K))) M=I
   20    CONTINUE
      L(K)=M
      D=D*A(M,K)
      IF(DABS(A(M,K)).LE.EE) GO TO 100   
      IF(M.EQ.K) GO TO 30
      D=-D
         DO 40 J=K,N
         S=A(M,J)
         A(M,J)=A(K,J)
   40    A(K,J)=S
   30    DO 50 I=KP,N
         S=A(I,K)/A(K,K)
         A(I,K)=S
            DO 50 J=KP,N
   50       A(I,J)=A(I,J)-S*A(K,J)
   10 CONTINUE
      D=D*A(N,N)
      IF(DABS(A(N,N)).LE.EE) GO TO 100
      IF(IOPT.EQ.1) RETURN
 
   60 DO 70 K=1,NM
      KP=K+1
      M=L(K)
      S=B(M)
      B(M)=B(K)
      B(K)=S
         DO 70 I=KP,N
   70    B(I)=B(I)-A(I,K)*S
 
      B(N)=B(N)/A(N,N)
      DO 80 I=1,NM
      K=N-I
      KP=K+1
      S=0.D0
         DO 90 J=KP,N
   90    S=S+A(K,J)*B(J)
   80 B(K)=(B(K)-S)/A(K,K)
      RETURN
  100 IREG=0
      RETURN
      END
 

prog.f95:30.5:

     END
     1
Error: END DO statement expected at (1)
prog.f95:37.3:

   SUBROUTINE HALLO3(A,X)
   1
Error: Unclassifiable statement at (1)
prog.f95:38.16:

   PARAMETER(n=6)
                1
Error: PARAMETER attribute of 'n' conflicts with PARAMETER attribute at (1)
prog.f95:39.13:

   INTEGER::i,j,k,IOPT,IREG,ND,MON(n,n)
             1
Error: Symbol 'i' at (1) already has basic type of INTEGER
prog.f95:40.22:

    DOUBLE PRECISION A(n,n),AI(n,n+1),B(n),X(n,n),D
                      1
Error: Symbol 'a' at (1) already has basic type of REAL
prog.f95:42.11:

   DO i=1,n
           1
prog.f95:25.13:

     DO i=2,n
             2
Error: Variable 'i' at (1) cannot be redefined inside loop beginning at (2)
prog.f95:45.16:

      MON(i,j)=1
                1
Error: Unexpected STATEMENT FUNCTION statement at (1)
prog.f95:47.16:

      MON(i,j)=0
                1
Error: Unexpected STATEMENT FUNCTION statement at (1)
prog.f95:49.20:

      AI(i,j)=A(i,j)
                    1
Error: Unexpected STATEMENT FUNCTION statement at (1)
prog.f95:56.12:

    DO i=1,n
            1
prog.f95:25.13:

     DO i=2,n
             2
Error: Variable 'i' at (1) cannot be redefined inside loop beginning at (2)
prog.f95:57.5:

     AI(i,n+1)=MON(i,k)
     1
Error: Unclassifiable statement at (1)
prog.f95:60.12:

    DO i=1,n
            1
prog.f95:25.13:

     DO i=2,n
             2
Error: Variable 'i' at (1) cannot be redefined inside loop beginning at (2)
prog.f95:64.2:

  END
  1
Error: END DO statement expected at (1)
prog.f95:69.2:

  SUBROUTINE DGAUSS(ND,N,A,B,IOPT,D,IREG)
  1
Error: Unclassifiable statement at (1)
prog.f95:70.21:

   DOUBLE PRECISION A,B,D,S,EE
                     1
Error: Symbol 'a' at (1) already has basic type of REAL
prog.f95:71.17:

      DIMENSION A(ND,N),B(ND),L(ND)
                 1
Error: Duplicate DIMENSION attribute specified at (1)
prog.f95:81.21:

         DO 20 I=KP,N
                     1
prog.f95:25.13:

     DO i=2,n
             2
Error: Variable 'i' at (1) cannot be redefined inside loop beginning at (2)
prog.f95:84.12:

      L(K)=M
            1
Error: Unexpected STATEMENT FUNCTION statement at (1)
prog.f95:92.17:

   40    A(K,J)=S
                 1
Warning: Obsolescent feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
prog.f95:93.21:

   30    DO 50 I=KP,N
                     1
prog.f95:25.13:

     DO i=2,n
             2
Error: Variable 'i' at (1) cannot be redefined inside loop beginning at (2)
prog.f95:96.24:

            DO 50 J=KP,N
                        1
Warning: Obsolescent feature: Shared DO termination label 50 at (1)
prog.f95:97.34:

   50       A(I,J)=A(I,J)-S*A(K,J)
                                  1
Warning: Obsolescent feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
prog.f95:101.26:

      IF(IOPT.EQ.1) RETURN
                          1
Error: GNU Extension: RETURN statement in main program at (1)
prog.f95:107.11:

      B(M)=B(K)
           1
Error: Statement function at (1) is recursive
prog.f95:108.12:

      B(K)=S
            1
Error: Unexpected STATEMENT FUNCTION statement at (1)
prog.f95:109.21:

         DO 70 I=KP,N
                     1
prog.f95:25.13:

     DO i=2,n
             2
Error: Variable 'i' at (1) cannot be redefined inside loop beginning at (2)
prog.f95:109.21:

         DO 70 I=KP,N
                     1
Warning: Obsolescent feature: Shared DO termination label 70 at (1)
prog.f95:110.14:

   70    B(I)=B(I)-A(I,K)*S
              1
Error: Statement function at (1) is recursive
prog.f95:112.11:

      B(N)=B(N)/A(N,N)
           1
Error: Statement function at (1) is recursive
Fatal Error: Error count reached limit of 25.

Standard output is empty