LAPACK  3.6.1
LAPACK: Linear Algebra PACKage
subroutine stest ( integer  LEN,
double precision, dimension(len)  SCOMP,
double precision, dimension(len)  STRUE,
double precision, dimension(len)  SSIZE,
double precision  SFAC 
)

Definition at line 564 of file zblat1.f.

564 * ********************************* STEST **************************
565 *
566 * THIS SUBR COMPARES ARRAYS SCOMP() AND STRUE() OF LENGTH LEN TO
567 * SEE IF THE TERM BY TERM DIFFERENCES, MULTIPLIED BY SFAC, ARE
568 * NEGLIGIBLE.
569 *
570 * C. L. LAWSON, JPL, 1974 DEC 10
571 *
572 * .. Parameters ..
573  INTEGER nout
574  DOUBLE PRECISION zero
575  parameter (nout=6, zero=0.0d0)
576 * .. Scalar Arguments ..
577  DOUBLE PRECISION sfac
578  INTEGER len
579 * .. Array Arguments ..
580  DOUBLE PRECISION scomp(len), ssize(len), strue(len)
581 * .. Scalars in Common ..
582  INTEGER icase, incx, incy, mode, n
583  LOGICAL pass
584 * .. Local Scalars ..
585  DOUBLE PRECISION sd
586  INTEGER i
587 * .. External Functions ..
588  DOUBLE PRECISION sdiff
589  EXTERNAL sdiff
590 * .. Intrinsic Functions ..
591  INTRINSIC abs
592 * .. Common blocks ..
593  COMMON /combla/icase, n, incx, incy, mode, pass
594 * .. Executable Statements ..
595 *
596  DO 40 i = 1, len
597  sd = scomp(i) - strue(i)
598  IF (abs(sfac*sd) .LE. abs(ssize(i))*epsilon(zero))
599  + GO TO 40
600 *
601 * HERE SCOMP(I) IS NOT CLOSE TO STRUE(I).
602 *
603  IF ( .NOT. pass) GO TO 20
604 * PRINT FAIL MESSAGE AND HEADER.
605  pass = .false.
606  WRITE (nout,99999)
607  WRITE (nout,99998)
608  20 WRITE (nout,99997) icase, n, incx, incy, mode, i, scomp(i),
609  + strue(i), sd, ssize(i)
610  40 CONTINUE
611  RETURN
612 *
613 99999 FORMAT (' FAIL')
614 99998 FORMAT (/' CASE N INCX INCY MODE I ',
615  + ' COMP(I) TRUE(I) DIFFERENCE',
616  + ' SIZE(I)',/1x)
617 99997 FORMAT (1x,i4,i3,3i5,i3,2d36.8,2d12.4)
real function sdiff(SA, SB)
Definition: cblat1.f:645