1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PROGRAM DBLAT1 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* Test program for the DOUBLE PRECISION Level 1 BLAS. 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* Based upon the original BLAS test routine together with: 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* F06EAF Example Program Text 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Parameters .. 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER NOUT 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PARAMETER (NOUT=6) 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalars in Common .. 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ICASE, INCX, INCY, MODE, N 10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LOGICAL PASS 11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Scalars .. 12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SFAC 13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER IC 14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. External Subroutines .. 15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath EXTERNAL CHECK0, CHECK1, CHECK2, CHECK3, HEADER 16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Common blocks .. 17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS 18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Data statements .. 19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA SFAC/9.765625D-4/ 20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Executable Statements .. 21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,99999) 22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 20 IC = 1, 10 23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ICASE = IC 24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL HEADER 25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Initialize PASS, INCX, INCY, and MODE for a new case. .. 27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. the value 9999 for INCX, INCY or MODE will appear in the .. 28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. detailed output, if any, for cases that do not involve .. 29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. these parameters .. 30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PASS = .TRUE. 32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INCX = 9999 33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INCY = 9999 34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MODE = 9999 35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF (ICASE.EQ.3) THEN 36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL CHECK0(SFAC) 37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE IF (ICASE.EQ.7 .OR. ICASE.EQ.8 .OR. ICASE.EQ.9 .OR. 38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + ICASE.EQ.10) THEN 39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL CHECK1(SFAC) 40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE IF (ICASE.EQ.1 .OR. ICASE.EQ.2 .OR. ICASE.EQ.5 .OR. 41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + ICASE.EQ.6) THEN 42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL CHECK2(SFAC) 43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE IF (ICASE.EQ.4) THEN 44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL CHECK3(SFAC) 45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END IF 46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* -- Print 47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF (PASS) WRITE (NOUT,99998) 48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 20 CONTINUE 49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STOP 50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath99999 FORMAT (' Real BLAS Test Program Results',/1X) 52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath99998 FORMAT (' ----- PASS -----') 53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END 54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SUBROUTINE HEADER 55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Parameters .. 56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER NOUT 57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PARAMETER (NOUT=6) 58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalars in Common .. 59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ICASE, INCX, INCY, MODE, N 60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LOGICAL PASS 61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Arrays .. 62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CHARACTER*6 L(10) 63c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Common blocks .. 64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS 65c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Data statements .. 66c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA L(1)/' DDOT '/ 67c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA L(2)/'DAXPY '/ 68c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA L(3)/'DROTG '/ 69c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA L(4)/' DROT '/ 70c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA L(5)/'DCOPY '/ 71c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA L(6)/'DSWAP '/ 72c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA L(7)/'DNRM2 '/ 73c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA L(8)/'DASUM '/ 74c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA L(9)/'DSCAL '/ 75c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA L(10)/'IDAMAX'/ 76c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Executable Statements .. 77c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,99999) ICASE, L(ICASE) 78c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath RETURN 79c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 80c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath99999 FORMAT (/' Test of subprogram number',I3,12X,A6) 81c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END 82c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SUBROUTINE CHECK0(SFAC) 83c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Parameters .. 84c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER NOUT 85c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PARAMETER (NOUT=6) 86c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalar Arguments .. 87c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SFAC 88c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalars in Common .. 89c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ICASE, INCX, INCY, MODE, N 90c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LOGICAL PASS 91c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Scalars .. 92c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION D12, SA, SB, SC, SS 93c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER K 94c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Arrays .. 95c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION DA1(8), DATRUE(8), DB1(8), DBTRUE(8), DC1(8), 96c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + DS1(8) 97c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. External Subroutines .. 98c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath EXTERNAL DROTG, STEST1 99c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Common blocks .. 100c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS 101c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Data statements .. 102c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DA1/0.3D0, 0.4D0, -0.3D0, -0.4D0, -0.3D0, 0.0D0, 103c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 1.0D0/ 104c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DB1/0.4D0, 0.3D0, 0.4D0, 0.3D0, -0.4D0, 0.0D0, 105c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 1.0D0, 0.0D0/ 106c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DC1/0.6D0, 0.8D0, -0.6D0, 0.8D0, 0.6D0, 1.0D0, 107c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 1.0D0/ 108c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DS1/0.8D0, 0.6D0, 0.8D0, -0.6D0, 0.8D0, 0.0D0, 109c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 1.0D0, 0.0D0/ 110c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DATRUE/0.5D0, 0.5D0, 0.5D0, -0.5D0, -0.5D0, 111c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 1.0D0, 1.0D0/ 112c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DBTRUE/0.0D0, 0.6D0, 0.0D0, -0.6D0, 0.0D0, 113c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 1.0D0, 0.0D0/ 114c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA D12/4096.0D0/ 115c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Executable Statements .. 116c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 117c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* Compute true values which cannot be prestored 118c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* in decimal notation 119c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 120c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DBTRUE(1) = 1.0D0/0.6D0 121c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DBTRUE(3) = -1.0D0/0.6D0 122c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DBTRUE(5) = 1.0D0/0.6D0 123c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 124c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 20 K = 1, 8 125c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Set N=K for identification in output if any .. 126c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath N = K 127c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF (ICASE.EQ.3) THEN 128c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. DROTG .. 129c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF (K.GT.8) GO TO 40 130c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SA = DA1(K) 131c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SB = DB1(K) 132c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL DROTG(SA,SB,SC,SS) 133c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST1(SA,DATRUE(K),DATRUE(K),SFAC) 134c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST1(SB,DBTRUE(K),DBTRUE(K),SFAC) 135c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST1(SC,DC1(K),DC1(K),SFAC) 136c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST1(SS,DS1(K),DS1(K),SFAC) 137c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE 138c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,*) ' Shouldn''t be here in CHECK0' 139c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STOP 140c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END IF 141c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 20 CONTINUE 142c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 40 RETURN 143c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END 144c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SUBROUTINE CHECK1(SFAC) 145c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Parameters .. 146c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER NOUT 147c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PARAMETER (NOUT=6) 148c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalar Arguments .. 149c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SFAC 150c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalars in Common .. 151c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ICASE, INCX, INCY, MODE, N 152c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LOGICAL PASS 153c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Scalars .. 154c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER I, LEN, NP1 155c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Arrays .. 156c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION DTRUE1(5), DTRUE3(5), DTRUE5(8,5,2), DV(8,5,2), 157c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + SA(10), STEMP(1), STRUE(8), SX(8) 158c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ITRUE2(5) 159c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. External Functions .. 160c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION DASUM, DNRM2 161c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER IDAMAX 162c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath EXTERNAL DASUM, DNRM2, IDAMAX 163c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. External Subroutines .. 164c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath EXTERNAL ITEST1, DSCAL, STEST, STEST1 165c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Intrinsic Functions .. 166c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTRINSIC MAX 167c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Common blocks .. 168c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS 169c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Data statements .. 170c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA SA/0.3D0, -1.0D0, 0.0D0, 1.0D0, 0.3D0, 0.3D0, 171c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.3D0, 0.3D0, 0.3D0, 0.3D0/ 172c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DV/0.1D0, 2.0D0, 2.0D0, 2.0D0, 2.0D0, 2.0D0, 173c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 2.0D0, 2.0D0, 0.3D0, 3.0D0, 3.0D0, 3.0D0, 3.0D0, 174c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 3.0D0, 3.0D0, 3.0D0, 0.3D0, -0.4D0, 4.0D0, 175c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 4.0D0, 4.0D0, 4.0D0, 4.0D0, 4.0D0, 0.2D0, 176c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.6D0, 0.3D0, 5.0D0, 5.0D0, 5.0D0, 5.0D0, 177c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 5.0D0, 0.1D0, -0.3D0, 0.5D0, -0.1D0, 6.0D0, 178c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 6.0D0, 6.0D0, 6.0D0, 0.1D0, 8.0D0, 8.0D0, 8.0D0, 179c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 8.0D0, 8.0D0, 8.0D0, 8.0D0, 0.3D0, 9.0D0, 9.0D0, 180c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 9.0D0, 9.0D0, 9.0D0, 9.0D0, 9.0D0, 0.3D0, 2.0D0, 181c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.4D0, 2.0D0, 2.0D0, 2.0D0, 2.0D0, 2.0D0, 182c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.2D0, 3.0D0, -0.6D0, 5.0D0, 0.3D0, 2.0D0, 183c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 2.0D0, 2.0D0, 0.1D0, 4.0D0, -0.3D0, 6.0D0, 184c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.5D0, 7.0D0, -0.1D0, 3.0D0/ 185c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DTRUE1/0.0D0, 0.3D0, 0.5D0, 0.7D0, 0.6D0/ 186c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DTRUE3/0.0D0, 0.3D0, 0.7D0, 1.1D0, 1.0D0/ 187c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DTRUE5/0.10D0, 2.0D0, 2.0D0, 2.0D0, 2.0D0, 188c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 2.0D0, 2.0D0, 2.0D0, -0.3D0, 3.0D0, 3.0D0, 189c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 3.0D0, 3.0D0, 3.0D0, 3.0D0, 3.0D0, 0.0D0, 0.0D0, 190c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 4.0D0, 4.0D0, 4.0D0, 4.0D0, 4.0D0, 4.0D0, 191c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.20D0, -0.60D0, 0.30D0, 5.0D0, 5.0D0, 5.0D0, 192c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 5.0D0, 5.0D0, 0.03D0, -0.09D0, 0.15D0, -0.03D0, 193c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 6.0D0, 6.0D0, 6.0D0, 6.0D0, 0.10D0, 8.0D0, 194c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 8.0D0, 8.0D0, 8.0D0, 8.0D0, 8.0D0, 8.0D0, 195c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.09D0, 9.0D0, 9.0D0, 9.0D0, 9.0D0, 9.0D0, 196c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 9.0D0, 9.0D0, 0.09D0, 2.0D0, -0.12D0, 2.0D0, 197c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 2.0D0, 2.0D0, 2.0D0, 2.0D0, 0.06D0, 3.0D0, 198c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.18D0, 5.0D0, 0.09D0, 2.0D0, 2.0D0, 2.0D0, 199c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.03D0, 4.0D0, -0.09D0, 6.0D0, -0.15D0, 7.0D0, 200c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.03D0, 3.0D0/ 201c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA ITRUE2/0, 1, 2, 2, 3/ 202c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Executable Statements .. 203c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 80 INCX = 1, 2 204c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 60 NP1 = 1, 5 205c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath N = NP1 - 1 206c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LEN = 2*MAX(N,1) 207c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Set vector arguments .. 208c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 20 I = 1, LEN 209c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SX(I) = DV(I,NP1,INCX) 210c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 20 CONTINUE 211c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 212c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF (ICASE.EQ.7) THEN 213c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. DNRM2 .. 214c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STEMP(1) = DTRUE1(NP1) 215c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST1(DNRM2(N,SX,INCX),STEMP,STEMP,SFAC) 216c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE IF (ICASE.EQ.8) THEN 217c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. DASUM .. 218c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STEMP(1) = DTRUE3(NP1) 219c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST1(DASUM(N,SX,INCX),STEMP,STEMP,SFAC) 220c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE IF (ICASE.EQ.9) THEN 221c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. DSCAL .. 222c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL DSCAL(N,SA((INCX-1)*5+NP1),SX,INCX) 223c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 40 I = 1, LEN 224c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STRUE(I) = DTRUE5(I,NP1,INCX) 225c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 40 CONTINUE 226c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST(LEN,SX,STRUE,STRUE,SFAC) 227c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE IF (ICASE.EQ.10) THEN 228c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. IDAMAX .. 229c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL ITEST1(IDAMAX(N,SX,INCX),ITRUE2(NP1)) 230c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE 231c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,*) ' Shouldn''t be here in CHECK1' 232c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STOP 233c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END IF 234c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 60 CONTINUE 235c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 80 CONTINUE 236c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath RETURN 237c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END 238c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SUBROUTINE CHECK2(SFAC) 239c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Parameters .. 240c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER NOUT 241c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PARAMETER (NOUT=6) 242c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalar Arguments .. 243c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SFAC 244c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalars in Common .. 245c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ICASE, INCX, INCY, MODE, N 246c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LOGICAL PASS 247c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Scalars .. 248c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SA, SC, SS 249c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER I, J, KI, KN, KSIZE, LENX, LENY, MX, MY 250c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Arrays .. 251c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION DT10X(7,4,4), DT10Y(7,4,4), DT7(4,4), 252c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + DT8(7,4,4), DT9X(7,4,4), DT9Y(7,4,4), DX1(7), 253c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + DY1(7), SSIZE1(4), SSIZE2(14,2), STX(7), STY(7), 254c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + SX(7), SY(7) 255c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER INCXS(4), INCYS(4), LENS(4,2), NS(4) 256c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. External Functions .. 257c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION DDOT 258c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath EXTERNAL DDOT 259c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. External Subroutines .. 260c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath EXTERNAL DAXPY, DCOPY, DSWAP, STEST, STEST1 261c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Intrinsic Functions .. 262c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTRINSIC ABS, MIN 263c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Common blocks .. 264c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS 265c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Data statements .. 266c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA SA/0.3D0/ 267c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA INCXS/1, 2, -2, -1/ 268c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA INCYS/1, -2, 1, -2/ 269c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA LENS/1, 1, 2, 4, 1, 1, 3, 7/ 270c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA NS/0, 1, 2, 4/ 271c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DX1/0.6D0, 0.1D0, -0.5D0, 0.8D0, 0.9D0, -0.3D0, 272c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.4D0/ 273c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DY1/0.5D0, -0.9D0, 0.3D0, 0.7D0, -0.6D0, 0.2D0, 274c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.8D0/ 275c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA SC, SS/0.8D0, 0.6D0/ 276c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DT7/0.0D0, 0.30D0, 0.21D0, 0.62D0, 0.0D0, 277c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.30D0, -0.07D0, 0.85D0, 0.0D0, 0.30D0, -0.79D0, 278c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.74D0, 0.0D0, 0.30D0, 0.33D0, 1.27D0/ 279c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DT8/0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 280c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.68D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 281c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.68D0, -0.87D0, 0.0D0, 0.0D0, 282c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.68D0, -0.87D0, 0.15D0, 283c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.94D0, 0.0D0, 0.0D0, 0.0D0, 0.5D0, 0.0D0, 284c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.68D0, 285c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 286c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.35D0, -0.9D0, 0.48D0, 0.0D0, 0.0D0, 0.0D0, 287c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.38D0, -0.9D0, 0.57D0, 0.7D0, -0.75D0, 288c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.2D0, 0.98D0, 0.5D0, 0.0D0, 0.0D0, 0.0D0, 289c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.68D0, 0.0D0, 0.0D0, 290c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.35D0, -0.72D0, 291c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.38D0, 292c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.63D0, 0.15D0, 0.88D0, 0.0D0, 0.0D0, 0.0D0, 293c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 294c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.68D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 295c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.68D0, -0.9D0, 0.33D0, 0.0D0, 0.0D0, 296c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.68D0, -0.9D0, 0.33D0, 0.7D0, 297c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.75D0, 0.2D0, 1.04D0/ 298c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DT9X/0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 299c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.78D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 300c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.78D0, -0.46D0, 0.0D0, 0.0D0, 301c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.78D0, -0.46D0, -0.22D0, 302c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 1.06D0, 0.0D0, 0.0D0, 0.0D0, 0.6D0, 0.0D0, 303c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.78D0, 304c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 305c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.66D0, 0.1D0, -0.1D0, 0.0D0, 0.0D0, 0.0D0, 306c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.96D0, 0.1D0, -0.76D0, 0.8D0, 0.90D0, 307c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.3D0, -0.02D0, 0.6D0, 0.0D0, 0.0D0, 0.0D0, 308c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.78D0, 0.0D0, 0.0D0, 309c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, -0.06D0, 0.1D0, 310c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.1D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.90D0, 311c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.1D0, -0.22D0, 0.8D0, 0.18D0, -0.3D0, -0.02D0, 312c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 313c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.78D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 314c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.78D0, 0.26D0, 0.0D0, 0.0D0, 0.0D0, 315c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.78D0, 0.26D0, -0.76D0, 1.12D0, 316c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0/ 317c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DT9Y/0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 318c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.04D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 319c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.04D0, -0.78D0, 0.0D0, 0.0D0, 320c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.04D0, -0.78D0, 0.54D0, 321c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.08D0, 0.0D0, 0.0D0, 0.0D0, 0.5D0, 0.0D0, 322c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.04D0, 323c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.7D0, 324c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.9D0, -0.12D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 325c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.64D0, -0.9D0, -0.30D0, 0.7D0, -0.18D0, 0.2D0, 326c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.28D0, 0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 327c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.04D0, 0.0D0, 0.0D0, 0.0D0, 328c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.7D0, -1.08D0, 0.0D0, 329c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.64D0, -1.26D0, 330c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.54D0, 0.20D0, 0.0D0, 0.0D0, 0.0D0, 0.5D0, 331c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 332c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.04D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 333c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.04D0, -0.9D0, 0.18D0, 0.0D0, 0.0D0, 334c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.04D0, -0.9D0, 0.18D0, 0.7D0, 335c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.18D0, 0.2D0, 0.16D0/ 336c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DT10X/0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 337c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 338c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.5D0, -0.9D0, 0.0D0, 0.0D0, 0.0D0, 339c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.5D0, -0.9D0, 0.3D0, 0.7D0, 340c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.6D0, 0.0D0, 0.0D0, 0.0D0, 341c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.5D0, 0.0D0, 0.0D0, 0.0D0, 342c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.3D0, 0.1D0, 0.5D0, 0.0D0, 343c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.8D0, 0.1D0, -0.6D0, 344c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.8D0, 0.3D0, -0.3D0, 0.5D0, 0.6D0, 0.0D0, 345c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.5D0, 0.0D0, 346c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, -0.9D0, 347c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.1D0, 0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.7D0, 348c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.1D0, 0.3D0, 0.8D0, -0.9D0, -0.3D0, 0.5D0, 349c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 350c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 351c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.5D0, 0.3D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 352c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.5D0, 0.3D0, -0.6D0, 0.8D0, 0.0D0, 0.0D0, 353c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0/ 354c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DT10Y/0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 355c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 356c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.6D0, 0.1D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 357c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.6D0, 0.1D0, -0.5D0, 0.8D0, 0.0D0, 358c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 359c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 360c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, -0.5D0, -0.9D0, 0.6D0, 0.0D0, 361c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, -0.4D0, -0.9D0, 0.9D0, 362c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.7D0, -0.5D0, 0.2D0, 0.6D0, 0.5D0, 0.0D0, 363c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.6D0, 0.0D0, 364c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, -0.5D0, 365c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 366c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.4D0, 0.9D0, -0.5D0, 0.6D0, 0.0D0, 0.0D0, 367c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 368c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 369c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.6D0, -0.9D0, 0.1D0, 0.0D0, 0.0D0, 370c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.6D0, -0.9D0, 0.1D0, 0.7D0, 371c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.5D0, 0.2D0, 0.8D0/ 372c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA SSIZE1/0.0D0, 0.3D0, 1.6D0, 3.2D0/ 373c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA SSIZE2/0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 374c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 375c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0, 376c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 1.17D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0, 377c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 1.17D0, 1.17D0, 1.17D0/ 378c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Executable Statements .. 379c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 380c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 120 KI = 1, 4 381c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INCX = INCXS(KI) 382c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INCY = INCYS(KI) 383c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MX = ABS(INCX) 384c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MY = ABS(INCY) 385c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 386c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 100 KN = 1, 4 387c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath N = NS(KN) 388c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath KSIZE = MIN(2,KN) 389c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LENX = LENS(KN,MX) 390c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LENY = LENS(KN,MY) 391c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Initialize all argument arrays .. 392c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 20 I = 1, 7 393c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SX(I) = DX1(I) 394c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SY(I) = DY1(I) 395c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 20 CONTINUE 396c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 397c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF (ICASE.EQ.1) THEN 398c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. DDOT .. 399c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST1(DDOT(N,SX,INCX,SY,INCY),DT7(KN,KI),SSIZE1(KN) 400c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + ,SFAC) 401c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE IF (ICASE.EQ.2) THEN 402c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. DAXPY .. 403c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL DAXPY(N,SA,SX,INCX,SY,INCY) 404c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 40 J = 1, LENY 405c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STY(J) = DT8(J,KN,KI) 406c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 40 CONTINUE 407c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST(LENY,SY,STY,SSIZE2(1,KSIZE),SFAC) 408c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE IF (ICASE.EQ.5) THEN 409c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. DCOPY .. 410c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 60 I = 1, 7 411c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STY(I) = DT10Y(I,KN,KI) 412c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 60 CONTINUE 413c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL DCOPY(N,SX,INCX,SY,INCY) 414c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST(LENY,SY,STY,SSIZE2(1,1),1.0D0) 415c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE IF (ICASE.EQ.6) THEN 416c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. DSWAP .. 417c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL DSWAP(N,SX,INCX,SY,INCY) 418c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 80 I = 1, 7 419c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STX(I) = DT10X(I,KN,KI) 420c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STY(I) = DT10Y(I,KN,KI) 421c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 80 CONTINUE 422c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST(LENX,SX,STX,SSIZE2(1,1),1.0D0) 423c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST(LENY,SY,STY,SSIZE2(1,1),1.0D0) 424c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE 425c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,*) ' Shouldn''t be here in CHECK2' 426c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STOP 427c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END IF 428c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 100 CONTINUE 429c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 120 CONTINUE 430c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath RETURN 431c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END 432c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SUBROUTINE CHECK3(SFAC) 433c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Parameters .. 434c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER NOUT 435c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PARAMETER (NOUT=6) 436c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalar Arguments .. 437c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SFAC 438c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalars in Common .. 439c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ICASE, INCX, INCY, MODE, N 440c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LOGICAL PASS 441c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Scalars .. 442c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SA, SC, SS 443c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER I, K, KI, KN, KSIZE, LENX, LENY, MX, MY 444c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Arrays .. 445c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION COPYX(5), COPYY(5), DT9X(7,4,4), DT9Y(7,4,4), 446c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + DX1(7), DY1(7), MWPC(11), MWPS(11), MWPSTX(5), 447c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + MWPSTY(5), MWPTX(11,5), MWPTY(11,5), MWPX(5), 448c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + MWPY(5), SSIZE2(14,2), STX(7), STY(7), SX(7), 449c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + SY(7) 450c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER INCXS(4), INCYS(4), LENS(4,2), MWPINX(11), 451c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + MWPINY(11), MWPN(11), NS(4) 452c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. External Subroutines .. 453c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath EXTERNAL DROT, STEST 454c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Intrinsic Functions .. 455c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTRINSIC ABS, MIN 456c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Common blocks .. 457c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS 458c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Data statements .. 459c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA SA/0.3D0/ 460c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA INCXS/1, 2, -2, -1/ 461c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA INCYS/1, -2, 1, -2/ 462c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA LENS/1, 1, 2, 4, 1, 1, 3, 7/ 463c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA NS/0, 1, 2, 4/ 464c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DX1/0.6D0, 0.1D0, -0.5D0, 0.8D0, 0.9D0, -0.3D0, 465c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.4D0/ 466c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DY1/0.5D0, -0.9D0, 0.3D0, 0.7D0, -0.6D0, 0.2D0, 467c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.8D0/ 468c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA SC, SS/0.8D0, 0.6D0/ 469c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DT9X/0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 470c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.78D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 471c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.78D0, -0.46D0, 0.0D0, 0.0D0, 472c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.78D0, -0.46D0, -0.22D0, 473c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 1.06D0, 0.0D0, 0.0D0, 0.0D0, 0.6D0, 0.0D0, 474c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.78D0, 475c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 476c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.66D0, 0.1D0, -0.1D0, 0.0D0, 0.0D0, 0.0D0, 477c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.96D0, 0.1D0, -0.76D0, 0.8D0, 0.90D0, 478c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.3D0, -0.02D0, 0.6D0, 0.0D0, 0.0D0, 0.0D0, 479c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.78D0, 0.0D0, 0.0D0, 480c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, -0.06D0, 0.1D0, 481c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.1D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.90D0, 482c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.1D0, -0.22D0, 0.8D0, 0.18D0, -0.3D0, -0.02D0, 483c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 484c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.78D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 485c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.78D0, 0.26D0, 0.0D0, 0.0D0, 0.0D0, 486c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.78D0, 0.26D0, -0.76D0, 1.12D0, 487c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0/ 488c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA DT9Y/0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 489c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.04D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 490c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.04D0, -0.78D0, 0.0D0, 0.0D0, 491c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.04D0, -0.78D0, 0.54D0, 492c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.08D0, 0.0D0, 0.0D0, 0.0D0, 0.5D0, 0.0D0, 493c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.04D0, 494c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.7D0, 495c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.9D0, -0.12D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 496c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.64D0, -0.9D0, -0.30D0, 0.7D0, -0.18D0, 0.2D0, 497c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.28D0, 0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 498c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.04D0, 0.0D0, 0.0D0, 0.0D0, 499c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.7D0, -1.08D0, 0.0D0, 500c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.64D0, -1.26D0, 501c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.54D0, 0.20D0, 0.0D0, 0.0D0, 0.0D0, 0.5D0, 502c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 503c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.04D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 504c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.04D0, -0.9D0, 0.18D0, 0.0D0, 0.0D0, 505c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.04D0, -0.9D0, 0.18D0, 0.7D0, 506c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + -0.18D0, 0.2D0, 0.16D0/ 507c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DATA SSIZE2/0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 508c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 509c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 0.0D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0, 510c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 1.17D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0, 511c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + 1.17D0, 1.17D0, 1.17D0/ 512c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Executable Statements .. 513c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 514c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 60 KI = 1, 4 515c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INCX = INCXS(KI) 516c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INCY = INCYS(KI) 517c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MX = ABS(INCX) 518c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MY = ABS(INCY) 519c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 520c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 40 KN = 1, 4 521c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath N = NS(KN) 522c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath KSIZE = MIN(2,KN) 523c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LENX = LENS(KN,MX) 524c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LENY = LENS(KN,MY) 525c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 526c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF (ICASE.EQ.4) THEN 527c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. DROT .. 528c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 20 I = 1, 7 529c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SX(I) = DX1(I) 530c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SY(I) = DY1(I) 531c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STX(I) = DT9X(I,KN,KI) 532c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STY(I) = DT9Y(I,KN,KI) 533c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 20 CONTINUE 534c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL DROT(N,SX,INCX,SY,INCY,SC,SS) 535c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST(LENX,SX,STX,SSIZE2(1,KSIZE),SFAC) 536c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST(LENY,SY,STY,SSIZE2(1,KSIZE),SFAC) 537c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ELSE 538c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,*) ' Shouldn''t be here in CHECK3' 539c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STOP 540c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END IF 541c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 40 CONTINUE 542c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 60 CONTINUE 543c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 544c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPC(1) = 1 545c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 80 I = 2, 11 546c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPC(I) = 0 547c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 80 CONTINUE 548c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPS(1) = 0 549c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 100 I = 2, 6 550c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPS(I) = 1 551c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 100 CONTINUE 552c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 120 I = 7, 11 553c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPS(I) = -1 554c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 120 CONTINUE 555c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(1) = 1 556c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(2) = 1 557c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(3) = 1 558c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(4) = -1 559c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(5) = 1 560c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(6) = -1 561c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(7) = 1 562c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(8) = 1 563c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(9) = -1 564c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(10) = 1 565c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINX(11) = -1 566c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(1) = 1 567c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(2) = 1 568c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(3) = -1 569c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(4) = -1 570c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(5) = 2 571c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(6) = 1 572c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(7) = 1 573c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(8) = -1 574c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(9) = -1 575c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(10) = 2 576c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPINY(11) = 1 577c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 140 I = 1, 11 578c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPN(I) = 5 579c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 140 CONTINUE 580c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPN(5) = 3 581c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPN(10) = 3 582c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 160 I = 1, 5 583c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPX(I) = I 584c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPY(I) = I 585c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(1,I) = I 586c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(1,I) = I 587c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(2,I) = I 588c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(2,I) = -I 589c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(3,I) = 6 - I 590c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(3,I) = I - 6 591c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(4,I) = I 592c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(4,I) = -I 593c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(6,I) = 6 - I 594c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(6,I) = I - 6 595c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(7,I) = -I 596c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(7,I) = I 597c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(8,I) = I - 6 598c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(8,I) = 6 - I 599c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(9,I) = -I 600c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(9,I) = I 601c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(11,I) = I - 6 602c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(11,I) = 6 - I 603c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 160 CONTINUE 604c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(5,1) = 1 605c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(5,2) = 3 606c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(5,3) = 5 607c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(5,4) = 4 608c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(5,5) = 5 609c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(5,1) = -1 610c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(5,2) = 2 611c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(5,3) = -2 612c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(5,4) = 4 613c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(5,5) = -3 614c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(10,1) = -1 615c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(10,2) = -3 616c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(10,3) = -5 617c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(10,4) = 4 618c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTX(10,5) = 5 619c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(10,1) = 1 620c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(10,2) = 2 621c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(10,3) = 2 622c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(10,4) = 4 623c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPTY(10,5) = 3 624c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 200 I = 1, 11 625c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INCX = MWPINX(I) 626c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INCY = MWPINY(I) 627c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 180 K = 1, 5 628c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath COPYX(K) = MWPX(K) 629c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath COPYY(K) = MWPY(K) 630c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPSTX(K) = MWPTX(I,K) 631c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MWPSTY(K) = MWPTY(I,K) 632c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 180 CONTINUE 633c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL DROT(MWPN(I),COPYX,INCX,COPYY,INCY,MWPC(I),MWPS(I)) 634c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST(5,COPYX,MWPSTX,MWPSTX,SFAC) 635c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST(5,COPYY,MWPSTY,MWPSTY,SFAC) 636c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 200 CONTINUE 637c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath RETURN 638c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END 639c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SUBROUTINE STEST(LEN,SCOMP,STRUE,SSIZE,SFAC) 640c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* ********************************* STEST ************************** 641c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 642c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* THIS SUBR COMPARES ARRAYS SCOMP() AND STRUE() OF LENGTH LEN TO 643c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* SEE IF THE TERM BY TERM DIFFERENCES, MULTIPLIED BY SFAC, ARE 644c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* NEGLIGIBLE. 645c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 646c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* C. L. LAWSON, JPL, 1974 DEC 10 647c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 648c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Parameters .. 649c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER NOUT 650c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PARAMETER (NOUT=6) 651c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalar Arguments .. 652c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SFAC 653c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER LEN 654c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Array Arguments .. 655c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SCOMP(LEN), SSIZE(LEN), STRUE(LEN) 656c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalars in Common .. 657c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ICASE, INCX, INCY, MODE, N 658c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LOGICAL PASS 659c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Scalars .. 660c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SD 661c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER I 662c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. External Functions .. 663c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SDIFF 664c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath EXTERNAL SDIFF 665c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Intrinsic Functions .. 666c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTRINSIC ABS 667c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Common blocks .. 668c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS 669c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Executable Statements .. 670c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 671c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DO 40 I = 1, LEN 672c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SD = SCOMP(I) - STRUE(I) 673c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF (SDIFF(ABS(SSIZE(I))+ABS(SFAC*SD),ABS(SSIZE(I))).EQ.0.0D0) 674c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + GO TO 40 675c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 676c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* HERE SCOMP(I) IS NOT CLOSE TO STRUE(I). 677c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 678c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF ( .NOT. PASS) GO TO 20 679c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* PRINT FAIL MESSAGE AND HEADER. 680c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PASS = .FALSE. 681c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,99999) 682c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,99998) 683c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 20 WRITE (NOUT,99997) ICASE, N, INCX, INCY, MODE, I, SCOMP(I), 684c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + STRUE(I), SD, SSIZE(I) 685c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 40 CONTINUE 686c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath RETURN 687c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 688c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath99999 FORMAT (' FAIL') 689c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath99998 FORMAT (/' CASE N INCX INCY MODE I ', 690c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + ' COMP(I) TRUE(I) DIFFERENCE', 691c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + ' SIZE(I)',/1X) 692c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath99997 FORMAT (1X,I4,I3,3I5,I3,2D36.8,2D12.4) 693c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END 694c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SUBROUTINE STEST1(SCOMP1,STRUE1,SSIZE,SFAC) 695c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* ************************* STEST1 ***************************** 696c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 697c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* THIS IS AN INTERFACE SUBROUTINE TO ACCOMODATE THE FORTRAN 698c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* REQUIREMENT THAT WHEN A DUMMY ARGUMENT IS AN ARRAY, THE 699c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* ACTUAL ARGUMENT MUST ALSO BE AN ARRAY OR AN ARRAY ELEMENT. 700c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 701c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* C.L. LAWSON, JPL, 1978 DEC 6 702c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 703c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalar Arguments .. 704c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SCOMP1, SFAC, STRUE1 705c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Array Arguments .. 706c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SSIZE(*) 707c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Arrays .. 708c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SCOMP(1), STRUE(1) 709c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. External Subroutines .. 710c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath EXTERNAL STEST 711c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Executable Statements .. 712c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 713c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SCOMP(1) = SCOMP1 714c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath STRUE(1) = STRUE1 715c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL STEST(1,SCOMP,STRUE,SSIZE,SFAC) 716c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 717c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath RETURN 718c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END 719c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION FUNCTION SDIFF(SA,SB) 720c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* ********************************* SDIFF ************************** 721c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* COMPUTES DIFFERENCE OF TWO NUMBERS. C. L. LAWSON, JPL 1974 FEB 15 722c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 723c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalar Arguments .. 724c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DOUBLE PRECISION SA, SB 725c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Executable Statements .. 726c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SDIFF = SA - SB 727c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath RETURN 728c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END 729c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SUBROUTINE ITEST1(ICOMP,ITRUE) 730c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* ********************************* ITEST1 ************************* 731c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 732c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* THIS SUBROUTINE COMPARES THE VARIABLES ICOMP AND ITRUE FOR 733c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* EQUALITY. 734c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* C. L. LAWSON, JPL, 1974 DEC 10 735c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 736c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Parameters .. 737c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER NOUT 738c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PARAMETER (NOUT=6) 739c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalar Arguments .. 740c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ICOMP, ITRUE 741c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Scalars in Common .. 742c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ICASE, INCX, INCY, MODE, N 743c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath LOGICAL PASS 744c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Local Scalars .. 745c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INTEGER ID 746c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Common blocks .. 747c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS 748c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* .. Executable Statements .. 749c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 750c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF (ICOMP.EQ.ITRUE) GO TO 40 751c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 752c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* HERE ICOMP IS NOT EQUAL TO ITRUE. 753c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 754c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath IF ( .NOT. PASS) GO TO 20 755c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* PRINT FAIL MESSAGE AND HEADER. 756c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath PASS = .FALSE. 757c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,99999) 758c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,99998) 759c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 20 ID = ICOMP - ITRUE 760c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath WRITE (NOUT,99997) ICASE, N, INCX, INCY, MODE, ICOMP, ITRUE, ID 761c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 40 CONTINUE 762c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath RETURN 763c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath* 764c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath99999 FORMAT (' FAIL') 765c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath99998 FORMAT (/' CASE N INCX INCY MODE ', 766c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + ' COMP TRUE DIFFERENCE', 767c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath + /1X) 768c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath99997 FORMAT (1X,I4,I3,3I5,2I36,I12) 769c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath END 770