LAPACK  3.9.0
LAPACK: Linear Algebra PACKage

◆ testdsdot()

subroutine testdsdot ( real  SCOMP,
real  STRUE,
real  SSIZE,
real  SFAC 
)

Definition at line 944 of file dblat1.f.

944 * ********************************* STEST **************************
945 *
946 * THIS SUBR COMPARES ARRAYS SCOMP() AND STRUE() OF LENGTH LEN TO
947 * SEE IF THE TERM BY TERM DIFFERENCES, MULTIPLIED BY SFAC, ARE
948 * NEGLIGIBLE.
949 *
950 * C. L. LAWSON, JPL, 1974 DEC 10
951 *
952 * .. Parameters ..
953  INTEGER NOUT
954  REAL ZERO
955  parameter(nout=6, zero=0.0e0)
956 * .. Scalar Arguments ..
957  REAL SFAC, SCOMP, SSIZE, STRUE
958 * .. Scalars in Common ..
959  INTEGER ICASE, INCX, INCY, N
960  LOGICAL PASS
961 * .. Local Scalars ..
962  REAL SD
963 * .. Intrinsic Functions ..
964  INTRINSIC abs
965 * .. Common blocks ..
966  COMMON /combla/icase, n, incx, incy, pass
967 * .. Executable Statements ..
968 *
969  sd = scomp - strue
970  IF (abs(sfac*sd) .LE. abs(ssize) * epsilon(zero))
971  + GO TO 40
972 *
973 * HERE SCOMP(I) IS NOT CLOSE TO STRUE(I).
974 *
975  IF ( .NOT. pass) GO TO 20
976 * PRINT FAIL MESSAGE AND HEADER.
977  pass = .false.
978  WRITE (nout,99999)
979  WRITE (nout,99998)
980  20 WRITE (nout,99997) icase, n, incx, incy, scomp,
981  + strue, sd, ssize
982  40 CONTINUE
983  RETURN
984 *
985 99999 FORMAT (' FAIL')
986 99998 FORMAT (/' CASE N INCX INCY ',
987  + ' COMP(I) TRUE(I) DIFFERENCE',
988  + ' SIZE(I)',/1x)
989 99997 FORMAT (1x,i4,i3,1i5,i3,2e36.8,2e12.4)
Here is the caller graph for this function: