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