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.

cerrvx.f 42 kB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014101510161017101810191020102110221023102410251026102710281029103010311032103310341035103610371038103910401041104210431044104510461047104810491050105110521053105410551056
  1. *> \brief \b CERRVX
  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 CERRVX( PATH, NUNIT )
  12. *
  13. * .. Scalar Arguments ..
  14. * CHARACTER*3 PATH
  15. * INTEGER NUNIT
  16. * ..
  17. *
  18. *
  19. *> \par Purpose:
  20. * =============
  21. *>
  22. *> \verbatim
  23. *>
  24. *> CERRVX tests the error exits for the COMPLEX driver routines
  25. *> for solving linear systems of equations.
  26. *> \endverbatim
  27. *
  28. * Arguments:
  29. * ==========
  30. *
  31. *> \param[in] PATH
  32. *> \verbatim
  33. *> PATH is CHARACTER*3
  34. *> The LAPACK path name for the routines to be tested.
  35. *> \endverbatim
  36. *>
  37. *> \param[in] NUNIT
  38. *> \verbatim
  39. *> NUNIT is INTEGER
  40. *> The unit number for output.
  41. *> \endverbatim
  42. *
  43. * Authors:
  44. * ========
  45. *
  46. *> \author Univ. of Tennessee
  47. *> \author Univ. of California Berkeley
  48. *> \author Univ. of Colorado Denver
  49. *> \author NAG Ltd.
  50. *
  51. *> \ingroup complex_lin
  52. *
  53. * =====================================================================
  54. SUBROUTINE CERRVX( PATH, NUNIT )
  55. *
  56. * -- LAPACK test routine --
  57. * -- LAPACK is a software package provided by Univ. of Tennessee, --
  58. * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
  59. *
  60. * .. Scalar Arguments ..
  61. CHARACTER*3 PATH
  62. INTEGER NUNIT
  63. * ..
  64. *
  65. * =====================================================================
  66. *
  67. * .. Parameters ..
  68. INTEGER NMAX
  69. PARAMETER ( NMAX = 4 )
  70. * ..
  71. * .. Local Scalars ..
  72. CHARACTER EQ
  73. CHARACTER*2 C2
  74. INTEGER I, INFO, J
  75. REAL RCOND
  76. * ..
  77. * .. Local Arrays ..
  78. INTEGER IP( NMAX )
  79. REAL C( NMAX ), R( NMAX ), R1( NMAX ), R2( NMAX ),
  80. $ RF( NMAX ), RW( NMAX )
  81. COMPLEX A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
  82. $ E( NMAX ), W( 2*NMAX ), X( NMAX )
  83. * ..
  84. * .. External Functions ..
  85. LOGICAL LSAMEN
  86. EXTERNAL LSAMEN
  87. * ..
  88. * .. External Subroutines ..
  89. EXTERNAL CGBSV, CGBSVX, CGESV, CGESVX, CGTSV, CGTSVX,
  90. $ CHESV, CHESV_RK, CHESV_ROOK, CHESVX, CHKXER,
  91. $ CHPSV, CHPSVX, CPBSV, CPBSVX, CPOSV, CPOSVX,
  92. $ CPPSV, CPPSVX, CPTSV, CPTSVX, CSPSV, CSPSVX,
  93. $ CSYSV, CSYSV_AA, CSYSV_RK, CSYSV_ROOK,
  94. $ CSYSVX, CSYSV_AA_2STAGE
  95. * ..
  96. * .. Scalars in Common ..
  97. LOGICAL LERR, OK
  98. CHARACTER*32 SRNAMT
  99. INTEGER INFOT, NOUT
  100. * ..
  101. * .. Common blocks ..
  102. COMMON / INFOC / INFOT, NOUT, OK, LERR
  103. COMMON / SRNAMC / SRNAMT
  104. * ..
  105. * .. Intrinsic Functions ..
  106. INTRINSIC CMPLX, REAL
  107. * ..
  108. * .. Executable Statements ..
  109. *
  110. NOUT = NUNIT
  111. WRITE( NOUT, FMT = * )
  112. C2 = PATH( 2: 3 )
  113. *
  114. * Set the variables to innocuous values.
  115. *
  116. DO 20 J = 1, NMAX
  117. DO 10 I = 1, NMAX
  118. A( I, J ) = CMPLX( 1. / REAL( I+J ), -1. / REAL( I+J ) )
  119. AF( I, J ) = CMPLX( 1. / REAL( I+J ), -1. / REAL( I+J ) )
  120. 10 CONTINUE
  121. B( J ) = 0.E+0
  122. E( J ) = 0.E+0
  123. R1( J ) = 0.E+0
  124. R2( J ) = 0.E+0
  125. W( J ) = 0.E+0
  126. X( J ) = 0.E+0
  127. C( J ) = 0.E+0
  128. R( J ) = 0.E+0
  129. IP( J ) = J
  130. 20 CONTINUE
  131. EQ = ' '
  132. OK = .TRUE.
  133. *
  134. IF( LSAMEN( 2, C2, 'GE' ) ) THEN
  135. *
  136. * CGESV
  137. *
  138. SRNAMT = 'CGESV '
  139. INFOT = 1
  140. CALL CGESV( -1, 0, A, 1, IP, B, 1, INFO )
  141. CALL CHKXER( 'CGESV ', INFOT, NOUT, LERR, OK )
  142. INFOT = 2
  143. CALL CGESV( 0, -1, A, 1, IP, B, 1, INFO )
  144. CALL CHKXER( 'CGESV ', INFOT, NOUT, LERR, OK )
  145. INFOT = 4
  146. CALL CGESV( 2, 1, A, 1, IP, B, 2, INFO )
  147. CALL CHKXER( 'CGESV ', INFOT, NOUT, LERR, OK )
  148. INFOT = 7
  149. CALL CGESV( 2, 1, A, 2, IP, B, 1, INFO )
  150. CALL CHKXER( 'CGESV ', INFOT, NOUT, LERR, OK )
  151. *
  152. * CGESVX
  153. *
  154. SRNAMT = 'CGESVX'
  155. INFOT = 1
  156. CALL CGESVX( '/', 'N', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
  157. $ X, 1, RCOND, R1, R2, W, RW, INFO )
  158. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  159. INFOT = 2
  160. CALL CGESVX( 'N', '/', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
  161. $ X, 1, RCOND, R1, R2, W, RW, INFO )
  162. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  163. INFOT = 3
  164. CALL CGESVX( 'N', 'N', -1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
  165. $ X, 1, RCOND, R1, R2, W, RW, INFO )
  166. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  167. INFOT = 4
  168. CALL CGESVX( 'N', 'N', 0, -1, A, 1, AF, 1, IP, EQ, R, C, B, 1,
  169. $ X, 1, RCOND, R1, R2, W, RW, INFO )
  170. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  171. INFOT = 6
  172. CALL CGESVX( 'N', 'N', 2, 1, A, 1, AF, 2, IP, EQ, R, C, B, 2,
  173. $ X, 2, RCOND, R1, R2, W, RW, INFO )
  174. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  175. INFOT = 8
  176. CALL CGESVX( 'N', 'N', 2, 1, A, 2, AF, 1, IP, EQ, R, C, B, 2,
  177. $ X, 2, RCOND, R1, R2, W, RW, INFO )
  178. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  179. INFOT = 10
  180. EQ = '/'
  181. CALL CGESVX( 'F', 'N', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
  182. $ X, 1, RCOND, R1, R2, W, RW, INFO )
  183. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  184. INFOT = 11
  185. EQ = 'R'
  186. CALL CGESVX( 'F', 'N', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
  187. $ X, 1, RCOND, R1, R2, W, RW, INFO )
  188. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  189. INFOT = 12
  190. EQ = 'C'
  191. CALL CGESVX( 'F', 'N', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
  192. $ X, 1, RCOND, R1, R2, W, RW, INFO )
  193. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  194. INFOT = 14
  195. CALL CGESVX( 'N', 'N', 2, 1, A, 2, AF, 2, IP, EQ, R, C, B, 1,
  196. $ X, 2, RCOND, R1, R2, W, RW, INFO )
  197. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  198. INFOT = 16
  199. CALL CGESVX( 'N', 'N', 2, 1, A, 2, AF, 2, IP, EQ, R, C, B, 2,
  200. $ X, 1, RCOND, R1, R2, W, RW, INFO )
  201. CALL CHKXER( 'CGESVX', INFOT, NOUT, LERR, OK )
  202. *
  203. ELSE IF( LSAMEN( 2, C2, 'GB' ) ) THEN
  204. *
  205. * CGBSV
  206. *
  207. SRNAMT = 'CGBSV '
  208. INFOT = 1
  209. CALL CGBSV( -1, 0, 0, 0, A, 1, IP, B, 1, INFO )
  210. CALL CHKXER( 'CGBSV ', INFOT, NOUT, LERR, OK )
  211. INFOT = 2
  212. CALL CGBSV( 1, -1, 0, 0, A, 1, IP, B, 1, INFO )
  213. CALL CHKXER( 'CGBSV ', INFOT, NOUT, LERR, OK )
  214. INFOT = 3
  215. CALL CGBSV( 1, 0, -1, 0, A, 1, IP, B, 1, INFO )
  216. CALL CHKXER( 'CGBSV ', INFOT, NOUT, LERR, OK )
  217. INFOT = 4
  218. CALL CGBSV( 0, 0, 0, -1, A, 1, IP, B, 1, INFO )
  219. CALL CHKXER( 'CGBSV ', INFOT, NOUT, LERR, OK )
  220. INFOT = 6
  221. CALL CGBSV( 1, 1, 1, 0, A, 3, IP, B, 1, INFO )
  222. CALL CHKXER( 'CGBSV ', INFOT, NOUT, LERR, OK )
  223. INFOT = 9
  224. CALL CGBSV( 2, 0, 0, 0, A, 1, IP, B, 1, INFO )
  225. CALL CHKXER( 'CGBSV ', INFOT, NOUT, LERR, OK )
  226. *
  227. * CGBSVX
  228. *
  229. SRNAMT = 'CGBSVX'
  230. INFOT = 1
  231. CALL CGBSVX( '/', 'N', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
  232. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  233. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  234. INFOT = 2
  235. CALL CGBSVX( 'N', '/', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
  236. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  237. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  238. INFOT = 3
  239. CALL CGBSVX( 'N', 'N', -1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
  240. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  241. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  242. INFOT = 4
  243. CALL CGBSVX( 'N', 'N', 1, -1, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
  244. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  245. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  246. INFOT = 5
  247. CALL CGBSVX( 'N', 'N', 1, 0, -1, 0, A, 1, AF, 1, IP, EQ, R, C,
  248. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  249. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  250. INFOT = 6
  251. CALL CGBSVX( 'N', 'N', 0, 0, 0, -1, A, 1, AF, 1, IP, EQ, R, C,
  252. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  253. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  254. INFOT = 8
  255. CALL CGBSVX( 'N', 'N', 1, 1, 1, 0, A, 2, AF, 4, IP, EQ, R, C,
  256. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  257. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  258. INFOT = 10
  259. CALL CGBSVX( 'N', 'N', 1, 1, 1, 0, A, 3, AF, 3, IP, EQ, R, C,
  260. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  261. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  262. INFOT = 12
  263. EQ = '/'
  264. CALL CGBSVX( 'F', 'N', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
  265. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  266. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  267. INFOT = 13
  268. EQ = 'R'
  269. CALL CGBSVX( 'F', 'N', 1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
  270. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  271. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  272. INFOT = 14
  273. EQ = 'C'
  274. CALL CGBSVX( 'F', 'N', 1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
  275. $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  276. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  277. INFOT = 16
  278. CALL CGBSVX( 'N', 'N', 2, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
  279. $ B, 1, X, 2, RCOND, R1, R2, W, RW, INFO )
  280. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  281. INFOT = 18
  282. CALL CGBSVX( 'N', 'N', 2, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
  283. $ B, 2, X, 1, RCOND, R1, R2, W, RW, INFO )
  284. CALL CHKXER( 'CGBSVX', INFOT, NOUT, LERR, OK )
  285. *
  286. ELSE IF( LSAMEN( 2, C2, 'GT' ) ) THEN
  287. *
  288. * CGTSV
  289. *
  290. SRNAMT = 'CGTSV '
  291. INFOT = 1
  292. CALL CGTSV( -1, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
  293. $ INFO )
  294. CALL CHKXER( 'CGTSV ', INFOT, NOUT, LERR, OK )
  295. INFOT = 2
  296. CALL CGTSV( 0, -1, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
  297. $ INFO )
  298. CALL CHKXER( 'CGTSV ', INFOT, NOUT, LERR, OK )
  299. INFOT = 7
  300. CALL CGTSV( 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1, INFO )
  301. CALL CHKXER( 'CGTSV ', INFOT, NOUT, LERR, OK )
  302. *
  303. * CGTSVX
  304. *
  305. SRNAMT = 'CGTSVX'
  306. INFOT = 1
  307. CALL CGTSVX( '/', 'N', 0, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
  308. $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
  309. $ IP, B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  310. CALL CHKXER( 'CGTSVX', INFOT, NOUT, LERR, OK )
  311. INFOT = 2
  312. CALL CGTSVX( 'N', '/', 0, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
  313. $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
  314. $ IP, B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  315. CALL CHKXER( 'CGTSVX', INFOT, NOUT, LERR, OK )
  316. INFOT = 3
  317. CALL CGTSVX( 'N', 'N', -1, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
  318. $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
  319. $ IP, B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  320. CALL CHKXER( 'CGTSVX', INFOT, NOUT, LERR, OK )
  321. INFOT = 4
  322. CALL CGTSVX( 'N', 'N', 0, -1, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
  323. $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
  324. $ IP, B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
  325. CALL CHKXER( 'CGTSVX', INFOT, NOUT, LERR, OK )
  326. INFOT = 14
  327. CALL CGTSVX( 'N', 'N', 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
  328. $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
  329. $ IP, B, 1, X, 2, RCOND, R1, R2, W, RW, INFO )
  330. CALL CHKXER( 'CGTSVX', INFOT, NOUT, LERR, OK )
  331. INFOT = 16
  332. CALL CGTSVX( 'N', 'N', 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
  333. $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
  334. $ IP, B, 2, X, 1, RCOND, R1, R2, W, RW, INFO )
  335. CALL CHKXER( 'CGTSVX', INFOT, NOUT, LERR, OK )
  336. *
  337. ELSE IF( LSAMEN( 2, C2, 'PO' ) ) THEN
  338. *
  339. * CPOSV
  340. *
  341. SRNAMT = 'CPOSV '
  342. INFOT = 1
  343. CALL CPOSV( '/', 0, 0, A, 1, B, 1, INFO )
  344. CALL CHKXER( 'CPOSV ', INFOT, NOUT, LERR, OK )
  345. INFOT = 2
  346. CALL CPOSV( 'U', -1, 0, A, 1, B, 1, INFO )
  347. CALL CHKXER( 'CPOSV ', INFOT, NOUT, LERR, OK )
  348. INFOT = 3
  349. CALL CPOSV( 'U', 0, -1, A, 1, B, 1, INFO )
  350. CALL CHKXER( 'CPOSV ', INFOT, NOUT, LERR, OK )
  351. INFOT = 5
  352. CALL CPOSV( 'U', 2, 0, A, 1, B, 2, INFO )
  353. CALL CHKXER( 'CPOSV ', INFOT, NOUT, LERR, OK )
  354. INFOT = 7
  355. CALL CPOSV( 'U', 2, 0, A, 2, B, 1, INFO )
  356. CALL CHKXER( 'CPOSV ', INFOT, NOUT, LERR, OK )
  357. *
  358. * CPOSVX
  359. *
  360. SRNAMT = 'CPOSVX'
  361. INFOT = 1
  362. CALL CPOSVX( '/', 'U', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
  363. $ RCOND, R1, R2, W, RW, INFO )
  364. CALL CHKXER( 'CPOSVX', INFOT, NOUT, LERR, OK )
  365. INFOT = 2
  366. CALL CPOSVX( 'N', '/', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
  367. $ RCOND, R1, R2, W, RW, INFO )
  368. CALL CHKXER( 'CPOSVX', INFOT, NOUT, LERR, OK )
  369. INFOT = 3
  370. CALL CPOSVX( 'N', 'U', -1, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
  371. $ RCOND, R1, R2, W, RW, INFO )
  372. CALL CHKXER( 'CPOSVX', INFOT, NOUT, LERR, OK )
  373. INFOT = 4
  374. CALL CPOSVX( 'N', 'U', 0, -1, A, 1, AF, 1, EQ, C, B, 1, X, 1,
  375. $ RCOND, R1, R2, W, RW, INFO )
  376. CALL CHKXER( 'CPOSVX', INFOT, NOUT, LERR, OK )
  377. INFOT = 6
  378. CALL CPOSVX( 'N', 'U', 2, 0, A, 1, AF, 2, EQ, C, B, 2, X, 2,
  379. $ RCOND, R1, R2, W, RW, INFO )
  380. CALL CHKXER( 'CPOSVX', INFOT, NOUT, LERR, OK )
  381. INFOT = 8
  382. CALL CPOSVX( 'N', 'U', 2, 0, A, 2, AF, 1, EQ, C, B, 2, X, 2,
  383. $ RCOND, R1, R2, W, RW, INFO )
  384. CALL CHKXER( 'CPOSVX', INFOT, NOUT, LERR, OK )
  385. INFOT = 9
  386. EQ = '/'
  387. CALL CPOSVX( 'F', 'U', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
  388. $ RCOND, R1, R2, W, RW, INFO )
  389. CALL CHKXER( 'CPOSVX', INFOT, NOUT, LERR, OK )
  390. INFOT = 10
  391. EQ = 'Y'
  392. CALL CPOSVX( 'F', 'U', 1, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
  393. $ RCOND, R1, R2, W, RW, INFO )
  394. CALL CHKXER( 'CPOSVX', INFOT, NOUT, LERR, OK )
  395. INFOT = 12
  396. CALL CPOSVX( 'N', 'U', 2, 0, A, 2, AF, 2, EQ, C, B, 1, X, 2,
  397. $ RCOND, R1, R2, W, RW, INFO )
  398. CALL CHKXER( 'CPOSVX', INFOT, NOUT, LERR, OK )
  399. INFOT = 14
  400. CALL CPOSVX( 'N', 'U', 2, 0, A, 2, AF, 2, EQ, C, B, 2, X, 1,
  401. $ RCOND, R1, R2, W, RW, INFO )
  402. CALL CHKXER( 'CPOSVX', INFOT, NOUT, LERR, OK )
  403. *
  404. ELSE IF( LSAMEN( 2, C2, 'PP' ) ) THEN
  405. *
  406. * CPPSV
  407. *
  408. SRNAMT = 'CPPSV '
  409. INFOT = 1
  410. CALL CPPSV( '/', 0, 0, A, B, 1, INFO )
  411. CALL CHKXER( 'CPPSV ', INFOT, NOUT, LERR, OK )
  412. INFOT = 2
  413. CALL CPPSV( 'U', -1, 0, A, B, 1, INFO )
  414. CALL CHKXER( 'CPPSV ', INFOT, NOUT, LERR, OK )
  415. INFOT = 3
  416. CALL CPPSV( 'U', 0, -1, A, B, 1, INFO )
  417. CALL CHKXER( 'CPPSV ', INFOT, NOUT, LERR, OK )
  418. INFOT = 6
  419. CALL CPPSV( 'U', 2, 0, A, B, 1, INFO )
  420. CALL CHKXER( 'CPPSV ', INFOT, NOUT, LERR, OK )
  421. *
  422. * CPPSVX
  423. *
  424. SRNAMT = 'CPPSVX'
  425. INFOT = 1
  426. CALL CPPSVX( '/', 'U', 0, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
  427. $ R1, R2, W, RW, INFO )
  428. CALL CHKXER( 'CPPSVX', INFOT, NOUT, LERR, OK )
  429. INFOT = 2
  430. CALL CPPSVX( 'N', '/', 0, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
  431. $ R1, R2, W, RW, INFO )
  432. CALL CHKXER( 'CPPSVX', INFOT, NOUT, LERR, OK )
  433. INFOT = 3
  434. CALL CPPSVX( 'N', 'U', -1, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
  435. $ R1, R2, W, RW, INFO )
  436. CALL CHKXER( 'CPPSVX', INFOT, NOUT, LERR, OK )
  437. INFOT = 4
  438. CALL CPPSVX( 'N', 'U', 0, -1, A, AF, EQ, C, B, 1, X, 1, RCOND,
  439. $ R1, R2, W, RW, INFO )
  440. CALL CHKXER( 'CPPSVX', INFOT, NOUT, LERR, OK )
  441. INFOT = 7
  442. EQ = '/'
  443. CALL CPPSVX( 'F', 'U', 0, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
  444. $ R1, R2, W, RW, INFO )
  445. CALL CHKXER( 'CPPSVX', INFOT, NOUT, LERR, OK )
  446. INFOT = 8
  447. EQ = 'Y'
  448. CALL CPPSVX( 'F', 'U', 1, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
  449. $ R1, R2, W, RW, INFO )
  450. CALL CHKXER( 'CPPSVX', INFOT, NOUT, LERR, OK )
  451. INFOT = 10
  452. CALL CPPSVX( 'N', 'U', 2, 0, A, AF, EQ, C, B, 1, X, 2, RCOND,
  453. $ R1, R2, W, RW, INFO )
  454. CALL CHKXER( 'CPPSVX', INFOT, NOUT, LERR, OK )
  455. INFOT = 12
  456. CALL CPPSVX( 'N', 'U', 2, 0, A, AF, EQ, C, B, 2, X, 1, RCOND,
  457. $ R1, R2, W, RW, INFO )
  458. CALL CHKXER( 'CPPSVX', INFOT, NOUT, LERR, OK )
  459. *
  460. ELSE IF( LSAMEN( 2, C2, 'PB' ) ) THEN
  461. *
  462. * CPBSV
  463. *
  464. SRNAMT = 'CPBSV '
  465. INFOT = 1
  466. CALL CPBSV( '/', 0, 0, 0, A, 1, B, 1, INFO )
  467. CALL CHKXER( 'CPBSV ', INFOT, NOUT, LERR, OK )
  468. INFOT = 2
  469. CALL CPBSV( 'U', -1, 0, 0, A, 1, B, 1, INFO )
  470. CALL CHKXER( 'CPBSV ', INFOT, NOUT, LERR, OK )
  471. INFOT = 3
  472. CALL CPBSV( 'U', 1, -1, 0, A, 1, B, 1, INFO )
  473. CALL CHKXER( 'CPBSV ', INFOT, NOUT, LERR, OK )
  474. INFOT = 4
  475. CALL CPBSV( 'U', 0, 0, -1, A, 1, B, 1, INFO )
  476. CALL CHKXER( 'CPBSV ', INFOT, NOUT, LERR, OK )
  477. INFOT = 6
  478. CALL CPBSV( 'U', 1, 1, 0, A, 1, B, 2, INFO )
  479. CALL CHKXER( 'CPBSV ', INFOT, NOUT, LERR, OK )
  480. INFOT = 8
  481. CALL CPBSV( 'U', 2, 0, 0, A, 1, B, 1, INFO )
  482. CALL CHKXER( 'CPBSV ', INFOT, NOUT, LERR, OK )
  483. *
  484. * CPBSVX
  485. *
  486. SRNAMT = 'CPBSVX'
  487. INFOT = 1
  488. CALL CPBSVX( '/', 'U', 0, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
  489. $ RCOND, R1, R2, W, RW, INFO )
  490. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  491. INFOT = 2
  492. CALL CPBSVX( 'N', '/', 0, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
  493. $ RCOND, R1, R2, W, RW, INFO )
  494. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  495. INFOT = 3
  496. CALL CPBSVX( 'N', 'U', -1, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X,
  497. $ 1, RCOND, R1, R2, W, RW, INFO )
  498. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  499. INFOT = 4
  500. CALL CPBSVX( 'N', 'U', 1, -1, 0, A, 1, AF, 1, EQ, C, B, 1, X,
  501. $ 1, RCOND, R1, R2, W, RW, INFO )
  502. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  503. INFOT = 5
  504. CALL CPBSVX( 'N', 'U', 0, 0, -1, A, 1, AF, 1, EQ, C, B, 1, X,
  505. $ 1, RCOND, R1, R2, W, RW, INFO )
  506. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  507. INFOT = 7
  508. CALL CPBSVX( 'N', 'U', 1, 1, 0, A, 1, AF, 2, EQ, C, B, 2, X, 2,
  509. $ RCOND, R1, R2, W, RW, INFO )
  510. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  511. INFOT = 9
  512. CALL CPBSVX( 'N', 'U', 1, 1, 0, A, 2, AF, 1, EQ, C, B, 2, X, 2,
  513. $ RCOND, R1, R2, W, RW, INFO )
  514. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  515. INFOT = 10
  516. EQ = '/'
  517. CALL CPBSVX( 'F', 'U', 0, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
  518. $ RCOND, R1, R2, W, RW, INFO )
  519. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  520. INFOT = 11
  521. EQ = 'Y'
  522. CALL CPBSVX( 'F', 'U', 1, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
  523. $ RCOND, R1, R2, W, RW, INFO )
  524. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  525. INFOT = 13
  526. CALL CPBSVX( 'N', 'U', 2, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 2,
  527. $ RCOND, R1, R2, W, RW, INFO )
  528. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  529. INFOT = 15
  530. CALL CPBSVX( 'N', 'U', 2, 0, 0, A, 1, AF, 1, EQ, C, B, 2, X, 1,
  531. $ RCOND, R1, R2, W, RW, INFO )
  532. CALL CHKXER( 'CPBSVX', INFOT, NOUT, LERR, OK )
  533. *
  534. ELSE IF( LSAMEN( 2, C2, 'PT' ) ) THEN
  535. *
  536. * CPTSV
  537. *
  538. SRNAMT = 'CPTSV '
  539. INFOT = 1
  540. CALL CPTSV( -1, 0, R, A( 1, 1 ), B, 1, INFO )
  541. CALL CHKXER( 'CPTSV ', INFOT, NOUT, LERR, OK )
  542. INFOT = 2
  543. CALL CPTSV( 0, -1, R, A( 1, 1 ), B, 1, INFO )
  544. CALL CHKXER( 'CPTSV ', INFOT, NOUT, LERR, OK )
  545. INFOT = 6
  546. CALL CPTSV( 2, 0, R, A( 1, 1 ), B, 1, INFO )
  547. CALL CHKXER( 'CPTSV ', INFOT, NOUT, LERR, OK )
  548. *
  549. * CPTSVX
  550. *
  551. SRNAMT = 'CPTSVX'
  552. INFOT = 1
  553. CALL CPTSVX( '/', 0, 0, R, A( 1, 1 ), RF, AF( 1, 1 ), B, 1, X,
  554. $ 1, RCOND, R1, R2, W, RW, INFO )
  555. CALL CHKXER( 'CPTSVX', INFOT, NOUT, LERR, OK )
  556. INFOT = 2
  557. CALL CPTSVX( 'N', -1, 0, R, A( 1, 1 ), RF, AF( 1, 1 ), B, 1, X,
  558. $ 1, RCOND, R1, R2, W, RW, INFO )
  559. CALL CHKXER( 'CPTSVX', INFOT, NOUT, LERR, OK )
  560. INFOT = 3
  561. CALL CPTSVX( 'N', 0, -1, R, A( 1, 1 ), RF, AF( 1, 1 ), B, 1, X,
  562. $ 1, RCOND, R1, R2, W, RW, INFO )
  563. CALL CHKXER( 'CPTSVX', INFOT, NOUT, LERR, OK )
  564. INFOT = 9
  565. CALL CPTSVX( 'N', 2, 0, R, A( 1, 1 ), RF, AF( 1, 1 ), B, 1, X,
  566. $ 2, RCOND, R1, R2, W, RW, INFO )
  567. CALL CHKXER( 'CPTSVX', INFOT, NOUT, LERR, OK )
  568. INFOT = 11
  569. CALL CPTSVX( 'N', 2, 0, R, A( 1, 1 ), RF, AF( 1, 1 ), B, 2, X,
  570. $ 1, RCOND, R1, R2, W, RW, INFO )
  571. CALL CHKXER( 'CPTSVX', INFOT, NOUT, LERR, OK )
  572. *
  573. ELSE IF( LSAMEN( 2, C2, 'HE' ) ) THEN
  574. *
  575. * CHESV
  576. *
  577. SRNAMT = 'CHESV '
  578. INFOT = 1
  579. CALL CHESV( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
  580. CALL CHKXER( 'CHESV ', INFOT, NOUT, LERR, OK )
  581. INFOT = 2
  582. CALL CHESV( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
  583. CALL CHKXER( 'CHESV ', INFOT, NOUT, LERR, OK )
  584. INFOT = 3
  585. CALL CHESV( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
  586. CALL CHKXER( 'CHESV ', INFOT, NOUT, LERR, OK )
  587. INFOT = 5
  588. CALL CHESV( 'U', 2, 0, A, 1, IP, B, 2, W, 1, INFO )
  589. CALL CHKXER( 'CHESV ', INFOT, NOUT, LERR, OK )
  590. INFOT = 8
  591. CALL CHESV( 'U', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
  592. CALL CHKXER( 'CHESV ', INFOT, NOUT, LERR, OK )
  593. INFOT = 10
  594. CALL CHESV( 'U', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
  595. CALL CHKXER( 'CHESV ', INFOT, NOUT, LERR, OK )
  596. INFOT = 10
  597. CALL CHESV( 'U', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
  598. CALL CHKXER( 'CHESV ', INFOT, NOUT, LERR, OK )
  599. *
  600. * CHESVX
  601. *
  602. SRNAMT = 'CHESVX'
  603. INFOT = 1
  604. CALL CHESVX( '/', 'U', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1,
  605. $ RCOND, R1, R2, W, 1, RW, INFO )
  606. CALL CHKXER( 'CHESVX', INFOT, NOUT, LERR, OK )
  607. INFOT = 2
  608. CALL CHESVX( 'N', '/', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1,
  609. $ RCOND, R1, R2, W, 1, RW, INFO )
  610. CALL CHKXER( 'CHESVX', INFOT, NOUT, LERR, OK )
  611. INFOT = 3
  612. CALL CHESVX( 'N', 'U', -1, 0, A, 1, AF, 1, IP, B, 1, X, 1,
  613. $ RCOND, R1, R2, W, 1, RW, INFO )
  614. CALL CHKXER( 'CHESVX', INFOT, NOUT, LERR, OK )
  615. INFOT = 4
  616. CALL CHESVX( 'N', 'U', 0, -1, A, 1, AF, 1, IP, B, 1, X, 1,
  617. $ RCOND, R1, R2, W, 1, RW, INFO )
  618. CALL CHKXER( 'CHESVX', INFOT, NOUT, LERR, OK )
  619. INFOT = 6
  620. CALL CHESVX( 'N', 'U', 2, 0, A, 1, AF, 2, IP, B, 2, X, 2,
  621. $ RCOND, R1, R2, W, 4, RW, INFO )
  622. CALL CHKXER( 'CHESVX', INFOT, NOUT, LERR, OK )
  623. INFOT = 8
  624. CALL CHESVX( 'N', 'U', 2, 0, A, 2, AF, 1, IP, B, 2, X, 2,
  625. $ RCOND, R1, R2, W, 4, RW, INFO )
  626. CALL CHKXER( 'CHESVX', INFOT, NOUT, LERR, OK )
  627. INFOT = 11
  628. CALL CHESVX( 'N', 'U', 2, 0, A, 2, AF, 2, IP, B, 1, X, 2,
  629. $ RCOND, R1, R2, W, 4, RW, INFO )
  630. CALL CHKXER( 'CHESVX', INFOT, NOUT, LERR, OK )
  631. INFOT = 13
  632. CALL CHESVX( 'N', 'U', 2, 0, A, 2, AF, 2, IP, B, 2, X, 1,
  633. $ RCOND, R1, R2, W, 4, RW, INFO )
  634. CALL CHKXER( 'CHESVX', INFOT, NOUT, LERR, OK )
  635. INFOT = 18
  636. CALL CHESVX( 'N', 'U', 2, 0, A, 2, AF, 2, IP, B, 2, X, 2,
  637. $ RCOND, R1, R2, W, 3, RW, INFO )
  638. CALL CHKXER( 'CHESVX', INFOT, NOUT, LERR, OK )
  639. *
  640. ELSE IF( LSAMEN( 2, C2, 'HR' ) ) THEN
  641. *
  642. * CHESV_ROOK
  643. *
  644. SRNAMT = 'CHESV_ROOK'
  645. INFOT = 1
  646. CALL CHESV_ROOK( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
  647. CALL CHKXER( 'CHESV_ROOK', INFOT, NOUT, LERR, OK )
  648. INFOT = 2
  649. CALL CHESV_ROOK( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
  650. CALL CHKXER( 'CHESV_ROOK', INFOT, NOUT, LERR, OK )
  651. INFOT = 3
  652. CALL CHESV_ROOK( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
  653. CALL CHKXER( 'CHESV_ROOK', INFOT, NOUT, LERR, OK )
  654. INFOT = 5
  655. CALL CHESV_ROOK( 'U', 2, 0, A, 1, IP, B, 2, W, 1, INFO )
  656. CALL CHKXER( 'CHESV_ROOK', INFOT, NOUT, LERR, OK )
  657. INFOT = 8
  658. CALL CHESV_ROOK( 'U', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
  659. CALL CHKXER( 'CHESV_ROOK', INFOT, NOUT, LERR, OK )
  660. INFOT = 10
  661. CALL CHESV_ROOK( 'U', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
  662. CALL CHKXER( 'CHESV_ROOK', INFOT, NOUT, LERR, OK )
  663. INFOT = 10
  664. CALL CHESV_ROOK( 'U', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
  665. CALL CHKXER( 'CHESV_ROOK', INFOT, NOUT, LERR, OK )
  666. *
  667. ELSE IF( LSAMEN( 2, C2, 'HK' ) ) THEN
  668. *
  669. * CHESV_RK
  670. *
  671. * Test error exits of the driver that uses factorization
  672. * of a symmetric indefinite matrix with rook
  673. * (bounded Bunch-Kaufman) pivoting with the new storage
  674. * format for factors L ( or U) and D.
  675. *
  676. * L (or U) is stored in A, diagonal of D is stored on the
  677. * diagonal of A, subdiagonal of D is stored in a separate array E.
  678. *
  679. SRNAMT = 'CHESV_RK'
  680. INFOT = 1
  681. CALL CHESV_RK( '/', 0, 0, A, 1, E, IP, B, 1, W, 1, INFO )
  682. CALL CHKXER( 'CHESV_RK', INFOT, NOUT, LERR, OK )
  683. INFOT = 2
  684. CALL CHESV_RK( 'U', -1, 0, A, 1, E, IP, B, 1, W, 1, INFO )
  685. CALL CHKXER( 'CHESV_RK', INFOT, NOUT, LERR, OK )
  686. INFOT = 3
  687. CALL CHESV_RK( 'U', 0, -1, A, 1, E, IP, B, 1, W, 1, INFO )
  688. CALL CHKXER( 'CHESV_RK', INFOT, NOUT, LERR, OK )
  689. INFOT = 5
  690. CALL CHESV_RK( 'U', 2, 0, A, 1, E, IP, B, 2, W, 1, INFO )
  691. CALL CHKXER( 'CHESV_RK', INFOT, NOUT, LERR, OK )
  692. INFOT = 9
  693. CALL CHESV_RK( 'U', 2, 0, A, 2, E, IP, B, 1, W, 1, INFO )
  694. CALL CHKXER( 'CHESV_RK', INFOT, NOUT, LERR, OK )
  695. INFOT = 11
  696. CALL CHESV_RK( 'U', 0, 0, A, 1, E, IP, B, 1, W, 0, INFO )
  697. CALL CHKXER( 'CHESV_RK', INFOT, NOUT, LERR, OK )
  698. INFOT = 11
  699. CALL CHESV_RK( 'U', 0, 0, A, 1, E, IP, B, 1, W, -2, INFO )
  700. CALL CHKXER( 'CHESV_RK', INFOT, NOUT, LERR, OK )
  701. *
  702. ELSE IF( LSAMEN( 2, C2, 'HA' ) ) THEN
  703. *
  704. * CHESV_AASEN
  705. *
  706. SRNAMT = 'CHESV_AA'
  707. INFOT = 1
  708. CALL CHESV_AA( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
  709. CALL CHKXER( 'CHESV_AA', INFOT, NOUT, LERR, OK )
  710. INFOT = 2
  711. CALL CHESV_AA( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
  712. CALL CHKXER( 'CHESV_AA', INFOT, NOUT, LERR, OK )
  713. INFOT = 3
  714. CALL CHESV_AA( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
  715. CALL CHKXER( 'CHESV_AA', INFOT, NOUT, LERR, OK )
  716. INFOT = 5
  717. CALL CHESV_AA( 'U', 2, 0, A, 1, IP, B, 2, W, 1, INFO )
  718. CALL CHKXER( 'CHESV_AA', INFOT, NOUT, LERR, OK )
  719. INFOT = 8
  720. CALL CHESV_AA( 'U', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
  721. CALL CHKXER( 'CHESV_AA', INFOT, NOUT, LERR, OK )
  722. INFOT = 10
  723. CALL CHESV_AA( 'U', 3, 1, A, 3, IP, B, 3, W, 6, INFO )
  724. CALL CHKXER( 'CHESV_AA', INFOT, NOUT, LERR, OK )
  725. *
  726. ELSE IF( LSAMEN( 2, C2, 'H2' ) ) THEN
  727. *
  728. * CHESV_AASEN_2STAGE
  729. *
  730. SRNAMT = 'CHESV_AA_2STAGE'
  731. INFOT = 1
  732. CALL CHESV_AA_2STAGE( '/', 0, 0, A, 1, A, 1, IP, IP, B, 1,
  733. $ W, 1, INFO )
  734. CALL CHKXER( 'CHESV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  735. INFOT = 2
  736. CALL CHESV_AA_2STAGE( 'U', -1, 0, A, 1, A, 1, IP, IP, B, 1,
  737. $ W, 1, INFO )
  738. CALL CHKXER( 'CHESV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  739. INFOT = 3
  740. CALL CHESV_AA_2STAGE( 'U', 0, -1, A, 1, A, 1, IP, IP, B, 1,
  741. $ W, 1, INFO )
  742. CALL CHKXER( 'CHESV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  743. INFOT = 5
  744. CALL CHESV_AA_2STAGE( 'U', 2, 1, A, 1, A, 1, IP, IP, B, 1,
  745. $ W, 1, INFO )
  746. CALL CHKXER( 'CHESV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  747. INFOT = 7
  748. CALL CHESV_AA_2STAGE( 'U', 2, 1, A, 2, A, 1, IP, IP, B, 2,
  749. $ W, 1, INFO )
  750. CALL CHKXER( 'CHESV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  751. INFOT = 11
  752. CALL CHESV_AA_2STAGE( 'U', 2, 1, A, 2, A, 8, IP, IP, B, 1,
  753. $ W, 1, INFO )
  754. CALL CHKXER( 'CHESV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  755. INFOT = 13
  756. CALL CHESV_AA_2STAGE( 'U', 2, 1, A, 2, A, 8, IP, IP, B, 2,
  757. $ W, 1, INFO )
  758. CALL CHKXER( 'CHESV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  759. *
  760. ELSE IF( LSAMEN( 2, C2, 'SA' ) ) THEN
  761. *
  762. * CSYSV_AASEN
  763. *
  764. SRNAMT = 'CSYSV_AA'
  765. INFOT = 1
  766. CALL CSYSV_AA( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
  767. CALL CHKXER( 'CSYSV_AA', INFOT, NOUT, LERR, OK )
  768. INFOT = 2
  769. CALL CSYSV_AA( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
  770. CALL CHKXER( 'CSYSV_AA', INFOT, NOUT, LERR, OK )
  771. INFOT = 3
  772. CALL CSYSV_AA( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
  773. CALL CHKXER( 'CSYSV_AA', INFOT, NOUT, LERR, OK )
  774. INFOT = 5
  775. CALL CSYSV_AA( 'U', 2, 0, A, 1, IP, B, 2, W, 1, INFO )
  776. CALL CHKXER( 'CSYSV_AA', INFOT, NOUT, LERR, OK )
  777. INFOT = 8
  778. CALL CSYSV_AA( 'U', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
  779. CALL CHKXER( 'CSYSV_AA', INFOT, NOUT, LERR, OK )
  780. INFOT = 10
  781. CALL CSYSV_AA( 'U', 3, 1, A, 3, IP, B, 3, W, 6, INFO )
  782. CALL CHKXER( 'CSYSV_AA', INFOT, NOUT, LERR, OK )
  783. *
  784. ELSE IF( LSAMEN( 2, C2, 'S2' ) ) THEN
  785. *
  786. * CSYSV_AASEN_2STAGE
  787. *
  788. SRNAMT = 'CSYSV_AA_2STAGE'
  789. INFOT = 1
  790. CALL CSYSV_AA_2STAGE( '/', 0, 0, A, 1, A, 1, IP, IP, B, 1,
  791. $ W, 1, INFO )
  792. CALL CHKXER( 'CSYSV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  793. INFOT = 2
  794. CALL CSYSV_AA_2STAGE( 'U', -1, 0, A, 1, A, 1, IP, IP, B, 1,
  795. $ W, 1, INFO )
  796. CALL CHKXER( 'CSYSV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  797. INFOT = 3
  798. CALL CSYSV_AA_2STAGE( 'U', 0, -1, A, 1, A, 1, IP, IP, B, 1,
  799. $ W, 1, INFO )
  800. CALL CHKXER( 'CSYSV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  801. INFOT = 5
  802. CALL CSYSV_AA_2STAGE( 'U', 2, 1, A, 1, A, 1, IP, IP, B, 1,
  803. $ W, 1, INFO )
  804. CALL CHKXER( 'CSYSV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  805. INFOT = 7
  806. CALL CSYSV_AA_2STAGE( 'U', 2, 1, A, 2, A, 1, IP, IP, B, 2,
  807. $ W, 1, INFO )
  808. CALL CHKXER( 'CSYSV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  809. INFOT = 11
  810. CALL CSYSV_AA_2STAGE( 'U', 2, 1, A, 2, A, 8, IP, IP, B, 1,
  811. $ W, 1, INFO )
  812. CALL CHKXER( 'CSYSV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  813. INFOT = 13
  814. CALL CSYSV_AA_2STAGE( 'U', 2, 1, A, 2, A, 8, IP, IP, B, 2,
  815. $ W, 1, INFO )
  816. CALL CHKXER( 'CSYSV_AA_2STAGE', INFOT, NOUT, LERR, OK )
  817. *
  818. ELSE IF( LSAMEN( 2, C2, 'HP' ) ) THEN
  819. *
  820. * CHPSV
  821. *
  822. SRNAMT = 'CHPSV '
  823. INFOT = 1
  824. CALL CHPSV( '/', 0, 0, A, IP, B, 1, INFO )
  825. CALL CHKXER( 'CHPSV ', INFOT, NOUT, LERR, OK )
  826. INFOT = 2
  827. CALL CHPSV( 'U', -1, 0, A, IP, B, 1, INFO )
  828. CALL CHKXER( 'CHPSV ', INFOT, NOUT, LERR, OK )
  829. INFOT = 3
  830. CALL CHPSV( 'U', 0, -1, A, IP, B, 1, INFO )
  831. CALL CHKXER( 'CHPSV ', INFOT, NOUT, LERR, OK )
  832. INFOT = 7
  833. CALL CHPSV( 'U', 2, 0, A, IP, B, 1, INFO )
  834. CALL CHKXER( 'CHPSV ', INFOT, NOUT, LERR, OK )
  835. *
  836. * CHPSVX
  837. *
  838. SRNAMT = 'CHPSVX'
  839. INFOT = 1
  840. CALL CHPSVX( '/', 'U', 0, 0, A, AF, IP, B, 1, X, 1, RCOND, R1,
  841. $ R2, W, RW, INFO )
  842. CALL CHKXER( 'CHPSVX', INFOT, NOUT, LERR, OK )
  843. INFOT = 2
  844. CALL CHPSVX( 'N', '/', 0, 0, A, AF, IP, B, 1, X, 1, RCOND, R1,
  845. $ R2, W, RW, INFO )
  846. CALL CHKXER( 'CHPSVX', INFOT, NOUT, LERR, OK )
  847. INFOT = 3
  848. CALL CHPSVX( 'N', 'U', -1, 0, A, AF, IP, B, 1, X, 1, RCOND, R1,
  849. $ R2, W, RW, INFO )
  850. CALL CHKXER( 'CHPSVX', INFOT, NOUT, LERR, OK )
  851. INFOT = 4
  852. CALL CHPSVX( 'N', 'U', 0, -1, A, AF, IP, B, 1, X, 1, RCOND, R1,
  853. $ R2, W, RW, INFO )
  854. CALL CHKXER( 'CHPSVX', INFOT, NOUT, LERR, OK )
  855. INFOT = 9
  856. CALL CHPSVX( 'N', 'U', 2, 0, A, AF, IP, B, 1, X, 2, RCOND, R1,
  857. $ R2, W, RW, INFO )
  858. CALL CHKXER( 'CHPSVX', INFOT, NOUT, LERR, OK )
  859. INFOT = 11
  860. CALL CHPSVX( 'N', 'U', 2, 0, A, AF, IP, B, 2, X, 1, RCOND, R1,
  861. $ R2, W, RW, INFO )
  862. CALL CHKXER( 'CHPSVX', INFOT, NOUT, LERR, OK )
  863. *
  864. ELSE IF( LSAMEN( 2, C2, 'SY' ) ) THEN
  865. *
  866. * CSYSV
  867. *
  868. SRNAMT = 'CSYSV '
  869. INFOT = 1
  870. CALL CSYSV( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
  871. CALL CHKXER( 'CSYSV ', INFOT, NOUT, LERR, OK )
  872. INFOT = 2
  873. CALL CSYSV( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
  874. CALL CHKXER( 'CSYSV ', INFOT, NOUT, LERR, OK )
  875. INFOT = 3
  876. CALL CSYSV( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
  877. CALL CHKXER( 'CSYSV ', INFOT, NOUT, LERR, OK )
  878. INFOT = 5
  879. CALL CSYSV( 'U', 2, 0, A, 1, IP, B, 2, W, 1, INFO )
  880. CALL CHKXER( 'CSYSV ', INFOT, NOUT, LERR, OK )
  881. INFOT = 8
  882. CALL CSYSV( 'U', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
  883. CALL CHKXER( 'CSYSV ', INFOT, NOUT, LERR, OK )
  884. INFOT = 10
  885. CALL CSYSV( 'U', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
  886. CALL CHKXER( 'CSYSV ', INFOT, NOUT, LERR, OK )
  887. INFOT = 10
  888. CALL CSYSV( 'U', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
  889. CALL CHKXER( 'CSYSV ', INFOT, NOUT, LERR, OK )
  890. *
  891. * CSYSVX
  892. *
  893. SRNAMT = 'CSYSVX'
  894. INFOT = 1
  895. CALL CSYSVX( '/', 'U', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1,
  896. $ RCOND, R1, R2, W, 1, RW, INFO )
  897. CALL CHKXER( 'CSYSVX', INFOT, NOUT, LERR, OK )
  898. INFOT = 2
  899. CALL CSYSVX( 'N', '/', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1,
  900. $ RCOND, R1, R2, W, 1, RW, INFO )
  901. CALL CHKXER( 'CSYSVX', INFOT, NOUT, LERR, OK )
  902. INFOT = 3
  903. CALL CSYSVX( 'N', 'U', -1, 0, A, 1, AF, 1, IP, B, 1, X, 1,
  904. $ RCOND, R1, R2, W, 1, RW, INFO )
  905. CALL CHKXER( 'CSYSVX', INFOT, NOUT, LERR, OK )
  906. INFOT = 4
  907. CALL CSYSVX( 'N', 'U', 0, -1, A, 1, AF, 1, IP, B, 1, X, 1,
  908. $ RCOND, R1, R2, W, 1, RW, INFO )
  909. CALL CHKXER( 'CSYSVX', INFOT, NOUT, LERR, OK )
  910. INFOT = 6
  911. CALL CSYSVX( 'N', 'U', 2, 0, A, 1, AF, 2, IP, B, 2, X, 2,
  912. $ RCOND, R1, R2, W, 4, RW, INFO )
  913. CALL CHKXER( 'CSYSVX', INFOT, NOUT, LERR, OK )
  914. INFOT = 8
  915. CALL CSYSVX( 'N', 'U', 2, 0, A, 2, AF, 1, IP, B, 2, X, 2,
  916. $ RCOND, R1, R2, W, 4, RW, INFO )
  917. CALL CHKXER( 'CSYSVX', INFOT, NOUT, LERR, OK )
  918. INFOT = 11
  919. CALL CSYSVX( 'N', 'U', 2, 0, A, 2, AF, 2, IP, B, 1, X, 2,
  920. $ RCOND, R1, R2, W, 4, RW, INFO )
  921. CALL CHKXER( 'CSYSVX', INFOT, NOUT, LERR, OK )
  922. INFOT = 13
  923. CALL CSYSVX( 'N', 'U', 2, 0, A, 2, AF, 2, IP, B, 2, X, 1,
  924. $ RCOND, R1, R2, W, 4, RW, INFO )
  925. CALL CHKXER( 'CSYSVX', INFOT, NOUT, LERR, OK )
  926. INFOT = 18
  927. CALL CSYSVX( 'N', 'U', 2, 0, A, 2, AF, 2, IP, B, 2, X, 2,
  928. $ RCOND, R1, R2, W, 3, RW, INFO )
  929. CALL CHKXER( 'CSYSVX', INFOT, NOUT, LERR, OK )
  930. *
  931. ELSE IF( LSAMEN( 2, C2, 'SR' ) ) THEN
  932. *
  933. * CSYSV_ROOK
  934. *
  935. SRNAMT = 'CSYSV_ROOK'
  936. INFOT = 1
  937. CALL CSYSV_ROOK( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
  938. CALL CHKXER( 'CSYSV_ROOK', INFOT, NOUT, LERR, OK )
  939. INFOT = 2
  940. CALL CSYSV_ROOK( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
  941. CALL CHKXER( 'CSYSV_ROOK', INFOT, NOUT, LERR, OK )
  942. INFOT = 3
  943. CALL CSYSV_ROOK( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
  944. CALL CHKXER( 'CSYSV_ROOK', INFOT, NOUT, LERR, OK )
  945. INFOT = 5
  946. CALL CSYSV_ROOK( 'U', 2, 0, A, 1, IP, B, 2, W, 1, INFO )
  947. CALL CHKXER( 'CSYSV_ROOK', INFOT, NOUT, LERR, OK )
  948. INFOT = 8
  949. CALL CSYSV_ROOK( 'U', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
  950. CALL CHKXER( 'CSYSV_ROOK', INFOT, NOUT, LERR, OK )
  951. INFOT = 10
  952. CALL CSYSV_ROOK( 'U', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
  953. CALL CHKXER( 'CSYSV_ROOK', INFOT, NOUT, LERR, OK )
  954. INFOT = 10
  955. CALL CSYSV_ROOK( 'U', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
  956. CALL CHKXER( 'CSYSV_ROOK', INFOT, NOUT, LERR, OK )
  957. *
  958. ELSE IF( LSAMEN( 2, C2, 'SK' ) ) THEN
  959. *
  960. * CSYSV_RK
  961. *
  962. * Test error exits of the driver that uses factorization
  963. * of a symmetric indefinite matrix with rook
  964. * (bounded Bunch-Kaufman) pivoting with the new storage
  965. * format for factors L ( or U) and D.
  966. *
  967. * L (or U) is stored in A, diagonal of D is stored on the
  968. * diagonal of A, subdiagonal of D is stored in a separate array E.
  969. *
  970. SRNAMT = 'CSYSV_RK'
  971. INFOT = 1
  972. CALL CSYSV_RK( '/', 0, 0, A, 1, E, IP, B, 1, W, 1, INFO )
  973. CALL CHKXER( 'CSYSV_RK', INFOT, NOUT, LERR, OK )
  974. INFOT = 2
  975. CALL CSYSV_RK( 'U', -1, 0, A, 1, E, IP, B, 1, W, 1, INFO )
  976. CALL CHKXER( 'CSYSV_RK', INFOT, NOUT, LERR, OK )
  977. INFOT = 3
  978. CALL CSYSV_RK( 'U', 0, -1, A, 1, E, IP, B, 1, W, 1, INFO )
  979. CALL CHKXER( 'CSYSV_RK', INFOT, NOUT, LERR, OK )
  980. INFOT = 5
  981. CALL CSYSV_RK( 'U', 2, 0, A, 1, E, IP, B, 2, W, 1, INFO )
  982. CALL CHKXER( 'CSYSV_RK', INFOT, NOUT, LERR, OK )
  983. INFOT = 9
  984. CALL CSYSV_RK( 'U', 2, 0, A, 2, E, IP, B, 1, W, 1, INFO )
  985. CALL CHKXER( 'CSYSV_RK', INFOT, NOUT, LERR, OK )
  986. INFOT = 11
  987. CALL CSYSV_RK( 'U', 0, 0, A, 1, E, IP, B, 1, W, 0, INFO )
  988. CALL CHKXER( 'CSYSV_RK', INFOT, NOUT, LERR, OK )
  989. INFOT = 11
  990. CALL CSYSV_RK( 'U', 0, 0, A, 1, E, IP, B, 1, W, -2, INFO )
  991. CALL CHKXER( 'CSYSV_RK', INFOT, NOUT, LERR, OK )
  992. *
  993. ELSE IF( LSAMEN( 2, C2, 'SP' ) ) THEN
  994. *
  995. * CSPSV
  996. *
  997. SRNAMT = 'CSPSV '
  998. INFOT = 1
  999. CALL CSPSV( '/', 0, 0, A, IP, B, 1, INFO )
  1000. CALL CHKXER( 'CSPSV ', INFOT, NOUT, LERR, OK )
  1001. INFOT = 2
  1002. CALL CSPSV( 'U', -1, 0, A, IP, B, 1, INFO )
  1003. CALL CHKXER( 'CSPSV ', INFOT, NOUT, LERR, OK )
  1004. INFOT = 3
  1005. CALL CSPSV( 'U', 0, -1, A, IP, B, 1, INFO )
  1006. CALL CHKXER( 'CSPSV ', INFOT, NOUT, LERR, OK )
  1007. INFOT = 7
  1008. CALL CSPSV( 'U', 2, 0, A, IP, B, 1, INFO )
  1009. CALL CHKXER( 'CSPSV ', INFOT, NOUT, LERR, OK )
  1010. *
  1011. * CSPSVX
  1012. *
  1013. SRNAMT = 'CSPSVX'
  1014. INFOT = 1
  1015. CALL CSPSVX( '/', 'U', 0, 0, A, AF, IP, B, 1, X, 1, RCOND, R1,
  1016. $ R2, W, RW, INFO )
  1017. CALL CHKXER( 'CSPSVX', INFOT, NOUT, LERR, OK )
  1018. INFOT = 2
  1019. CALL CSPSVX( 'N', '/', 0, 0, A, AF, IP, B, 1, X, 1, RCOND, R1,
  1020. $ R2, W, RW, INFO )
  1021. CALL CHKXER( 'CSPSVX', INFOT, NOUT, LERR, OK )
  1022. INFOT = 3
  1023. CALL CSPSVX( 'N', 'U', -1, 0, A, AF, IP, B, 1, X, 1, RCOND, R1,
  1024. $ R2, W, RW, INFO )
  1025. CALL CHKXER( 'CSPSVX', INFOT, NOUT, LERR, OK )
  1026. INFOT = 4
  1027. CALL CSPSVX( 'N', 'U', 0, -1, A, AF, IP, B, 1, X, 1, RCOND, R1,
  1028. $ R2, W, RW, INFO )
  1029. CALL CHKXER( 'CSPSVX', INFOT, NOUT, LERR, OK )
  1030. INFOT = 9
  1031. CALL CSPSVX( 'N', 'U', 2, 0, A, AF, IP, B, 1, X, 2, RCOND, R1,
  1032. $ R2, W, RW, INFO )
  1033. CALL CHKXER( 'CSPSVX', INFOT, NOUT, LERR, OK )
  1034. INFOT = 11
  1035. CALL CSPSVX( 'N', 'U', 2, 0, A, AF, IP, B, 2, X, 1, RCOND, R1,
  1036. $ R2, W, RW, INFO )
  1037. CALL CHKXER( 'CSPSVX', INFOT, NOUT, LERR, OK )
  1038. END IF
  1039. *
  1040. * Print a summary line.
  1041. *
  1042. IF( OK ) THEN
  1043. WRITE( NOUT, FMT = 9999 )PATH
  1044. ELSE
  1045. WRITE( NOUT, FMT = 9998 )PATH
  1046. END IF
  1047. *
  1048. 9999 FORMAT( 1X, A3, ' drivers passed the tests of the error exits' )
  1049. 9998 FORMAT( ' *** ', A3, ' drivers failed the tests of the error ',
  1050. $ 'exits ***' )
  1051. *
  1052. RETURN
  1053. *
  1054. * End of CERRVX
  1055. *
  1056. END