You can not select more than 25 topics Topics must start with a chinese character,a letter or number, can include dashes ('-') and can be up to 35 characters long.

sdrvst.f 105 kB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014101510161017101810191020102110221023102410251026102710281029103010311032103310341035103610371038103910401041104210431044104510461047104810491050105110521053105410551056105710581059106010611062106310641065106610671068106910701071107210731074107510761077107810791080108110821083108410851086108710881089109010911092109310941095109610971098109911001101110211031104110511061107110811091110111111121113111411151116111711181119112011211122112311241125112611271128112911301131113211331134113511361137113811391140114111421143114411451146114711481149115011511152115311541155115611571158115911601161116211631164116511661167116811691170117111721173117411751176117711781179118011811182118311841185118611871188118911901191119211931194119511961197119811991200120112021203120412051206120712081209121012111212121312141215121612171218121912201221122212231224122512261227122812291230123112321233123412351236123712381239124012411242124312441245124612471248124912501251125212531254125512561257125812591260126112621263126412651266126712681269127012711272127312741275127612771278127912801281128212831284128512861287128812891290129112921293129412951296129712981299130013011302130313041305130613071308130913101311131213131314131513161317131813191320132113221323132413251326132713281329133013311332133313341335133613371338133913401341134213431344134513461347134813491350135113521353135413551356135713581359136013611362136313641365136613671368136913701371137213731374137513761377137813791380138113821383138413851386138713881389139013911392139313941395139613971398139914001401140214031404140514061407140814091410141114121413141414151416141714181419142014211422142314241425142614271428142914301431143214331434143514361437143814391440144114421443144414451446144714481449145014511452145314541455145614571458145914601461146214631464146514661467146814691470147114721473147414751476147714781479148014811482148314841485148614871488148914901491149214931494149514961497149814991500150115021503150415051506150715081509151015111512151315141515151615171518151915201521152215231524152515261527152815291530153115321533153415351536153715381539154015411542154315441545154615471548154915501551155215531554155515561557155815591560156115621563156415651566156715681569157015711572157315741575157615771578157915801581158215831584158515861587158815891590159115921593159415951596159715981599160016011602160316041605160616071608160916101611161216131614161516161617161816191620162116221623162416251626162716281629163016311632163316341635163616371638163916401641164216431644164516461647164816491650165116521653165416551656165716581659166016611662166316641665166616671668166916701671167216731674167516761677167816791680168116821683168416851686168716881689169016911692169316941695169616971698169917001701170217031704170517061707170817091710171117121713171417151716171717181719172017211722172317241725172617271728172917301731173217331734173517361737173817391740174117421743174417451746174717481749175017511752175317541755175617571758175917601761176217631764176517661767176817691770177117721773177417751776177717781779178017811782178317841785178617871788178917901791179217931794179517961797179817991800180118021803180418051806180718081809181018111812181318141815181618171818181918201821182218231824182518261827182818291830183118321833183418351836183718381839184018411842184318441845184618471848184918501851185218531854185518561857185818591860186118621863186418651866186718681869187018711872187318741875187618771878187918801881188218831884188518861887188818891890189118921893189418951896189718981899190019011902190319041905190619071908190919101911191219131914191519161917191819191920192119221923192419251926192719281929193019311932193319341935193619371938193919401941194219431944194519461947194819491950195119521953195419551956195719581959196019611962196319641965196619671968196919701971197219731974197519761977197819791980198119821983198419851986198719881989199019911992199319941995199619971998199920002001200220032004200520062007200820092010201120122013201420152016201720182019202020212022202320242025202620272028202920302031203220332034203520362037203820392040204120422043204420452046204720482049205020512052205320542055205620572058205920602061206220632064206520662067206820692070207120722073207420752076207720782079208020812082208320842085208620872088208920902091209220932094209520962097209820992100210121022103210421052106210721082109211021112112211321142115211621172118211921202121212221232124212521262127212821292130213121322133213421352136213721382139214021412142214321442145214621472148214921502151215221532154215521562157215821592160216121622163216421652166216721682169217021712172217321742175217621772178217921802181218221832184218521862187218821892190219121922193219421952196219721982199220022012202220322042205220622072208220922102211221222132214221522162217221822192220222122222223222422252226222722282229223022312232223322342235223622372238223922402241224222432244224522462247224822492250225122522253225422552256225722582259226022612262226322642265226622672268226922702271227222732274227522762277227822792280228122822283228422852286228722882289229022912292229322942295229622972298229923002301230223032304230523062307230823092310231123122313231423152316231723182319232023212322232323242325232623272328232923302331233223332334233523362337233823392340234123422343234423452346234723482349235023512352235323542355235623572358235923602361236223632364236523662367236823692370237123722373237423752376237723782379238023812382238323842385238623872388238923902391239223932394239523962397239823992400240124022403240424052406240724082409241024112412241324142415241624172418241924202421242224232424242524262427242824292430243124322433243424352436243724382439244024412442244324442445244624472448244924502451245224532454245524562457245824592460246124622463246424652466246724682469247024712472247324742475247624772478247924802481248224832484248524862487248824892490249124922493249424952496249724982499250025012502250325042505250625072508250925102511251225132514251525162517251825192520252125222523252425252526252725282529253025312532253325342535253625372538253925402541254225432544254525462547254825492550255125522553255425552556255725582559256025612562256325642565256625672568256925702571257225732574257525762577257825792580258125822583258425852586258725882589259025912592259325942595259625972598259926002601260226032604260526062607260826092610261126122613261426152616261726182619262026212622262326242625262626272628262926302631263226332634263526362637263826392640264126422643264426452646264726482649265026512652265326542655265626572658265926602661266226632664266526662667266826692670267126722673267426752676267726782679268026812682268326842685268626872688268926902691269226932694269526962697269826992700270127022703270427052706270727082709271027112712271327142715271627172718271927202721272227232724272527262727272827292730273127322733273427352736273727382739274027412742274327442745274627472748274927502751275227532754275527562757275827592760276127622763276427652766276727682769277027712772277327742775277627772778277927802781278227832784278527862787278827892790279127922793279427952796279727982799280028012802280328042805280628072808280928102811281228132814281528162817281828192820282128222823282428252826282728282829283028312832283328342835283628372838283928402841284228432844284528462847
  1. *> \brief \b SDRVST
  2. *
  3. * =========== DOCUMENTATION ===========
  4. *
  5. * Online html documentation available at
  6. * http://www.netlib.org/lapack/explore-html/
  7. *
  8. * Definition:
  9. * ===========
  10. *
  11. * SUBROUTINE SDRVST( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
  12. * NOUNIT, A, LDA, D1, D2, D3, D4, EVEIGS, WA1,
  13. * WA2, WA3, U, LDU, V, TAU, Z, WORK, LWORK,
  14. * IWORK, LIWORK, RESULT, INFO )
  15. *
  16. * .. Scalar Arguments ..
  17. * INTEGER INFO, LDA, LDU, LIWORK, LWORK, NOUNIT, NSIZES,
  18. * $ NTYPES
  19. * REAL THRESH
  20. * ..
  21. * .. Array Arguments ..
  22. * LOGICAL DOTYPE( * )
  23. * INTEGER ISEED( 4 ), IWORK( * ), NN( * )
  24. * REAL A( LDA, * ), D1( * ), D2( * ), D3( * ),
  25. * $ D4( * ), EVEIGS( * ), RESULT( * ), TAU( * ),
  26. * $ U( LDU, * ), V( LDU, * ), WA1( * ), WA2( * ),
  27. * $ WA3( * ), WORK( * ), Z( LDU, * )
  28. * ..
  29. *
  30. *
  31. *> \par Purpose:
  32. * =============
  33. *>
  34. *> \verbatim
  35. *>
  36. *> SDRVST checks the symmetric eigenvalue problem drivers.
  37. *>
  38. *> SSTEV computes all eigenvalues and, optionally,
  39. *> eigenvectors of a real symmetric tridiagonal matrix.
  40. *>
  41. *> SSTEVX computes selected eigenvalues and, optionally,
  42. *> eigenvectors of a real symmetric tridiagonal matrix.
  43. *>
  44. *> SSTEVR computes selected eigenvalues and, optionally,
  45. *> eigenvectors of a real symmetric tridiagonal matrix
  46. *> using the Relatively Robust Representation where it can.
  47. *>
  48. *> SSYEV computes all eigenvalues and, optionally,
  49. *> eigenvectors of a real symmetric matrix.
  50. *>
  51. *> SSYEVX computes selected eigenvalues and, optionally,
  52. *> eigenvectors of a real symmetric matrix.
  53. *>
  54. *> SSYEVR computes selected eigenvalues and, optionally,
  55. *> eigenvectors of a real symmetric matrix
  56. *> using the Relatively Robust Representation where it can.
  57. *>
  58. *> SSPEV computes all eigenvalues and, optionally,
  59. *> eigenvectors of a real symmetric matrix in packed
  60. *> storage.
  61. *>
  62. *> SSPEVX computes selected eigenvalues and, optionally,
  63. *> eigenvectors of a real symmetric matrix in packed
  64. *> storage.
  65. *>
  66. *> SSBEV computes all eigenvalues and, optionally,
  67. *> eigenvectors of a real symmetric band matrix.
  68. *>
  69. *> SSBEVX computes selected eigenvalues and, optionally,
  70. *> eigenvectors of a real symmetric band matrix.
  71. *>
  72. *> SSYEVD computes all eigenvalues and, optionally,
  73. *> eigenvectors of a real symmetric matrix using
  74. *> a divide and conquer algorithm.
  75. *>
  76. *> SSPEVD computes all eigenvalues and, optionally,
  77. *> eigenvectors of a real symmetric matrix in packed
  78. *> storage, using a divide and conquer algorithm.
  79. *>
  80. *> SSBEVD computes all eigenvalues and, optionally,
  81. *> eigenvectors of a real symmetric band matrix,
  82. *> using a divide and conquer algorithm.
  83. *>
  84. *> When SDRVST is called, a number of matrix "sizes" ("n's") and a
  85. *> number of matrix "types" are specified. For each size ("n")
  86. *> and each type of matrix, one matrix will be generated and used
  87. *> to test the appropriate drivers. For each matrix and each
  88. *> driver routine called, the following tests will be performed:
  89. *>
  90. *> (1) | A - Z D Z' | / ( |A| n ulp )
  91. *>
  92. *> (2) | I - Z Z' | / ( n ulp )
  93. *>
  94. *> (3) | D1 - D2 | / ( |D1| ulp )
  95. *>
  96. *> where Z is the matrix of eigenvectors returned when the
  97. *> eigenvector option is given and D1 and D2 are the eigenvalues
  98. *> returned with and without the eigenvector option.
  99. *>
  100. *> The "sizes" are specified by an array NN(1:NSIZES); the value of
  101. *> each element NN(j) specifies one size.
  102. *> The "types" are specified by a logical array DOTYPE( 1:NTYPES );
  103. *> if DOTYPE(j) is .TRUE., then matrix type "j" will be generated.
  104. *> Currently, the list of possible types is:
  105. *>
  106. *> (1) The zero matrix.
  107. *> (2) The identity matrix.
  108. *>
  109. *> (3) A diagonal matrix with evenly spaced eigenvalues
  110. *> 1, ..., ULP and random signs.
  111. *> (ULP = (first number larger than 1) - 1 )
  112. *> (4) A diagonal matrix with geometrically spaced eigenvalues
  113. *> 1, ..., ULP and random signs.
  114. *> (5) A diagonal matrix with "clustered" eigenvalues
  115. *> 1, ULP, ..., ULP and random signs.
  116. *>
  117. *> (6) Same as (4), but multiplied by SQRT( overflow threshold )
  118. *> (7) Same as (4), but multiplied by SQRT( underflow threshold )
  119. *>
  120. *> (8) A matrix of the form U' D U, where U is orthogonal and
  121. *> D has evenly spaced entries 1, ..., ULP with random signs
  122. *> on the diagonal.
  123. *>
  124. *> (9) A matrix of the form U' D U, where U is orthogonal and
  125. *> D has geometrically spaced entries 1, ..., ULP with random
  126. *> signs on the diagonal.
  127. *>
  128. *> (10) A matrix of the form U' D U, where U is orthogonal and
  129. *> D has "clustered" entries 1, ULP,..., ULP with random
  130. *> signs on the diagonal.
  131. *>
  132. *> (11) Same as (8), but multiplied by SQRT( overflow threshold )
  133. *> (12) Same as (8), but multiplied by SQRT( underflow threshold )
  134. *>
  135. *> (13) Symmetric matrix with random entries chosen from (-1,1).
  136. *> (14) Same as (13), but multiplied by SQRT( overflow threshold )
  137. *> (15) Same as (13), but multiplied by SQRT( underflow threshold )
  138. *> (16) A band matrix with half bandwidth randomly chosen between
  139. *> 0 and N-1, with evenly spaced eigenvalues 1, ..., ULP
  140. *> with random signs.
  141. *> (17) Same as (16), but multiplied by SQRT( overflow threshold )
  142. *> (18) Same as (16), but multiplied by SQRT( underflow threshold )
  143. *> \endverbatim
  144. *
  145. * Arguments:
  146. * ==========
  147. *
  148. *> \verbatim
  149. *> NSIZES INTEGER
  150. *> The number of sizes of matrices to use. If it is zero,
  151. *> SDRVST does nothing. It must be at least zero.
  152. *> Not modified.
  153. *>
  154. *> NN INTEGER array, dimension (NSIZES)
  155. *> An array containing the sizes to be used for the matrices.
  156. *> Zero values will be skipped. The values must be at least
  157. *> zero.
  158. *> Not modified.
  159. *>
  160. *> NTYPES INTEGER
  161. *> The number of elements in DOTYPE. If it is zero, SDRVST
  162. *> does nothing. It must be at least zero. If it is MAXTYP+1
  163. *> and NSIZES is 1, then an additional type, MAXTYP+1 is
  164. *> defined, which is to use whatever matrix is in A. This
  165. *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and
  166. *> DOTYPE(MAXTYP+1) is .TRUE. .
  167. *> Not modified.
  168. *>
  169. *> DOTYPE LOGICAL array, dimension (NTYPES)
  170. *> If DOTYPE(j) is .TRUE., then for each size in NN a
  171. *> matrix of that size and of type j will be generated.
  172. *> If NTYPES is smaller than the maximum number of types
  173. *> defined (PARAMETER MAXTYP), then types NTYPES+1 through
  174. *> MAXTYP will not be generated. If NTYPES is larger
  175. *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES)
  176. *> will be ignored.
  177. *> Not modified.
  178. *>
  179. *> ISEED INTEGER array, dimension (4)
  180. *> On entry ISEED specifies the seed of the random number
  181. *> generator. The array elements should be between 0 and 4095;
  182. *> if not they will be reduced mod 4096. Also, ISEED(4) must
  183. *> be odd. The random number generator uses a linear
  184. *> congruential sequence limited to small integers, and so
  185. *> should produce machine independent random numbers. The
  186. *> values of ISEED are changed on exit, and can be used in the
  187. *> next call to SDRVST to continue the same random number
  188. *> sequence.
  189. *> Modified.
  190. *>
  191. *> THRESH REAL
  192. *> A test will count as "failed" if the "error", computed as
  193. *> described above, exceeds THRESH. Note that the error
  194. *> is scaled to be O(1), so THRESH should be a reasonably
  195. *> small multiple of 1, e.g., 10 or 100. In particular,
  196. *> it should not depend on the precision (single vs. double)
  197. *> or the size of the matrix. It must be at least zero.
  198. *> Not modified.
  199. *>
  200. *> NOUNIT INTEGER
  201. *> The FORTRAN unit number for printing out error messages
  202. *> (e.g., if a routine returns IINFO not equal to 0.)
  203. *> Not modified.
  204. *>
  205. *> A REAL array, dimension (LDA , max(NN))
  206. *> Used to hold the matrix whose eigenvalues are to be
  207. *> computed. On exit, A contains the last matrix actually
  208. *> used.
  209. *> Modified.
  210. *>
  211. *> LDA INTEGER
  212. *> The leading dimension of A. It must be at
  213. *> least 1 and at least max( NN ).
  214. *> Not modified.
  215. *>
  216. *> D1 REAL array, dimension (max(NN))
  217. *> The eigenvalues of A, as computed by SSTEQR simlutaneously
  218. *> with Z. On exit, the eigenvalues in D1 correspond with the
  219. *> matrix in A.
  220. *> Modified.
  221. *>
  222. *> D2 REAL array, dimension (max(NN))
  223. *> The eigenvalues of A, as computed by SSTEQR if Z is not
  224. *> computed. On exit, the eigenvalues in D2 correspond with
  225. *> the matrix in A.
  226. *> Modified.
  227. *>
  228. *> D3 REAL array, dimension (max(NN))
  229. *> The eigenvalues of A, as computed by SSTERF. On exit, the
  230. *> eigenvalues in D3 correspond with the matrix in A.
  231. *> Modified.
  232. *>
  233. *> D4 REAL array, dimension
  234. *>
  235. *> EVEIGS REAL array, dimension (max(NN))
  236. *> The eigenvalues as computed by SSTEV('N', ... )
  237. *> (I reserve the right to change this to the output of
  238. *> whichever algorithm computes the most accurate eigenvalues).
  239. *>
  240. *> WA1 REAL array, dimension
  241. *>
  242. *> WA2 REAL array, dimension
  243. *>
  244. *> WA3 REAL array, dimension
  245. *>
  246. *> U REAL array, dimension (LDU, max(NN))
  247. *> The orthogonal matrix computed by SSYTRD + SORGTR.
  248. *> Modified.
  249. *>
  250. *> LDU INTEGER
  251. *> The leading dimension of U, Z, and V. It must be at
  252. *> least 1 and at least max( NN ).
  253. *> Not modified.
  254. *>
  255. *> V REAL array, dimension (LDU, max(NN))
  256. *> The Housholder vectors computed by SSYTRD in reducing A to
  257. *> tridiagonal form.
  258. *> Modified.
  259. *>
  260. *> TAU REAL array, dimension (max(NN))
  261. *> The Householder factors computed by SSYTRD in reducing A
  262. *> to tridiagonal form.
  263. *> Modified.
  264. *>
  265. *> Z REAL array, dimension (LDU, max(NN))
  266. *> The orthogonal matrix of eigenvectors computed by SSTEQR,
  267. *> SPTEQR, and SSTEIN.
  268. *> Modified.
  269. *>
  270. *> WORK REAL array, dimension (LWORK)
  271. *> Workspace.
  272. *> Modified.
  273. *>
  274. *> LWORK INTEGER
  275. *> The number of entries in WORK. This must be at least
  276. *> 1 + 4 * Nmax + 2 * Nmax * lg Nmax + 4 * Nmax**2
  277. *> where Nmax = max( NN(j), 2 ) and lg = log base 2.
  278. *> Not modified.
  279. *>
  280. *> IWORK INTEGER array,
  281. *> dimension (6 + 6*Nmax + 5 * Nmax * lg Nmax )
  282. *> where Nmax = max( NN(j), 2 ) and lg = log base 2.
  283. *> Workspace.
  284. *> Modified.
  285. *>
  286. *> RESULT REAL array, dimension (105)
  287. *> The values computed by the tests described above.
  288. *> The values are currently limited to 1/ulp, to avoid
  289. *> overflow.
  290. *> Modified.
  291. *>
  292. *> INFO INTEGER
  293. *> If 0, then everything ran OK.
  294. *> -1: NSIZES < 0
  295. *> -2: Some NN(j) < 0
  296. *> -3: NTYPES < 0
  297. *> -5: THRESH < 0
  298. *> -9: LDA < 1 or LDA < NMAX, where NMAX is max( NN(j) ).
  299. *> -16: LDU < 1 or LDU < NMAX.
  300. *> -21: LWORK too small.
  301. *> If SLATMR, SLATMS, SSYTRD, SORGTR, SSTEQR, SSTERF,
  302. *> or SORMTR returns an error code, the
  303. *> absolute value of it is returned.
  304. *> Modified.
  305. *>
  306. *>-----------------------------------------------------------------------
  307. *>
  308. *> Some Local Variables and Parameters:
  309. *> ---- ----- --------- --- ----------
  310. *> ZERO, ONE Real 0 and 1.
  311. *> MAXTYP The number of types defined.
  312. *> NTEST The number of tests performed, or which can
  313. *> be performed so far, for the current matrix.
  314. *> NTESTT The total number of tests performed so far.
  315. *> NMAX Largest value in NN.
  316. *> NMATS The number of matrices generated so far.
  317. *> NERRS The number of tests which have exceeded THRESH
  318. *> so far (computed by SLAFTS).
  319. *> COND, IMODE Values to be passed to the matrix generators.
  320. *> ANORM Norm of A; passed to matrix generators.
  321. *>
  322. *> OVFL, UNFL Overflow and underflow thresholds.
  323. *> ULP, ULPINV Finest relative precision and its inverse.
  324. *> RTOVFL, RTUNFL Square roots of the previous 2 values.
  325. *> The following four arrays decode JTYPE:
  326. *> KTYPE(j) The general type (1-10) for type "j".
  327. *> KMODE(j) The MODE value to be passed to the matrix
  328. *> generator for type "j".
  329. *> KMAGN(j) The order of magnitude ( O(1),
  330. *> O(overflow^(1/2) ), O(underflow^(1/2) )
  331. *>
  332. *> The tests performed are: Routine tested
  333. *> 1= | A - U S U' | / ( |A| n ulp ) SSTEV('V', ... )
  334. *> 2= | I - U U' | / ( n ulp ) SSTEV('V', ... )
  335. *> 3= |D(with Z) - D(w/o Z)| / (|D| ulp) SSTEV('N', ... )
  336. *> 4= | A - U S U' | / ( |A| n ulp ) SSTEVX('V','A', ... )
  337. *> 5= | I - U U' | / ( n ulp ) SSTEVX('V','A', ... )
  338. *> 6= |D(with Z) - EVEIGS| / (|D| ulp) SSTEVX('N','A', ... )
  339. *> 7= | A - U S U' | / ( |A| n ulp ) SSTEVR('V','A', ... )
  340. *> 8= | I - U U' | / ( n ulp ) SSTEVR('V','A', ... )
  341. *> 9= |D(with Z) - EVEIGS| / (|D| ulp) SSTEVR('N','A', ... )
  342. *> 10= | A - U S U' | / ( |A| n ulp ) SSTEVX('V','I', ... )
  343. *> 11= | I - U U' | / ( n ulp ) SSTEVX('V','I', ... )
  344. *> 12= |D(with Z) - D(w/o Z)| / (|D| ulp) SSTEVX('N','I', ... )
  345. *> 13= | A - U S U' | / ( |A| n ulp ) SSTEVX('V','V', ... )
  346. *> 14= | I - U U' | / ( n ulp ) SSTEVX('V','V', ... )
  347. *> 15= |D(with Z) - D(w/o Z)| / (|D| ulp) SSTEVX('N','V', ... )
  348. *> 16= | A - U S U' | / ( |A| n ulp ) SSTEVD('V', ... )
  349. *> 17= | I - U U' | / ( n ulp ) SSTEVD('V', ... )
  350. *> 18= |D(with Z) - EVEIGS| / (|D| ulp) SSTEVD('N', ... )
  351. *> 19= | A - U S U' | / ( |A| n ulp ) SSTEVR('V','I', ... )
  352. *> 20= | I - U U' | / ( n ulp ) SSTEVR('V','I', ... )
  353. *> 21= |D(with Z) - D(w/o Z)| / (|D| ulp) SSTEVR('N','I', ... )
  354. *> 22= | A - U S U' | / ( |A| n ulp ) SSTEVR('V','V', ... )
  355. *> 23= | I - U U' | / ( n ulp ) SSTEVR('V','V', ... )
  356. *> 24= |D(with Z) - D(w/o Z)| / (|D| ulp) SSTEVR('N','V', ... )
  357. *>
  358. *> 25= | A - U S U' | / ( |A| n ulp ) SSYEV('L','V', ... )
  359. *> 26= | I - U U' | / ( n ulp ) SSYEV('L','V', ... )
  360. *> 27= |D(with Z) - D(w/o Z)| / (|D| ulp) SSYEV('L','N', ... )
  361. *> 28= | A - U S U' | / ( |A| n ulp ) SSYEVX('L','V','A', ... )
  362. *> 29= | I - U U' | / ( n ulp ) SSYEVX('L','V','A', ... )
  363. *> 30= |D(with Z) - D(w/o Z)| / (|D| ulp) SSYEVX('L','N','A', ... )
  364. *> 31= | A - U S U' | / ( |A| n ulp ) SSYEVX('L','V','I', ... )
  365. *> 32= | I - U U' | / ( n ulp ) SSYEVX('L','V','I', ... )
  366. *> 33= |D(with Z) - D(w/o Z)| / (|D| ulp) SSYEVX('L','N','I', ... )
  367. *> 34= | A - U S U' | / ( |A| n ulp ) SSYEVX('L','V','V', ... )
  368. *> 35= | I - U U' | / ( n ulp ) SSYEVX('L','V','V', ... )
  369. *> 36= |D(with Z) - D(w/o Z)| / (|D| ulp) SSYEVX('L','N','V', ... )
  370. *> 37= | A - U S U' | / ( |A| n ulp ) SSPEV('L','V', ... )
  371. *> 38= | I - U U' | / ( n ulp ) SSPEV('L','V', ... )
  372. *> 39= |D(with Z) - D(w/o Z)| / (|D| ulp) SSPEV('L','N', ... )
  373. *> 40= | A - U S U' | / ( |A| n ulp ) SSPEVX('L','V','A', ... )
  374. *> 41= | I - U U' | / ( n ulp ) SSPEVX('L','V','A', ... )
  375. *> 42= |D(with Z) - D(w/o Z)| / (|D| ulp) SSPEVX('L','N','A', ... )
  376. *> 43= | A - U S U' | / ( |A| n ulp ) SSPEVX('L','V','I', ... )
  377. *> 44= | I - U U' | / ( n ulp ) SSPEVX('L','V','I', ... )
  378. *> 45= |D(with Z) - D(w/o Z)| / (|D| ulp) SSPEVX('L','N','I', ... )
  379. *> 46= | A - U S U' | / ( |A| n ulp ) SSPEVX('L','V','V', ... )
  380. *> 47= | I - U U' | / ( n ulp ) SSPEVX('L','V','V', ... )
  381. *> 48= |D(with Z) - D(w/o Z)| / (|D| ulp) SSPEVX('L','N','V', ... )
  382. *> 49= | A - U S U' | / ( |A| n ulp ) SSBEV('L','V', ... )
  383. *> 50= | I - U U' | / ( n ulp ) SSBEV('L','V', ... )
  384. *> 51= |D(with Z) - D(w/o Z)| / (|D| ulp) SSBEV('L','N', ... )
  385. *> 52= | A - U S U' | / ( |A| n ulp ) SSBEVX('L','V','A', ... )
  386. *> 53= | I - U U' | / ( n ulp ) SSBEVX('L','V','A', ... )
  387. *> 54= |D(with Z) - D(w/o Z)| / (|D| ulp) SSBEVX('L','N','A', ... )
  388. *> 55= | A - U S U' | / ( |A| n ulp ) SSBEVX('L','V','I', ... )
  389. *> 56= | I - U U' | / ( n ulp ) SSBEVX('L','V','I', ... )
  390. *> 57= |D(with Z) - D(w/o Z)| / (|D| ulp) SSBEVX('L','N','I', ... )
  391. *> 58= | A - U S U' | / ( |A| n ulp ) SSBEVX('L','V','V', ... )
  392. *> 59= | I - U U' | / ( n ulp ) SSBEVX('L','V','V', ... )
  393. *> 60= |D(with Z) - D(w/o Z)| / (|D| ulp) SSBEVX('L','N','V', ... )
  394. *> 61= | A - U S U' | / ( |A| n ulp ) SSYEVD('L','V', ... )
  395. *> 62= | I - U U' | / ( n ulp ) SSYEVD('L','V', ... )
  396. *> 63= |D(with Z) - D(w/o Z)| / (|D| ulp) SSYEVD('L','N', ... )
  397. *> 64= | A - U S U' | / ( |A| n ulp ) SSPEVD('L','V', ... )
  398. *> 65= | I - U U' | / ( n ulp ) SSPEVD('L','V', ... )
  399. *> 66= |D(with Z) - D(w/o Z)| / (|D| ulp) SSPEVD('L','N', ... )
  400. *> 67= | A - U S U' | / ( |A| n ulp ) SSBEVD('L','V', ... )
  401. *> 68= | I - U U' | / ( n ulp ) SSBEVD('L','V', ... )
  402. *> 69= |D(with Z) - D(w/o Z)| / (|D| ulp) SSBEVD('L','N', ... )
  403. *> 70= | A - U S U' | / ( |A| n ulp ) SSYEVR('L','V','A', ... )
  404. *> 71= | I - U U' | / ( n ulp ) SSYEVR('L','V','A', ... )
  405. *> 72= |D(with Z) - D(w/o Z)| / (|D| ulp) SSYEVR('L','N','A', ... )
  406. *> 73= | A - U S U' | / ( |A| n ulp ) SSYEVR('L','V','I', ... )
  407. *> 74= | I - U U' | / ( n ulp ) SSYEVR('L','V','I', ... )
  408. *> 75= |D(with Z) - D(w/o Z)| / (|D| ulp) SSYEVR('L','N','I', ... )
  409. *> 76= | A - U S U' | / ( |A| n ulp ) SSYEVR('L','V','V', ... )
  410. *> 77= | I - U U' | / ( n ulp ) SSYEVR('L','V','V', ... )
  411. *> 78= |D(with Z) - D(w/o Z)| / (|D| ulp) SSYEVR('L','N','V', ... )
  412. *>
  413. *> Tests 25 through 78 are repeated (as tests 79 through 132)
  414. *> with UPLO='U'
  415. *>
  416. *> To be added in 1999
  417. *>
  418. *> 79= | A - U S U' | / ( |A| n ulp ) SSPEVR('L','V','A', ... )
  419. *> 80= | I - U U' | / ( n ulp ) SSPEVR('L','V','A', ... )
  420. *> 81= |D(with Z) - D(w/o Z)| / (|D| ulp) SSPEVR('L','N','A', ... )
  421. *> 82= | A - U S U' | / ( |A| n ulp ) SSPEVR('L','V','I', ... )
  422. *> 83= | I - U U' | / ( n ulp ) SSPEVR('L','V','I', ... )
  423. *> 84= |D(with Z) - D(w/o Z)| / (|D| ulp) SSPEVR('L','N','I', ... )
  424. *> 85= | A - U S U' | / ( |A| n ulp ) SSPEVR('L','V','V', ... )
  425. *> 86= | I - U U' | / ( n ulp ) SSPEVR('L','V','V', ... )
  426. *> 87= |D(with Z) - D(w/o Z)| / (|D| ulp) SSPEVR('L','N','V', ... )
  427. *> 88= | A - U S U' | / ( |A| n ulp ) SSBEVR('L','V','A', ... )
  428. *> 89= | I - U U' | / ( n ulp ) SSBEVR('L','V','A', ... )
  429. *> 90= |D(with Z) - D(w/o Z)| / (|D| ulp) SSBEVR('L','N','A', ... )
  430. *> 91= | A - U S U' | / ( |A| n ulp ) SSBEVR('L','V','I', ... )
  431. *> 92= | I - U U' | / ( n ulp ) SSBEVR('L','V','I', ... )
  432. *> 93= |D(with Z) - D(w/o Z)| / (|D| ulp) SSBEVR('L','N','I', ... )
  433. *> 94= | A - U S U' | / ( |A| n ulp ) SSBEVR('L','V','V', ... )
  434. *> 95= | I - U U' | / ( n ulp ) SSBEVR('L','V','V', ... )
  435. *> 96= |D(with Z) - D(w/o Z)| / (|D| ulp) SSBEVR('L','N','V', ... )
  436. *> \endverbatim
  437. *
  438. * Authors:
  439. * ========
  440. *
  441. *> \author Univ. of Tennessee
  442. *> \author Univ. of California Berkeley
  443. *> \author Univ. of Colorado Denver
  444. *> \author NAG Ltd.
  445. *
  446. *> \ingroup single_eig
  447. *
  448. * =====================================================================
  449. SUBROUTINE SDRVST( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
  450. $ NOUNIT, A, LDA, D1, D2, D3, D4, EVEIGS, WA1,
  451. $ WA2, WA3, U, LDU, V, TAU, Z, WORK, LWORK,
  452. $ IWORK, LIWORK, RESULT, INFO )
  453. *
  454. * -- LAPACK test routine --
  455. * -- LAPACK is a software package provided by Univ. of Tennessee, --
  456. * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
  457. *
  458. * .. Scalar Arguments ..
  459. INTEGER INFO, LDA, LDU, LIWORK, LWORK, NOUNIT, NSIZES,
  460. $ NTYPES
  461. REAL THRESH
  462. * ..
  463. * .. Array Arguments ..
  464. LOGICAL DOTYPE( * )
  465. INTEGER ISEED( 4 ), IWORK( * ), NN( * )
  466. REAL A( LDA, * ), D1( * ), D2( * ), D3( * ),
  467. $ D4( * ), EVEIGS( * ), RESULT( * ), TAU( * ),
  468. $ U( LDU, * ), V( LDU, * ), WA1( * ), WA2( * ),
  469. $ WA3( * ), WORK( * ), Z( LDU, * )
  470. * ..
  471. *
  472. * =====================================================================
  473. *
  474. * .. Parameters ..
  475. REAL ZERO, ONE, TWO, TEN
  476. PARAMETER ( ZERO = 0.0E0, ONE = 1.0E0, TWO = 2.0E0,
  477. $ TEN = 10.0E0 )
  478. REAL HALF
  479. PARAMETER ( HALF = 0.5E0 )
  480. INTEGER MAXTYP
  481. PARAMETER ( MAXTYP = 18 )
  482. * ..
  483. * .. Local Scalars ..
  484. LOGICAL BADNN
  485. CHARACTER UPLO
  486. INTEGER I, IDIAG, IHBW, IINFO, IL, IMODE, INDX, IROW,
  487. $ ITEMP, ITYPE, IU, IUPLO, J, J1, J2, JCOL,
  488. $ JSIZE, JTYPE, KD, LGN, LIWEDC, LWEDC, M, M2,
  489. $ M3, MTYPES, N, NERRS, NMATS, NMAX, NTEST,
  490. $ NTESTT
  491. REAL ABSTOL, ANINV, ANORM, COND, OVFL, RTOVFL,
  492. $ RTUNFL, TEMP1, TEMP2, TEMP3, ULP, ULPINV, UNFL,
  493. $ VL, VU
  494. * ..
  495. * .. Local Arrays ..
  496. INTEGER IDUMMA( 1 ), IOLDSD( 4 ), ISEED2( 4 ),
  497. $ ISEED3( 4 ), KMAGN( MAXTYP ), KMODE( MAXTYP ),
  498. $ KTYPE( MAXTYP )
  499. * ..
  500. * .. External Functions ..
  501. REAL SLAMCH, SLARND, SSXT1
  502. EXTERNAL SLAMCH, SLARND, SSXT1
  503. * ..
  504. * .. External Subroutines ..
  505. EXTERNAL ALASVM, SLABAD, SLACPY, SLAFTS, SLASET, SLATMR,
  506. $ SLATMS, SSBEV, SSBEVD, SSBEVX, SSPEV, SSPEVD,
  507. $ SSPEVX, SSTEV, SSTEVD, SSTEVR, SSTEVX, SSTT21,
  508. $ SSTT22, SSYEV, SSYEVD, SSYEVR, SSYEVX, SSYT21,
  509. $ SSYT22, XERBLA
  510. * ..
  511. * .. Scalars in Common ..
  512. CHARACTER*32 SRNAMT
  513. * ..
  514. * .. Common blocks ..
  515. COMMON / SRNAMC / SRNAMT
  516. * ..
  517. * .. Intrinsic Functions ..
  518. INTRINSIC ABS, INT, LOG, MAX, MIN, REAL, SQRT
  519. * ..
  520. * .. Data statements ..
  521. DATA KTYPE / 1, 2, 5*4, 5*5, 3*8, 3*9 /
  522. DATA KMAGN / 2*1, 1, 1, 1, 2, 3, 1, 1, 1, 2, 3, 1,
  523. $ 2, 3, 1, 2, 3 /
  524. DATA KMODE / 2*0, 4, 3, 1, 4, 4, 4, 3, 1, 4, 4, 0,
  525. $ 0, 0, 4, 4, 4 /
  526. * ..
  527. * .. Executable Statements ..
  528. *
  529. * Keep ftrnchek happy
  530. *
  531. VL = ZERO
  532. VU = ZERO
  533. *
  534. * 1) Check for errors
  535. *
  536. NTESTT = 0
  537. INFO = 0
  538. *
  539. BADNN = .FALSE.
  540. NMAX = 1
  541. DO 10 J = 1, NSIZES
  542. NMAX = MAX( NMAX, NN( J ) )
  543. IF( NN( J ).LT.0 )
  544. $ BADNN = .TRUE.
  545. 10 CONTINUE
  546. *
  547. * Check for errors
  548. *
  549. IF( NSIZES.LT.0 ) THEN
  550. INFO = -1
  551. ELSE IF( BADNN ) THEN
  552. INFO = -2
  553. ELSE IF( NTYPES.LT.0 ) THEN
  554. INFO = -3
  555. ELSE IF( LDA.LT.NMAX ) THEN
  556. INFO = -9
  557. ELSE IF( LDU.LT.NMAX ) THEN
  558. INFO = -16
  559. ELSE IF( 2*MAX( 2, NMAX )**2.GT.LWORK ) THEN
  560. INFO = -21
  561. END IF
  562. *
  563. IF( INFO.NE.0 ) THEN
  564. CALL XERBLA( 'SDRVST', -INFO )
  565. RETURN
  566. END IF
  567. *
  568. * Quick return if nothing to do
  569. *
  570. IF( NSIZES.EQ.0 .OR. NTYPES.EQ.0 )
  571. $ RETURN
  572. *
  573. * More Important constants
  574. *
  575. UNFL = SLAMCH( 'Safe minimum' )
  576. OVFL = SLAMCH( 'Overflow' )
  577. CALL SLABAD( UNFL, OVFL )
  578. ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' )
  579. ULPINV = ONE / ULP
  580. RTUNFL = SQRT( UNFL )
  581. RTOVFL = SQRT( OVFL )
  582. *
  583. * Loop over sizes, types
  584. *
  585. DO 20 I = 1, 4
  586. ISEED2( I ) = ISEED( I )
  587. ISEED3( I ) = ISEED( I )
  588. 20 CONTINUE
  589. *
  590. NERRS = 0
  591. NMATS = 0
  592. *
  593. *
  594. DO 1740 JSIZE = 1, NSIZES
  595. N = NN( JSIZE )
  596. IF( N.GT.0 ) THEN
  597. LGN = INT( LOG( REAL( N ) ) / LOG( TWO ) )
  598. IF( 2**LGN.LT.N )
  599. $ LGN = LGN + 1
  600. IF( 2**LGN.LT.N )
  601. $ LGN = LGN + 1
  602. LWEDC = 1 + 4*N + 2*N*LGN + 4*N**2
  603. c LIWEDC = 6 + 6*N + 5*N*LGN
  604. LIWEDC = 3 + 5*N
  605. ELSE
  606. LWEDC = 9
  607. c LIWEDC = 12
  608. LIWEDC = 8
  609. END IF
  610. ANINV = ONE / REAL( MAX( 1, N ) )
  611. *
  612. IF( NSIZES.NE.1 ) THEN
  613. MTYPES = MIN( MAXTYP, NTYPES )
  614. ELSE
  615. MTYPES = MIN( MAXTYP+1, NTYPES )
  616. END IF
  617. *
  618. DO 1730 JTYPE = 1, MTYPES
  619. *
  620. IF( .NOT.DOTYPE( JTYPE ) )
  621. $ GO TO 1730
  622. NMATS = NMATS + 1
  623. NTEST = 0
  624. *
  625. DO 30 J = 1, 4
  626. IOLDSD( J ) = ISEED( J )
  627. 30 CONTINUE
  628. *
  629. * 2) Compute "A"
  630. *
  631. * Control parameters:
  632. *
  633. * KMAGN KMODE KTYPE
  634. * =1 O(1) clustered 1 zero
  635. * =2 large clustered 2 identity
  636. * =3 small exponential (none)
  637. * =4 arithmetic diagonal, (w/ eigenvalues)
  638. * =5 random log symmetric, w/ eigenvalues
  639. * =6 random (none)
  640. * =7 random diagonal
  641. * =8 random symmetric
  642. * =9 band symmetric, w/ eigenvalues
  643. *
  644. IF( MTYPES.GT.MAXTYP )
  645. $ GO TO 110
  646. *
  647. ITYPE = KTYPE( JTYPE )
  648. IMODE = KMODE( JTYPE )
  649. *
  650. * Compute norm
  651. *
  652. GO TO ( 40, 50, 60 )KMAGN( JTYPE )
  653. *
  654. 40 CONTINUE
  655. ANORM = ONE
  656. GO TO 70
  657. *
  658. 50 CONTINUE
  659. ANORM = ( RTOVFL*ULP )*ANINV
  660. GO TO 70
  661. *
  662. 60 CONTINUE
  663. ANORM = RTUNFL*N*ULPINV
  664. GO TO 70
  665. *
  666. 70 CONTINUE
  667. *
  668. CALL SLASET( 'Full', LDA, N, ZERO, ZERO, A, LDA )
  669. IINFO = 0
  670. COND = ULPINV
  671. *
  672. * Special Matrices -- Identity & Jordan block
  673. *
  674. * Zero
  675. *
  676. IF( ITYPE.EQ.1 ) THEN
  677. IINFO = 0
  678. *
  679. ELSE IF( ITYPE.EQ.2 ) THEN
  680. *
  681. * Identity
  682. *
  683. DO 80 JCOL = 1, N
  684. A( JCOL, JCOL ) = ANORM
  685. 80 CONTINUE
  686. *
  687. ELSE IF( ITYPE.EQ.4 ) THEN
  688. *
  689. * Diagonal Matrix, [Eigen]values Specified
  690. *
  691. CALL SLATMS( N, N, 'S', ISEED, 'S', WORK, IMODE, COND,
  692. $ ANORM, 0, 0, 'N', A, LDA, WORK( N+1 ),
  693. $ IINFO )
  694. *
  695. ELSE IF( ITYPE.EQ.5 ) THEN
  696. *
  697. * Symmetric, eigenvalues specified
  698. *
  699. CALL SLATMS( N, N, 'S', ISEED, 'S', WORK, IMODE, COND,
  700. $ ANORM, N, N, 'N', A, LDA, WORK( N+1 ),
  701. $ IINFO )
  702. *
  703. ELSE IF( ITYPE.EQ.7 ) THEN
  704. *
  705. * Diagonal, random eigenvalues
  706. *
  707. IDUMMA( 1 ) = 1
  708. CALL SLATMR( N, N, 'S', ISEED, 'S', WORK, 6, ONE, ONE,
  709. $ 'T', 'N', WORK( N+1 ), 1, ONE,
  710. $ WORK( 2*N+1 ), 1, ONE, 'N', IDUMMA, 0, 0,
  711. $ ZERO, ANORM, 'NO', A, LDA, IWORK, IINFO )
  712. *
  713. ELSE IF( ITYPE.EQ.8 ) THEN
  714. *
  715. * Symmetric, random eigenvalues
  716. *
  717. IDUMMA( 1 ) = 1
  718. CALL SLATMR( N, N, 'S', ISEED, 'S', WORK, 6, ONE, ONE,
  719. $ 'T', 'N', WORK( N+1 ), 1, ONE,
  720. $ WORK( 2*N+1 ), 1, ONE, 'N', IDUMMA, N, N,
  721. $ ZERO, ANORM, 'NO', A, LDA, IWORK, IINFO )
  722. *
  723. ELSE IF( ITYPE.EQ.9 ) THEN
  724. *
  725. * Symmetric banded, eigenvalues specified
  726. *
  727. IHBW = INT( ( N-1 )*SLARND( 1, ISEED3 ) )
  728. CALL SLATMS( N, N, 'S', ISEED, 'S', WORK, IMODE, COND,
  729. $ ANORM, IHBW, IHBW, 'Z', U, LDU, WORK( N+1 ),
  730. $ IINFO )
  731. *
  732. * Store as dense matrix for most routines.
  733. *
  734. CALL SLASET( 'Full', LDA, N, ZERO, ZERO, A, LDA )
  735. DO 100 IDIAG = -IHBW, IHBW
  736. IROW = IHBW - IDIAG + 1
  737. J1 = MAX( 1, IDIAG+1 )
  738. J2 = MIN( N, N+IDIAG )
  739. DO 90 J = J1, J2
  740. I = J - IDIAG
  741. A( I, J ) = U( IROW, J )
  742. 90 CONTINUE
  743. 100 CONTINUE
  744. ELSE
  745. IINFO = 1
  746. END IF
  747. *
  748. IF( IINFO.NE.0 ) THEN
  749. WRITE( NOUNIT, FMT = 9999 )'Generator', IINFO, N, JTYPE,
  750. $ IOLDSD
  751. INFO = ABS( IINFO )
  752. RETURN
  753. END IF
  754. *
  755. 110 CONTINUE
  756. *
  757. ABSTOL = UNFL + UNFL
  758. IF( N.LE.1 ) THEN
  759. IL = 1
  760. IU = N
  761. ELSE
  762. IL = 1 + INT( ( N-1 )*SLARND( 1, ISEED2 ) )
  763. IU = 1 + INT( ( N-1 )*SLARND( 1, ISEED2 ) )
  764. IF( IL.GT.IU ) THEN
  765. ITEMP = IL
  766. IL = IU
  767. IU = ITEMP
  768. END IF
  769. END IF
  770. *
  771. * 3) If matrix is tridiagonal, call SSTEV and SSTEVX.
  772. *
  773. IF( JTYPE.LE.7 ) THEN
  774. NTEST = 1
  775. DO 120 I = 1, N
  776. D1( I ) = REAL( A( I, I ) )
  777. 120 CONTINUE
  778. DO 130 I = 1, N - 1
  779. D2( I ) = REAL( A( I+1, I ) )
  780. 130 CONTINUE
  781. SRNAMT = 'SSTEV'
  782. CALL SSTEV( 'V', N, D1, D2, Z, LDU, WORK, IINFO )
  783. IF( IINFO.NE.0 ) THEN
  784. WRITE( NOUNIT, FMT = 9999 )'SSTEV(V)', IINFO, N,
  785. $ JTYPE, IOLDSD
  786. INFO = ABS( IINFO )
  787. IF( IINFO.LT.0 ) THEN
  788. RETURN
  789. ELSE
  790. RESULT( 1 ) = ULPINV
  791. RESULT( 2 ) = ULPINV
  792. RESULT( 3 ) = ULPINV
  793. GO TO 180
  794. END IF
  795. END IF
  796. *
  797. * Do tests 1 and 2.
  798. *
  799. DO 140 I = 1, N
  800. D3( I ) = REAL( A( I, I ) )
  801. 140 CONTINUE
  802. DO 150 I = 1, N - 1
  803. D4( I ) = REAL( A( I+1, I ) )
  804. 150 CONTINUE
  805. CALL SSTT21( N, 0, D3, D4, D1, D2, Z, LDU, WORK,
  806. $ RESULT( 1 ) )
  807. *
  808. NTEST = 3
  809. DO 160 I = 1, N - 1
  810. D4( I ) = REAL( A( I+1, I ) )
  811. 160 CONTINUE
  812. SRNAMT = 'SSTEV'
  813. CALL SSTEV( 'N', N, D3, D4, Z, LDU, WORK, IINFO )
  814. IF( IINFO.NE.0 ) THEN
  815. WRITE( NOUNIT, FMT = 9999 )'SSTEV(N)', IINFO, N,
  816. $ JTYPE, IOLDSD
  817. INFO = ABS( IINFO )
  818. IF( IINFO.LT.0 ) THEN
  819. RETURN
  820. ELSE
  821. RESULT( 3 ) = ULPINV
  822. GO TO 180
  823. END IF
  824. END IF
  825. *
  826. * Do test 3.
  827. *
  828. TEMP1 = ZERO
  829. TEMP2 = ZERO
  830. DO 170 J = 1, N
  831. TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
  832. TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
  833. 170 CONTINUE
  834. RESULT( 3 ) = TEMP2 / MAX( UNFL,
  835. $ ULP*MAX( TEMP1, TEMP2 ) )
  836. *
  837. 180 CONTINUE
  838. *
  839. NTEST = 4
  840. DO 190 I = 1, N
  841. EVEIGS( I ) = D3( I )
  842. D1( I ) = REAL( A( I, I ) )
  843. 190 CONTINUE
  844. DO 200 I = 1, N - 1
  845. D2( I ) = REAL( A( I+1, I ) )
  846. 200 CONTINUE
  847. SRNAMT = 'SSTEVX'
  848. CALL SSTEVX( 'V', 'A', N, D1, D2, VL, VU, IL, IU, ABSTOL,
  849. $ M, WA1, Z, LDU, WORK, IWORK, IWORK( 5*N+1 ),
  850. $ IINFO )
  851. IF( IINFO.NE.0 ) THEN
  852. WRITE( NOUNIT, FMT = 9999 )'SSTEVX(V,A)', IINFO, N,
  853. $ JTYPE, IOLDSD
  854. INFO = ABS( IINFO )
  855. IF( IINFO.LT.0 ) THEN
  856. RETURN
  857. ELSE
  858. RESULT( 4 ) = ULPINV
  859. RESULT( 5 ) = ULPINV
  860. RESULT( 6 ) = ULPINV
  861. GO TO 250
  862. END IF
  863. END IF
  864. IF( N.GT.0 ) THEN
  865. TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
  866. ELSE
  867. TEMP3 = ZERO
  868. END IF
  869. *
  870. * Do tests 4 and 5.
  871. *
  872. DO 210 I = 1, N
  873. D3( I ) = REAL( A( I, I ) )
  874. 210 CONTINUE
  875. DO 220 I = 1, N - 1
  876. D4( I ) = REAL( A( I+1, I ) )
  877. 220 CONTINUE
  878. CALL SSTT21( N, 0, D3, D4, WA1, D2, Z, LDU, WORK,
  879. $ RESULT( 4 ) )
  880. *
  881. NTEST = 6
  882. DO 230 I = 1, N - 1
  883. D4( I ) = REAL( A( I+1, I ) )
  884. 230 CONTINUE
  885. SRNAMT = 'SSTEVX'
  886. CALL SSTEVX( 'N', 'A', N, D3, D4, VL, VU, IL, IU, ABSTOL,
  887. $ M2, WA2, Z, LDU, WORK, IWORK,
  888. $ IWORK( 5*N+1 ), IINFO )
  889. IF( IINFO.NE.0 ) THEN
  890. WRITE( NOUNIT, FMT = 9999 )'SSTEVX(N,A)', IINFO, N,
  891. $ JTYPE, IOLDSD
  892. INFO = ABS( IINFO )
  893. IF( IINFO.LT.0 ) THEN
  894. RETURN
  895. ELSE
  896. RESULT( 6 ) = ULPINV
  897. GO TO 250
  898. END IF
  899. END IF
  900. *
  901. * Do test 6.
  902. *
  903. TEMP1 = ZERO
  904. TEMP2 = ZERO
  905. DO 240 J = 1, N
  906. TEMP1 = MAX( TEMP1, ABS( WA2( J ) ),
  907. $ ABS( EVEIGS( J ) ) )
  908. TEMP2 = MAX( TEMP2, ABS( WA2( J )-EVEIGS( J ) ) )
  909. 240 CONTINUE
  910. RESULT( 6 ) = TEMP2 / MAX( UNFL,
  911. $ ULP*MAX( TEMP1, TEMP2 ) )
  912. *
  913. 250 CONTINUE
  914. *
  915. NTEST = 7
  916. DO 260 I = 1, N
  917. D1( I ) = REAL( A( I, I ) )
  918. 260 CONTINUE
  919. DO 270 I = 1, N - 1
  920. D2( I ) = REAL( A( I+1, I ) )
  921. 270 CONTINUE
  922. SRNAMT = 'SSTEVR'
  923. CALL SSTEVR( 'V', 'A', N, D1, D2, VL, VU, IL, IU, ABSTOL,
  924. $ M, WA1, Z, LDU, IWORK, WORK, LWORK,
  925. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  926. IF( IINFO.NE.0 ) THEN
  927. WRITE( NOUNIT, FMT = 9999 )'SSTEVR(V,A)', IINFO, N,
  928. $ JTYPE, IOLDSD
  929. INFO = ABS( IINFO )
  930. IF( IINFO.LT.0 ) THEN
  931. RETURN
  932. ELSE
  933. RESULT( 7 ) = ULPINV
  934. RESULT( 8 ) = ULPINV
  935. GO TO 320
  936. END IF
  937. END IF
  938. IF( N.GT.0 ) THEN
  939. TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
  940. ELSE
  941. TEMP3 = ZERO
  942. END IF
  943. *
  944. * Do tests 7 and 8.
  945. *
  946. DO 280 I = 1, N
  947. D3( I ) = REAL( A( I, I ) )
  948. 280 CONTINUE
  949. DO 290 I = 1, N - 1
  950. D4( I ) = REAL( A( I+1, I ) )
  951. 290 CONTINUE
  952. CALL SSTT21( N, 0, D3, D4, WA1, D2, Z, LDU, WORK,
  953. $ RESULT( 7 ) )
  954. *
  955. NTEST = 9
  956. DO 300 I = 1, N - 1
  957. D4( I ) = REAL( A( I+1, I ) )
  958. 300 CONTINUE
  959. SRNAMT = 'SSTEVR'
  960. CALL SSTEVR( 'N', 'A', N, D3, D4, VL, VU, IL, IU, ABSTOL,
  961. $ M2, WA2, Z, LDU, IWORK, WORK, LWORK,
  962. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  963. IF( IINFO.NE.0 ) THEN
  964. WRITE( NOUNIT, FMT = 9999 )'SSTEVR(N,A)', IINFO, N,
  965. $ JTYPE, IOLDSD
  966. INFO = ABS( IINFO )
  967. IF( IINFO.LT.0 ) THEN
  968. RETURN
  969. ELSE
  970. RESULT( 9 ) = ULPINV
  971. GO TO 320
  972. END IF
  973. END IF
  974. *
  975. * Do test 9.
  976. *
  977. TEMP1 = ZERO
  978. TEMP2 = ZERO
  979. DO 310 J = 1, N
  980. TEMP1 = MAX( TEMP1, ABS( WA2( J ) ),
  981. $ ABS( EVEIGS( J ) ) )
  982. TEMP2 = MAX( TEMP2, ABS( WA2( J )-EVEIGS( J ) ) )
  983. 310 CONTINUE
  984. RESULT( 9 ) = TEMP2 / MAX( UNFL,
  985. $ ULP*MAX( TEMP1, TEMP2 ) )
  986. *
  987. 320 CONTINUE
  988. *
  989. *
  990. NTEST = 10
  991. DO 330 I = 1, N
  992. D1( I ) = REAL( A( I, I ) )
  993. 330 CONTINUE
  994. DO 340 I = 1, N - 1
  995. D2( I ) = REAL( A( I+1, I ) )
  996. 340 CONTINUE
  997. SRNAMT = 'SSTEVX'
  998. CALL SSTEVX( 'V', 'I', N, D1, D2, VL, VU, IL, IU, ABSTOL,
  999. $ M2, WA2, Z, LDU, WORK, IWORK,
  1000. $ IWORK( 5*N+1 ), IINFO )
  1001. IF( IINFO.NE.0 ) THEN
  1002. WRITE( NOUNIT, FMT = 9999 )'SSTEVX(V,I)', IINFO, N,
  1003. $ JTYPE, IOLDSD
  1004. INFO = ABS( IINFO )
  1005. IF( IINFO.LT.0 ) THEN
  1006. RETURN
  1007. ELSE
  1008. RESULT( 10 ) = ULPINV
  1009. RESULT( 11 ) = ULPINV
  1010. RESULT( 12 ) = ULPINV
  1011. GO TO 380
  1012. END IF
  1013. END IF
  1014. *
  1015. * Do tests 10 and 11.
  1016. *
  1017. DO 350 I = 1, N
  1018. D3( I ) = REAL( A( I, I ) )
  1019. 350 CONTINUE
  1020. DO 360 I = 1, N - 1
  1021. D4( I ) = REAL( A( I+1, I ) )
  1022. 360 CONTINUE
  1023. CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
  1024. $ MAX( 1, M2 ), RESULT( 10 ) )
  1025. *
  1026. *
  1027. NTEST = 12
  1028. DO 370 I = 1, N - 1
  1029. D4( I ) = REAL( A( I+1, I ) )
  1030. 370 CONTINUE
  1031. SRNAMT = 'SSTEVX'
  1032. CALL SSTEVX( 'N', 'I', N, D3, D4, VL, VU, IL, IU, ABSTOL,
  1033. $ M3, WA3, Z, LDU, WORK, IWORK,
  1034. $ IWORK( 5*N+1 ), IINFO )
  1035. IF( IINFO.NE.0 ) THEN
  1036. WRITE( NOUNIT, FMT = 9999 )'SSTEVX(N,I)', IINFO, N,
  1037. $ JTYPE, IOLDSD
  1038. INFO = ABS( IINFO )
  1039. IF( IINFO.LT.0 ) THEN
  1040. RETURN
  1041. ELSE
  1042. RESULT( 12 ) = ULPINV
  1043. GO TO 380
  1044. END IF
  1045. END IF
  1046. *
  1047. * Do test 12.
  1048. *
  1049. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  1050. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  1051. RESULT( 12 ) = ( TEMP1+TEMP2 ) / MAX( UNFL, ULP*TEMP3 )
  1052. *
  1053. 380 CONTINUE
  1054. *
  1055. NTEST = 12
  1056. IF( N.GT.0 ) THEN
  1057. IF( IL.NE.1 ) THEN
  1058. VL = WA1( IL ) - MAX( HALF*
  1059. $ ( WA1( IL )-WA1( IL-1 ) ), TEN*ULP*TEMP3,
  1060. $ TEN*RTUNFL )
  1061. ELSE
  1062. VL = WA1( 1 ) - MAX( HALF*( WA1( N )-WA1( 1 ) ),
  1063. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1064. END IF
  1065. IF( IU.NE.N ) THEN
  1066. VU = WA1( IU ) + MAX( HALF*
  1067. $ ( WA1( IU+1 )-WA1( IU ) ), TEN*ULP*TEMP3,
  1068. $ TEN*RTUNFL )
  1069. ELSE
  1070. VU = WA1( N ) + MAX( HALF*( WA1( N )-WA1( 1 ) ),
  1071. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1072. END IF
  1073. ELSE
  1074. VL = ZERO
  1075. VU = ONE
  1076. END IF
  1077. *
  1078. DO 390 I = 1, N
  1079. D1( I ) = REAL( A( I, I ) )
  1080. 390 CONTINUE
  1081. DO 400 I = 1, N - 1
  1082. D2( I ) = REAL( A( I+1, I ) )
  1083. 400 CONTINUE
  1084. SRNAMT = 'SSTEVX'
  1085. CALL SSTEVX( 'V', 'V', N, D1, D2, VL, VU, IL, IU, ABSTOL,
  1086. $ M2, WA2, Z, LDU, WORK, IWORK,
  1087. $ IWORK( 5*N+1 ), IINFO )
  1088. IF( IINFO.NE.0 ) THEN
  1089. WRITE( NOUNIT, FMT = 9999 )'SSTEVX(V,V)', IINFO, N,
  1090. $ JTYPE, IOLDSD
  1091. INFO = ABS( IINFO )
  1092. IF( IINFO.LT.0 ) THEN
  1093. RETURN
  1094. ELSE
  1095. RESULT( 13 ) = ULPINV
  1096. RESULT( 14 ) = ULPINV
  1097. RESULT( 15 ) = ULPINV
  1098. GO TO 440
  1099. END IF
  1100. END IF
  1101. *
  1102. IF( M2.EQ.0 .AND. N.GT.0 ) THEN
  1103. RESULT( 13 ) = ULPINV
  1104. RESULT( 14 ) = ULPINV
  1105. RESULT( 15 ) = ULPINV
  1106. GO TO 440
  1107. END IF
  1108. *
  1109. * Do tests 13 and 14.
  1110. *
  1111. DO 410 I = 1, N
  1112. D3( I ) = REAL( A( I, I ) )
  1113. 410 CONTINUE
  1114. DO 420 I = 1, N - 1
  1115. D4( I ) = REAL( A( I+1, I ) )
  1116. 420 CONTINUE
  1117. CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
  1118. $ MAX( 1, M2 ), RESULT( 13 ) )
  1119. *
  1120. NTEST = 15
  1121. DO 430 I = 1, N - 1
  1122. D4( I ) = REAL( A( I+1, I ) )
  1123. 430 CONTINUE
  1124. SRNAMT = 'SSTEVX'
  1125. CALL SSTEVX( 'N', 'V', N, D3, D4, VL, VU, IL, IU, ABSTOL,
  1126. $ M3, WA3, Z, LDU, WORK, IWORK,
  1127. $ IWORK( 5*N+1 ), IINFO )
  1128. IF( IINFO.NE.0 ) THEN
  1129. WRITE( NOUNIT, FMT = 9999 )'SSTEVX(N,V)', IINFO, N,
  1130. $ JTYPE, IOLDSD
  1131. INFO = ABS( IINFO )
  1132. IF( IINFO.LT.0 ) THEN
  1133. RETURN
  1134. ELSE
  1135. RESULT( 15 ) = ULPINV
  1136. GO TO 440
  1137. END IF
  1138. END IF
  1139. *
  1140. * Do test 15.
  1141. *
  1142. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  1143. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  1144. RESULT( 15 ) = ( TEMP1+TEMP2 ) / MAX( UNFL, TEMP3*ULP )
  1145. *
  1146. 440 CONTINUE
  1147. *
  1148. NTEST = 16
  1149. DO 450 I = 1, N
  1150. D1( I ) = REAL( A( I, I ) )
  1151. 450 CONTINUE
  1152. DO 460 I = 1, N - 1
  1153. D2( I ) = REAL( A( I+1, I ) )
  1154. 460 CONTINUE
  1155. SRNAMT = 'SSTEVD'
  1156. CALL SSTEVD( 'V', N, D1, D2, Z, LDU, WORK, LWEDC, IWORK,
  1157. $ LIWEDC, IINFO )
  1158. IF( IINFO.NE.0 ) THEN
  1159. WRITE( NOUNIT, FMT = 9999 )'SSTEVD(V)', IINFO, N,
  1160. $ JTYPE, IOLDSD
  1161. INFO = ABS( IINFO )
  1162. IF( IINFO.LT.0 ) THEN
  1163. RETURN
  1164. ELSE
  1165. RESULT( 16 ) = ULPINV
  1166. RESULT( 17 ) = ULPINV
  1167. RESULT( 18 ) = ULPINV
  1168. GO TO 510
  1169. END IF
  1170. END IF
  1171. *
  1172. * Do tests 16 and 17.
  1173. *
  1174. DO 470 I = 1, N
  1175. D3( I ) = REAL( A( I, I ) )
  1176. 470 CONTINUE
  1177. DO 480 I = 1, N - 1
  1178. D4( I ) = REAL( A( I+1, I ) )
  1179. 480 CONTINUE
  1180. CALL SSTT21( N, 0, D3, D4, D1, D2, Z, LDU, WORK,
  1181. $ RESULT( 16 ) )
  1182. *
  1183. NTEST = 18
  1184. DO 490 I = 1, N - 1
  1185. D4( I ) = REAL( A( I+1, I ) )
  1186. 490 CONTINUE
  1187. SRNAMT = 'SSTEVD'
  1188. CALL SSTEVD( 'N', N, D3, D4, Z, LDU, WORK, LWEDC, IWORK,
  1189. $ LIWEDC, IINFO )
  1190. IF( IINFO.NE.0 ) THEN
  1191. WRITE( NOUNIT, FMT = 9999 )'SSTEVD(N)', IINFO, N,
  1192. $ JTYPE, IOLDSD
  1193. INFO = ABS( IINFO )
  1194. IF( IINFO.LT.0 ) THEN
  1195. RETURN
  1196. ELSE
  1197. RESULT( 18 ) = ULPINV
  1198. GO TO 510
  1199. END IF
  1200. END IF
  1201. *
  1202. * Do test 18.
  1203. *
  1204. TEMP1 = ZERO
  1205. TEMP2 = ZERO
  1206. DO 500 J = 1, N
  1207. TEMP1 = MAX( TEMP1, ABS( EVEIGS( J ) ),
  1208. $ ABS( D3( J ) ) )
  1209. TEMP2 = MAX( TEMP2, ABS( EVEIGS( J )-D3( J ) ) )
  1210. 500 CONTINUE
  1211. RESULT( 18 ) = TEMP2 / MAX( UNFL,
  1212. $ ULP*MAX( TEMP1, TEMP2 ) )
  1213. *
  1214. 510 CONTINUE
  1215. *
  1216. NTEST = 19
  1217. DO 520 I = 1, N
  1218. D1( I ) = REAL( A( I, I ) )
  1219. 520 CONTINUE
  1220. DO 530 I = 1, N - 1
  1221. D2( I ) = REAL( A( I+1, I ) )
  1222. 530 CONTINUE
  1223. SRNAMT = 'SSTEVR'
  1224. CALL SSTEVR( 'V', 'I', N, D1, D2, VL, VU, IL, IU, ABSTOL,
  1225. $ M2, WA2, Z, LDU, IWORK, WORK, LWORK,
  1226. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  1227. IF( IINFO.NE.0 ) THEN
  1228. WRITE( NOUNIT, FMT = 9999 )'SSTEVR(V,I)', IINFO, N,
  1229. $ JTYPE, IOLDSD
  1230. INFO = ABS( IINFO )
  1231. IF( IINFO.LT.0 ) THEN
  1232. RETURN
  1233. ELSE
  1234. RESULT( 19 ) = ULPINV
  1235. RESULT( 20 ) = ULPINV
  1236. RESULT( 21 ) = ULPINV
  1237. GO TO 570
  1238. END IF
  1239. END IF
  1240. *
  1241. * DO tests 19 and 20.
  1242. *
  1243. DO 540 I = 1, N
  1244. D3( I ) = REAL( A( I, I ) )
  1245. 540 CONTINUE
  1246. DO 550 I = 1, N - 1
  1247. D4( I ) = REAL( A( I+1, I ) )
  1248. 550 CONTINUE
  1249. CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
  1250. $ MAX( 1, M2 ), RESULT( 19 ) )
  1251. *
  1252. *
  1253. NTEST = 21
  1254. DO 560 I = 1, N - 1
  1255. D4( I ) = REAL( A( I+1, I ) )
  1256. 560 CONTINUE
  1257. SRNAMT = 'SSTEVR'
  1258. CALL SSTEVR( 'N', 'I', N, D3, D4, VL, VU, IL, IU, ABSTOL,
  1259. $ M3, WA3, Z, LDU, IWORK, WORK, LWORK,
  1260. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  1261. IF( IINFO.NE.0 ) THEN
  1262. WRITE( NOUNIT, FMT = 9999 )'SSTEVR(N,I)', IINFO, N,
  1263. $ JTYPE, IOLDSD
  1264. INFO = ABS( IINFO )
  1265. IF( IINFO.LT.0 ) THEN
  1266. RETURN
  1267. ELSE
  1268. RESULT( 21 ) = ULPINV
  1269. GO TO 570
  1270. END IF
  1271. END IF
  1272. *
  1273. * Do test 21.
  1274. *
  1275. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  1276. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  1277. RESULT( 21 ) = ( TEMP1+TEMP2 ) / MAX( UNFL, ULP*TEMP3 )
  1278. *
  1279. 570 CONTINUE
  1280. *
  1281. NTEST = 21
  1282. IF( N.GT.0 ) THEN
  1283. IF( IL.NE.1 ) THEN
  1284. VL = WA1( IL ) - MAX( HALF*
  1285. $ ( WA1( IL )-WA1( IL-1 ) ), TEN*ULP*TEMP3,
  1286. $ TEN*RTUNFL )
  1287. ELSE
  1288. VL = WA1( 1 ) - MAX( HALF*( WA1( N )-WA1( 1 ) ),
  1289. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1290. END IF
  1291. IF( IU.NE.N ) THEN
  1292. VU = WA1( IU ) + MAX( HALF*
  1293. $ ( WA1( IU+1 )-WA1( IU ) ), TEN*ULP*TEMP3,
  1294. $ TEN*RTUNFL )
  1295. ELSE
  1296. VU = WA1( N ) + MAX( HALF*( WA1( N )-WA1( 1 ) ),
  1297. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1298. END IF
  1299. ELSE
  1300. VL = ZERO
  1301. VU = ONE
  1302. END IF
  1303. *
  1304. DO 580 I = 1, N
  1305. D1( I ) = REAL( A( I, I ) )
  1306. 580 CONTINUE
  1307. DO 590 I = 1, N - 1
  1308. D2( I ) = REAL( A( I+1, I ) )
  1309. 590 CONTINUE
  1310. SRNAMT = 'SSTEVR'
  1311. CALL SSTEVR( 'V', 'V', N, D1, D2, VL, VU, IL, IU, ABSTOL,
  1312. $ M2, WA2, Z, LDU, IWORK, WORK, LWORK,
  1313. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  1314. IF( IINFO.NE.0 ) THEN
  1315. WRITE( NOUNIT, FMT = 9999 )'SSTEVR(V,V)', IINFO, N,
  1316. $ JTYPE, IOLDSD
  1317. INFO = ABS( IINFO )
  1318. IF( IINFO.LT.0 ) THEN
  1319. RETURN
  1320. ELSE
  1321. RESULT( 22 ) = ULPINV
  1322. RESULT( 23 ) = ULPINV
  1323. RESULT( 24 ) = ULPINV
  1324. GO TO 630
  1325. END IF
  1326. END IF
  1327. *
  1328. IF( M2.EQ.0 .AND. N.GT.0 ) THEN
  1329. RESULT( 22 ) = ULPINV
  1330. RESULT( 23 ) = ULPINV
  1331. RESULT( 24 ) = ULPINV
  1332. GO TO 630
  1333. END IF
  1334. *
  1335. * Do tests 22 and 23.
  1336. *
  1337. DO 600 I = 1, N
  1338. D3( I ) = REAL( A( I, I ) )
  1339. 600 CONTINUE
  1340. DO 610 I = 1, N - 1
  1341. D4( I ) = REAL( A( I+1, I ) )
  1342. 610 CONTINUE
  1343. CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
  1344. $ MAX( 1, M2 ), RESULT( 22 ) )
  1345. *
  1346. NTEST = 24
  1347. DO 620 I = 1, N - 1
  1348. D4( I ) = REAL( A( I+1, I ) )
  1349. 620 CONTINUE
  1350. SRNAMT = 'SSTEVR'
  1351. CALL SSTEVR( 'N', 'V', N, D3, D4, VL, VU, IL, IU, ABSTOL,
  1352. $ M3, WA3, Z, LDU, IWORK, WORK, LWORK,
  1353. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  1354. IF( IINFO.NE.0 ) THEN
  1355. WRITE( NOUNIT, FMT = 9999 )'SSTEVR(N,V)', IINFO, N,
  1356. $ JTYPE, IOLDSD
  1357. INFO = ABS( IINFO )
  1358. IF( IINFO.LT.0 ) THEN
  1359. RETURN
  1360. ELSE
  1361. RESULT( 24 ) = ULPINV
  1362. GO TO 630
  1363. END IF
  1364. END IF
  1365. *
  1366. * Do test 24.
  1367. *
  1368. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  1369. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  1370. RESULT( 24 ) = ( TEMP1+TEMP2 ) / MAX( UNFL, TEMP3*ULP )
  1371. *
  1372. 630 CONTINUE
  1373. *
  1374. *
  1375. *
  1376. ELSE
  1377. *
  1378. DO 640 I = 1, 24
  1379. RESULT( I ) = ZERO
  1380. 640 CONTINUE
  1381. NTEST = 24
  1382. END IF
  1383. *
  1384. * Perform remaining tests storing upper or lower triangular
  1385. * part of matrix.
  1386. *
  1387. DO 1720 IUPLO = 0, 1
  1388. IF( IUPLO.EQ.0 ) THEN
  1389. UPLO = 'L'
  1390. ELSE
  1391. UPLO = 'U'
  1392. END IF
  1393. *
  1394. * 4) Call SSYEV and SSYEVX.
  1395. *
  1396. CALL SLACPY( ' ', N, N, A, LDA, V, LDU )
  1397. *
  1398. NTEST = NTEST + 1
  1399. SRNAMT = 'SSYEV'
  1400. CALL SSYEV( 'V', UPLO, N, A, LDU, D1, WORK, LWORK,
  1401. $ IINFO )
  1402. IF( IINFO.NE.0 ) THEN
  1403. WRITE( NOUNIT, FMT = 9999 )'SSYEV(V,' // UPLO // ')',
  1404. $ IINFO, N, JTYPE, IOLDSD
  1405. INFO = ABS( IINFO )
  1406. IF( IINFO.LT.0 ) THEN
  1407. RETURN
  1408. ELSE
  1409. RESULT( NTEST ) = ULPINV
  1410. RESULT( NTEST+1 ) = ULPINV
  1411. RESULT( NTEST+2 ) = ULPINV
  1412. GO TO 660
  1413. END IF
  1414. END IF
  1415. *
  1416. * Do tests 25 and 26 (or +54)
  1417. *
  1418. CALL SSYT21( 1, UPLO, N, 0, V, LDU, D1, D2, A, LDU, Z,
  1419. $ LDU, TAU, WORK, RESULT( NTEST ) )
  1420. *
  1421. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  1422. *
  1423. NTEST = NTEST + 2
  1424. SRNAMT = 'SSYEV'
  1425. CALL SSYEV( 'N', UPLO, N, A, LDU, D3, WORK, LWORK,
  1426. $ IINFO )
  1427. IF( IINFO.NE.0 ) THEN
  1428. WRITE( NOUNIT, FMT = 9999 )'SSYEV(N,' // UPLO // ')',
  1429. $ IINFO, N, JTYPE, IOLDSD
  1430. INFO = ABS( IINFO )
  1431. IF( IINFO.LT.0 ) THEN
  1432. RETURN
  1433. ELSE
  1434. RESULT( NTEST ) = ULPINV
  1435. GO TO 660
  1436. END IF
  1437. END IF
  1438. *
  1439. * Do test 27 (or +54)
  1440. *
  1441. TEMP1 = ZERO
  1442. TEMP2 = ZERO
  1443. DO 650 J = 1, N
  1444. TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
  1445. TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
  1446. 650 CONTINUE
  1447. RESULT( NTEST ) = TEMP2 / MAX( UNFL,
  1448. $ ULP*MAX( TEMP1, TEMP2 ) )
  1449. *
  1450. 660 CONTINUE
  1451. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  1452. *
  1453. NTEST = NTEST + 1
  1454. *
  1455. IF( N.GT.0 ) THEN
  1456. TEMP3 = MAX( ABS( D1( 1 ) ), ABS( D1( N ) ) )
  1457. IF( IL.NE.1 ) THEN
  1458. VL = D1( IL ) - MAX( HALF*( D1( IL )-D1( IL-1 ) ),
  1459. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1460. ELSE IF( N.GT.0 ) THEN
  1461. VL = D1( 1 ) - MAX( HALF*( D1( N )-D1( 1 ) ),
  1462. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1463. END IF
  1464. IF( IU.NE.N ) THEN
  1465. VU = D1( IU ) + MAX( HALF*( D1( IU+1 )-D1( IU ) ),
  1466. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1467. ELSE IF( N.GT.0 ) THEN
  1468. VU = D1( N ) + MAX( HALF*( D1( N )-D1( 1 ) ),
  1469. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1470. END IF
  1471. ELSE
  1472. TEMP3 = ZERO
  1473. VL = ZERO
  1474. VU = ONE
  1475. END IF
  1476. *
  1477. SRNAMT = 'SSYEVX'
  1478. CALL SSYEVX( 'V', 'A', UPLO, N, A, LDU, VL, VU, IL, IU,
  1479. $ ABSTOL, M, WA1, Z, LDU, WORK, LWORK, IWORK,
  1480. $ IWORK( 5*N+1 ), IINFO )
  1481. IF( IINFO.NE.0 ) THEN
  1482. WRITE( NOUNIT, FMT = 9999 )'SSYEVX(V,A,' // UPLO //
  1483. $ ')', IINFO, N, JTYPE, IOLDSD
  1484. INFO = ABS( IINFO )
  1485. IF( IINFO.LT.0 ) THEN
  1486. RETURN
  1487. ELSE
  1488. RESULT( NTEST ) = ULPINV
  1489. RESULT( NTEST+1 ) = ULPINV
  1490. RESULT( NTEST+2 ) = ULPINV
  1491. GO TO 680
  1492. END IF
  1493. END IF
  1494. *
  1495. * Do tests 28 and 29 (or +54)
  1496. *
  1497. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  1498. *
  1499. CALL SSYT21( 1, UPLO, N, 0, A, LDU, D1, D2, Z, LDU, V,
  1500. $ LDU, TAU, WORK, RESULT( NTEST ) )
  1501. *
  1502. NTEST = NTEST + 2
  1503. SRNAMT = 'SSYEVX'
  1504. CALL SSYEVX( 'N', 'A', UPLO, N, A, LDU, VL, VU, IL, IU,
  1505. $ ABSTOL, M2, WA2, Z, LDU, WORK, LWORK, IWORK,
  1506. $ IWORK( 5*N+1 ), IINFO )
  1507. IF( IINFO.NE.0 ) THEN
  1508. WRITE( NOUNIT, FMT = 9999 )'SSYEVX(N,A,' // UPLO //
  1509. $ ')', IINFO, N, JTYPE, IOLDSD
  1510. INFO = ABS( IINFO )
  1511. IF( IINFO.LT.0 ) THEN
  1512. RETURN
  1513. ELSE
  1514. RESULT( NTEST ) = ULPINV
  1515. GO TO 680
  1516. END IF
  1517. END IF
  1518. *
  1519. * Do test 30 (or +54)
  1520. *
  1521. TEMP1 = ZERO
  1522. TEMP2 = ZERO
  1523. DO 670 J = 1, N
  1524. TEMP1 = MAX( TEMP1, ABS( WA1( J ) ), ABS( WA2( J ) ) )
  1525. TEMP2 = MAX( TEMP2, ABS( WA1( J )-WA2( J ) ) )
  1526. 670 CONTINUE
  1527. RESULT( NTEST ) = TEMP2 / MAX( UNFL,
  1528. $ ULP*MAX( TEMP1, TEMP2 ) )
  1529. *
  1530. 680 CONTINUE
  1531. *
  1532. NTEST = NTEST + 1
  1533. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  1534. SRNAMT = 'SSYEVX'
  1535. CALL SSYEVX( 'V', 'I', UPLO, N, A, LDU, VL, VU, IL, IU,
  1536. $ ABSTOL, M2, WA2, Z, LDU, WORK, LWORK, IWORK,
  1537. $ IWORK( 5*N+1 ), IINFO )
  1538. IF( IINFO.NE.0 ) THEN
  1539. WRITE( NOUNIT, FMT = 9999 )'SSYEVX(V,I,' // UPLO //
  1540. $ ')', IINFO, N, JTYPE, IOLDSD
  1541. INFO = ABS( IINFO )
  1542. IF( IINFO.LT.0 ) THEN
  1543. RETURN
  1544. ELSE
  1545. RESULT( NTEST ) = ULPINV
  1546. RESULT( NTEST+1 ) = ULPINV
  1547. RESULT( NTEST+2 ) = ULPINV
  1548. GO TO 690
  1549. END IF
  1550. END IF
  1551. *
  1552. * Do tests 31 and 32 (or +54)
  1553. *
  1554. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  1555. *
  1556. CALL SSYT22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
  1557. $ V, LDU, TAU, WORK, RESULT( NTEST ) )
  1558. *
  1559. NTEST = NTEST + 2
  1560. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  1561. SRNAMT = 'SSYEVX'
  1562. CALL SSYEVX( 'N', 'I', UPLO, N, A, LDU, VL, VU, IL, IU,
  1563. $ ABSTOL, M3, WA3, Z, LDU, WORK, LWORK, IWORK,
  1564. $ IWORK( 5*N+1 ), IINFO )
  1565. IF( IINFO.NE.0 ) THEN
  1566. WRITE( NOUNIT, FMT = 9999 )'SSYEVX(N,I,' // UPLO //
  1567. $ ')', IINFO, N, JTYPE, IOLDSD
  1568. INFO = ABS( IINFO )
  1569. IF( IINFO.LT.0 ) THEN
  1570. RETURN
  1571. ELSE
  1572. RESULT( NTEST ) = ULPINV
  1573. GO TO 690
  1574. END IF
  1575. END IF
  1576. *
  1577. * Do test 33 (or +54)
  1578. *
  1579. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  1580. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  1581. RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
  1582. $ MAX( UNFL, ULP*TEMP3 )
  1583. 690 CONTINUE
  1584. *
  1585. NTEST = NTEST + 1
  1586. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  1587. SRNAMT = 'SSYEVX'
  1588. CALL SSYEVX( 'V', 'V', UPLO, N, A, LDU, VL, VU, IL, IU,
  1589. $ ABSTOL, M2, WA2, Z, LDU, WORK, LWORK, IWORK,
  1590. $ IWORK( 5*N+1 ), IINFO )
  1591. IF( IINFO.NE.0 ) THEN
  1592. WRITE( NOUNIT, FMT = 9999 )'SSYEVX(V,V,' // UPLO //
  1593. $ ')', IINFO, N, JTYPE, IOLDSD
  1594. INFO = ABS( IINFO )
  1595. IF( IINFO.LT.0 ) THEN
  1596. RETURN
  1597. ELSE
  1598. RESULT( NTEST ) = ULPINV
  1599. RESULT( NTEST+1 ) = ULPINV
  1600. RESULT( NTEST+2 ) = ULPINV
  1601. GO TO 700
  1602. END IF
  1603. END IF
  1604. *
  1605. * Do tests 34 and 35 (or +54)
  1606. *
  1607. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  1608. *
  1609. CALL SSYT22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
  1610. $ V, LDU, TAU, WORK, RESULT( NTEST ) )
  1611. *
  1612. NTEST = NTEST + 2
  1613. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  1614. SRNAMT = 'SSYEVX'
  1615. CALL SSYEVX( 'N', 'V', UPLO, N, A, LDU, VL, VU, IL, IU,
  1616. $ ABSTOL, M3, WA3, Z, LDU, WORK, LWORK, IWORK,
  1617. $ IWORK( 5*N+1 ), IINFO )
  1618. IF( IINFO.NE.0 ) THEN
  1619. WRITE( NOUNIT, FMT = 9999 )'SSYEVX(N,V,' // UPLO //
  1620. $ ')', IINFO, N, JTYPE, IOLDSD
  1621. INFO = ABS( IINFO )
  1622. IF( IINFO.LT.0 ) THEN
  1623. RETURN
  1624. ELSE
  1625. RESULT( NTEST ) = ULPINV
  1626. GO TO 700
  1627. END IF
  1628. END IF
  1629. *
  1630. IF( M3.EQ.0 .AND. N.GT.0 ) THEN
  1631. RESULT( NTEST ) = ULPINV
  1632. GO TO 700
  1633. END IF
  1634. *
  1635. * Do test 36 (or +54)
  1636. *
  1637. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  1638. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  1639. IF( N.GT.0 ) THEN
  1640. TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
  1641. ELSE
  1642. TEMP3 = ZERO
  1643. END IF
  1644. RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
  1645. $ MAX( UNFL, TEMP3*ULP )
  1646. *
  1647. 700 CONTINUE
  1648. *
  1649. * 5) Call SSPEV and SSPEVX.
  1650. *
  1651. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  1652. *
  1653. * Load array WORK with the upper or lower triangular
  1654. * part of the matrix in packed form.
  1655. *
  1656. IF( IUPLO.EQ.1 ) THEN
  1657. INDX = 1
  1658. DO 720 J = 1, N
  1659. DO 710 I = 1, J
  1660. WORK( INDX ) = A( I, J )
  1661. INDX = INDX + 1
  1662. 710 CONTINUE
  1663. 720 CONTINUE
  1664. ELSE
  1665. INDX = 1
  1666. DO 740 J = 1, N
  1667. DO 730 I = J, N
  1668. WORK( INDX ) = A( I, J )
  1669. INDX = INDX + 1
  1670. 730 CONTINUE
  1671. 740 CONTINUE
  1672. END IF
  1673. *
  1674. NTEST = NTEST + 1
  1675. SRNAMT = 'SSPEV'
  1676. CALL SSPEV( 'V', UPLO, N, WORK, D1, Z, LDU, V, IINFO )
  1677. IF( IINFO.NE.0 ) THEN
  1678. WRITE( NOUNIT, FMT = 9999 )'SSPEV(V,' // UPLO // ')',
  1679. $ IINFO, N, JTYPE, IOLDSD
  1680. INFO = ABS( IINFO )
  1681. IF( IINFO.LT.0 ) THEN
  1682. RETURN
  1683. ELSE
  1684. RESULT( NTEST ) = ULPINV
  1685. RESULT( NTEST+1 ) = ULPINV
  1686. RESULT( NTEST+2 ) = ULPINV
  1687. GO TO 800
  1688. END IF
  1689. END IF
  1690. *
  1691. * Do tests 37 and 38 (or +54)
  1692. *
  1693. CALL SSYT21( 1, UPLO, N, 0, A, LDA, D1, D2, Z, LDU, V,
  1694. $ LDU, TAU, WORK, RESULT( NTEST ) )
  1695. *
  1696. IF( IUPLO.EQ.1 ) THEN
  1697. INDX = 1
  1698. DO 760 J = 1, N
  1699. DO 750 I = 1, J
  1700. WORK( INDX ) = A( I, J )
  1701. INDX = INDX + 1
  1702. 750 CONTINUE
  1703. 760 CONTINUE
  1704. ELSE
  1705. INDX = 1
  1706. DO 780 J = 1, N
  1707. DO 770 I = J, N
  1708. WORK( INDX ) = A( I, J )
  1709. INDX = INDX + 1
  1710. 770 CONTINUE
  1711. 780 CONTINUE
  1712. END IF
  1713. *
  1714. NTEST = NTEST + 2
  1715. SRNAMT = 'SSPEV'
  1716. CALL SSPEV( 'N', UPLO, N, WORK, D3, Z, LDU, V, IINFO )
  1717. IF( IINFO.NE.0 ) THEN
  1718. WRITE( NOUNIT, FMT = 9999 )'SSPEV(N,' // UPLO // ')',
  1719. $ IINFO, N, JTYPE, IOLDSD
  1720. INFO = ABS( IINFO )
  1721. IF( IINFO.LT.0 ) THEN
  1722. RETURN
  1723. ELSE
  1724. RESULT( NTEST ) = ULPINV
  1725. GO TO 800
  1726. END IF
  1727. END IF
  1728. *
  1729. * Do test 39 (or +54)
  1730. *
  1731. TEMP1 = ZERO
  1732. TEMP2 = ZERO
  1733. DO 790 J = 1, N
  1734. TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
  1735. TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
  1736. 790 CONTINUE
  1737. RESULT( NTEST ) = TEMP2 / MAX( UNFL,
  1738. $ ULP*MAX( TEMP1, TEMP2 ) )
  1739. *
  1740. * Load array WORK with the upper or lower triangular part
  1741. * of the matrix in packed form.
  1742. *
  1743. 800 CONTINUE
  1744. IF( IUPLO.EQ.1 ) THEN
  1745. INDX = 1
  1746. DO 820 J = 1, N
  1747. DO 810 I = 1, J
  1748. WORK( INDX ) = A( I, J )
  1749. INDX = INDX + 1
  1750. 810 CONTINUE
  1751. 820 CONTINUE
  1752. ELSE
  1753. INDX = 1
  1754. DO 840 J = 1, N
  1755. DO 830 I = J, N
  1756. WORK( INDX ) = A( I, J )
  1757. INDX = INDX + 1
  1758. 830 CONTINUE
  1759. 840 CONTINUE
  1760. END IF
  1761. *
  1762. NTEST = NTEST + 1
  1763. *
  1764. IF( N.GT.0 ) THEN
  1765. TEMP3 = MAX( ABS( D1( 1 ) ), ABS( D1( N ) ) )
  1766. IF( IL.NE.1 ) THEN
  1767. VL = D1( IL ) - MAX( HALF*( D1( IL )-D1( IL-1 ) ),
  1768. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1769. ELSE IF( N.GT.0 ) THEN
  1770. VL = D1( 1 ) - MAX( HALF*( D1( N )-D1( 1 ) ),
  1771. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1772. END IF
  1773. IF( IU.NE.N ) THEN
  1774. VU = D1( IU ) + MAX( HALF*( D1( IU+1 )-D1( IU ) ),
  1775. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1776. ELSE IF( N.GT.0 ) THEN
  1777. VU = D1( N ) + MAX( HALF*( D1( N )-D1( 1 ) ),
  1778. $ TEN*ULP*TEMP3, TEN*RTUNFL )
  1779. END IF
  1780. ELSE
  1781. TEMP3 = ZERO
  1782. VL = ZERO
  1783. VU = ONE
  1784. END IF
  1785. *
  1786. SRNAMT = 'SSPEVX'
  1787. CALL SSPEVX( 'V', 'A', UPLO, N, WORK, VL, VU, IL, IU,
  1788. $ ABSTOL, M, WA1, Z, LDU, V, IWORK,
  1789. $ IWORK( 5*N+1 ), IINFO )
  1790. IF( IINFO.NE.0 ) THEN
  1791. WRITE( NOUNIT, FMT = 9999 )'SSPEVX(V,A,' // UPLO //
  1792. $ ')', IINFO, N, JTYPE, IOLDSD
  1793. INFO = ABS( IINFO )
  1794. IF( IINFO.LT.0 ) THEN
  1795. RETURN
  1796. ELSE
  1797. RESULT( NTEST ) = ULPINV
  1798. RESULT( NTEST+1 ) = ULPINV
  1799. RESULT( NTEST+2 ) = ULPINV
  1800. GO TO 900
  1801. END IF
  1802. END IF
  1803. *
  1804. * Do tests 40 and 41 (or +54)
  1805. *
  1806. CALL SSYT21( 1, UPLO, N, 0, A, LDU, WA1, D2, Z, LDU, V,
  1807. $ LDU, TAU, WORK, RESULT( NTEST ) )
  1808. *
  1809. NTEST = NTEST + 2
  1810. *
  1811. IF( IUPLO.EQ.1 ) THEN
  1812. INDX = 1
  1813. DO 860 J = 1, N
  1814. DO 850 I = 1, J
  1815. WORK( INDX ) = A( I, J )
  1816. INDX = INDX + 1
  1817. 850 CONTINUE
  1818. 860 CONTINUE
  1819. ELSE
  1820. INDX = 1
  1821. DO 880 J = 1, N
  1822. DO 870 I = J, N
  1823. WORK( INDX ) = A( I, J )
  1824. INDX = INDX + 1
  1825. 870 CONTINUE
  1826. 880 CONTINUE
  1827. END IF
  1828. *
  1829. SRNAMT = 'SSPEVX'
  1830. CALL SSPEVX( 'N', 'A', UPLO, N, WORK, VL, VU, IL, IU,
  1831. $ ABSTOL, M2, WA2, Z, LDU, V, IWORK,
  1832. $ IWORK( 5*N+1 ), IINFO )
  1833. IF( IINFO.NE.0 ) THEN
  1834. WRITE( NOUNIT, FMT = 9999 )'SSPEVX(N,A,' // UPLO //
  1835. $ ')', IINFO, N, JTYPE, IOLDSD
  1836. INFO = ABS( IINFO )
  1837. IF( IINFO.LT.0 ) THEN
  1838. RETURN
  1839. ELSE
  1840. RESULT( NTEST ) = ULPINV
  1841. GO TO 900
  1842. END IF
  1843. END IF
  1844. *
  1845. * Do test 42 (or +54)
  1846. *
  1847. TEMP1 = ZERO
  1848. TEMP2 = ZERO
  1849. DO 890 J = 1, N
  1850. TEMP1 = MAX( TEMP1, ABS( WA1( J ) ), ABS( WA2( J ) ) )
  1851. TEMP2 = MAX( TEMP2, ABS( WA1( J )-WA2( J ) ) )
  1852. 890 CONTINUE
  1853. RESULT( NTEST ) = TEMP2 / MAX( UNFL,
  1854. $ ULP*MAX( TEMP1, TEMP2 ) )
  1855. *
  1856. 900 CONTINUE
  1857. IF( IUPLO.EQ.1 ) THEN
  1858. INDX = 1
  1859. DO 920 J = 1, N
  1860. DO 910 I = 1, J
  1861. WORK( INDX ) = A( I, J )
  1862. INDX = INDX + 1
  1863. 910 CONTINUE
  1864. 920 CONTINUE
  1865. ELSE
  1866. INDX = 1
  1867. DO 940 J = 1, N
  1868. DO 930 I = J, N
  1869. WORK( INDX ) = A( I, J )
  1870. INDX = INDX + 1
  1871. 930 CONTINUE
  1872. 940 CONTINUE
  1873. END IF
  1874. *
  1875. NTEST = NTEST + 1
  1876. *
  1877. SRNAMT = 'SSPEVX'
  1878. CALL SSPEVX( 'V', 'I', UPLO, N, WORK, VL, VU, IL, IU,
  1879. $ ABSTOL, M2, WA2, Z, LDU, V, IWORK,
  1880. $ IWORK( 5*N+1 ), IINFO )
  1881. IF( IINFO.NE.0 ) THEN
  1882. WRITE( NOUNIT, FMT = 9999 )'SSPEVX(V,I,' // UPLO //
  1883. $ ')', IINFO, N, JTYPE, IOLDSD
  1884. INFO = ABS( IINFO )
  1885. IF( IINFO.LT.0 ) THEN
  1886. RETURN
  1887. ELSE
  1888. RESULT( NTEST ) = ULPINV
  1889. RESULT( NTEST+1 ) = ULPINV
  1890. RESULT( NTEST+2 ) = ULPINV
  1891. GO TO 990
  1892. END IF
  1893. END IF
  1894. *
  1895. * Do tests 43 and 44 (or +54)
  1896. *
  1897. CALL SSYT22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
  1898. $ V, LDU, TAU, WORK, RESULT( NTEST ) )
  1899. *
  1900. NTEST = NTEST + 2
  1901. *
  1902. IF( IUPLO.EQ.1 ) THEN
  1903. INDX = 1
  1904. DO 960 J = 1, N
  1905. DO 950 I = 1, J
  1906. WORK( INDX ) = A( I, J )
  1907. INDX = INDX + 1
  1908. 950 CONTINUE
  1909. 960 CONTINUE
  1910. ELSE
  1911. INDX = 1
  1912. DO 980 J = 1, N
  1913. DO 970 I = J, N
  1914. WORK( INDX ) = A( I, J )
  1915. INDX = INDX + 1
  1916. 970 CONTINUE
  1917. 980 CONTINUE
  1918. END IF
  1919. *
  1920. SRNAMT = 'SSPEVX'
  1921. CALL SSPEVX( 'N', 'I', UPLO, N, WORK, VL, VU, IL, IU,
  1922. $ ABSTOL, M3, WA3, Z, LDU, V, IWORK,
  1923. $ IWORK( 5*N+1 ), IINFO )
  1924. IF( IINFO.NE.0 ) THEN
  1925. WRITE( NOUNIT, FMT = 9999 )'SSPEVX(N,I,' // UPLO //
  1926. $ ')', IINFO, N, JTYPE, IOLDSD
  1927. INFO = ABS( IINFO )
  1928. IF( IINFO.LT.0 ) THEN
  1929. RETURN
  1930. ELSE
  1931. RESULT( NTEST ) = ULPINV
  1932. GO TO 990
  1933. END IF
  1934. END IF
  1935. *
  1936. IF( M3.EQ.0 .AND. N.GT.0 ) THEN
  1937. RESULT( NTEST ) = ULPINV
  1938. GO TO 990
  1939. END IF
  1940. *
  1941. * Do test 45 (or +54)
  1942. *
  1943. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  1944. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  1945. IF( N.GT.0 ) THEN
  1946. TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
  1947. ELSE
  1948. TEMP3 = ZERO
  1949. END IF
  1950. RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
  1951. $ MAX( UNFL, TEMP3*ULP )
  1952. *
  1953. 990 CONTINUE
  1954. IF( IUPLO.EQ.1 ) THEN
  1955. INDX = 1
  1956. DO 1010 J = 1, N
  1957. DO 1000 I = 1, J
  1958. WORK( INDX ) = A( I, J )
  1959. INDX = INDX + 1
  1960. 1000 CONTINUE
  1961. 1010 CONTINUE
  1962. ELSE
  1963. INDX = 1
  1964. DO 1030 J = 1, N
  1965. DO 1020 I = J, N
  1966. WORK( INDX ) = A( I, J )
  1967. INDX = INDX + 1
  1968. 1020 CONTINUE
  1969. 1030 CONTINUE
  1970. END IF
  1971. *
  1972. NTEST = NTEST + 1
  1973. *
  1974. SRNAMT = 'SSPEVX'
  1975. CALL SSPEVX( 'V', 'V', UPLO, N, WORK, VL, VU, IL, IU,
  1976. $ ABSTOL, M2, WA2, Z, LDU, V, IWORK,
  1977. $ IWORK( 5*N+1 ), IINFO )
  1978. IF( IINFO.NE.0 ) THEN
  1979. WRITE( NOUNIT, FMT = 9999 )'SSPEVX(V,V,' // UPLO //
  1980. $ ')', IINFO, N, JTYPE, IOLDSD
  1981. INFO = ABS( IINFO )
  1982. IF( IINFO.LT.0 ) THEN
  1983. RETURN
  1984. ELSE
  1985. RESULT( NTEST ) = ULPINV
  1986. RESULT( NTEST+1 ) = ULPINV
  1987. RESULT( NTEST+2 ) = ULPINV
  1988. GO TO 1080
  1989. END IF
  1990. END IF
  1991. *
  1992. * Do tests 46 and 47 (or +54)
  1993. *
  1994. CALL SSYT22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
  1995. $ V, LDU, TAU, WORK, RESULT( NTEST ) )
  1996. *
  1997. NTEST = NTEST + 2
  1998. *
  1999. IF( IUPLO.EQ.1 ) THEN
  2000. INDX = 1
  2001. DO 1050 J = 1, N
  2002. DO 1040 I = 1, J
  2003. WORK( INDX ) = A( I, J )
  2004. INDX = INDX + 1
  2005. 1040 CONTINUE
  2006. 1050 CONTINUE
  2007. ELSE
  2008. INDX = 1
  2009. DO 1070 J = 1, N
  2010. DO 1060 I = J, N
  2011. WORK( INDX ) = A( I, J )
  2012. INDX = INDX + 1
  2013. 1060 CONTINUE
  2014. 1070 CONTINUE
  2015. END IF
  2016. *
  2017. SRNAMT = 'SSPEVX'
  2018. CALL SSPEVX( 'N', 'V', UPLO, N, WORK, VL, VU, IL, IU,
  2019. $ ABSTOL, M3, WA3, Z, LDU, V, IWORK,
  2020. $ IWORK( 5*N+1 ), IINFO )
  2021. IF( IINFO.NE.0 ) THEN
  2022. WRITE( NOUNIT, FMT = 9999 )'SSPEVX(N,V,' // UPLO //
  2023. $ ')', IINFO, N, JTYPE, IOLDSD
  2024. INFO = ABS( IINFO )
  2025. IF( IINFO.LT.0 ) THEN
  2026. RETURN
  2027. ELSE
  2028. RESULT( NTEST ) = ULPINV
  2029. GO TO 1080
  2030. END IF
  2031. END IF
  2032. *
  2033. IF( M3.EQ.0 .AND. N.GT.0 ) THEN
  2034. RESULT( NTEST ) = ULPINV
  2035. GO TO 1080
  2036. END IF
  2037. *
  2038. * Do test 48 (or +54)
  2039. *
  2040. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  2041. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  2042. IF( N.GT.0 ) THEN
  2043. TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
  2044. ELSE
  2045. TEMP3 = ZERO
  2046. END IF
  2047. RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
  2048. $ MAX( UNFL, TEMP3*ULP )
  2049. *
  2050. 1080 CONTINUE
  2051. *
  2052. * 6) Call SSBEV and SSBEVX.
  2053. *
  2054. IF( JTYPE.LE.7 ) THEN
  2055. KD = 1
  2056. ELSE IF( JTYPE.GE.8 .AND. JTYPE.LE.15 ) THEN
  2057. KD = MAX( N-1, 0 )
  2058. ELSE
  2059. KD = IHBW
  2060. END IF
  2061. *
  2062. * Load array V with the upper or lower triangular part
  2063. * of the matrix in band form.
  2064. *
  2065. IF( IUPLO.EQ.1 ) THEN
  2066. DO 1100 J = 1, N
  2067. DO 1090 I = MAX( 1, J-KD ), J
  2068. V( KD+1+I-J, J ) = A( I, J )
  2069. 1090 CONTINUE
  2070. 1100 CONTINUE
  2071. ELSE
  2072. DO 1120 J = 1, N
  2073. DO 1110 I = J, MIN( N, J+KD )
  2074. V( 1+I-J, J ) = A( I, J )
  2075. 1110 CONTINUE
  2076. 1120 CONTINUE
  2077. END IF
  2078. *
  2079. NTEST = NTEST + 1
  2080. SRNAMT = 'SSBEV'
  2081. CALL SSBEV( 'V', UPLO, N, KD, V, LDU, D1, Z, LDU, WORK,
  2082. $ IINFO )
  2083. IF( IINFO.NE.0 ) THEN
  2084. WRITE( NOUNIT, FMT = 9999 )'SSBEV(V,' // UPLO // ')',
  2085. $ IINFO, N, JTYPE, IOLDSD
  2086. INFO = ABS( IINFO )
  2087. IF( IINFO.LT.0 ) THEN
  2088. RETURN
  2089. ELSE
  2090. RESULT( NTEST ) = ULPINV
  2091. RESULT( NTEST+1 ) = ULPINV
  2092. RESULT( NTEST+2 ) = ULPINV
  2093. GO TO 1180
  2094. END IF
  2095. END IF
  2096. *
  2097. * Do tests 49 and 50 (or ... )
  2098. *
  2099. CALL SSYT21( 1, UPLO, N, 0, A, LDA, D1, D2, Z, LDU, V,
  2100. $ LDU, TAU, WORK, RESULT( NTEST ) )
  2101. *
  2102. IF( IUPLO.EQ.1 ) THEN
  2103. DO 1140 J = 1, N
  2104. DO 1130 I = MAX( 1, J-KD ), J
  2105. V( KD+1+I-J, J ) = A( I, J )
  2106. 1130 CONTINUE
  2107. 1140 CONTINUE
  2108. ELSE
  2109. DO 1160 J = 1, N
  2110. DO 1150 I = J, MIN( N, J+KD )
  2111. V( 1+I-J, J ) = A( I, J )
  2112. 1150 CONTINUE
  2113. 1160 CONTINUE
  2114. END IF
  2115. *
  2116. NTEST = NTEST + 2
  2117. SRNAMT = 'SSBEV'
  2118. CALL SSBEV( 'N', UPLO, N, KD, V, LDU, D3, Z, LDU, WORK,
  2119. $ IINFO )
  2120. IF( IINFO.NE.0 ) THEN
  2121. WRITE( NOUNIT, FMT = 9999 )'SSBEV(N,' // UPLO // ')',
  2122. $ IINFO, N, JTYPE, IOLDSD
  2123. INFO = ABS( IINFO )
  2124. IF( IINFO.LT.0 ) THEN
  2125. RETURN
  2126. ELSE
  2127. RESULT( NTEST ) = ULPINV
  2128. GO TO 1180
  2129. END IF
  2130. END IF
  2131. *
  2132. * Do test 51 (or +54)
  2133. *
  2134. TEMP1 = ZERO
  2135. TEMP2 = ZERO
  2136. DO 1170 J = 1, N
  2137. TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
  2138. TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
  2139. 1170 CONTINUE
  2140. RESULT( NTEST ) = TEMP2 / MAX( UNFL,
  2141. $ ULP*MAX( TEMP1, TEMP2 ) )
  2142. *
  2143. * Load array V with the upper or lower triangular part
  2144. * of the matrix in band form.
  2145. *
  2146. 1180 CONTINUE
  2147. IF( IUPLO.EQ.1 ) THEN
  2148. DO 1200 J = 1, N
  2149. DO 1190 I = MAX( 1, J-KD ), J
  2150. V( KD+1+I-J, J ) = A( I, J )
  2151. 1190 CONTINUE
  2152. 1200 CONTINUE
  2153. ELSE
  2154. DO 1220 J = 1, N
  2155. DO 1210 I = J, MIN( N, J+KD )
  2156. V( 1+I-J, J ) = A( I, J )
  2157. 1210 CONTINUE
  2158. 1220 CONTINUE
  2159. END IF
  2160. *
  2161. NTEST = NTEST + 1
  2162. SRNAMT = 'SSBEVX'
  2163. CALL SSBEVX( 'V', 'A', UPLO, N, KD, V, LDU, U, LDU, VL,
  2164. $ VU, IL, IU, ABSTOL, M, WA2, Z, LDU, WORK,
  2165. $ IWORK, IWORK( 5*N+1 ), IINFO )
  2166. IF( IINFO.NE.0 ) THEN
  2167. WRITE( NOUNIT, FMT = 9999 )'SSBEVX(V,A,' // UPLO //
  2168. $ ')', IINFO, N, JTYPE, IOLDSD
  2169. INFO = ABS( IINFO )
  2170. IF( IINFO.LT.0 ) THEN
  2171. RETURN
  2172. ELSE
  2173. RESULT( NTEST ) = ULPINV
  2174. RESULT( NTEST+1 ) = ULPINV
  2175. RESULT( NTEST+2 ) = ULPINV
  2176. GO TO 1280
  2177. END IF
  2178. END IF
  2179. *
  2180. * Do tests 52 and 53 (or +54)
  2181. *
  2182. CALL SSYT21( 1, UPLO, N, 0, A, LDU, WA2, D2, Z, LDU, V,
  2183. $ LDU, TAU, WORK, RESULT( NTEST ) )
  2184. *
  2185. NTEST = NTEST + 2
  2186. *
  2187. IF( IUPLO.EQ.1 ) THEN
  2188. DO 1240 J = 1, N
  2189. DO 1230 I = MAX( 1, J-KD ), J
  2190. V( KD+1+I-J, J ) = A( I, J )
  2191. 1230 CONTINUE
  2192. 1240 CONTINUE
  2193. ELSE
  2194. DO 1260 J = 1, N
  2195. DO 1250 I = J, MIN( N, J+KD )
  2196. V( 1+I-J, J ) = A( I, J )
  2197. 1250 CONTINUE
  2198. 1260 CONTINUE
  2199. END IF
  2200. *
  2201. SRNAMT = 'SSBEVX'
  2202. CALL SSBEVX( 'N', 'A', UPLO, N, KD, V, LDU, U, LDU, VL,
  2203. $ VU, IL, IU, ABSTOL, M3, WA3, Z, LDU, WORK,
  2204. $ IWORK, IWORK( 5*N+1 ), IINFO )
  2205. IF( IINFO.NE.0 ) THEN
  2206. WRITE( NOUNIT, FMT = 9999 )'SSBEVX(N,A,' // UPLO //
  2207. $ ')', IINFO, N, JTYPE, IOLDSD
  2208. INFO = ABS( IINFO )
  2209. IF( IINFO.LT.0 ) THEN
  2210. RETURN
  2211. ELSE
  2212. RESULT( NTEST ) = ULPINV
  2213. GO TO 1280
  2214. END IF
  2215. END IF
  2216. *
  2217. * Do test 54 (or +54)
  2218. *
  2219. TEMP1 = ZERO
  2220. TEMP2 = ZERO
  2221. DO 1270 J = 1, N
  2222. TEMP1 = MAX( TEMP1, ABS( WA2( J ) ), ABS( WA3( J ) ) )
  2223. TEMP2 = MAX( TEMP2, ABS( WA2( J )-WA3( J ) ) )
  2224. 1270 CONTINUE
  2225. RESULT( NTEST ) = TEMP2 / MAX( UNFL,
  2226. $ ULP*MAX( TEMP1, TEMP2 ) )
  2227. *
  2228. 1280 CONTINUE
  2229. NTEST = NTEST + 1
  2230. IF( IUPLO.EQ.1 ) THEN
  2231. DO 1300 J = 1, N
  2232. DO 1290 I = MAX( 1, J-KD ), J
  2233. V( KD+1+I-J, J ) = A( I, J )
  2234. 1290 CONTINUE
  2235. 1300 CONTINUE
  2236. ELSE
  2237. DO 1320 J = 1, N
  2238. DO 1310 I = J, MIN( N, J+KD )
  2239. V( 1+I-J, J ) = A( I, J )
  2240. 1310 CONTINUE
  2241. 1320 CONTINUE
  2242. END IF
  2243. *
  2244. SRNAMT = 'SSBEVX'
  2245. CALL SSBEVX( 'V', 'I', UPLO, N, KD, V, LDU, U, LDU, VL,
  2246. $ VU, IL, IU, ABSTOL, M2, WA2, Z, LDU, WORK,
  2247. $ IWORK, IWORK( 5*N+1 ), IINFO )
  2248. IF( IINFO.NE.0 ) THEN
  2249. WRITE( NOUNIT, FMT = 9999 )'SSBEVX(V,I,' // UPLO //
  2250. $ ')', IINFO, N, JTYPE, IOLDSD
  2251. INFO = ABS( IINFO )
  2252. IF( IINFO.LT.0 ) THEN
  2253. RETURN
  2254. ELSE
  2255. RESULT( NTEST ) = ULPINV
  2256. RESULT( NTEST+1 ) = ULPINV
  2257. RESULT( NTEST+2 ) = ULPINV
  2258. GO TO 1370
  2259. END IF
  2260. END IF
  2261. *
  2262. * Do tests 55 and 56 (or +54)
  2263. *
  2264. CALL SSYT22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
  2265. $ V, LDU, TAU, WORK, RESULT( NTEST ) )
  2266. *
  2267. NTEST = NTEST + 2
  2268. *
  2269. IF( IUPLO.EQ.1 ) THEN
  2270. DO 1340 J = 1, N
  2271. DO 1330 I = MAX( 1, J-KD ), J
  2272. V( KD+1+I-J, J ) = A( I, J )
  2273. 1330 CONTINUE
  2274. 1340 CONTINUE
  2275. ELSE
  2276. DO 1360 J = 1, N
  2277. DO 1350 I = J, MIN( N, J+KD )
  2278. V( 1+I-J, J ) = A( I, J )
  2279. 1350 CONTINUE
  2280. 1360 CONTINUE
  2281. END IF
  2282. *
  2283. SRNAMT = 'SSBEVX'
  2284. CALL SSBEVX( 'N', 'I', UPLO, N, KD, V, LDU, U, LDU, VL,
  2285. $ VU, IL, IU, ABSTOL, M3, WA3, Z, LDU, WORK,
  2286. $ IWORK, IWORK( 5*N+1 ), IINFO )
  2287. IF( IINFO.NE.0 ) THEN
  2288. WRITE( NOUNIT, FMT = 9999 )'SSBEVX(N,I,' // UPLO //
  2289. $ ')', IINFO, N, JTYPE, IOLDSD
  2290. INFO = ABS( IINFO )
  2291. IF( IINFO.LT.0 ) THEN
  2292. RETURN
  2293. ELSE
  2294. RESULT( NTEST ) = ULPINV
  2295. GO TO 1370
  2296. END IF
  2297. END IF
  2298. *
  2299. * Do test 57 (or +54)
  2300. *
  2301. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  2302. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  2303. IF( N.GT.0 ) THEN
  2304. TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
  2305. ELSE
  2306. TEMP3 = ZERO
  2307. END IF
  2308. RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
  2309. $ MAX( UNFL, TEMP3*ULP )
  2310. *
  2311. 1370 CONTINUE
  2312. NTEST = NTEST + 1
  2313. IF( IUPLO.EQ.1 ) THEN
  2314. DO 1390 J = 1, N
  2315. DO 1380 I = MAX( 1, J-KD ), J
  2316. V( KD+1+I-J, J ) = A( I, J )
  2317. 1380 CONTINUE
  2318. 1390 CONTINUE
  2319. ELSE
  2320. DO 1410 J = 1, N
  2321. DO 1400 I = J, MIN( N, J+KD )
  2322. V( 1+I-J, J ) = A( I, J )
  2323. 1400 CONTINUE
  2324. 1410 CONTINUE
  2325. END IF
  2326. *
  2327. SRNAMT = 'SSBEVX'
  2328. CALL SSBEVX( 'V', 'V', UPLO, N, KD, V, LDU, U, LDU, VL,
  2329. $ VU, IL, IU, ABSTOL, M2, WA2, Z, LDU, WORK,
  2330. $ IWORK, IWORK( 5*N+1 ), IINFO )
  2331. IF( IINFO.NE.0 ) THEN
  2332. WRITE( NOUNIT, FMT = 9999 )'SSBEVX(V,V,' // UPLO //
  2333. $ ')', IINFO, N, JTYPE, IOLDSD
  2334. INFO = ABS( IINFO )
  2335. IF( IINFO.LT.0 ) THEN
  2336. RETURN
  2337. ELSE
  2338. RESULT( NTEST ) = ULPINV
  2339. RESULT( NTEST+1 ) = ULPINV
  2340. RESULT( NTEST+2 ) = ULPINV
  2341. GO TO 1460
  2342. END IF
  2343. END IF
  2344. *
  2345. * Do tests 58 and 59 (or +54)
  2346. *
  2347. CALL SSYT22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
  2348. $ V, LDU, TAU, WORK, RESULT( NTEST ) )
  2349. *
  2350. NTEST = NTEST + 2
  2351. *
  2352. IF( IUPLO.EQ.1 ) THEN
  2353. DO 1430 J = 1, N
  2354. DO 1420 I = MAX( 1, J-KD ), J
  2355. V( KD+1+I-J, J ) = A( I, J )
  2356. 1420 CONTINUE
  2357. 1430 CONTINUE
  2358. ELSE
  2359. DO 1450 J = 1, N
  2360. DO 1440 I = J, MIN( N, J+KD )
  2361. V( 1+I-J, J ) = A( I, J )
  2362. 1440 CONTINUE
  2363. 1450 CONTINUE
  2364. END IF
  2365. *
  2366. SRNAMT = 'SSBEVX'
  2367. CALL SSBEVX( 'N', 'V', UPLO, N, KD, V, LDU, U, LDU, VL,
  2368. $ VU, IL, IU, ABSTOL, M3, WA3, Z, LDU, WORK,
  2369. $ IWORK, IWORK( 5*N+1 ), IINFO )
  2370. IF( IINFO.NE.0 ) THEN
  2371. WRITE( NOUNIT, FMT = 9999 )'SSBEVX(N,V,' // UPLO //
  2372. $ ')', IINFO, N, JTYPE, IOLDSD
  2373. INFO = ABS( IINFO )
  2374. IF( IINFO.LT.0 ) THEN
  2375. RETURN
  2376. ELSE
  2377. RESULT( NTEST ) = ULPINV
  2378. GO TO 1460
  2379. END IF
  2380. END IF
  2381. *
  2382. IF( M3.EQ.0 .AND. N.GT.0 ) THEN
  2383. RESULT( NTEST ) = ULPINV
  2384. GO TO 1460
  2385. END IF
  2386. *
  2387. * Do test 60 (or +54)
  2388. *
  2389. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  2390. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  2391. IF( N.GT.0 ) THEN
  2392. TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
  2393. ELSE
  2394. TEMP3 = ZERO
  2395. END IF
  2396. RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
  2397. $ MAX( UNFL, TEMP3*ULP )
  2398. *
  2399. 1460 CONTINUE
  2400. *
  2401. * 7) Call SSYEVD
  2402. *
  2403. CALL SLACPY( ' ', N, N, A, LDA, V, LDU )
  2404. *
  2405. NTEST = NTEST + 1
  2406. SRNAMT = 'SSYEVD'
  2407. CALL SSYEVD( 'V', UPLO, N, A, LDU, D1, WORK, LWEDC,
  2408. $ IWORK, LIWEDC, IINFO )
  2409. IF( IINFO.NE.0 ) THEN
  2410. WRITE( NOUNIT, FMT = 9999 )'SSYEVD(V,' // UPLO //
  2411. $ ')', IINFO, N, JTYPE, IOLDSD
  2412. INFO = ABS( IINFO )
  2413. IF( IINFO.LT.0 ) THEN
  2414. RETURN
  2415. ELSE
  2416. RESULT( NTEST ) = ULPINV
  2417. RESULT( NTEST+1 ) = ULPINV
  2418. RESULT( NTEST+2 ) = ULPINV
  2419. GO TO 1480
  2420. END IF
  2421. END IF
  2422. *
  2423. * Do tests 61 and 62 (or +54)
  2424. *
  2425. CALL SSYT21( 1, UPLO, N, 0, V, LDU, D1, D2, A, LDU, Z,
  2426. $ LDU, TAU, WORK, RESULT( NTEST ) )
  2427. *
  2428. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  2429. *
  2430. NTEST = NTEST + 2
  2431. SRNAMT = 'SSYEVD'
  2432. CALL SSYEVD( 'N', UPLO, N, A, LDU, D3, WORK, LWEDC,
  2433. $ IWORK, LIWEDC, IINFO )
  2434. IF( IINFO.NE.0 ) THEN
  2435. WRITE( NOUNIT, FMT = 9999 )'SSYEVD(N,' // UPLO //
  2436. $ ')', IINFO, N, JTYPE, IOLDSD
  2437. INFO = ABS( IINFO )
  2438. IF( IINFO.LT.0 ) THEN
  2439. RETURN
  2440. ELSE
  2441. RESULT( NTEST ) = ULPINV
  2442. GO TO 1480
  2443. END IF
  2444. END IF
  2445. *
  2446. * Do test 63 (or +54)
  2447. *
  2448. TEMP1 = ZERO
  2449. TEMP2 = ZERO
  2450. DO 1470 J = 1, N
  2451. TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
  2452. TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
  2453. 1470 CONTINUE
  2454. RESULT( NTEST ) = TEMP2 / MAX( UNFL,
  2455. $ ULP*MAX( TEMP1, TEMP2 ) )
  2456. *
  2457. 1480 CONTINUE
  2458. *
  2459. * 8) Call SSPEVD.
  2460. *
  2461. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  2462. *
  2463. * Load array WORK with the upper or lower triangular
  2464. * part of the matrix in packed form.
  2465. *
  2466. IF( IUPLO.EQ.1 ) THEN
  2467. INDX = 1
  2468. DO 1500 J = 1, N
  2469. DO 1490 I = 1, J
  2470. WORK( INDX ) = A( I, J )
  2471. INDX = INDX + 1
  2472. 1490 CONTINUE
  2473. 1500 CONTINUE
  2474. ELSE
  2475. INDX = 1
  2476. DO 1520 J = 1, N
  2477. DO 1510 I = J, N
  2478. WORK( INDX ) = A( I, J )
  2479. INDX = INDX + 1
  2480. 1510 CONTINUE
  2481. 1520 CONTINUE
  2482. END IF
  2483. *
  2484. NTEST = NTEST + 1
  2485. SRNAMT = 'SSPEVD'
  2486. CALL SSPEVD( 'V', UPLO, N, WORK, D1, Z, LDU,
  2487. $ WORK( INDX ), LWEDC-INDX+1, IWORK, LIWEDC,
  2488. $ IINFO )
  2489. IF( IINFO.NE.0 ) THEN
  2490. WRITE( NOUNIT, FMT = 9999 )'SSPEVD(V,' // UPLO //
  2491. $ ')', IINFO, N, JTYPE, IOLDSD
  2492. INFO = ABS( IINFO )
  2493. IF( IINFO.LT.0 ) THEN
  2494. RETURN
  2495. ELSE
  2496. RESULT( NTEST ) = ULPINV
  2497. RESULT( NTEST+1 ) = ULPINV
  2498. RESULT( NTEST+2 ) = ULPINV
  2499. GO TO 1580
  2500. END IF
  2501. END IF
  2502. *
  2503. * Do tests 64 and 65 (or +54)
  2504. *
  2505. CALL SSYT21( 1, UPLO, N, 0, A, LDA, D1, D2, Z, LDU, V,
  2506. $ LDU, TAU, WORK, RESULT( NTEST ) )
  2507. *
  2508. IF( IUPLO.EQ.1 ) THEN
  2509. INDX = 1
  2510. DO 1540 J = 1, N
  2511. DO 1530 I = 1, J
  2512. *
  2513. WORK( INDX ) = A( I, J )
  2514. INDX = INDX + 1
  2515. 1530 CONTINUE
  2516. 1540 CONTINUE
  2517. ELSE
  2518. INDX = 1
  2519. DO 1560 J = 1, N
  2520. DO 1550 I = J, N
  2521. WORK( INDX ) = A( I, J )
  2522. INDX = INDX + 1
  2523. 1550 CONTINUE
  2524. 1560 CONTINUE
  2525. END IF
  2526. *
  2527. NTEST = NTEST + 2
  2528. SRNAMT = 'SSPEVD'
  2529. CALL SSPEVD( 'N', UPLO, N, WORK, D3, Z, LDU,
  2530. $ WORK( INDX ), LWEDC-INDX+1, IWORK, LIWEDC,
  2531. $ IINFO )
  2532. IF( IINFO.NE.0 ) THEN
  2533. WRITE( NOUNIT, FMT = 9999 )'SSPEVD(N,' // UPLO //
  2534. $ ')', IINFO, N, JTYPE, IOLDSD
  2535. INFO = ABS( IINFO )
  2536. IF( IINFO.LT.0 ) THEN
  2537. RETURN
  2538. ELSE
  2539. RESULT( NTEST ) = ULPINV
  2540. GO TO 1580
  2541. END IF
  2542. END IF
  2543. *
  2544. * Do test 66 (or +54)
  2545. *
  2546. TEMP1 = ZERO
  2547. TEMP2 = ZERO
  2548. DO 1570 J = 1, N
  2549. TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
  2550. TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
  2551. 1570 CONTINUE
  2552. RESULT( NTEST ) = TEMP2 / MAX( UNFL,
  2553. $ ULP*MAX( TEMP1, TEMP2 ) )
  2554. 1580 CONTINUE
  2555. *
  2556. * 9) Call SSBEVD.
  2557. *
  2558. IF( JTYPE.LE.7 ) THEN
  2559. KD = 1
  2560. ELSE IF( JTYPE.GE.8 .AND. JTYPE.LE.15 ) THEN
  2561. KD = MAX( N-1, 0 )
  2562. ELSE
  2563. KD = IHBW
  2564. END IF
  2565. *
  2566. * Load array V with the upper or lower triangular part
  2567. * of the matrix in band form.
  2568. *
  2569. IF( IUPLO.EQ.1 ) THEN
  2570. DO 1600 J = 1, N
  2571. DO 1590 I = MAX( 1, J-KD ), J
  2572. V( KD+1+I-J, J ) = A( I, J )
  2573. 1590 CONTINUE
  2574. 1600 CONTINUE
  2575. ELSE
  2576. DO 1620 J = 1, N
  2577. DO 1610 I = J, MIN( N, J+KD )
  2578. V( 1+I-J, J ) = A( I, J )
  2579. 1610 CONTINUE
  2580. 1620 CONTINUE
  2581. END IF
  2582. *
  2583. NTEST = NTEST + 1
  2584. SRNAMT = 'SSBEVD'
  2585. CALL SSBEVD( 'V', UPLO, N, KD, V, LDU, D1, Z, LDU, WORK,
  2586. $ LWEDC, IWORK, LIWEDC, IINFO )
  2587. IF( IINFO.NE.0 ) THEN
  2588. WRITE( NOUNIT, FMT = 9999 )'SSBEVD(V,' // UPLO //
  2589. $ ')', IINFO, N, JTYPE, IOLDSD
  2590. INFO = ABS( IINFO )
  2591. IF( IINFO.LT.0 ) THEN
  2592. RETURN
  2593. ELSE
  2594. RESULT( NTEST ) = ULPINV
  2595. RESULT( NTEST+1 ) = ULPINV
  2596. RESULT( NTEST+2 ) = ULPINV
  2597. GO TO 1680
  2598. END IF
  2599. END IF
  2600. *
  2601. * Do tests 67 and 68 (or +54)
  2602. *
  2603. CALL SSYT21( 1, UPLO, N, 0, A, LDA, D1, D2, Z, LDU, V,
  2604. $ LDU, TAU, WORK, RESULT( NTEST ) )
  2605. *
  2606. IF( IUPLO.EQ.1 ) THEN
  2607. DO 1640 J = 1, N
  2608. DO 1630 I = MAX( 1, J-KD ), J
  2609. V( KD+1+I-J, J ) = A( I, J )
  2610. 1630 CONTINUE
  2611. 1640 CONTINUE
  2612. ELSE
  2613. DO 1660 J = 1, N
  2614. DO 1650 I = J, MIN( N, J+KD )
  2615. V( 1+I-J, J ) = A( I, J )
  2616. 1650 CONTINUE
  2617. 1660 CONTINUE
  2618. END IF
  2619. *
  2620. NTEST = NTEST + 2
  2621. SRNAMT = 'SSBEVD'
  2622. CALL SSBEVD( 'N', UPLO, N, KD, V, LDU, D3, Z, LDU, WORK,
  2623. $ LWEDC, IWORK, LIWEDC, IINFO )
  2624. IF( IINFO.NE.0 ) THEN
  2625. WRITE( NOUNIT, FMT = 9999 )'SSBEVD(N,' // UPLO //
  2626. $ ')', IINFO, N, JTYPE, IOLDSD
  2627. INFO = ABS( IINFO )
  2628. IF( IINFO.LT.0 ) THEN
  2629. RETURN
  2630. ELSE
  2631. RESULT( NTEST ) = ULPINV
  2632. GO TO 1680
  2633. END IF
  2634. END IF
  2635. *
  2636. * Do test 69 (or +54)
  2637. *
  2638. TEMP1 = ZERO
  2639. TEMP2 = ZERO
  2640. DO 1670 J = 1, N
  2641. TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
  2642. TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
  2643. 1670 CONTINUE
  2644. RESULT( NTEST ) = TEMP2 / MAX( UNFL,
  2645. $ ULP*MAX( TEMP1, TEMP2 ) )
  2646. *
  2647. 1680 CONTINUE
  2648. *
  2649. *
  2650. CALL SLACPY( ' ', N, N, A, LDA, V, LDU )
  2651. NTEST = NTEST + 1
  2652. SRNAMT = 'SSYEVR'
  2653. CALL SSYEVR( 'V', 'A', UPLO, N, A, LDU, VL, VU, IL, IU,
  2654. $ ABSTOL, M, WA1, Z, LDU, IWORK, WORK, LWORK,
  2655. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  2656. IF( IINFO.NE.0 ) THEN
  2657. WRITE( NOUNIT, FMT = 9999 )'SSYEVR(V,A,' // UPLO //
  2658. $ ')', IINFO, N, JTYPE, IOLDSD
  2659. INFO = ABS( IINFO )
  2660. IF( IINFO.LT.0 ) THEN
  2661. RETURN
  2662. ELSE
  2663. RESULT( NTEST ) = ULPINV
  2664. RESULT( NTEST+1 ) = ULPINV
  2665. RESULT( NTEST+2 ) = ULPINV
  2666. GO TO 1700
  2667. END IF
  2668. END IF
  2669. *
  2670. * Do tests 70 and 71 (or ... )
  2671. *
  2672. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  2673. *
  2674. CALL SSYT21( 1, UPLO, N, 0, A, LDU, WA1, D2, Z, LDU, V,
  2675. $ LDU, TAU, WORK, RESULT( NTEST ) )
  2676. *
  2677. NTEST = NTEST + 2
  2678. SRNAMT = 'SSYEVR'
  2679. CALL SSYEVR( 'N', 'A', UPLO, N, A, LDU, VL, VU, IL, IU,
  2680. $ ABSTOL, M2, WA2, Z, LDU, IWORK, WORK, LWORK,
  2681. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  2682. IF( IINFO.NE.0 ) THEN
  2683. WRITE( NOUNIT, FMT = 9999 )'SSYEVR(N,A,' // UPLO //
  2684. $ ')', IINFO, N, JTYPE, IOLDSD
  2685. INFO = ABS( IINFO )
  2686. IF( IINFO.LT.0 ) THEN
  2687. RETURN
  2688. ELSE
  2689. RESULT( NTEST ) = ULPINV
  2690. GO TO 1700
  2691. END IF
  2692. END IF
  2693. *
  2694. * Do test 72 (or ... )
  2695. *
  2696. TEMP1 = ZERO
  2697. TEMP2 = ZERO
  2698. DO 1690 J = 1, N
  2699. TEMP1 = MAX( TEMP1, ABS( WA1( J ) ), ABS( WA2( J ) ) )
  2700. TEMP2 = MAX( TEMP2, ABS( WA1( J )-WA2( J ) ) )
  2701. 1690 CONTINUE
  2702. RESULT( NTEST ) = TEMP2 / MAX( UNFL,
  2703. $ ULP*MAX( TEMP1, TEMP2 ) )
  2704. *
  2705. 1700 CONTINUE
  2706. *
  2707. NTEST = NTEST + 1
  2708. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  2709. SRNAMT = 'SSYEVR'
  2710. CALL SSYEVR( 'V', 'I', UPLO, N, A, LDU, VL, VU, IL, IU,
  2711. $ ABSTOL, M2, WA2, Z, LDU, IWORK, WORK, LWORK,
  2712. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  2713. IF( IINFO.NE.0 ) THEN
  2714. WRITE( NOUNIT, FMT = 9999 )'SSYEVR(V,I,' // UPLO //
  2715. $ ')', IINFO, N, JTYPE, IOLDSD
  2716. INFO = ABS( IINFO )
  2717. IF( IINFO.LT.0 ) THEN
  2718. RETURN
  2719. ELSE
  2720. RESULT( NTEST ) = ULPINV
  2721. RESULT( NTEST+1 ) = ULPINV
  2722. RESULT( NTEST+2 ) = ULPINV
  2723. GO TO 1710
  2724. END IF
  2725. END IF
  2726. *
  2727. * Do tests 73 and 74 (or +54)
  2728. *
  2729. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  2730. *
  2731. CALL SSYT22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
  2732. $ V, LDU, TAU, WORK, RESULT( NTEST ) )
  2733. *
  2734. NTEST = NTEST + 2
  2735. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  2736. SRNAMT = 'SSYEVR'
  2737. CALL SSYEVR( 'N', 'I', UPLO, N, A, LDU, VL, VU, IL, IU,
  2738. $ ABSTOL, M3, WA3, Z, LDU, IWORK, WORK, LWORK,
  2739. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  2740. IF( IINFO.NE.0 ) THEN
  2741. WRITE( NOUNIT, FMT = 9999 )'SSYEVR(N,I,' // UPLO //
  2742. $ ')', IINFO, N, JTYPE, IOLDSD
  2743. INFO = ABS( IINFO )
  2744. IF( IINFO.LT.0 ) THEN
  2745. RETURN
  2746. ELSE
  2747. RESULT( NTEST ) = ULPINV
  2748. GO TO 1710
  2749. END IF
  2750. END IF
  2751. *
  2752. * Do test 75 (or +54)
  2753. *
  2754. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  2755. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  2756. RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
  2757. $ MAX( UNFL, ULP*TEMP3 )
  2758. 1710 CONTINUE
  2759. *
  2760. NTEST = NTEST + 1
  2761. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  2762. SRNAMT = 'SSYEVR'
  2763. CALL SSYEVR( 'V', 'V', UPLO, N, A, LDU, VL, VU, IL, IU,
  2764. $ ABSTOL, M2, WA2, Z, LDU, IWORK, WORK, LWORK,
  2765. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  2766. IF( IINFO.NE.0 ) THEN
  2767. WRITE( NOUNIT, FMT = 9999 )'SSYEVR(V,V,' // UPLO //
  2768. $ ')', IINFO, N, JTYPE, IOLDSD
  2769. INFO = ABS( IINFO )
  2770. IF( IINFO.LT.0 ) THEN
  2771. RETURN
  2772. ELSE
  2773. RESULT( NTEST ) = ULPINV
  2774. RESULT( NTEST+1 ) = ULPINV
  2775. RESULT( NTEST+2 ) = ULPINV
  2776. GO TO 700
  2777. END IF
  2778. END IF
  2779. *
  2780. * Do tests 76 and 77 (or +54)
  2781. *
  2782. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  2783. *
  2784. CALL SSYT22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
  2785. $ V, LDU, TAU, WORK, RESULT( NTEST ) )
  2786. *
  2787. NTEST = NTEST + 2
  2788. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  2789. SRNAMT = 'SSYEVR'
  2790. CALL SSYEVR( 'N', 'V', UPLO, N, A, LDU, VL, VU, IL, IU,
  2791. $ ABSTOL, M3, WA3, Z, LDU, IWORK, WORK, LWORK,
  2792. $ IWORK(2*N+1), LIWORK-2*N, IINFO )
  2793. IF( IINFO.NE.0 ) THEN
  2794. WRITE( NOUNIT, FMT = 9999 )'SSYEVR(N,V,' // UPLO //
  2795. $ ')', IINFO, N, JTYPE, IOLDSD
  2796. INFO = ABS( IINFO )
  2797. IF( IINFO.LT.0 ) THEN
  2798. RETURN
  2799. ELSE
  2800. RESULT( NTEST ) = ULPINV
  2801. GO TO 700
  2802. END IF
  2803. END IF
  2804. *
  2805. IF( M3.EQ.0 .AND. N.GT.0 ) THEN
  2806. RESULT( NTEST ) = ULPINV
  2807. GO TO 700
  2808. END IF
  2809. *
  2810. * Do test 78 (or +54)
  2811. *
  2812. TEMP1 = SSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
  2813. TEMP2 = SSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
  2814. IF( N.GT.0 ) THEN
  2815. TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
  2816. ELSE
  2817. TEMP3 = ZERO
  2818. END IF
  2819. RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
  2820. $ MAX( UNFL, TEMP3*ULP )
  2821. *
  2822. CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
  2823. *
  2824. 1720 CONTINUE
  2825. *
  2826. * End of Loop -- Check for RESULT(j) > THRESH
  2827. *
  2828. NTESTT = NTESTT + NTEST
  2829. *
  2830. CALL SLAFTS( 'SST', N, N, JTYPE, NTEST, RESULT, IOLDSD,
  2831. $ THRESH, NOUNIT, NERRS )
  2832. *
  2833. 1730 CONTINUE
  2834. 1740 CONTINUE
  2835. *
  2836. * Summary
  2837. *
  2838. CALL ALASVM( 'SST', NOUNIT, NERRS, NTESTT, 0 )
  2839. *
  2840. 9999 FORMAT( ' SDRVST: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
  2841. $ I6, ', JTYPE=', I6, ', ISEED=(', 3( I5, ',' ), I5, ')' )
  2842. *
  2843. RETURN
  2844. *
  2845. * End of SDRVST
  2846. *
  2847. END