LAPACK  3.9.1
LAPACK: Linear Algebra PACKage

◆ testdsdot()

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

Definition at line 975 of file dblat1.f.

976 * ********************************* STEST **************************
977 *
978 * THIS SUBR COMPARES ARRAYS SCOMP() AND STRUE() OF LENGTH LEN TO
979 * SEE IF THE TERM BY TERM DIFFERENCES, MULTIPLIED BY SFAC, ARE
980 * NEGLIGIBLE.
981 *
982 * C. L. LAWSON, JPL, 1974 DEC 10
983 *
984 * .. Parameters ..
985  INTEGER NOUT
986  REAL ZERO
987  parameter(nout=6, zero=0.0e0)
988 * .. Scalar Arguments ..
989  REAL SFAC, SCOMP, SSIZE, STRUE
990 * .. Scalars in Common ..
991  INTEGER ICASE, INCX, INCY, N
992  LOGICAL PASS
993 * .. Local Scalars ..
994  REAL SD
995 * .. Intrinsic Functions ..
996  INTRINSIC abs
997 * .. Common blocks ..
998  COMMON /combla/icase, n, incx, incy, pass
999 * .. Executable Statements ..
1000 *
1001  sd = scomp - strue
1002  IF (abs(sfac*sd) .LE. abs(ssize) * epsilon(zero))
1003  + GO TO 40
1004 *
1005 * HERE SCOMP(I) IS NOT CLOSE TO STRUE(I).
1006 *
1007  IF ( .NOT. pass) GO TO 20
1008 * PRINT FAIL MESSAGE AND HEADER.
1009  pass = .false.
1010  WRITE (nout,99999)
1011  WRITE (nout,99998)
1012  20 WRITE (nout,99997) icase, n, incx, incy, scomp,
1013  + strue, sd, ssize
1014  40 CONTINUE
1015  RETURN
1016 *
1017 99999 FORMAT (' FAIL')
1018 99998 FORMAT (/' CASE N INCX INCY ',
1019  + ' COMP(I) TRUE(I) DIFFERENCE',
1020  + ' SIZE(I)',/1x)
1021 99997 FORMAT (1x,i4,i3,1i5,i3,2e36.8,2e12.4)
Here is the caller graph for this function: