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