1002
 1003
 1004
 1005
 1006
 1007
 1008
 1009
 1010
 1011      INTEGER          NOUT
 1012      REAL             ZERO
 1013      parameter(nout=6, zero=0.0e0)
 1014
 1015      REAL             SFAC, SCOMP, SSIZE, STRUE
 1016
 1017      INTEGER          ICASE, INCX, INCY, N
 1018      LOGICAL          PASS
 1019
 1020      REAL             SD
 1021
 1022      INTRINSIC        abs
 1023
 1024      COMMON           /combla/icase, n, incx, incy, pass
 1025
 1026
 1027         sd = scomp - strue
 1028         IF (abs(sfac*sd) .LE. abs(ssize) * epsilon(zero))
 1029     +       GO TO 40
 1030
 1031
 1032
 1033         IF ( .NOT. pass) GO TO 20
 1034
 1035         pass = .false.
 1036         WRITE (nout,99999)
 1037         WRITE (nout,99998)
 1038   20    WRITE (nout,99997) icase, n, incx, incy, scomp,
 1039     +     strue, sd, ssize
 1040   40 CONTINUE
 1041      RETURN
 1042
 104399999 FORMAT ('                                       FAIL')
 104499998 FORMAT (/' CASE  N INCX INCY                           ',
 1045     +       ' COMP(I)                             TRUE(I)  DIFFERENCE',
 1046     +       '     SIZE(I)',/1x)
 104799997 FORMAT (1x,i4,i3,1i5,i3,2e36.8,2e12.4)
 1048
 1049
 1050