1
        2
        3
        4
        5
        6
        7
        8
        9
       10
       11
       12
       13
       14
       15
       16
       17
       18
       19
       20
       21
       22
       23
       24
       25
       26
       27
       28
       29
       30
       31
       32
       33
       34
       35
       36
       37
       38
       39
       40
       41
       42
       43
       44
       45
       46
       47
       48
       49
       50
       51
       52
       53
       54
       55
       56
       57
       58
       59
       60
       61
       62
       63
       64
       65
       66
       67
       68
       69
       70
       71
       72
       73
       74
       75
       76
       77
       78
       79
       80
       81
       82
       83
       84
       85
       86
       87
       88
       89
       90
       91
       92
       93
       94
       95
       96
       97
       98
       99
      100
      101
      102
      103
      104
      105
      106
      107
      108
      109
      110
      111
      112
      113
      114
      115
      116
      117
      118
      119
      120
      121
      122
      123
      124
      125
      126
      127
      128
      129
      130
      131
      132
      133
      134
      135
      136
      137
      138
      139
      140
      141
      142
      143
      144
      145
      146
      147
      148
      149
      150
      151
      152
      153
      154
      155
      156
      157
      158
      159
      160
      161
      162
      163
      164
      165
      166
      167
      168
      169
      170
      171
      172
      173
      174
      175
      176
      177
      178
      179
      180
      181
      182
      183
      184
      185
      186
      187
      188
      189
      190
      191
      192
      193
      194
      195
      196
      197
      198
      199
      200
      201
      202
      203
      204
      205
      206
      207
      208
      209
      210
      211
      212
      213
      214
      215
      216
      217
      218
      219
      220
      221
      222
      223
      224
      225
      226
      227
      228
      229
      230
      231
      232
      233
      234
      235
      236
      237
      238
      239
      240
      241
      242
      243
      244
      245
      246
      247
      248
      249
      250
      251
      252
      253
      254
      255
      256
      257
      258
      259
      260
      261
      262
      263
      264
      265
      266
      267
      268
      269
      270
      271
      272
      273
      274
      275
      276
      277
      278
      279
      280
      281
      282
      283
      284
      285
      286
      287
      288
      289
      290
      291
      292
      293
      294
      295
      296
      297
      298
      299
      300
      301
      302
      303
      304
      305
      306
      307
      308
      309
      310
      311
      312
      313
      314
      315
      316
      317
      318
      319
      320
      321
      322
      323
      324
      325
      326
      327
      328
      329
      330
      331
      332
      333
      334
      335
      336
      337
      338
      339
      340
      341
      342
      343
      344
      345
      346
      347
      348
      349
      350
      351
      352
      353
      354
      355
      356
      357
      358
      359
      360
      361
      362
      363
      364
      365
      366
      367
      368
      369
      370
      371
      372
      373
      374
      375
      376
      377
      378
      379
      380
      381
      382
      383
      384
      385
      386
      387
      388
      389
      390
      391
      392
      393
      394
      395
      396
      397
      398
      399
      400
      401
      402
      403
      404
      405
      406
      407
      408
      409
      410
      411
      412
      413
      414
      415
      416
      417
      418
      419
      420
      421
      422
      423
      424
      425
      426
      427
      428
      429
      430
      431
      432
      433
      434
      435
      436
      437
      438
      439
      440
      441
      442
      443
      444
      445
      446
      447
      448
      449
      450
      451
      452
      453
      454
      455
      456
      457
      458
      459
      460
      461
      462
      463
      464
      465
      466
      467
      468
      469
      470
      471
      472
      473
      474
      475
      476
      477
      478
      479
      480
      481
      482
      483
      484
      485
      486
      487
      488
      489
      490
      491
      492
      493
      494
      495
      496
      497
      498
      499
      500
      501
      502
      503
      504
      505
      506
      507
      508
      509
      510
      511
      512
      513
      514
      515
      516
      517
      518
      519
      520
      521
      522
      523
      524
      525
      526
      527
      528
      529
      530
      531
      532
      533
      534
      535
      536
      537
      538
      539
      540
      541
      542
      543
      544
      545
      546
      547
      548
      549
      550
      551
      552
      553
      554
      555
      556
      557
      558
      559
      560
      561
      562
      563
      564
      565
      566
      567
      568
      569
      570
      571
      572
      573
      574
      575
      576
      577
      578
      579
      580
      581
      582
      583
      584
      585
      586
      587
      588
      589
      590
      591
      592
      593
      594
      595
      596
      597
      598
      599
      600
      601
      602
      603
      604
      605
      606
      607
      608
      609
      610
      611
      612
      613
      614
      615
      616
      617
      618
      619
      620
      621
      622
      623
      624
      625
      626
      627
      628
      629
      630
      631
      632
      633
      634
      635
      636
      637
      638
      639
      640
      641
      642
      643
      644
      645
      646
      647
      648
      649
      650
      651
      652
      653
      654
      655
      656
      657
      658
      659
      660
      661
      662
      663
      664
      665
      666
      667
      668
      669
      670
      671
      672
      673
      674
      675
      676
      677
      678
      679
      680
      681
      682
      683
      684
      685
      686
      687
      688
      689
      690
      691
      692
      693
      694
      695
      696
      697
      698
      699
      700
      701
      702
      703
      704
      705
      706
      707
      708
      709
      710
      711
      712
      713
      714
      715
      716
      717
      718
      719
      720
      721
      722
      723
      724
      725
      726
      727
      728
      729
      730
      731
      732
      733
      734
      735
      736
      737
      738
      739
      740
      741
      742
      743
      744
      745
      746
      747
      748
      749
      750
      751
      752
      753
      754
      755
      756
      757
      758
      759
      760
      761
      762
      763
      764
      765
      766
      767
      768
      769
      770
      771
      772
      773
      774
      775
      776
      777
      778
      779
      780
      781
      782
      783
      784
      785
      786
      787
      788
      789
      790
      791
      792
      793
      794
      795
      796
      797
      798
      799
      800
      801
      802
      803
      804
      805
      806
      807
      808
      809
      810
      811
      812
      813
      814
      815
      816
      817
      818
      819
      820
      821
      822
      823
      824
      825
      826
      827
      828
      829
      830
      831
      832
      833
      834
      835
      836
      837
      838
      839
      840
      841
      842
      843
      844
      845
      846
      847
      848
      849
      850
      851
      852
      853
      854
      855
      856
      857
      858
      859
      860
      861
      862
      863
      864
      865
      866
      867
      868
      869
      870
      871
      872
      873
      874
      875
      876
      877
      878
      879
      880
      881
      882
      883
      884
      885
      886
      887
      888
      889
      890
      891
      892
      893
      894
      895
      896
      897
      898
      899
      900
      901
      902
      903
      904
      905
      906
      907
      908
      909
      910
      911
      912
      913
      914
      915
      916
      917
      918
      919
      920
      921
      922
      923
      924
      925
      926
      927
      928
      929
      930
      931
      932
      933
      934
      935
      936
      937
      938
      939
      940
      941
      942
      943
      944
      945
      946
      947
      948
      949
      950
      951
      952
      953
      954
      955
      956
      957
      958
      959
      960
      961
      962
      963
      964
      965
      966
      967
      968
      969
      970
      971
      972
      973
      974
      975
      976
      977
      978
      979
      980
      981
      982
      983
      984
      985
      986
      987
      988
      989
      990
      991
      992
      993
      994
      995
      996
      997
      998
      999
     1000
     1001
     1002
     1003
     1004
     1005
     1006
     1007
     1008
     1009
     1010
     1011
     1012
     1013
     1014
     1015
     1016
     1017
     1018
     1019
     1020
     1021
     1022
     1023
     1024
     1025
     1026
     1027
     1028
     1029
     1030
     1031
     1032
     1033
     1034
     1035
     1036
     1037
     1038
     1039
     1040
     1041
     1042
     1043
     1044
     1045
     1046
     1047
     1048
     1049
     1050
     1051
     1052
     1053
     1054
     1055
     1056
     1057
     1058
     1059
     1060
     1061
     1062
     1063
     1064
     1065
     1066
     1067
     1068
     1069
     1070
     1071
     1072
     1073
     1074
     1075
     1076
     1077
     1078
     1079
     1080
     1081
     1082
     1083
     1084
     1085
     1086
     1087
     1088
     1089
     1090
     1091
     1092
     1093
     1094
     1095
     1096
     1097
     1098
     1099
     1100
     1101
     1102
     1103
     1104
     1105
     1106
     1107
     1108
     1109
     1110
     1111
     1112
     1113
     1114
     1115
     1116
     1117
     1118
     1119
     1120
     1121
     1122
     1123
     1124
     1125
     1126
     1127
     1128
     1129
     1130
     1131
     1132
     1133
     1134
     1135
     1136
     1137
     1138
     1139
     1140
     1141
     1142
     1143
     1144
     1145
     1146
     1147
     1148
     1149
     1150
     1151
     1152
     1153
     1154
     1155
     1156
     1157
     1158
     1159
     1160
     1161
     1162
     1163
     1164
     1165
     1166
     1167
     1168
     1169
     1170
     1171
     1172
     1173
     1174
     1175
     1176
     1177
     1178
     1179
     1180
     1181
     1182
     1183
     1184
     1185
     1186
     1187
     1188
     1189
     1190
     1191
     1192
     1193
     1194
     1195
     1196
     1197
     1198
     1199
     1200
     1201
     1202
     1203
     1204
     1205
     1206
     1207
     1208
     1209
     1210
     1211
     1212
     1213
     1214
     1215
     1216
     1217
     1218
     1219
     1220
     1221
     1222
     1223
     1224
     1225
     1226
     1227
     1228
     1229
     1230
     1231
     1232
     1233
     1234
     1235
     1236
     1237
     1238
     1239
     1240
     1241
     1242
     1243
     1244
     1245
     1246
     1247
     1248
     1249
     1250
     1251
     1252
     1253
     1254
     1255
     1256
     1257
     1258
     1259
     1260
     1261
     1262
     1263
     1264
     1265
     1266
     1267
     1268
     1269
     1270
     1271
     1272
     1273
     1274
     1275
     1276
     1277
     1278
     1279
     1280
     1281
     1282
     1283
     1284
     1285
     1286
     1287
     1288
     1289
     1290
     1291
     1292
     1293
     1294
     1295
     1296
     1297
     1298
     1299
     1300
     1301
     1302
     1303
     1304
     1305
     1306
     1307
     1308
     1309
     1310
     1311
     1312
     1313
     1314
     1315
     1316
     1317
     1318
     1319
     1320
     1321
     1322
     1323
     1324
     1325
     1326
     1327
     1328
     1329
     1330
     1331
     1332
     1333
     1334
     1335
     1336
     1337
     1338
     1339
     1340
     1341
     1342
     1343
     1344
     1345
     1346
     1347
     1348
     1349
     1350
     1351
     1352
     1353
     1354
     1355
     1356
     1357
     1358
     1359
     1360
     1361
     1362
     1363
     1364
     1365
     1366
     1367
     1368
     1369
     1370
     1371
     1372
     1373
     1374
     1375
     1376
     1377
     1378
     1379
     1380
     1381
     1382
     1383
     1384
     1385
     1386
     1387
     1388
     1389
     1390
     1391
     1392
     1393
     1394
     1395
     1396
     1397
     1398
     1399
     1400
     1401
     1402
     1403
     1404
     1405
     1406
     1407
     1408
     1409
     1410
     1411
     1412
     1413
     1414
     1415
     1416
     1417
     1418
     1419
     1420
     1421
     1422
     1423
     1424
     1425
     1426
     1427
     1428
     1429
     1430
     1431
     1432
     1433
     1434
     1435
     1436
     1437
     1438
     1439
     1440
     1441
     1442
     1443
     1444
     1445
     1446
     1447
     1448
     1449
     1450
     1451
     1452
     1453
     1454
     1455
     1456
     1457
     1458
     1459
     1460
     1461
     1462
     1463
     1464
     1465
     1466
     1467
     1468
     1469
     1470
     1471
     1472
     1473
     1474
     1475
     1476
     1477
     1478
     1479
     1480
     1481
     1482
     1483
     1484
     1485
     1486
     1487
     1488
     1489
     1490
     1491
     1492
     1493
     1494
     1495
     1496
     1497
     1498
     1499
     1500
     1501
     1502
     1503
     1504
     1505
     1506
     1507
     1508
     1509
     1510
     1511
     1512
     1513
     1514
     1515
     1516
     1517
     1518
     1519
     1520
     1521
     1522
     1523
     1524
     1525
     1526
     1527
     1528
     1529
     1530
     1531
     1532
     1533
     1534
     1535
     1536
     1537
     1538
     1539
     1540
     1541
     1542
     1543
     1544
     1545
     1546
     1547
     1548
     1549
     1550
     1551
     1552
     1553
     1554
     1555
     1556
     1557
     1558
     1559
     1560
     1561
     1562
     1563
     1564
     1565
     1566
     1567
     1568
     1569
     1570
     1571
     1572
     1573
     1574
     1575
     1576
     1577
     1578
     1579
     1580
     1581
     1582
     1583
     1584
     1585
     1586
     1587
     1588
     1589
     1590
     1591
     1592
     1593
     1594
     1595
     1596
     1597
     1598
     1599
     1600
     1601
     1602
     1603
     1604
     1605
     1606
     1607
     1608
     1609
     1610
     1611
     1612
     1613
     1614
     1615
     1616
     1617
     1618
     1619
     1620
     1621
     1622
     1623
     1624
     1625
     1626
     1627
     1628
     1629
     1630
     1631
     1632
     1633
     1634
     1635
     1636
     1637
     1638
     1639
     1640
     1641
     1642
     1643
     1644
     1645
     1646
     1647
     1648
     1649
     1650
     1651
     1652
     1653
     1654
     1655
     1656
     1657
     1658
     1659
     1660
     1661
     1662
     1663
     1664
     1665
     1666
     1667
     1668
     1669
     1670
     1671
     1672
     1673
     1674
     1675
     1676
     1677
     1678
     1679
     1680
     1681
     1682
     1683
     1684
     1685
     1686
     1687
     1688
     1689
     1690
     1691
     1692
     1693
     1694
     1695
     1696
     1697
     1698
     1699
     1700
     1701
     1702
     1703
     1704
     1705
     1706
     1707
     1708
     1709
     1710
     1711
     1712
     1713
     1714
     1715
     1716
     1717
     1718
     1719
     1720
     1721
     1722
     1723
     1724
     1725
     1726
     1727
     1728
     1729
     1730
     1731
     1732
     1733
     1734
     1735
     1736
     1737
     1738
     1739
     1740
     1741
     1742
     1743
     1744
     1745
     1746
     1747
     1748
     1749
     1750
     1751
     1752
     1753
     1754
     1755
     1756
     1757
     1758
     1759
     1760
     1761
     1762
     1763
     1764
     1765
     1766
     1767
     1768
     1769
     1770
     1771
     1772
     1773
     1774
     1775
     1776
     1777
     1778
     1779
     1780
     1781
     1782
     1783
     1784
     1785
     1786
     1787
     1788
     1789
     1790
     1791
     1792
     1793
     1794
     1795
     1796
     1797
     1798
     1799
     1800
     1801
     1802
     1803
     1804
     1805
     1806
     1807
     1808
     1809
     1810
     1811
     1812
     1813
     1814
     1815
     1816
     1817
     1818
     1819
     1820
     1821
     1822
     1823
      SUBROUTINE DCHKSTNSIZESNNNTYPESDOTYPEISEEDTHRESH,
     $                   NOUNITALDAAPSDSED1D2D3D4D5,
     $                   WA1WA2WA3WRULDUVVPTAUZWORK,
     $                   LWORKIWORKLIWORKRESULTINFO )
      IMPLICIT NONE
*
*  -- LAPACK test routine (version 3.1) --
*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
*     November 2006
*
*     .. Scalar Arguments ..
      INTEGER            INFOLDALDULIWORKLWORKNOUNITNSIZES,
     $                   NTYPES
      DOUBLE PRECISION   THRESH
*     ..
*     .. Array Arguments ..
      LOGICAL            DOTYPE* )
      INTEGER            ISEED4 ), IWORK* ), NN* )
      DOUBLE PRECISION   ALDA* ), AP* ), D1* ), D2* ),
     $                   D3* ), D4* ), D5* ), RESULT* ),
     $                   SD* ), SE* ), TAU* ), ULDU* ),
     $                   VLDU* ), VP* ), WA1* ), WA2* ),
     $                   WA3* ), WORK* ), WR* ), ZLDU* )
*     ..
*
*  Purpose
*  =======
*
*  DCHKST  checks the symmetric eigenvalue problem routines.
*
*     DSYTRD factors A as  U S U' , where ' means transpose,
*     S is symmetric tridiagonal, and U is orthogonal.
*     DSYTRD can use either just the lower or just the upper triangle
*     of A; DCHKST checks both cases.
*     U is represented as a product of Householder
*     transformations, whose vectors are stored in the first
*     n-1 columns of V, and whose scale factors are in TAU.
*
*     DSPTRD does the same as DSYTRD, except that A and V are stored
*     in "packed" format.
*
*     DORGTR constructs the matrix U from the contents of V and TAU.
*
*     DOPGTR constructs the matrix U from the contents of VP and TAU.
*
*     DSTEQR factors S as  Z D1 Z' , where Z is the orthogonal
*     matrix of eigenvectors and D1 is a diagonal matrix with
*     the eigenvalues on the diagonal.  D2 is the matrix of
*     eigenvalues computed when Z is not computed.
*
*     DSTERF computes D3, the matrix of eigenvalues, by the
*     PWK method, which does not yield eigenvectors.
*
*     DPTEQR factors S as  Z4 D4 Z4' , for a
*     symmetric positive definite tridiagonal matrix.
*     D5 is the matrix of eigenvalues computed when Z is not
*     computed.
*
*     DSTEBZ computes selected eigenvalues.  WA1, WA2, and
*     WA3 will denote eigenvalues computed to high
*     absolute accuracy, with different range options.
*     WR will denote eigenvalues computed to high relative
*     accuracy.
*
*     DSTEIN computes Y, the eigenvectors of S, given the
*     eigenvalues.
*
*     DSTEDC factors S as Z D1 Z' , where Z is the orthogonal
*     matrix of eigenvectors and D1 is a diagonal matrix with
*     the eigenvalues on the diagonal ('I' option). It may also
*     update an input orthogonal matrix, usually the output
*     from DSYTRD/DORGTR or DSPTRD/DOPGTR ('V' option). It may
*     also just compute eigenvalues ('N' option).
*
*     DSTEMR factors S as Z D1 Z' , where Z is the orthogonal
*     matrix of eigenvectors and D1 is a diagonal matrix with
*     the eigenvalues on the diagonal ('I' option).  DSTEMR
*     uses the Relatively Robust Representation whenever possible.
*
*  When DCHKST is called, a number of matrix "sizes" ("n's") and a
*  number of matrix "types" are specified.  For each size ("n")
*  and each type of matrix, one matrix will be generated and used
*  to test the symmetric eigenroutines.  For each matrix, a number
*  of tests will be performed:
*
*  (1)     | A - V S V' | / ( |A| n ulp ) DSYTRD( UPLO='U', ... )
*
*  (2)     | I - UV' | / ( n ulp )        DORGTR( UPLO='U', ... )
*
*  (3)     | A - V S V' | / ( |A| n ulp ) DSYTRD( UPLO='L', ... )
*
*  (4)     | I - UV' | / ( n ulp )        DORGTR( UPLO='L', ... )
*
*  (5-8)   Same as 1-4, but for DSPTRD and DOPGTR.
*
*  (9)     | S - Z D Z' | / ( |S| n ulp ) DSTEQR('V',...)
*
*  (10)    | I - ZZ' | / ( n ulp )        DSTEQR('V',...)
*
*  (11)    | D1 - D2 | / ( |D1| ulp )        DSTEQR('N',...)
*
*  (12)    | D1 - D3 | / ( |D1| ulp )        DSTERF
*
*  (13)    0 if the true eigenvalues (computed by sturm count)
*          of S are within THRESH of
*          those in D1.  2*THRESH if they are not.  (Tested using
*          DSTECH)
*
*  For S positive definite,
*
*  (14)    | S - Z4 D4 Z4' | / ( |S| n ulp ) DPTEQR('V',...)
*
*  (15)    | I - Z4 Z4' | / ( n ulp )        DPTEQR('V',...)
*
*  (16)    | D4 - D5 | / ( 100 |D4| ulp )       DPTEQR('N',...)
*
*  When S is also diagonally dominant by the factor gamma < 1,
*
*  (17)    max | D4(i) - WR(i) | / ( |D4(i)| omega ) ,
*           i
*          omega = 2 (2n-1) ULP (1 + 8 gamma**2) / (1 - gamma)**4
*                                               DSTEBZ( 'A', 'E', ...)
*
*  (18)    | WA1 - D3 | / ( |D3| ulp )          DSTEBZ( 'A', 'E', ...)
*
*  (19)    ( max { min | WA2(i)-WA3(j) | } +
*             i     j
*            max { min | WA3(i)-WA2(j) | } ) / ( |D3| ulp )
*             i     j
*                                               DSTEBZ( 'I', 'E', ...)
*
*  (20)    | S - Y WA1 Y' | / ( |S| n ulp )  DSTEBZ, SSTEIN
*
*  (21)    | I - Y Y' | / ( n ulp )          DSTEBZ, SSTEIN
*
*  (22)    | S - Z D Z' | / ( |S| n ulp )    DSTEDC('I')
*
*  (23)    | I - ZZ' | / ( n ulp )           DSTEDC('I')
*
*  (24)    | S - Z D Z' | / ( |S| n ulp )    DSTEDC('V')
*
*  (25)    | I - ZZ' | / ( n ulp )           DSTEDC('V')
*
*  (26)    | D1 - D2 | / ( |D1| ulp )           DSTEDC('V') and
*                                               DSTEDC('N')
*
*  Test 27 is disabled at the moment because DSTEMR does not
*  guarantee high relatvie accuracy.
*
*  (27)    max | D6(i) - WR(i) | / ( |D6(i)| omega ) ,
*           i
*          omega = 2 (2n-1) ULP (1 + 8 gamma**2) / (1 - gamma)**4
*                                               DSTEMR('V', 'A')
*
*  (28)    max | D6(i) - WR(i) | / ( |D6(i)| omega ) ,
*           i
*          omega = 2 (2n-1) ULP (1 + 8 gamma**2) / (1 - gamma)**4
*                                               DSTEMR('V', 'I')
*
*  Tests 29 through 34 are disable at present because DSTEMR
*  does not handle partial specturm requests.
*
*  (29)    | S - Z D Z' | / ( |S| n ulp )    DSTEMR('V', 'I')
*
*  (30)    | I - ZZ' | / ( n ulp )           DSTEMR('V', 'I')
*
*  (31)    ( max { min | WA2(i)-WA3(j) | } +
*             i     j
*            max { min | WA3(i)-WA2(j) | } ) / ( |D3| ulp )
*             i     j
*          DSTEMR('N', 'I') vs. SSTEMR('V', 'I')
*
*  (32)    | S - Z D Z' | / ( |S| n ulp )    DSTEMR('V', 'V')
*
*  (33)    | I - ZZ' | / ( n ulp )           DSTEMR('V', 'V')
*
*  (34)    ( max { min | WA2(i)-WA3(j) | } +
*             i     j
*            max { min | WA3(i)-WA2(j) | } ) / ( |D3| ulp )
*             i     j
*          DSTEMR('N', 'V') vs. SSTEMR('V', 'V')
*
*  (35)    | S - Z D Z' | / ( |S| n ulp )    DSTEMR('V', 'A')
*
*  (36)    | I - ZZ' | / ( n ulp )           DSTEMR('V', 'A')
*
*  (37)    ( max { min | WA2(i)-WA3(j) | } +
*             i     j
*            max { min | WA3(i)-WA2(j) | } ) / ( |D3| ulp )
*             i     j
*          DSTEMR('N', 'A') vs. SSTEMR('V', 'A')
*
*  The "sizes" are specified by an array NN(1:NSIZES); the value of
*  each element NN(j) specifies one size.
*  The "types" are specified by a logical array DOTYPE( 1:NTYPES );
*  if DOTYPE(j) is .TRUE., then matrix type "j" will be generated.
*  Currently, the list of possible types is:
*
*  (1)  The zero matrix.
*  (2)  The identity matrix.
*
*  (3)  A diagonal matrix with evenly spaced entries
*       1, ..., ULP  and random signs.
*       (ULP = (first number larger than 1) - 1 )
*  (4)  A diagonal matrix with geometrically spaced entries
*       1, ..., ULP  and random signs.
*  (5)  A diagonal matrix with "clustered" entries 1, ULP, ..., ULP
*       and random signs.
*
*  (6)  Same as (4), but multiplied by SQRT( overflow threshold )
*  (7)  Same as (4), but multiplied by SQRT( underflow threshold )
*
*  (8)  A matrix of the form  U' D U, where U is orthogonal and
*       D has evenly spaced entries 1, ..., ULP with random signs
*       on the diagonal.
*
*  (9)  A matrix of the form  U' D U, where U is orthogonal and
*       D has geometrically spaced entries 1, ..., ULP with random
*       signs on the diagonal.
*
*  (10) A matrix of the form  U' D U, where U is orthogonal and
*       D has "clustered" entries 1, ULP,..., ULP with random
*       signs on the diagonal.
*
*  (11) Same as (8), but multiplied by SQRT( overflow threshold )
*  (12) Same as (8), but multiplied by SQRT( underflow threshold )
*
*  (13) Symmetric matrix with random entries chosen from (-1,1).
*  (14) Same as (13), but multiplied by SQRT( overflow threshold )
*  (15) Same as (13), but multiplied by SQRT( underflow threshold )
*  (16) Same as (8), but diagonal elements are all positive.
*  (17) Same as (9), but diagonal elements are all positive.
*  (18) Same as (10), but diagonal elements are all positive.
*  (19) Same as (16), but multiplied by SQRT( overflow threshold )
*  (20) Same as (16), but multiplied by SQRT( underflow threshold )
*  (21) A diagonally dominant tridiagonal matrix with geometrically
*       spaced diagonal entries 1, ..., ULP.
*
*  Arguments
*  =========
*
*  NSIZES  (input) INTEGER
*          The number of sizes of matrices to use.  If it is zero,
*          DCHKST does nothing.  It must be at least zero.
*
*  NN      (input) INTEGER array, dimension (NSIZES)
*          An array containing the sizes to be used for the matrices.
*          Zero values will be skipped.  The values must be at least
*          zero.
*
*  NTYPES  (input) INTEGER
*          The number of elements in DOTYPE.   If it is zero, DCHKST
*          does nothing.  It must be at least zero.  If it is MAXTYP+1
*          and NSIZES is 1, then an additional type, MAXTYP+1 is
*          defined, which is to use whatever matrix is in A.  This
*          is only useful if DOTYPE(1:MAXTYP) is .FALSE. and
*          DOTYPE(MAXTYP+1) is .TRUE. .
*
*  DOTYPE  (input) LOGICAL array, dimension (NTYPES)
*          If DOTYPE(j) is .TRUE., then for each size in NN a
*          matrix of that size and of type j will be generated.
*          If NTYPES is smaller than the maximum number of types
*          defined (PARAMETER MAXTYP), then types NTYPES+1 through
*          MAXTYP will not be generated.  If NTYPES is larger
*          than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES)
*          will be ignored.
*
*  ISEED   (input/output) INTEGER array, dimension (4)
*          On entry ISEED specifies the seed of the random number
*          generator. The array elements should be between 0 and 4095;
*          if not they will be reduced mod 4096.  Also, ISEED(4) must
*          be odd.  The random number generator uses a linear
*          congruential sequence limited to small integers, and so
*          should produce machine independent random numbers. The
*          values of ISEED are changed on exit, and can be used in the
*          next call to DCHKST to continue the same random number
*          sequence.
*
*  THRESH  (input) DOUBLE PRECISION
*          A test will count as "failed" if the "error", computed as
*          described above, exceeds THRESH.  Note that the error
*          is scaled to be O(1), so THRESH should be a reasonably
*          small multiple of 1, e.g., 10 or 100.  In particular,
*          it should not depend on the precision (single vs. double)
*          or the size of the matrix.  It must be at least zero.
*
*  NOUNIT  (input) INTEGER
*          The FORTRAN unit number for printing out error messages
*          (e.g., if a routine returns IINFO not equal to 0.)
*
*  A       (input/workspace/output) DOUBLE PRECISION array of
*                                  dimension ( LDA , max(NN) )
*          Used to hold the matrix whose eigenvalues are to be
*          computed.  On exit, A contains the last matrix actually
*          used.
*
*  LDA     (input) INTEGER
*          The leading dimension of A.  It must be at
*          least 1 and at least max( NN ).
*
*  AP      (workspace) DOUBLE PRECISION array of
*                      dimension( max(NN)*max(NN+1)/2 )
*          The matrix A stored in packed format.
*
*  SD      (workspace/output) DOUBLE PRECISION array of
*                             dimension( max(NN) )
*          The diagonal of the tridiagonal matrix computed by DSYTRD.
*          On exit, SD and SE contain the tridiagonal form of the
*          matrix in A.
*
*  SE      (workspace/output) DOUBLE PRECISION array of
*                             dimension( max(NN) )
*          The off-diagonal of the tridiagonal matrix computed by
*          DSYTRD.  On exit, SD and SE contain the tridiagonal form of
*          the matrix in A.
*
*  D1      (workspace/output) DOUBLE PRECISION array of
*                             dimension( max(NN) )
*          The eigenvalues of A, as computed by DSTEQR simlutaneously
*          with Z.  On exit, the eigenvalues in D1 correspond with the
*          matrix in A.
*
*  D2      (workspace/output) DOUBLE PRECISION array of
*                             dimension( max(NN) )
*          The eigenvalues of A, as computed by DSTEQR if Z is not
*          computed.  On exit, the eigenvalues in D2 correspond with
*          the matrix in A.
*
*  D3      (workspace/output) DOUBLE PRECISION array of
*                             dimension( max(NN) )
*          The eigenvalues of A, as computed by DSTERF.  On exit, the
*          eigenvalues in D3 correspond with the matrix in A.
*
*  U       (workspace/output) DOUBLE PRECISION array of
*                             dimension( LDU, max(NN) ).
*          The orthogonal matrix computed by DSYTRD + DORGTR.
*
*  LDU     (input) INTEGER
*          The leading dimension of U, Z, and V.  It must be at least 1
*          and at least max( NN ).
*
*  V       (workspace/output) DOUBLE PRECISION array of
*                             dimension( LDU, max(NN) ).
*          The Housholder vectors computed by DSYTRD in reducing A to
*          tridiagonal form.  The vectors computed with UPLO='U' are
*          in the upper triangle, and the vectors computed with UPLO='L'
*          are in the lower triangle.  (As described in DSYTRD, the
*          sub- and superdiagonal are not set to 1, although the
*          true Householder vector has a 1 in that position.  The
*          routines that use V, such as DORGTR, set those entries to
*          1 before using them, and then restore them later.)
*
*  VP      (workspace) DOUBLE PRECISION array of
*                      dimension( max(NN)*max(NN+1)/2 )
*          The matrix V stored in packed format.
*
*  TAU     (workspace/output) DOUBLE PRECISION array of
*                             dimension( max(NN) )
*          The Householder factors computed by DSYTRD in reducing A
*          to tridiagonal form.
*
*  Z       (workspace/output) DOUBLE PRECISION array of
*                             dimension( LDU, max(NN) ).
*          The orthogonal matrix of eigenvectors computed by DSTEQR,
*          DPTEQR, and DSTEIN.
*
*  WORK    (workspace/output) DOUBLE PRECISION array of
*                      dimension( LWORK )
*
*  LWORK   (input) INTEGER
*          The number of entries in WORK.  This must be at least
*          1 + 4 * Nmax + 2 * Nmax * lg Nmax + 3 * Nmax**2
*          where Nmax = max( NN(j), 2 ) and lg = log base 2.
*
*  IWORK   (workspace/output) INTEGER array,
*             dimension (6 + 6*Nmax + 5 * Nmax * lg Nmax )
*          where Nmax = max( NN(j), 2 ) and lg = log base 2.
*          Workspace.
*
*  RESULT  (output) DOUBLE PRECISION array, dimension (26)
*          The values computed by the tests described above.
*          The values are currently limited to 1/ulp, to avoid
*          overflow.
*
*  INFO    (output) INTEGER
*          If 0, then everything ran OK.
*           -1: NSIZES < 0
*           -2: Some NN(j) < 0
*           -3: NTYPES < 0
*           -5: THRESH < 0
*           -9: LDA < 1 or LDA < NMAX, where NMAX is max( NN(j) ).
*          -23: LDU < 1 or LDU < NMAX.
*          -29: LWORK too small.
*          If  DLATMR, SLATMS, DSYTRD, DORGTR, DSTEQR, SSTERF,
*              or DORMC2 returns an error code, the
*              absolute value of it is returned.
*
*-----------------------------------------------------------------------
*
*       Some Local Variables and Parameters:
*       ---- ----- --------- --- ----------
*       ZERO, ONE       Real 0 and 1.
*       MAXTYP          The number of types defined.
*       NTEST           The number of tests performed, or which can
*                       be performed so far, for the current matrix.
*       NTESTT          The total number of tests performed so far.
*       NBLOCK          Blocksize as returned by ENVIR.
*       NMAX            Largest value in NN.
*       NMATS           The number of matrices generated so far.
*       NERRS           The number of tests which have exceeded THRESH
*                       so far.
*       COND, IMODE     Values to be passed to the matrix generators.
*       ANORM           Norm of A; passed to matrix generators.
*
*       OVFL, UNFL      Overflow and underflow thresholds.
*       ULP, ULPINV     Finest relative precision and its inverse.
*       RTOVFL, RTUNFL  Square roots of the previous 2 values.
*               The following four arrays decode JTYPE:
*       KTYPE(j)        The general type (1-10) for type "j".
*       KMODE(j)        The MODE value to be passed to the matrix
*                       generator for type "j".
*       KMAGN(j)        The order of magnitude ( O(1),
*                       O(overflow^(1/2) ), O(underflow^(1/2) )
*
*  =====================================================================
*
*     .. Parameters ..
      DOUBLE PRECISION   ZEROONETWOEIGHTTENHUN
      PARAMETER          ( ZERO = 0.0D0ONE = 1.0D0TWO = 2.0D0,
     $                   EIGHT = 8.0D0TEN = 10.0D0HUN = 100.0D0 )
      DOUBLE PRECISION   HALF
      PARAMETER          ( HALF = ONE / TWO )
      INTEGER            MAXTYP
      PARAMETER          ( MAXTYP = 21 )
      LOGICAL            SRANGE
      PARAMETER          ( SRANGE = .FALSE. )
      LOGICAL            SREL
      PARAMETER          ( SREL = .FALSE. )
*     ..
*     .. Local Scalars ..
      LOGICAL            BADNNTRYRAC
      INTEGER            IIINFOILIMODEITEMPITYPEIUJJC,
     $                   JRJSIZEJTYPELGNLIWEDCLOG2UILWEDC,
     $                   MM2M3MTYPESNNAPNBLOCKNERRS,
     $                   NMATSNMAXNSPLITNTESTNTESTT
      DOUBLE PRECISION   ABSTOLANINVANORMCONDOVFLRTOVFL,
     $                   RTUNFLTEMP1TEMP2TEMP3TEMP4ULP,
     $                   ULPINVUNFLVLVU
*     ..
*     .. Local Arrays ..
      INTEGER            IDUMMA1 ), IOLDSD4 ), ISEED24 ),
     $                   KMAGNMAXTYP ), KMODEMAXTYP ),
     $                   KTYPEMAXTYP )
      DOUBLE PRECISION   DUMMA1 )
*     ..
*     .. External Functions ..
      INTEGER            ILAENV
      DOUBLE PRECISION   DLAMCHDLARNDDSXT1
      EXTERNAL           ILAENVDLAMCHDLARNDDSXT1
*     ..
*     .. External Subroutines ..
      EXTERNAL           DCOPYDLABADDLACPYDLASETDLASUMDLATMR,
     $                   DLATMSDOPGTRDORGTRDPTEQRDSPT21DSPTRD,
     $                   DSTEBZDSTECHDSTEDCDSTEMRDSTEINDSTEQR,
     $                   DSTERFDSTT21DSTT22DSYT21DSYTRDXERBLA
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          ABSDBLEINTLOGMAXMINSQRT
*     ..
*     .. Data statements ..
      DATA               KTYPE / 1244444555558,
     $                   889999910 /
      DATA               KMAGN / 1111123111231,
     $                   23111231 /
      DATA               KMODE / 0043144431440,
     $                   00431443 /
*     ..
*     .. Executable Statements ..
*
*     Keep ftnchek happy
      IDUMMA1 ) = 1
*
*     Check for errors
*
      NTESTT = 0
      INFO = 0
*
*     Important constants
*
      BADNN = .FALSE.
      TRYRAC = .TRUE.
      NMAX = 1
      DO 10 J = 1NSIZES
         NMAX = MAXNMAXNNJ ) )
         IFNNJ ).LT.0 )
     $      BADNN = .TRUE.
   10 CONTINUE
*
      NBLOCK = ILAENV1'DSYTRD''L'NMAX-1-1-1 )
      NBLOCK = MINNMAXMAX1NBLOCK ) )
*
*     Check for errors
*
      IFNSIZES.LT.0 ) THEN
         INFO = -1
      ELSE IFBADNN ) THEN
         INFO = -2
      ELSE IFNTYPES.LT.0 ) THEN
         INFO = -3
      ELSE IFLDA.LT.NMAX ) THEN
         INFO = -9
      ELSE IFLDU.LT.NMAX ) THEN
         INFO = -23
      ELSE IF2*MAX2NMAX )**2.GT.LWORK ) THEN
         INFO = -29
      END IF
*
      IFINFO.NE.0 ) THEN
         CALL XERBLA'DCHKST'-INFO )
         RETURN
      END IF
*
*     Quick return if possible
*
      IFNSIZES.EQ.0 .OR. NTYPES.EQ.0 )
     $   RETURN
*
*     More Important constants
*
      UNFL = DLAMCH'Safe minimum' )
      OVFL = ONE / UNFL
      CALL DLABADUNFLOVFL )
      ULP = DLAMCH'Epsilon' )*DLAMCH'Base' )
      ULPINV = ONE / ULP
      LOG2UI = INTLOGULPINV ) / LOGTWO ) )
      RTUNFL = SQRTUNFL )
      RTOVFL = SQRTOVFL )
*
*     Loop over sizes, types
*
      DO 20 I = 14
         ISEED2I ) = ISEEDI )
   20 CONTINUE
      NERRS = 0
      NMATS = 0
*
      DO 310 JSIZE = 1NSIZES
         N = NNJSIZE )
         IFN.GT.0 ) THEN
            LGN = INTLOGDBLEN ) ) / LOGTWO ) )
            IF2**LGN.LT.N )
     $         LGN = LGN + 1
            IF2**LGN.LT.N )
     $         LGN = LGN + 1
            LWEDC = 1 + 4*N + 2*N*LGN + 3*N**2
            LIWEDC = 6 + 6*N + 5*N*LGN
         ELSE
            LWEDC = 8
            LIWEDC = 12
         END IF
         NAP = ( N*N+1 ) ) / 2
         ANINV = ONE / DBLEMAX1N ) )
*
         IFNSIZES.NE.1 ) THEN
            MTYPES = MINMAXTYPNTYPES )
         ELSE
            MTYPES = MINMAXTYP+1NTYPES )
         END IF
*
         DO 300 JTYPE = 1MTYPES
            IF.NOT.DOTYPEJTYPE ) )
     $         GO TO 300
            NMATS = NMATS + 1
            NTEST = 0
*
            DO 30 J = 14
               IOLDSDJ ) = ISEEDJ )
   30       CONTINUE
*
*           Compute "A"
*
*           Control parameters:
*
*               KMAGN  KMODE        KTYPE
*           =1  O(1)   clustered 1  zero
*           =2  large  clustered 2  identity
*           =3  small  exponential  (none)
*           =4         arithmetic   diagonal, (w/ eigenvalues)
*           =5         random log   symmetric, w/ eigenvalues
*           =6         random       (none)
*           =7                      random diagonal
*           =8                      random symmetric
*           =9                      positive definite
*           =10                     diagonally dominant tridiagonal
*
            IFMTYPES.GT.MAXTYP )
     $         GO TO 100
*
            ITYPE = KTYPEJTYPE )
            IMODE = KMODEJTYPE )
*
*           Compute norm
*
            GO TO ( 405060 )KMAGNJTYPE )
*
   40       CONTINUE
            ANORM = ONE
            GO TO 70
*
   50       CONTINUE
            ANORM = ( RTOVFL*ULP )*ANINV
            GO TO 70
*
   60       CONTINUE
            ANORM = RTUNFL*N*ULPINV
            GO TO 70
*
   70       CONTINUE
*
            CALL DLASET'Full'LDANZEROZEROALDA )
            IINFO = 0
            IFJTYPE.LE.15 ) THEN
               COND = ULPINV
            ELSE
               COND = ULPINV*ANINV / TEN
            END IF
*
*           Special Matrices -- Identity & Jordan block
*
*              Zero
*
            IFITYPE.EQ.1 ) THEN
               IINFO = 0
*
            ELSE IFITYPE.EQ.2 ) THEN
*
*              Identity
*
               DO 80 JC = 1N
                  AJCJC ) = ANORM
   80          CONTINUE
*
            ELSE IFITYPE.EQ.4 ) THEN
*
*              Diagonal Matrix, [Eigen]values Specified
*
               CALL DLATMSNN'S'ISEED'S'WORKIMODECOND,
     $                      ANORM00'N'ALDAWORKN+1 ),
     $                      IINFO )
*
*
            ELSE IFITYPE.EQ.5 ) THEN
*
*              Symmetric, eigenvalues specified
*
               CALL DLATMSNN'S'ISEED'S'WORKIMODECOND,
     $                      ANORMNN'N'ALDAWORKN+1 ),
     $                      IINFO )
*
            ELSE IFITYPE.EQ.7 ) THEN
*
*              Diagonal, random eigenvalues
*
               CALL DLATMRNN'S'ISEED'S'WORK6ONEONE,
     $                      'T''N'WORKN+1 ), 1ONE,
     $                      WORK2*N+1 ), 1ONE'N'IDUMMA00,
     $                      ZEROANORM'NO'ALDAIWORKIINFO )
*
            ELSE IFITYPE.EQ.8 ) THEN
*
*              Symmetric, random eigenvalues
*
               CALL DLATMRNN'S'ISEED'S'WORK6ONEONE,
     $                      'T''N'WORKN+1 ), 1ONE,
     $                      WORK2*N+1 ), 1ONE'N'IDUMMANN,
     $                      ZEROANORM'NO'ALDAIWORKIINFO )
*
            ELSE IFITYPE.EQ.9 ) THEN
*
*              Positive definite, eigenvalues specified.
*
               CALL DLATMSNN'S'ISEED'P'WORKIMODECOND,
     $                      ANORMNN'N'ALDAWORKN+1 ),
     $                      IINFO )
*
            ELSE IFITYPE.EQ.10 ) THEN
*
*              Positive definite tridiagonal, eigenvalues specified.
*
               CALL DLATMSNN'S'ISEED'P'WORKIMODECOND,
     $                      ANORM11'N'ALDAWORKN+1 ),
     $                      IINFO )
               DO 90 I = 2N
                  TEMP1 = ABSAI-1I ) ) /
     $                    SQRTABSAI-1I-1 )*AII ) ) )
                  IFTEMP1.GT.HALF ) THEN
                     AI-1I ) = HALF*SQRTABSAI-1I-1 )*AI,
     $                             I ) ) )
                     AII-1 ) = AI-1I )
                  END IF
   90          CONTINUE
*
            ELSE
*
               IINFO = 1
            END IF
*
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'Generator', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               RETURN
            END IF
*
  100       CONTINUE
*
*           Call DSYTRD and DORGTR to compute S and U from
*           upper triangle.
*
            CALL DLACPY'U'NNALDAVLDU )
*
            NTEST = 1
            CALL DSYTRD'U'NVLDUSDSETAUWORKLWORK,
     $                   IINFO )
*
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSYTRD(U)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT1 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
            CALL DLACPY'U'NNVLDUULDU )
*
            NTEST = 2
            CALL DORGTR'U'NULDUTAUWORKLWORKIINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DORGTR(U)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT2 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Do tests 1 and 2
*
            CALL DSYT212'Upper'N1ALDASDSEULDUV,
     $                   LDUTAUWORKRESULT1 ) )
            CALL DSYT213'Upper'N1ALDASDSEULDUV,
     $                   LDUTAUWORKRESULT2 ) )
*
*           Call DSYTRD and DORGTR to compute S and U from
*           lower triangle, do tests.
*
            CALL DLACPY'L'NNALDAVLDU )
*
            NTEST = 3
            CALL DSYTRD'L'NVLDUSDSETAUWORKLWORK,
     $                   IINFO )
*
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSYTRD(L)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT3 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
            CALL DLACPY'L'NNVLDUULDU )
*
            NTEST = 4
            CALL DORGTR'L'NULDUTAUWORKLWORKIINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DORGTR(L)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT4 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
            CALL DSYT212'Lower'N1ALDASDSEULDUV,
     $                   LDUTAUWORKRESULT3 ) )
            CALL DSYT213'Lower'N1ALDASDSEULDUV,
     $                   LDUTAUWORKRESULT4 ) )
*
*           Store the upper triangle of A in AP
*
            I = 0
            DO 120 JC = 1N
               DO 110 JR = 1JC
                  I = I + 1
                  API ) = AJRJC )
  110          CONTINUE
  120       CONTINUE
*
*           Call DSPTRD and DOPGTR to compute S and U from AP
*
            CALL DCOPYNAPAP1VP1 )
*
            NTEST = 5
            CALL DSPTRD'U'NVPSDSETAUIINFO )
*
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSPTRD(U)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT5 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
            NTEST = 6
            CALL DOPGTR'U'NVPTAUULDUWORKIINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DOPGTR(U)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT6 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Do tests 5 and 6
*
            CALL DSPT212'Upper'N1APSDSEULDUVPTAU,
     $                   WORKRESULT5 ) )
            CALL DSPT213'Upper'N1APSDSEULDUVPTAU,
     $                   WORKRESULT6 ) )
*
*           Store the lower triangle of A in AP
*
            I = 0
            DO 140 JC = 1N
               DO 130 JR = JCN
                  I = I + 1
                  API ) = AJRJC )
  130          CONTINUE
  140       CONTINUE
*
*           Call DSPTRD and DOPGTR to compute S and U from AP
*
            CALL DCOPYNAPAP1VP1 )
*
            NTEST = 7
            CALL DSPTRD'L'NVPSDSETAUIINFO )
*
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSPTRD(L)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT7 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
            NTEST = 8
            CALL DOPGTR'L'NVPTAUULDUWORKIINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DOPGTR(L)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT8 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
            CALL DSPT212'Lower'N1APSDSEULDUVPTAU,
     $                   WORKRESULT7 ) )
            CALL DSPT213'Lower'N1APSDSEULDUVPTAU,
     $                   WORKRESULT8 ) )
*
*           Call DSTEQR to compute D1, D2, and Z, do tests.
*
*           Compute D1 and Z
*
            CALL DCOPYNSD1D11 )
            IFN.GT.0 )
     $         CALL DCOPYN-1SE1WORK1 )
            CALL DLASET'Full'NNZEROONEZLDU )
*
            NTEST = 9
            CALL DSTEQR'V'ND1WORKZLDUWORKN+1 ), IINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTEQR(V)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT9 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Compute D2
*
            CALL DCOPYNSD1D21 )
            IFN.GT.0 )
     $         CALL DCOPYN-1SE1WORK1 )
*
            NTEST = 11
            CALL DSTEQR'N'ND2WORKWORKN+1 ), LDU,
     $                   WORKN+1 ), IINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTEQR(N)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT11 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Compute D3 (using PWK method)
*
            CALL DCOPYNSD1D31 )
            IFN.GT.0 )
     $         CALL DCOPYN-1SE1WORK1 )
*
            NTEST = 12
            CALL DSTERFND3WORKIINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTERF', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT12 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Do Tests 9 and 10
*
            CALL DSTT21N0SDSED1DUMMAZLDUWORK,
     $                   RESULT9 ) )
*
*           Do Tests 11 and 12
*
            TEMP1 = ZERO
            TEMP2 = ZERO
            TEMP3 = ZERO
            TEMP4 = ZERO
*
            DO 150 J = 1N
               TEMP1 = MAXTEMP1ABSD1J ) ), ABSD2J ) ) )
               TEMP2 = MAXTEMP2ABSD1J )-D2J ) ) )
               TEMP3 = MAXTEMP3ABSD1J ) ), ABSD3J ) ) )
               TEMP4 = MAXTEMP4ABSD1J )-D3J ) ) )
  150       CONTINUE
*
            RESULT11 ) = TEMP2 / MAXUNFLULP*MAXTEMP1TEMP2 ) )
            RESULT12 ) = TEMP4 / MAXUNFLULP*MAXTEMP3TEMP4 ) )
*
*           Do Test 13 -- Sturm Sequence Test of Eigenvalues
*                         Go up by factors of two until it succeeds
*
            NTEST = 13
            TEMP1 = THRESH*HALF-ULP )
*
            DO 160 J = 0LOG2UI
               CALL DSTECHNSDSED1TEMP1WORKIINFO )
               IFIINFO.EQ.0 )
     $            GO TO 170
               TEMP1 = TEMP1*TWO
  160       CONTINUE
*
  170       CONTINUE
            RESULT13 ) = TEMP1
*
*           For positive definite matrices ( JTYPE.GT.15 ) call DPTEQR
*           and do tests 14, 15, and 16 .
*
            IFJTYPE.GT.15 ) THEN
*
*              Compute D4 and Z4
*
               CALL DCOPYNSD1D41 )
               IFN.GT.0 )
     $            CALL DCOPYN-1SE1WORK1 )
               CALL DLASET'Full'NNZEROONEZLDU )
*
               NTEST = 14
               CALL DPTEQR'V'ND4WORKZLDUWORKN+1 ),
     $                      IINFO )
               IFIINFO.NE.0 ) THEN
                  WRITENOUNIT, FMT = 9999 )'DPTEQR(V)', IINFO, N,
     $               JTYPE, IOLDSD
                  INFO = ABSIINFO )
                  IFIINFO.LT.0 ) THEN
                     RETURN
                  ELSE
                     RESULT14 ) = ULPINV
                     GO TO 280
                  END IF
               END IF
*
*              Do Tests 14 and 15
*
               CALL DSTT21N0SDSED4DUMMAZLDUWORK,
     $                      RESULT14 ) )
*
*              Compute D5
*
               CALL DCOPYNSD1D51 )
               IFN.GT.0 )
     $            CALL DCOPYN-1SE1WORK1 )
*
               NTEST = 16
               CALL DPTEQR'N'ND5WORKZLDUWORKN+1 ),
     $                      IINFO )
               IFIINFO.NE.0 ) THEN
                  WRITENOUNIT, FMT = 9999 )'DPTEQR(N)', IINFO, N,
     $               JTYPE, IOLDSD
                  INFO = ABSIINFO )
                  IFIINFO.LT.0 ) THEN
                     RETURN
                  ELSE
                     RESULT16 ) = ULPINV
                     GO TO 280
                  END IF
               END IF
*
*              Do Test 16
*
               TEMP1 = ZERO
               TEMP2 = ZERO
               DO 180 J = 1N
                  TEMP1 = MAXTEMP1ABSD4J ) ), ABSD5J ) ) )
                  TEMP2 = MAXTEMP2ABSD4J )-D5J ) ) )
  180          CONTINUE
*
               RESULT16 ) = TEMP2 / MAXUNFL,
     $                        HUN*ULP*MAXTEMP1TEMP2 ) )
            ELSE
               RESULT14 ) = ZERO
               RESULT15 ) = ZERO
               RESULT16 ) = ZERO
            END IF
*
*           Call DSTEBZ with different options and do tests 17-18.
*
*              If S is positive definite and diagonally dominant,
*              ask for all eigenvalues with high relative accuracy.
*
            VL = ZERO
            VU = ZERO
            IL = 0
            IU = 0
            IFJTYPE.EQ.21 ) THEN
               NTEST = 17
               ABSTOL = UNFL + UNFL
               CALL DSTEBZ'A''E'NVLVUILIUABSTOLSDSE,
     $                      MNSPLITWRIWORK1 ), IWORKN+1 ),
     $                      WORKIWORK2*N+1 ), IINFO )
               IFIINFO.NE.0 ) THEN
                  WRITENOUNIT, FMT = 9999 )'DSTEBZ(A,rel)', IINFO, N,
     $               JTYPE, IOLDSD
                  INFO = ABSIINFO )
                  IFIINFO.LT.0 ) THEN
                     RETURN
                  ELSE
                     RESULT17 ) = ULPINV
                     GO TO 280
                  END IF
               END IF
*
*              Do test 17
*
               TEMP2 = TWO*TWO*N-ONE )*ULP*ONE+EIGHT*HALF**2 ) /
     $                 ( ONE-HALF )**4
*
               TEMP1 = ZERO
               DO 190 J = 1N
                  TEMP1 = MAXTEMP1ABSD4J )-WRN-J+1 ) ) /
     $                    ( ABSTOL+ABSD4J ) ) ) )
  190          CONTINUE
*
               RESULT17 ) = TEMP1 / TEMP2
            ELSE
               RESULT17 ) = ZERO
            END IF
*
*           Now ask for all eigenvalues with high absolute accuracy.
*
            NTEST = 18
            ABSTOL = UNFL + UNFL
            CALL DSTEBZ'A''E'NVLVUILIUABSTOLSDSEM,
     $                   NSPLITWA1IWORK1 ), IWORKN+1 ), WORK,
     $                   IWORK2*N+1 ), IINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTEBZ(A)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT18 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Do test 18
*
            TEMP1 = ZERO
            TEMP2 = ZERO
            DO 200 J = 1N
               TEMP1 = MAXTEMP1ABSD3J ) ), ABSWA1J ) ) )
               TEMP2 = MAXTEMP2ABSD3J )-WA1J ) ) )
  200       CONTINUE
*
            RESULT18 ) = TEMP2 / MAXUNFLULP*MAXTEMP1TEMP2 ) )
*
*           Choose random values for IL and IU, and ask for the
*           IL-th through IU-th eigenvalues.
*
            NTEST = 19
            IFN.LE.1 ) THEN
               IL = 1
               IU = N
            ELSE
               IL = 1 + ( N-1 )*INTDLARND1ISEED2 ) )
               IU = 1 + ( N-1 )*INTDLARND1ISEED2 ) )
               IFIU.LT.IL ) THEN
                  ITEMP = IU
                  IU = IL
                  IL = ITEMP
               END IF
            END IF
*
            CALL DSTEBZ'I''E'NVLVUILIUABSTOLSDSE,
     $                   M2NSPLITWA2IWORK1 ), IWORKN+1 ),
     $                   WORKIWORK2*N+1 ), IINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTEBZ(I)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT19 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Determine the values VL and VU of the IL-th and IU-th
*           eigenvalues and ask for all eigenvalues in this range.
*
            IFN.GT.0 ) THEN
               IFIL.NE.1 ) THEN
                  VL = WA1IL ) - MAXHALF*WA1IL )-WA1IL-1 ) ),
     $                 ULP*ANORMTWO*RTUNFL )
               ELSE
                  VL = WA11 ) - MAXHALF*WA1N )-WA11 ) ),
     $                 ULP*ANORMTWO*RTUNFL )
               END IF
               IFIU.NE.N ) THEN
                  VU = WA1IU ) + MAXHALF*WA1IU+1 )-WA1IU ) ),
     $                 ULP*ANORMTWO*RTUNFL )
               ELSE
                  VU = WA1N ) + MAXHALF*WA1N )-WA11 ) ),
     $                 ULP*ANORMTWO*RTUNFL )
               END IF
            ELSE
               VL = ZERO
               VU = ONE
            END IF
*
            CALL DSTEBZ'V''E'NVLVUILIUABSTOLSDSE,
     $                   M3NSPLITWA3IWORK1 ), IWORKN+1 ),
     $                   WORKIWORK2*N+1 ), IINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTEBZ(V)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT19 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
            IFM3.EQ.0 .AND. N.NE.0 ) THEN
               RESULT19 ) = ULPINV
               GO TO 280
            END IF
*
*           Do test 19
*
            TEMP1 = DSXT11WA2M2WA3M3ABSTOLULPUNFL )
            TEMP2 = DSXT11WA3M3WA2M2ABSTOLULPUNFL )
            IFN.GT.0 ) THEN
               TEMP3 = MAXABSWA1N ) ), ABSWA11 ) ) )
            ELSE
               TEMP3 = ZERO
            END IF
*
            RESULT19 ) = ( TEMP1+TEMP2 ) / MAXUNFLTEMP3*ULP )
*
*           Call DSTEIN to compute eigenvectors corresponding to
*           eigenvalues in WA1.  (First call DSTEBZ again, to make sure
*           it returns these eigenvalues in the correct order.)
*
            NTEST = 21
            CALL DSTEBZ'A''B'NVLVUILIUABSTOLSDSEM,
     $                   NSPLITWA1IWORK1 ), IWORKN+1 ), WORK,
     $                   IWORK2*N+1 ), IINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTEBZ(A,B)', IINFO, N,
     $            JTYPE, IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT20 ) = ULPINV
                  RESULT21 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
            CALL DSTEINNSDSEMWA1IWORK1 ), IWORKN+1 ), Z,
     $                   LDUWORKIWORK2*N+1 ), IWORK3*N+1 ),
     $                   IINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTEIN', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT20 ) = ULPINV
                  RESULT21 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Do tests 20 and 21
*
            CALL DSTT21N0SDSEWA1DUMMAZLDUWORK,
     $                   RESULT20 ) )
*
*           Call DSTEDC(I) to compute D1 and Z, do tests.
*
*           Compute D1 and Z
*
            CALL DCOPYNSD1D11 )
            IFN.GT.0 )
     $         CALL DCOPYN-1SE1WORK1 )
            CALL DLASET'Full'NNZEROONEZLDU )
*
            NTEST = 22
            CALL DSTEDC'I'ND1WORKZLDUWORKN+1 ), LWEDC-N,
     $                   IWORKLIWEDCIINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTEDC(I)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT22 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Do Tests 22 and 23
*
            CALL DSTT21N0SDSED1DUMMAZLDUWORK,
     $                   RESULT22 ) )
*
*           Call DSTEDC(V) to compute D1 and Z, do tests.
*
*           Compute D1 and Z
*
            CALL DCOPYNSD1D11 )
            IFN.GT.0 )
     $         CALL DCOPYN-1SE1WORK1 )
            CALL DLASET'Full'NNZEROONEZLDU )
*
            NTEST = 24
            CALL DSTEDC'V'ND1WORKZLDUWORKN+1 ), LWEDC-N,
     $                   IWORKLIWEDCIINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTEDC(V)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT24 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Do Tests 24 and 25
*
            CALL DSTT21N0SDSED1DUMMAZLDUWORK,
     $                   RESULT24 ) )
*
*           Call DSTEDC(N) to compute D2, do tests.
*
*           Compute D2
*
            CALL DCOPYNSD1D21 )
            IFN.GT.0 )
     $         CALL DCOPYN-1SE1WORK1 )
            CALL DLASET'Full'NNZEROONEZLDU )
*
            NTEST = 26
            CALL DSTEDC'N'ND2WORKZLDUWORKN+1 ), LWEDC-N,
     $                   IWORKLIWEDCIINFO )
            IFIINFO.NE.0 ) THEN
               WRITENOUNIT, FMT = 9999 )'DSTEDC(N)', IINFO, N, JTYPE,
     $            IOLDSD
               INFO = ABSIINFO )
               IFIINFO.LT.0 ) THEN
                  RETURN
               ELSE
                  RESULT26 ) = ULPINV
                  GO TO 280
               END IF
            END IF
*
*           Do Test 26
*
            TEMP1 = ZERO
            TEMP2 = ZERO
*
            DO 210 J = 1N
               TEMP1 = MAXTEMP1ABSD1J ) ), ABSD2J ) ) )
               TEMP2 = MAXTEMP2ABSD1J )-D2J ) ) )
  210       CONTINUE
*
            RESULT26 ) = TEMP2 / MAXUNFLULP*MAXTEMP1TEMP2 ) )
*
*           Only test DSTEMR if IEEE compliant
*
            IFILAENV10'DSTEMR''VA'1000 ).EQ.1 .AND.
     $          ILAENV11'DSTEMR''VA'1000 ).EQ.1 ) THEN
*
*           Call DSTEMR, do test 27 (relative eigenvalue accuracy)
*
*              If S is positive definite and diagonally dominant,
*              ask for all eigenvalues with high relative accuracy.
*
               VL = ZERO
               VU = ZERO
               IL = 0
               IU = 0
               IFJTYPE.EQ.21 .AND. SREL ) THEN
                  NTEST = 27
                  ABSTOL = UNFL + UNFL
                  CALL DSTEMR'V''A'NSDSEVLVUILIU,
     $                         MWRZLDUNIWORK1 ), TRYRAC,
     $                         WORKLWORKIWORK2*N+1 ), LWORK-2*N,
     $                         IINFO )
                  IFIINFO.NE.0 ) THEN
                     WRITENOUNIT, FMT = 9999 )'DSTEMR(V,A,rel)',
     $                  IINFO, N, JTYPE, IOLDSD
                     INFO = ABSIINFO )
                     IFIINFO.LT.0 ) THEN
                        RETURN
                     ELSE
                        RESULT27 ) = ULPINV
                        GO TO 270
                     END IF
                  END IF
*
*              Do test 27
*
                  TEMP2 = TWO*TWO*N-ONE )*ULP*ONE+EIGHT*HALF**2 ) /
     $                    ( ONE-HALF )**4
*
                  TEMP1 = ZERO
                  DO 220 J = 1N
                     TEMP1 = MAXTEMP1ABSD4J )-WRN-J+1 ) ) /
     $                       ( ABSTOL+ABSD4J ) ) ) )
  220             CONTINUE
*
                  RESULT27 ) = TEMP1 / TEMP2
*
                  IL = 1 + ( N-1 )*INTDLARND1ISEED2 ) )
                  IU = 1 + ( N-1 )*INTDLARND1ISEED2 ) )
                  IFIU.LT.IL ) THEN
                     ITEMP = IU
                     IU = IL
                     IL = ITEMP
                  END IF
*
                  IFSRANGE ) THEN
                     NTEST = 28
                     ABSTOL = UNFL + UNFL
                     CALL DSTEMR'V''I'NSDSEVLVUILIU,
     $                            MWRZLDUNIWORK1 ), TRYRAC,
     $                            WORKLWORKIWORK2*N+1 ),
     $                            LWORK-2*NIINFO )
*
                     IFIINFO.NE.0 ) THEN
                        WRITENOUNIT, FMT = 9999 )'DSTEMR(V,I,rel)',
     $                     IINFO, N, JTYPE, IOLDSD
                        INFO = ABSIINFO )
                        IFIINFO.LT.0 ) THEN
                           RETURN
                        ELSE
                           RESULT28 ) = ULPINV
                           GO TO 270
                        END IF
                     END IF
*
*
*                 Do test 28
*
                     TEMP2 = TWO*TWO*N-ONE )*ULP*
     $                       ( ONE+EIGHT*HALF**2 ) / ( ONE-HALF )**4
*
                     TEMP1 = ZERO
                     DO 230 J = ILIU
                        TEMP1 = MAXTEMP1ABSWRJ-IL+1 )-D4N-J+
     $                          1 ) ) / ( ABSTOL+ABSWRJ-IL+1 ) ) ) )
  230                CONTINUE
*
                     RESULT28 ) = TEMP1 / TEMP2
                  ELSE
                     RESULT28 ) = ZERO
                  END IF
               ELSE
                  RESULT27 ) = ZERO
                  RESULT28 ) = ZERO
               END IF
*
*           Call DSTEMR(V,I) to compute D1 and Z, do tests.
*
*           Compute D1 and Z
*
               CALL DCOPYNSD1D51 )
               IFN.GT.0 )
     $            CALL DCOPYN-1SE1WORK1 )
               CALL DLASET'Full'NNZEROONEZLDU )
*
               IFSRANGE ) THEN
                  NTEST = 29
                  IL = 1 + ( N-1 )*INTDLARND1ISEED2 ) )
                  IU = 1 + ( N-1 )*INTDLARND1ISEED2 ) )
                  IFIU.LT.IL ) THEN
                     ITEMP = IU
                     IU = IL
                     IL = ITEMP
                  END IF
                  CALL DSTEMR'V''I'ND5WORKVLVUILIU,
     $                         MD1ZLDUNIWORK1 ), TRYRAC,
     $                         WORKN+1 ), LWORK-NIWORK2*N+1 ),
     $                         LIWORK-2*NIINFO )
                  IFIINFO.NE.0 ) THEN
                     WRITENOUNIT, FMT = 9999 )'DSTEMR(V,I)', IINFO,
     $                  N, JTYPE, IOLDSD
                     INFO = ABSIINFO )
                     IFIINFO.LT.0 ) THEN
                        RETURN
                     ELSE
                        RESULT29 ) = ULPINV
                        GO TO 280
                     END IF
                  END IF
*
*           Do Tests 29 and 30
*
                  CALL DSTT22NM0SDSED1DUMMAZLDUWORK,
     $                         MRESULT29 ) )
*
*           Call DSTEMR to compute D2, do tests.
*
*           Compute D2
*
                  CALL DCOPYNSD1D51 )
                  IFN.GT.0 )
     $               CALL DCOPYN-1SE1WORK1 )
*
                  NTEST = 31
                  CALL DSTEMR'N''I'ND5WORKVLVUILIU,
     $                         MD2ZLDUNIWORK1 ), TRYRAC,
     $                         WORKN+1 ), LWORK-NIWORK2*N+1 ),
     $                         LIWORK-2*NIINFO )
                  IFIINFO.NE.0 ) THEN
                     WRITENOUNIT, FMT = 9999 )'DSTEMR(N,I)', IINFO,
     $                  N, JTYPE, IOLDSD
                     INFO = ABSIINFO )
                     IFIINFO.LT.0 ) THEN
                        RETURN
                     ELSE
                        RESULT31 ) = ULPINV
                        GO TO 280
                     END IF
                  END IF
*
*           Do Test 31
*
                  TEMP1 = ZERO
                  TEMP2 = ZERO
*
                  DO 240 J = 1IU - IL + 1
                     TEMP1 = MAXTEMP1ABSD1J ) ),
     $                       ABSD2J ) ) )
                     TEMP2 = MAXTEMP2ABSD1J )-D2J ) ) )
  240             CONTINUE
*
                  RESULT31 ) = TEMP2 / MAXUNFL,
     $                           ULP*MAXTEMP1TEMP2 ) )
*
*
*           Call DSTEMR(V,V) to compute D1 and Z, do tests.
*
*           Compute D1 and Z
*
                  CALL DCOPYNSD1D51 )
                  IFN.GT.0 )
     $               CALL DCOPYN-1SE1WORK1 )
                  CALL DLASET'Full'NNZEROONEZLDU )
*
                  NTEST = 32
*
                  IFN.GT.0 ) THEN
                     IFIL.NE.1 ) THEN
                        VL = D2IL ) - MAXHALF*
     $                       ( D2IL )-D2IL-1 ) ), ULP*ANORM,
     $                       TWO*RTUNFL )
                     ELSE
                        VL = D21 ) - MAXHALF*D2N )-D21 ) ),
     $                       ULP*ANORMTWO*RTUNFL )
                     END IF
                     IFIU.NE.N ) THEN
                        VU = D2IU ) + MAXHALF*
     $                       ( D2IU+1 )-D2IU ) ), ULP*ANORM,
     $                       TWO*RTUNFL )
                     ELSE
                        VU = D2N ) + MAXHALF*D2N )-D21 ) ),
     $                       ULP*ANORMTWO*RTUNFL )
                     END IF
                  ELSE
                     VL = ZERO
                     VU = ONE
                  END IF
*
                  CALL DSTEMR'V''V'ND5WORKVLVUILIU,
     $                         MD1ZLDUNIWORK1 ), TRYRAC,
     $                         WORKN+1 ), LWORK-NIWORK2*N+1 ),
     $                         LIWORK-2*NIINFO )
                  IFIINFO.NE.0 ) THEN
                     WRITENOUNIT, FMT = 9999 )'DSTEMR(V,V)', IINFO,
     $                  N, JTYPE, IOLDSD
                     INFO = ABSIINFO )
                     IFIINFO.LT.0 ) THEN
                        RETURN
                     ELSE
                        RESULT32 ) = ULPINV
                        GO TO 280
                     END IF
                  END IF
*
*           Do Tests 32 and 33
*
                  CALL DSTT22NM0SDSED1DUMMAZLDUWORK,
     $                         MRESULT32 ) )
*
*           Call DSTEMR to compute D2, do tests.
*
*           Compute D2
*
                  CALL DCOPYNSD1D51 )
                  IFN.GT.0 )
     $               CALL DCOPYN-1SE1WORK1 )
*
                  NTEST = 34
                  CALL DSTEMR'N''V'ND5WORKVLVUILIU,
     $                         MD2ZLDUNIWORK1 ), TRYRAC,
     $                         WORKN+1 ), LWORK-NIWORK2*N+1 ),
     $                         LIWORK-2*NIINFO )
                  IFIINFO.NE.0 ) THEN
                     WRITENOUNIT, FMT = 9999 )'DSTEMR(N,V)', IINFO,
     $                  N, JTYPE, IOLDSD
                     INFO = ABSIINFO )
                     IFIINFO.LT.0 ) THEN
                        RETURN
                     ELSE
                        RESULT34 ) = ULPINV
                        GO TO 280
                     END IF
                  END IF
*
*           Do Test 34
*
                  TEMP1 = ZERO
                  TEMP2 = ZERO
*
                  DO 250 J = 1IU - IL + 1
                     TEMP1 = MAXTEMP1ABSD1J ) ),
     $                       ABSD2J ) ) )
                     TEMP2 = MAXTEMP2ABSD1J )-D2J ) ) )
  250             CONTINUE
*
                  RESULT34 ) = TEMP2 / MAXUNFL,
     $                           ULP*MAXTEMP1TEMP2 ) )
               ELSE
                  RESULT29 ) = ZERO
                  RESULT30 ) = ZERO
                  RESULT31 ) = ZERO
                  RESULT32 ) = ZERO
                  RESULT33 ) = ZERO
                  RESULT34 ) = ZERO
               END IF
*
*
*           Call DSTEMR(V,A) to compute D1 and Z, do tests.
*
*           Compute D1 and Z
*
               CALL DCOPYNSD1D51 )
               IFN.GT.0 )
     $            CALL DCOPYN-1SE1WORK1 )
*
               NTEST = 35
*
               CALL DSTEMR'V''A'ND5WORKVLVUILIU,
     $                      MD1ZLDUNIWORK1 ), TRYRAC,
     $                      WORKN+1 ), LWORK-NIWORK2*N+1 ),
     $                      LIWORK-2*NIINFO )
               IFIINFO.NE.0 ) THEN
                  WRITENOUNIT, FMT = 9999 )'DSTEMR(V,A)', IINFO, N,
     $               JTYPE, IOLDSD
                  INFO = ABSIINFO )
                  IFIINFO.LT.0 ) THEN
                     RETURN
                  ELSE
                     RESULT35 ) = ULPINV
                     GO TO 280
                  END IF
               END IF
*
*           Do Tests 35 and 36
*
               CALL DSTT22NM0SDSED1DUMMAZLDUWORKM,
     $                      RESULT35 ) )
*
*           Call DSTEMR to compute D2, do tests.
*
*           Compute D2
*
               CALL DCOPYNSD1D51 )
               IFN.GT.0 )
     $            CALL DCOPYN-1SE1WORK1 )
*
               NTEST = 37
               CALL DSTEMR'N''A'ND5WORKVLVUILIU,
     $                      MD2ZLDUNIWORK1 ), TRYRAC,
     $                      WORKN+1 ), LWORK-NIWORK2*N+1 ),
     $                      LIWORK-2*NIINFO )
               IFIINFO.NE.0 ) THEN
                  WRITENOUNIT, FMT = 9999 )'DSTEMR(N,A)', IINFO, N,
     $               JTYPE, IOLDSD
                  INFO = ABSIINFO )
                  IFIINFO.LT.0 ) THEN
                     RETURN
                  ELSE
                     RESULT37 ) = ULPINV
                     GO TO 280
                  END IF
               END IF
*
*           Do Test 34
*
               TEMP1 = ZERO
               TEMP2 = ZERO
*
               DO 260 J = 1N
                  TEMP1 = MAXTEMP1ABSD1J ) ), ABSD2J ) ) )
                  TEMP2 = MAXTEMP2ABSD1J )-D2J ) ) )
  260          CONTINUE
*
               RESULT37 ) = TEMP2 / MAXUNFL,
     $                        ULP*MAXTEMP1TEMP2 ) )
            END IF
  270       CONTINUE
  280       CONTINUE
            NTESTT = NTESTT + NTEST
*
*           End of Loop -- Check for RESULT(j) > THRESH
*
*
*           Print out tests which fail.
*
            DO 290 JR = 1NTEST
               IFRESULTJR ).GE.THRESH ) THEN
*
*                 If this is the first test to fail,
*                 print a header to the data file.
*
                  IFNERRS.EQ.0 ) THEN
                     WRITENOUNIT, FMT = 9998 )'DST'
                     WRITENOUNIT, FMT = 9997 )
                     WRITENOUNIT, FMT = 9996 )
                     WRITENOUNIT, FMT = 9995 )'Symmetric'
                     WRITENOUNIT, FMT = 9994 )
*
*                    Tests performed
*
                     WRITENOUNIT, FMT = 9988 )
                  END IF
                  NERRS = NERRS + 1
                  WRITENOUNIT, FMT = 9990 )N, IOLDSD, JTYPE, JR,
     $               RESULT( JR )
               END IF
  290       CONTINUE
  300    CONTINUE
  310 CONTINUE
*
*     Summary
*
      CALL DLASUM'DST'NOUNITNERRSNTESTT )
      RETURN
*
 9999 FORMAT( ' DCHKST: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
     $      I6, ', JTYPE=', I6, ', ISEED=(', 3( I5, ',' ), I5, ')' )
*
 9998 FORMAT( / 1X, A3, ' -- Real Symmetric eigenvalue problem' )
 9997 FORMAT( ' Matrix types (see DCHKST for details): ' )
*
 9996 FORMAT( / ' Special Matrices:',
     $      / '  1=Zero matrix.                        ',
     $      '  5=Diagonal: clustered entries.',
     $      / '  2=Identity matrix.                    ',
     $      '  6=Diagonal: large, evenly spaced.',
     $      / '  3=Diagonal: evenly spaced entries.    ',
     $      '  7=Diagonal: small, evenly spaced.',
     $      / '  4=Diagonal: geometr. spaced entries.' )
 9995 FORMAT( ' Dense ', A, ' Matrices:',
     $      / '  8=Evenly spaced eigenvals.            ',
     $      ' 12=Small, evenly spaced eigenvals.',
     $      / '  9=Geometrically spaced eigenvals.     ',
     $      ' 13=Matrix with random O(1) entries.',
     $      / ' 10=Clustered eigenvalues.              ',
     $      ' 14=Matrix with large random entries.',
     $      / ' 11=Large, evenly spaced eigenvals.     ',
     $      ' 15=Matrix with small random entries.' )
 9994 FORMAT( ' 16=Positive definite, evenly spaced eigenvalues',
     $      / ' 17=Positive definite, geometrically spaced eigenvlaues',
     $      / ' 18=Positive definite, clustered eigenvalues',
     $      / ' 19=Positive definite, small evenly spaced eigenvalues',
     $      / ' 20=Positive definite, large evenly spaced eigenvalues',
     $      / ' 21=Diagonally dominant tridiagonal, geometrically',
     $      ' spaced eigenvalues' )
*
 9993 FORMAT( / ' Tests performed:   ',
     $      '(S is Tridiag, D is diagonal, U and Z are ', A, ',', / 20X,
     $      A, ', W is a diagonal matrix of eigenvalues,', / 20X,
     $      ' V is U represented by Householder vectors, and', / 20X,
     $      ' Y is a matrix of eigenvectors of S.)',
     $      / ' DSYTRD, UPLO=''U'':', / '  1= | A - V S V', A1,
     $      ' | / ( |A| n ulp )     ', '  2= | I - U V', A1,
     $      ' | / ( n ulp )', / ' DSYTRD, UPLO=''L'':',
     $      / '  3= | A - V S V', A1, ' | / ( |A| n ulp )     ',
     $      '  4= | I - U V', A1, ' | / ( n ulp )' )
 9992 FORMAT( ' DSPTRD, UPLO=''U'':', / '  5= | A - V S V', A1,
     $      ' | / ( |A| n ulp )     ', '  6= | I - U V', A1,
     $      ' | / ( n ulp )', / ' DSPTRD, UPLO=''L'':',
     $      / '  7= | A - V S V', A1, ' | / ( |A| n ulp )     ',
     $      '  8= | I - U V', A1, ' | / ( n ulp )',
     $      / '  9= | S - Z D Z', A1, ' | / ( |S| n ulp )     ',
     $      ' 10= | I - Z Z', A1, ' | / ( n ulp )',
     $      / ' 11= |D(with Z) - D(w/o Z)| / (|D| ulp) ',
     $      ' 12= | D(PWK) - D(QR) | / (|D| ulp)',
     $      / ' 13=   Sturm sequence test on W         ' )
 9991 FORMAT( ' 14= | S - Z4 D4 Z4', A1, ' | / (|S| n ulp)',
     $      / ' 15= | I - Z4 Z4', A1, ' | / (n ulp ) ',
     $      ' 16= | D4 - D5 | / ( 100 |D4| ulp ) ',
     $      / ' 17= max | D4(i) - WR(i) | / ( |D4(i)| (2n-1) ulp )',
     $      / ' 18= | WA1 - D3 | / ( |D3| ulp )',
     $      / ' 19= max | WA2(i) - WA3(ii) | / ( |D3| ulp )',
     $      / ' 20= | S - Y WA1 Y', A1, ' | / ( |S| n ulp )',
     $      / ' 21= | I - Y Y', A1, ' | / ( n ulp )' )
 9990 FORMAT( ' N=', I5, ', seed=', 4( I4, ',' ), ' type ', I2,
     $      ', test(', I2, ')=', G10.3 )
 9989 FORMAT( ' 22= | S - Z D Z', A1, '| / ( |S| n ulp ) for DSTEDC(I)',
     $      / ' 23= | I - Z Z', A1, '| / ( n ulp )       for DSTEDC(I)',
     $      / ' 24= | S - Z D Z', A1, '| / ( |S| n ulp ) for DSTEDC(V)',
     $      / ' 25= | I - Z Z', A1, '| / ( n ulp )       for DSTEDC(V)',
     $      / ' 26= | D1(DSTEDC(V)) - D2(SSTEDC(N)) | / ( |D1| ulp )' )
*
 9988 FORMAT( / 'Test performed:  see DCHKST for details.', / )
*     End of DCHKST
*
      END