1 SUBROUTINE RESULT(NR,N,X,F,G,A,P,ITNCNT,IFLG,IPR) 2 IMPLICIT DOUBLE PRECISION (A-H,O-Z) 3C 4C PURPOSE 5C ------- 6C PRINT INFORMATION 7C 8C PARAMETERS 9C ---------- 10C NR --> ROW DIMENSION OF MATRIX 11C N --> DIMENSION OF PROBLEM 12C X(N) --> ITERATE X[K] 13C F --> FUNCTION VALUE AT X[K] 14C G(N) --> GRADIENT AT X[K] 15C A(N,N) --> HESSIAN AT X[K] 16C P(N) --> STEP TAKEN 17C ITNCNT --> ITERATION NUMBER K 18C IFLG --> FLAG CONTROLLING INFO TO PRINT 19C IPR --> DEVICE TO WHICH TO SEND OUTPUT 20C 21C***REVISION HISTORY (YYMMDD) 22C 000330 Modified array declarations. (JEC) 23C 24 DIMENSION X(N),G(N),P(N),A(NR,*) 25C PRINT ITERATION NUMBER 26 WRITE(IPR,903) ITNCNT 27 IF(IFLG.EQ.0) GO TO 120 28C 29C PRINT STEP 30 WRITE(IPR,907) 31 WRITE(IPR,905) (P(I),I=1,N) 32C 33C PRINT CURRENT ITERATE 34 120 CONTINUE 35 WRITE(IPR,904) 36 WRITE(IPR,905) (X(I),I=1,N) 37C 38C PRINT FUNCTION VALUE 39 WRITE(IPR,906) 40 WRITE(IPR,905) F 41C 42C PRINT GRADIENT 43 WRITE(IPR,908) 44 WRITE(IPR,905) (G(I),I=1,N) 45C 46C PRINT HESSIAN FROM ITERATION K 47 IF(IFLG.EQ.0) GO TO 140 48 WRITE(IPR,901) 49 DO 130 I=1,N 50 WRITE(IPR,900) I 51 WRITE(IPR,902) (A(I,J),J=1,I) 52 130 CONTINUE 53C 54 140 RETURN 55 900 FORMAT(15H RESULT ROW,I5) 56 901 FORMAT(29H RESULT HESSIAN AT X(K)) 57 902 FORMAT(14H RESULT ,5(2X,E20.13)) 58 903 FORMAT(/21H0RESULT ITERATE K=,I5) 59 904 FORMAT(18H RESULT X(K)) 60 905 FORMAT(22H RESULT ,5(2X,E20.13) ) 61 906 FORMAT(30H RESULT FUNCTION AT X(K)) 62 907 FORMAT(18H RESULT STEP) 63 908 FORMAT(30H RESULT GRADIENT AT X(K)) 64 END 65