LAPACK  3.10.1
LAPACK: Linear Algebra PACKage

◆ testdsdot()

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

Definition at line 996 of file dblat1.f.

997 * ********************************* STEST **************************
998 *
999 * THIS SUBR COMPARES ARRAYS SCOMP() AND STRUE() OF LENGTH LEN TO
1000 * SEE IF THE TERM BY TERM DIFFERENCES, MULTIPLIED BY SFAC, ARE
1001 * NEGLIGIBLE.
1002 *
1003 * C. L. LAWSON, JPL, 1974 DEC 10
1004 *
1005 * .. Parameters ..
1006  INTEGER NOUT
1007  REAL ZERO
1008  parameter(nout=6, zero=0.0e0)
1009 * .. Scalar Arguments ..
1010  REAL SFAC, SCOMP, SSIZE, STRUE
1011 * .. Scalars in Common ..
1012  INTEGER ICASE, INCX, INCY, N
1013  LOGICAL PASS
1014 * .. Local Scalars ..
1015  REAL SD
1016 * .. Intrinsic Functions ..
1017  INTRINSIC abs
1018 * .. Common blocks ..
1019  COMMON /combla/icase, n, incx, incy, pass
1020 * .. Executable Statements ..
1021 *
1022  sd = scomp - strue
1023  IF (abs(sfac*sd) .LE. abs(ssize) * epsilon(zero))
1024  + GO TO 40
1025 *
1026 * HERE SCOMP(I) IS NOT CLOSE TO STRUE(I).
1027 *
1028  IF ( .NOT. pass) GO TO 20
1029 * PRINT FAIL MESSAGE AND HEADER.
1030  pass = .false.
1031  WRITE (nout,99999)
1032  WRITE (nout,99998)
1033  20 WRITE (nout,99997) icase, n, incx, incy, scomp,
1034  + strue, sd, ssize
1035  40 CONTINUE
1036  RETURN
1037 *
1038 99999 FORMAT (' FAIL')
1039 99998 FORMAT (/' CASE N INCX INCY ',
1040  + ' COMP(I) TRUE(I) DIFFERENCE',
1041  + ' SIZE(I)',/1x)
1042 99997 FORMAT (1x,i4,i3,1i5,i3,2e36.8,2e12.4)
1043 *
1044 * End of TESTDSDOT
1045 *
Here is the caller graph for this function: