LAPACK 3.12.0
LAPACK: Linear Algebra PACKage
Loading...
Searching...
No Matches

◆ testdsdot()

subroutine testdsdot ( real  scomp,
real  strue,
real  ssize,
real  sfac 
)

Definition at line 1001 of file dblat1.f.

1002* ********************************* STEST **************************
1003*
1004* THIS SUBR COMPARES ARRAYS SCOMP() AND STRUE() OF LENGTH LEN TO
1005* SEE IF THE TERM BY TERM DIFFERENCES, MULTIPLIED BY SFAC, ARE
1006* NEGLIGIBLE.
1007*
1008* C. L. LAWSON, JPL, 1974 DEC 10
1009*
1010* .. Parameters ..
1011 INTEGER NOUT
1012 REAL ZERO
1013 parameter(nout=6, zero=0.0e0)
1014* .. Scalar Arguments ..
1015 REAL SFAC, SCOMP, SSIZE, STRUE
1016* .. Scalars in Common ..
1017 INTEGER ICASE, INCX, INCY, N
1018 LOGICAL PASS
1019* .. Local Scalars ..
1020 REAL SD
1021* .. Intrinsic Functions ..
1022 INTRINSIC abs
1023* .. Common blocks ..
1024 COMMON /combla/icase, n, incx, incy, pass
1025* .. Executable Statements ..
1026*
1027 sd = scomp - strue
1028 IF (abs(sfac*sd) .LE. abs(ssize) * epsilon(zero))
1029 + GO TO 40
1030*
1031* HERE SCOMP(I) IS NOT CLOSE TO STRUE(I).
1032*
1033 IF ( .NOT. pass) GO TO 20
1034* PRINT FAIL MESSAGE AND HEADER.
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* End of TESTDSDOT
1050*
Here is the caller graph for this function: