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.

cerrsyx.f 21 kB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609
  1. *> \brief \b CERRSYX
  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 CERRSY( 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. *> CERRSY tests the error exits for the COMPLEX routines
  25. *> for symmetric indefinite matrices.
  26. *>
  27. *> Note that this file is used only when the XBLAS are available,
  28. *> otherwise cerrsy.f defines this subroutine.
  29. *> \endverbatim
  30. *
  31. * Arguments:
  32. * ==========
  33. *
  34. *> \param[in] PATH
  35. *> \verbatim
  36. *> PATH is CHARACTER*3
  37. *> The LAPACK path name for the routines to be tested.
  38. *> \endverbatim
  39. *>
  40. *> \param[in] NUNIT
  41. *> \verbatim
  42. *> NUNIT is INTEGER
  43. *> The unit number for output.
  44. *> \endverbatim
  45. *
  46. * Authors:
  47. * ========
  48. *
  49. *> \author Univ. of Tennessee
  50. *> \author Univ. of California Berkeley
  51. *> \author Univ. of Colorado Denver
  52. *> \author NAG Ltd.
  53. *
  54. *> \ingroup complex_lin
  55. *
  56. * =====================================================================
  57. SUBROUTINE CERRSY( PATH, NUNIT )
  58. *
  59. * -- LAPACK test routine --
  60. * -- LAPACK is a software package provided by Univ. of Tennessee, --
  61. * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
  62. *
  63. * .. Scalar Arguments ..
  64. CHARACTER*3 PATH
  65. INTEGER NUNIT
  66. * ..
  67. *
  68. * =====================================================================
  69. *
  70. * .. Parameters ..
  71. INTEGER NMAX
  72. PARAMETER ( NMAX = 4 )
  73. * ..
  74. * .. Local Scalars ..
  75. CHARACTER EQ
  76. CHARACTER*2 C2
  77. INTEGER I, INFO, J, N_ERR_BNDS, NPARAMS
  78. REAL ANRM, RCOND, BERR
  79. * ..
  80. * .. Local Arrays ..
  81. INTEGER IP( NMAX )
  82. REAL R( NMAX ), R1( NMAX ), R2( NMAX ),
  83. $ S( NMAX ), ERR_BNDS_N( NMAX, 3 ),
  84. $ ERR_BNDS_C( NMAX, 3 ), PARAMS( 1 )
  85. COMPLEX A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
  86. $ E( NMAX), W( 2*NMAX ), X( NMAX )
  87. * ..
  88. * .. External Functions ..
  89. LOGICAL LSAMEN
  90. EXTERNAL LSAMEN
  91. * ..
  92. * .. External Subroutines ..
  93. EXTERNAL ALAESM, CHKXER, CSPCON, CSPRFS, CSPTRF, CSPTRI,
  94. $ CSPTRS, CSYCON, CSYRFS, CSYTF2, CSYTRF, CSYTRI,
  95. $ CSYTRI2, CSYTRS, CSYRFSX, CSYCON_ROOK,
  96. $ CSYTF2_ROOK, CSYTRF_ROOK, CSYTRI_ROOK,
  97. $ CSYTRS_ROOK
  98. * ..
  99. * .. Scalars in Common ..
  100. LOGICAL LERR, OK
  101. CHARACTER*32 SRNAMT
  102. INTEGER INFOT, NOUT
  103. * ..
  104. * .. Common blocks ..
  105. COMMON / INFOC / INFOT, NOUT, OK, LERR
  106. COMMON / SRNAMC / SRNAMT
  107. * ..
  108. * .. Intrinsic Functions ..
  109. INTRINSIC CMPLX, REAL
  110. * ..
  111. * .. Executable Statements ..
  112. *
  113. NOUT = NUNIT
  114. WRITE( NOUT, FMT = * )
  115. C2 = PATH( 2: 3 )
  116. *
  117. * Set the variables to innocuous values.
  118. *
  119. DO 20 J = 1, NMAX
  120. DO 10 I = 1, NMAX
  121. A( I, J ) = CMPLX( 1. / REAL( I+J ), -1. / REAL( I+J ) )
  122. AF( I, J ) = CMPLX( 1. / REAL( I+J ), -1. / REAL( I+J ) )
  123. 10 CONTINUE
  124. B( J ) = 0.E0
  125. E( J ) = 0.E0
  126. R1( J ) = 0.E0
  127. R2( J ) = 0.E0
  128. W( J ) = 0.E0
  129. X( J ) = 0.E0
  130. IP( J ) = J
  131. 20 CONTINUE
  132. ANRM = 1.0
  133. OK = .TRUE.
  134. IF( LSAMEN( 2, C2, 'SY' ) ) THEN
  135. *
  136. * Test error exits of the routines that use factorization
  137. * of a symmetric indefinite matrix with partial
  138. * (Bunch-Kaufman) diagonal pivoting method.
  139. *
  140. * CSYTRF
  141. *
  142. SRNAMT = 'CSYTRF'
  143. INFOT = 1
  144. CALL CSYTRF( '/', 0, A, 1, IP, W, 1, INFO )
  145. CALL CHKXER( 'CSYTRF', INFOT, NOUT, LERR, OK )
  146. INFOT = 2
  147. CALL CSYTRF( 'U', -1, A, 1, IP, W, 1, INFO )
  148. CALL CHKXER( 'CSYTRF', INFOT, NOUT, LERR, OK )
  149. INFOT = 4
  150. CALL CSYTRF( 'U', 2, A, 1, IP, W, 4, INFO )
  151. CALL CHKXER( 'CSYTRF', INFOT, NOUT, LERR, OK )
  152. INFOT = 7
  153. CALL CSYTRF( 'U', 0, A, 1, IP, W, 0, INFO )
  154. CALL CHKXER( 'CSYTRF', INFOT, NOUT, LERR, OK )
  155. INFOT = 7
  156. CALL CSYTRF( 'U', 0, A, 1, IP, W, -2, INFO )
  157. CALL CHKXER( 'CSYTRF', INFOT, NOUT, LERR, OK )
  158. *
  159. * CSYTF2
  160. *
  161. SRNAMT = 'CSYTF2'
  162. INFOT = 1
  163. CALL CSYTF2( '/', 0, A, 1, IP, INFO )
  164. CALL CHKXER( 'CSYTF2', INFOT, NOUT, LERR, OK )
  165. INFOT = 2
  166. CALL CSYTF2( 'U', -1, A, 1, IP, INFO )
  167. CALL CHKXER( 'CSYTF2', INFOT, NOUT, LERR, OK )
  168. INFOT = 4
  169. CALL CSYTF2( 'U', 2, A, 1, IP, INFO )
  170. CALL CHKXER( 'CSYTF2', INFOT, NOUT, LERR, OK )
  171. *
  172. * CSYTRI
  173. *
  174. SRNAMT = 'CSYTRI'
  175. INFOT = 1
  176. CALL CSYTRI( '/', 0, A, 1, IP, W, INFO )
  177. CALL CHKXER( 'CSYTRI', INFOT, NOUT, LERR, OK )
  178. INFOT = 2
  179. CALL CSYTRI( 'U', -1, A, 1, IP, W, INFO )
  180. CALL CHKXER( 'CSYTRI', INFOT, NOUT, LERR, OK )
  181. INFOT = 4
  182. CALL CSYTRI( 'U', 2, A, 1, IP, W, INFO )
  183. CALL CHKXER( 'CSYTRI', INFOT, NOUT, LERR, OK )
  184. *
  185. * CSYTRI2
  186. *
  187. SRNAMT = 'CSYTRI2'
  188. INFOT = 1
  189. CALL CSYTRI2( '/', 0, A, 1, IP, W, 1, INFO )
  190. CALL CHKXER( 'CSYTRI2', INFOT, NOUT, LERR, OK )
  191. INFOT = 2
  192. CALL CSYTRI2( 'U', -1, A, 1, IP, W, 1, INFO )
  193. CALL CHKXER( 'CSYTRI2', INFOT, NOUT, LERR, OK )
  194. INFOT = 4
  195. CALL CSYTRI2( 'U', 2, A, 1, IP, W, 1, INFO )
  196. CALL CHKXER( 'CSYTRI2', INFOT, NOUT, LERR, OK )
  197. *
  198. * CSYTRI2X
  199. *
  200. SRNAMT = 'CSYTRI2X'
  201. INFOT = 1
  202. CALL CSYTRI2X( '/', 0, A, 1, IP, W, 1, INFO )
  203. CALL CHKXER( 'CSYTRI2X', INFOT, NOUT, LERR, OK )
  204. INFOT = 2
  205. CALL CSYTRI2X( 'U', -1, A, 1, IP, W, 1, INFO )
  206. CALL CHKXER( 'CSYTRI2X', INFOT, NOUT, LERR, OK )
  207. INFOT = 4
  208. CALL CSYTRI2X( 'U', 2, A, 1, IP, W, 1, INFO )
  209. CALL CHKXER( 'CSYTRI2X', INFOT, NOUT, LERR, OK )
  210. *
  211. * CSYTRS
  212. *
  213. SRNAMT = 'CSYTRS'
  214. INFOT = 1
  215. CALL CSYTRS( '/', 0, 0, A, 1, IP, B, 1, INFO )
  216. CALL CHKXER( 'CSYTRS', INFOT, NOUT, LERR, OK )
  217. INFOT = 2
  218. CALL CSYTRS( 'U', -1, 0, A, 1, IP, B, 1, INFO )
  219. CALL CHKXER( 'CSYTRS', INFOT, NOUT, LERR, OK )
  220. INFOT = 3
  221. CALL CSYTRS( 'U', 0, -1, A, 1, IP, B, 1, INFO )
  222. CALL CHKXER( 'CSYTRS', INFOT, NOUT, LERR, OK )
  223. INFOT = 5
  224. CALL CSYTRS( 'U', 2, 1, A, 1, IP, B, 2, INFO )
  225. CALL CHKXER( 'CSYTRS', INFOT, NOUT, LERR, OK )
  226. INFOT = 8
  227. CALL CSYTRS( 'U', 2, 1, A, 2, IP, B, 1, INFO )
  228. CALL CHKXER( 'CSYTRS', INFOT, NOUT, LERR, OK )
  229. *
  230. * CSYRFS
  231. *
  232. SRNAMT = 'CSYRFS'
  233. INFOT = 1
  234. CALL CSYRFS( '/', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1, R1, R2, W,
  235. $ R, INFO )
  236. CALL CHKXER( 'CSYRFS', INFOT, NOUT, LERR, OK )
  237. INFOT = 2
  238. CALL CSYRFS( 'U', -1, 0, A, 1, AF, 1, IP, B, 1, X, 1, R1, R2,
  239. $ W, R, INFO )
  240. CALL CHKXER( 'CSYRFS', INFOT, NOUT, LERR, OK )
  241. INFOT = 3
  242. CALL CSYRFS( 'U', 0, -1, A, 1, AF, 1, IP, B, 1, X, 1, R1, R2,
  243. $ W, R, INFO )
  244. CALL CHKXER( 'CSYRFS', INFOT, NOUT, LERR, OK )
  245. INFOT = 5
  246. CALL CSYRFS( 'U', 2, 1, A, 1, AF, 2, IP, B, 2, X, 2, R1, R2, W,
  247. $ R, INFO )
  248. CALL CHKXER( 'CSYRFS', INFOT, NOUT, LERR, OK )
  249. INFOT = 7
  250. CALL CSYRFS( 'U', 2, 1, A, 2, AF, 1, IP, B, 2, X, 2, R1, R2, W,
  251. $ R, INFO )
  252. CALL CHKXER( 'CSYRFS', INFOT, NOUT, LERR, OK )
  253. INFOT = 10
  254. CALL CSYRFS( 'U', 2, 1, A, 2, AF, 2, IP, B, 1, X, 2, R1, R2, W,
  255. $ R, INFO )
  256. CALL CHKXER( 'CSYRFS', INFOT, NOUT, LERR, OK )
  257. INFOT = 12
  258. CALL CSYRFS( 'U', 2, 1, A, 2, AF, 2, IP, B, 2, X, 1, R1, R2, W,
  259. $ R, INFO )
  260. CALL CHKXER( 'CSYRFS', INFOT, NOUT, LERR, OK )
  261. *
  262. * CSYRFSX
  263. *
  264. N_ERR_BNDS = 3
  265. NPARAMS = 0
  266. SRNAMT = 'CSYRFSX'
  267. INFOT = 1
  268. CALL CSYRFSX( '/', EQ, 0, 0, A, 1, AF, 1, IP, S, B, 1, X, 1,
  269. $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
  270. $ PARAMS, W, R, INFO )
  271. CALL CHKXER( 'CSYRFSX', INFOT, NOUT, LERR, OK )
  272. INFOT = 2
  273. CALL CSYRFSX( 'U', EQ, -1, 0, A, 1, AF, 1, IP, S, B, 1, X, 1,
  274. $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
  275. $ PARAMS, W, R, INFO )
  276. CALL CHKXER( 'CSYRFSX', INFOT, NOUT, LERR, OK )
  277. EQ = 'N'
  278. INFOT = 3
  279. CALL CSYRFSX( 'U', EQ, -1, 0, A, 1, AF, 1, IP, S, B, 1, X, 1,
  280. $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
  281. $ PARAMS, W, R, INFO )
  282. CALL CHKXER( 'CSYRFSX', INFOT, NOUT, LERR, OK )
  283. INFOT = 4
  284. CALL CSYRFSX( 'U', EQ, 0, -1, A, 1, AF, 1, IP, S, B, 1, X, 1,
  285. $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
  286. $ PARAMS, W, R, INFO )
  287. CALL CHKXER( 'CSYRFSX', INFOT, NOUT, LERR, OK )
  288. INFOT = 6
  289. CALL CSYRFSX( 'U', EQ, 2, 1, A, 1, AF, 2, IP, S, B, 2, X, 2,
  290. $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
  291. $ PARAMS, W, R, INFO )
  292. CALL CHKXER( 'CSYRFSX', INFOT, NOUT, LERR, OK )
  293. INFOT = 8
  294. CALL CSYRFSX( 'U', EQ, 2, 1, A, 2, AF, 1, IP, S, B, 2, X, 2,
  295. $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
  296. $ PARAMS, W, R, INFO )
  297. CALL CHKXER( 'CSYRFSX', INFOT, NOUT, LERR, OK )
  298. INFOT = 12
  299. CALL CSYRFSX( 'U', EQ, 2, 1, A, 2, AF, 2, IP, S, B, 1, X, 2,
  300. $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
  301. $ PARAMS, W, R, INFO )
  302. CALL CHKXER( 'CSYRFSX', INFOT, NOUT, LERR, OK )
  303. INFOT = 14
  304. CALL CSYRFSX( 'U', EQ, 2, 1, A, 2, AF, 2, IP, S, B, 2, X, 1,
  305. $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
  306. $ PARAMS, W, R, INFO )
  307. CALL CHKXER( 'CSYRFSX', INFOT, NOUT, LERR, OK )
  308. *
  309. * CSYCON
  310. *
  311. SRNAMT = 'CSYCON'
  312. INFOT = 1
  313. CALL CSYCON( '/', 0, A, 1, IP, ANRM, RCOND, W, INFO )
  314. CALL CHKXER( 'CSYCON', INFOT, NOUT, LERR, OK )
  315. INFOT = 2
  316. CALL CSYCON( 'U', -1, A, 1, IP, ANRM, RCOND, W, INFO )
  317. CALL CHKXER( 'CSYCON', INFOT, NOUT, LERR, OK )
  318. INFOT = 4
  319. CALL CSYCON( 'U', 2, A, 1, IP, ANRM, RCOND, W, INFO )
  320. CALL CHKXER( 'CSYCON', INFOT, NOUT, LERR, OK )
  321. INFOT = 6
  322. CALL CSYCON( 'U', 1, A, 1, IP, -ANRM, RCOND, W, INFO )
  323. CALL CHKXER( 'CSYCON', INFOT, NOUT, LERR, OK )
  324. *
  325. ELSE IF( LSAMEN( 2, C2, 'SR' ) ) THEN
  326. *
  327. * Test error exits of the routines that use factorization
  328. * of a symmetric indefinite matrix with rook
  329. * (bounded Bunch-Kaufman) diagonal pivoting method.
  330. *
  331. * CSYTRF_ROOK
  332. *
  333. SRNAMT = 'CSYTRF_ROOK'
  334. INFOT = 1
  335. CALL CSYTRF_ROOK( '/', 0, A, 1, IP, W, 1, INFO )
  336. CALL CHKXER( 'CSYTRF_ROOK', INFOT, NOUT, LERR, OK )
  337. INFOT = 2
  338. CALL CSYTRF_ROOK( 'U', -1, A, 1, IP, W, 1, INFO )
  339. CALL CHKXER( 'CSYTRF_ROOK', INFOT, NOUT, LERR, OK )
  340. INFOT = 4
  341. CALL CSYTRF_ROOK( 'U', 2, A, 1, IP, W, 4, INFO )
  342. CALL CHKXER( 'CSYTRF_ROOK', INFOT, NOUT, LERR, OK )
  343. INFOT = 7
  344. CALL CSYTRF_ROOK( 'U', 0, A, 1, IP, W, 0, INFO )
  345. CALL CHKXER( 'CSYTRF_ROOK', INFOT, NOUT, LERR, OK )
  346. INFOT = 7
  347. CALL CSYTRF_ROOK( 'U', 0, A, 1, IP, W, -2, INFO )
  348. CALL CHKXER( 'CSYTRF_ROOK', INFOT, NOUT, LERR, OK )
  349. *
  350. * CSYTF2_ROOK
  351. *
  352. SRNAMT = 'CSYTF2_ROOK'
  353. INFOT = 1
  354. CALL CSYTF2_ROOK( '/', 0, A, 1, IP, INFO )
  355. CALL CHKXER( 'CSYTF2_ROOK', INFOT, NOUT, LERR, OK )
  356. INFOT = 2
  357. CALL CSYTF2_ROOK( 'U', -1, A, 1, IP, INFO )
  358. CALL CHKXER( 'CSYTF2_ROOK', INFOT, NOUT, LERR, OK )
  359. INFOT = 4
  360. CALL CSYTF2_ROOK( 'U', 2, A, 1, IP, INFO )
  361. CALL CHKXER( 'CSYTF2_ROOK', INFOT, NOUT, LERR, OK )
  362. *
  363. * CSYTRI_ROOK
  364. *
  365. SRNAMT = 'CSYTRI_ROOK'
  366. INFOT = 1
  367. CALL CSYTRI_ROOK( '/', 0, A, 1, IP, W, INFO )
  368. CALL CHKXER( 'CSYTRI_ROOK', INFOT, NOUT, LERR, OK )
  369. INFOT = 2
  370. CALL CSYTRI_ROOK( 'U', -1, A, 1, IP, W, INFO )
  371. CALL CHKXER( 'CSYTRI_ROOK', INFOT, NOUT, LERR, OK )
  372. INFOT = 4
  373. CALL CSYTRI_ROOK( 'U', 2, A, 1, IP, W, INFO )
  374. CALL CHKXER( 'CSYTRI_ROOK', INFOT, NOUT, LERR, OK )
  375. *
  376. * CSYTRS_ROOK
  377. *
  378. SRNAMT = 'CSYTRS_ROOK'
  379. INFOT = 1
  380. CALL CSYTRS_ROOK( '/', 0, 0, A, 1, IP, B, 1, INFO )
  381. CALL CHKXER( 'CSYTRS_ROOK', INFOT, NOUT, LERR, OK )
  382. INFOT = 2
  383. CALL CSYTRS_ROOK( 'U', -1, 0, A, 1, IP, B, 1, INFO )
  384. CALL CHKXER( 'CSYTRS_ROOK', INFOT, NOUT, LERR, OK )
  385. INFOT = 3
  386. CALL CSYTRS_ROOK( 'U', 0, -1, A, 1, IP, B, 1, INFO )
  387. CALL CHKXER( 'CSYTRS_ROOK', INFOT, NOUT, LERR, OK )
  388. INFOT = 5
  389. CALL CSYTRS_ROOK( 'U', 2, 1, A, 1, IP, B, 2, INFO )
  390. CALL CHKXER( 'CSYTRS_ROOK', INFOT, NOUT, LERR, OK )
  391. INFOT = 8
  392. CALL CSYTRS_ROOK( 'U', 2, 1, A, 2, IP, B, 1, INFO )
  393. CALL CHKXER( 'CSYTRS_ROOK', INFOT, NOUT, LERR, OK )
  394. *
  395. * CSYCON_ROOK
  396. *
  397. SRNAMT = 'CSYCON_ROOK'
  398. INFOT = 1
  399. CALL CSYCON_ROOK( '/', 0, A, 1, IP, ANRM, RCOND, W, INFO )
  400. CALL CHKXER( 'CSYCON_ROOK', INFOT, NOUT, LERR, OK )
  401. INFOT = 2
  402. CALL CSYCON_ROOK( 'U', -1, A, 1, IP, ANRM, RCOND, W, INFO )
  403. CALL CHKXER( 'CSYCON_ROOK', INFOT, NOUT, LERR, OK )
  404. INFOT = 4
  405. CALL CSYCON_ROOK( 'U', 2, A, 1, IP, ANRM, RCOND, W, INFO )
  406. CALL CHKXER( 'CSYCON_ROOK', INFOT, NOUT, LERR, OK )
  407. INFOT = 6
  408. CALL CSYCON_ROOK( 'U', 1, A, 1, IP, -ANRM, RCOND, W, INFO )
  409. CALL CHKXER( 'CSYCON_ROOK', INFOT, NOUT, LERR, OK )
  410. *
  411. ELSE IF( LSAMEN( 2, C2, 'SK' ) ) THEN
  412. *
  413. * Test error exits of the routines that use factorization
  414. * of a symmetric indefinite matrix with rook
  415. * (bounded Bunch-Kaufman) pivoting with the new storage
  416. * format for factors L ( or U) and D.
  417. *
  418. * L (or U) is stored in A, diagonal of D is stored on the
  419. * diagonal of A, subdiagonal of D is stored in a separate array E.
  420. *
  421. * CSYTRF_RK
  422. *
  423. SRNAMT = 'CSYTRF_RK'
  424. INFOT = 1
  425. CALL CSYTRF_RK( '/', 0, A, 1, E, IP, W, 1, INFO )
  426. CALL CHKXER( 'CSYTRF_RK', INFOT, NOUT, LERR, OK )
  427. INFOT = 2
  428. CALL CSYTRF_RK( 'U', -1, A, 1, E, IP, W, 1, INFO )
  429. CALL CHKXER( 'CSYTRF_RK', INFOT, NOUT, LERR, OK )
  430. INFOT = 4
  431. CALL CSYTRF_RK( 'U', 2, A, 1, E, IP, W, 4, INFO )
  432. CALL CHKXER( 'CSYTRF_RK', INFOT, NOUT, LERR, OK )
  433. INFOT = 8
  434. CALL CSYTRF_RK( 'U', 0, A, 1, E, IP, W, 0, INFO )
  435. CALL CHKXER( 'CSYTRF_RK', INFOT, NOUT, LERR, OK )
  436. INFOT = 8
  437. CALL CSYTRF_RK( 'U', 0, A, 1, E, IP, W, -2, INFO )
  438. CALL CHKXER( 'CSYTRF_RK', INFOT, NOUT, LERR, OK )
  439. *
  440. * CSYTF2_RK
  441. *
  442. SRNAMT = 'CSYTF2_RK'
  443. INFOT = 1
  444. CALL CSYTF2_RK( '/', 0, A, 1, E, IP, INFO )
  445. CALL CHKXER( 'CSYTF2_RK', INFOT, NOUT, LERR, OK )
  446. INFOT = 2
  447. CALL CSYTF2_RK( 'U', -1, A, 1, E, IP, INFO )
  448. CALL CHKXER( 'CSYTF2_RK', INFOT, NOUT, LERR, OK )
  449. INFOT = 4
  450. CALL CSYTF2_RK( 'U', 2, A, 1, E, IP, INFO )
  451. CALL CHKXER( 'CSYTF2_RK', INFOT, NOUT, LERR, OK )
  452. *
  453. * CSYTRI_3
  454. *
  455. SRNAMT = 'CSYTRI_3'
  456. INFOT = 1
  457. CALL CSYTRI_3( '/', 0, A, 1, E, IP, W, 1, INFO )
  458. CALL CHKXER( 'CSYTRI_3', INFOT, NOUT, LERR, OK )
  459. INFOT = 2
  460. CALL CSYTRI_3( 'U', -1, A, 1, E, IP, W, 1, INFO )
  461. CALL CHKXER( 'CSYTRI_3', INFOT, NOUT, LERR, OK )
  462. INFOT = 4
  463. CALL CSYTRI_3( 'U', 2, A, 1, E, IP, W, 1, INFO )
  464. CALL CHKXER( 'CSYTRI_3', INFOT, NOUT, LERR, OK )
  465. INFOT = 8
  466. CALL CSYTRI_3( 'U', 0, A, 1, E, IP, W, 0, INFO )
  467. CALL CHKXER( 'CSYTRI_3', INFOT, NOUT, LERR, OK )
  468. INFOT = 8
  469. CALL CSYTRI_3( 'U', 0, A, 1, E, IP, W, -2, INFO )
  470. CALL CHKXER( 'CSYTRI_3', INFOT, NOUT, LERR, OK )
  471. *
  472. * CSYTRI_3X
  473. *
  474. SRNAMT = 'CSYTRI_3X'
  475. INFOT = 1
  476. CALL CSYTRI_3X( '/', 0, A, 1, E, IP, W, 1, INFO )
  477. CALL CHKXER( 'CSYTRI_3X', INFOT, NOUT, LERR, OK )
  478. INFOT = 2
  479. CALL CSYTRI_3X( 'U', -1, A, 1, E, IP, W, 1, INFO )
  480. CALL CHKXER( 'CSYTRI_3X', INFOT, NOUT, LERR, OK )
  481. INFOT = 4
  482. CALL CSYTRI_3X( 'U', 2, A, 1, E, IP, W, 1, INFO )
  483. CALL CHKXER( 'CSYTRI_3X', INFOT, NOUT, LERR, OK )
  484. *
  485. * CSYTRS_3
  486. *
  487. SRNAMT = 'CSYTRS_3'
  488. INFOT = 1
  489. CALL CSYTRS_3( '/', 0, 0, A, 1, E, IP, B, 1, INFO )
  490. CALL CHKXER( 'CSYTRS_3', INFOT, NOUT, LERR, OK )
  491. INFOT = 2
  492. CALL CSYTRS_3( 'U', -1, 0, A, 1, E, IP, B, 1, INFO )
  493. CALL CHKXER( 'CSYTRS_3', INFOT, NOUT, LERR, OK )
  494. INFOT = 3
  495. CALL CSYTRS_3( 'U', 0, -1, A, 1, E, IP, B, 1, INFO )
  496. CALL CHKXER( 'CSYTRS_3', INFOT, NOUT, LERR, OK )
  497. INFOT = 5
  498. CALL CSYTRS_3( 'U', 2, 1, A, 1, E, IP, B, 2, INFO )
  499. CALL CHKXER( 'CSYTRS_3', INFOT, NOUT, LERR, OK )
  500. INFOT = 9
  501. CALL CSYTRS_3( 'U', 2, 1, A, 2, E, IP, B, 1, INFO )
  502. CALL CHKXER( 'CSYTRS_3', INFOT, NOUT, LERR, OK )
  503. *
  504. * CSYCON_3
  505. *
  506. SRNAMT = 'CSYCON_3'
  507. INFOT = 1
  508. CALL CSYCON_3( '/', 0, A, 1, E, IP, ANRM, RCOND, W, INFO )
  509. CALL CHKXER( 'CSYCON_3', INFOT, NOUT, LERR, OK )
  510. INFOT = 2
  511. CALL CSYCON_3( 'U', -1, A, 1, E, IP, ANRM, RCOND, W, INFO )
  512. CALL CHKXER( 'CSYCON_3', INFOT, NOUT, LERR, OK )
  513. INFOT = 4
  514. CALL CSYCON_3( 'U', 2, A, 1, E, IP, ANRM, RCOND, W, INFO )
  515. CALL CHKXER( 'CSYCON_3', INFOT, NOUT, LERR, OK )
  516. INFOT = 7
  517. CALL CSYCON_3( 'U', 1, A, 1, E, IP, -1.0E0, RCOND, W, INFO)
  518. CALL CHKXER( 'CSYCON_3', INFOT, NOUT, LERR, OK )
  519. *
  520. ELSE IF( LSAMEN( 2, C2, 'SP' ) ) THEN
  521. *
  522. * Test error exits of the routines that use factorization
  523. * of a symmetric indefinite packed matrix with partial
  524. * (Bunch-Kaufman) diagonal pivoting method.
  525. *
  526. * CSPTRF
  527. *
  528. SRNAMT = 'CSPTRF'
  529. INFOT = 1
  530. CALL CSPTRF( '/', 0, A, IP, INFO )
  531. CALL CHKXER( 'CSPTRF', INFOT, NOUT, LERR, OK )
  532. INFOT = 2
  533. CALL CSPTRF( 'U', -1, A, IP, INFO )
  534. CALL CHKXER( 'CSPTRF', INFOT, NOUT, LERR, OK )
  535. *
  536. * CSPTRI
  537. *
  538. SRNAMT = 'CSPTRI'
  539. INFOT = 1
  540. CALL CSPTRI( '/', 0, A, IP, W, INFO )
  541. CALL CHKXER( 'CSPTRI', INFOT, NOUT, LERR, OK )
  542. INFOT = 2
  543. CALL CSPTRI( 'U', -1, A, IP, W, INFO )
  544. CALL CHKXER( 'CSPTRI', INFOT, NOUT, LERR, OK )
  545. *
  546. * CSPTRS
  547. *
  548. SRNAMT = 'CSPTRS'
  549. INFOT = 1
  550. CALL CSPTRS( '/', 0, 0, A, IP, B, 1, INFO )
  551. CALL CHKXER( 'CSPTRS', INFOT, NOUT, LERR, OK )
  552. INFOT = 2
  553. CALL CSPTRS( 'U', -1, 0, A, IP, B, 1, INFO )
  554. CALL CHKXER( 'CSPTRS', INFOT, NOUT, LERR, OK )
  555. INFOT = 3
  556. CALL CSPTRS( 'U', 0, -1, A, IP, B, 1, INFO )
  557. CALL CHKXER( 'CSPTRS', INFOT, NOUT, LERR, OK )
  558. INFOT = 7
  559. CALL CSPTRS( 'U', 2, 1, A, IP, B, 1, INFO )
  560. CALL CHKXER( 'CSPTRS', INFOT, NOUT, LERR, OK )
  561. *
  562. * CSPRFS
  563. *
  564. SRNAMT = 'CSPRFS'
  565. INFOT = 1
  566. CALL CSPRFS( '/', 0, 0, A, AF, IP, B, 1, X, 1, R1, R2, W, R,
  567. $ INFO )
  568. CALL CHKXER( 'CSPRFS', INFOT, NOUT, LERR, OK )
  569. INFOT = 2
  570. CALL CSPRFS( 'U', -1, 0, A, AF, IP, B, 1, X, 1, R1, R2, W, R,
  571. $ INFO )
  572. CALL CHKXER( 'CSPRFS', INFOT, NOUT, LERR, OK )
  573. INFOT = 3
  574. CALL CSPRFS( 'U', 0, -1, A, AF, IP, B, 1, X, 1, R1, R2, W, R,
  575. $ INFO )
  576. CALL CHKXER( 'CSPRFS', INFOT, NOUT, LERR, OK )
  577. INFOT = 8
  578. CALL CSPRFS( 'U', 2, 1, A, AF, IP, B, 1, X, 2, R1, R2, W, R,
  579. $ INFO )
  580. CALL CHKXER( 'CSPRFS', INFOT, NOUT, LERR, OK )
  581. INFOT = 10
  582. CALL CSPRFS( 'U', 2, 1, A, AF, IP, B, 2, X, 1, R1, R2, W, R,
  583. $ INFO )
  584. CALL CHKXER( 'CSPRFS', INFOT, NOUT, LERR, OK )
  585. *
  586. * CSPCON
  587. *
  588. SRNAMT = 'CSPCON'
  589. INFOT = 1
  590. CALL CSPCON( '/', 0, A, IP, ANRM, RCOND, W, INFO )
  591. CALL CHKXER( 'CSPCON', INFOT, NOUT, LERR, OK )
  592. INFOT = 2
  593. CALL CSPCON( 'U', -1, A, IP, ANRM, RCOND, W, INFO )
  594. CALL CHKXER( 'CSPCON', INFOT, NOUT, LERR, OK )
  595. INFOT = 5
  596. CALL CSPCON( 'U', 1, A, IP, -ANRM, RCOND, W, INFO )
  597. CALL CHKXER( 'CSPCON', INFOT, NOUT, LERR, OK )
  598. END IF
  599. *
  600. * Print a summary line.
  601. *
  602. CALL ALAESM( PATH, OK, NOUT )
  603. *
  604. RETURN
  605. *
  606. * End of CERRSYX
  607. *
  608. END