LAPACK  3.10.0
LAPACK: Linear Algebra PACKage

◆ stest()

subroutine stest ( integer  LEN,
real, dimension(len)  SCOMP,
real, dimension(len)  STRUE,
real, dimension(len)  SSIZE,
real  SFAC 
)

Definition at line 608 of file cblat1.f.

609 * ********************************* STEST **************************
610 *
611 * THIS SUBR COMPARES ARRAYS SCOMP() AND STRUE() OF LENGTH LEN TO
612 * SEE IF THE TERM BY TERM DIFFERENCES, MULTIPLIED BY SFAC, ARE
613 * NEGLIGIBLE.
614 *
615 * C. L. LAWSON, JPL, 1974 DEC 10
616 *
617 * .. Parameters ..
618  INTEGER NOUT
619  REAL ZERO
620  parameter(nout=6, zero=0.0e0)
621 * .. Scalar Arguments ..
622  REAL SFAC
623  INTEGER LEN
624 * .. Array Arguments ..
625  REAL SCOMP(LEN), SSIZE(LEN), STRUE(LEN)
626 * .. Scalars in Common ..
627  INTEGER ICASE, INCX, INCY, MODE, N
628  LOGICAL PASS
629 * .. Local Scalars ..
630  REAL SD
631  INTEGER I
632 * .. External Functions ..
633  REAL SDIFF
634  EXTERNAL sdiff
635 * .. Intrinsic Functions ..
636  INTRINSIC abs
637 * .. Common blocks ..
638  COMMON /combla/icase, n, incx, incy, mode, pass
639 * .. Executable Statements ..
640 *
641  DO 40 i = 1, len
642  sd = scomp(i) - strue(i)
643  IF (abs(sfac*sd) .LE. abs(ssize(i))*epsilon(zero))
644  + GO TO 40
645 *
646 * HERE SCOMP(I) IS NOT CLOSE TO STRUE(I).
647 *
648  IF ( .NOT. pass) GO TO 20
649 * PRINT FAIL MESSAGE AND HEADER.
650  pass = .false.
651  WRITE (nout,99999)
652  WRITE (nout,99998)
653  20 WRITE (nout,99997) icase, n, incx, incy, mode, i, scomp(i),
654  + strue(i), sd, ssize(i)
655  40 CONTINUE
656  RETURN
657 *
658 99999 FORMAT (' FAIL')
659 99998 FORMAT (/' CASE N INCX INCY MODE I ',
660  + ' COMP(I) TRUE(I) DIFFERENCE',
661  + ' SIZE(I)',/1x)
662 99997 FORMAT (1x,i4,i3,3i5,i3,2e36.8,2e12.4)
663 *
664 * End of STEST
665 *
real function sdiff(SA, SB)
Definition: cblat1.f:696
Here is the caller graph for this function: