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.

zerred.f 19 kB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518
  1. *> \brief \b ZERRED
  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 ZERRED( 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. *> ZERRED tests the error exits for the eigenvalue driver routines for
  25. *> DOUBLE COMPLEX PRECISION matrices:
  26. *>
  27. *> PATH driver description
  28. *> ---- ------ -----------
  29. *> ZEV ZGEEV find eigenvalues/eigenvectors for nonsymmetric A
  30. *> ZES ZGEES find eigenvalues/Schur form for nonsymmetric A
  31. *> ZVX ZGEEVX ZGEEV + balancing and condition estimation
  32. *> ZSX ZGEESX ZGEES + balancing and condition estimation
  33. *> ZBD ZGESVD compute SVD of an M-by-N matrix A
  34. *> ZGESDD compute SVD of an M-by-N matrix A(by divide and
  35. *> conquer)
  36. *> ZGEJSV compute SVD of an M-by-N matrix A where M >= N
  37. *> ZGESVDX compute SVD of an M-by-N matrix A(by bisection
  38. *> and inverse iteration)
  39. *> \endverbatim
  40. *
  41. * Arguments:
  42. * ==========
  43. *
  44. *> \param[in] PATH
  45. *> \verbatim
  46. *> PATH is CHARACTER*3
  47. *> The LAPACK path name for the routines to be tested.
  48. *> \endverbatim
  49. *>
  50. *> \param[in] NUNIT
  51. *> \verbatim
  52. *> NUNIT is INTEGER
  53. *> The unit number for output.
  54. *> \endverbatim
  55. *
  56. * Authors:
  57. * ========
  58. *
  59. *> \author Univ. of Tennessee
  60. *> \author Univ. of California Berkeley
  61. *> \author Univ. of Colorado Denver
  62. *> \author NAG Ltd.
  63. *
  64. *> \date June 2016
  65. *
  66. *> \ingroup complex16_eig
  67. *
  68. * =====================================================================
  69. SUBROUTINE ZERRED( PATH, NUNIT )
  70. *
  71. * -- LAPACK test routine (version 3.7.0) --
  72. * -- LAPACK is a software package provided by Univ. of Tennessee, --
  73. * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
  74. * June 2016
  75. *
  76. * .. Scalar Arguments ..
  77. CHARACTER*3 PATH
  78. INTEGER NUNIT
  79. * ..
  80. *
  81. * =====================================================================
  82. *
  83. * .. Parameters ..
  84. INTEGER NMAX, LW
  85. PARAMETER ( NMAX = 4, LW = 5*NMAX )
  86. DOUBLE PRECISION ONE, ZERO
  87. PARAMETER ( ONE = 1.0D0, ZERO = 0.0D0 )
  88. * ..
  89. * .. Local Scalars ..
  90. CHARACTER*2 C2
  91. INTEGER I, IHI, ILO, INFO, J, NS, NT, SDIM
  92. DOUBLE PRECISION ABNRM
  93. * ..
  94. * .. Local Arrays ..
  95. LOGICAL B( NMAX )
  96. INTEGER IW( 4*NMAX )
  97. DOUBLE PRECISION R1( NMAX ), R2( NMAX ), RW( LW ), S( NMAX )
  98. COMPLEX*16 A( NMAX, NMAX ), U( NMAX, NMAX ),
  99. $ VL( NMAX, NMAX ), VR( NMAX, NMAX ),
  100. $ VT( NMAX, NMAX ), W( 10*NMAX ), X( NMAX )
  101. * ..
  102. * .. External Subroutines ..
  103. EXTERNAL CHKXER, ZGEES, ZGEESX, ZGEEV, ZGEEVX, ZGESVJ,
  104. $ ZGESDD, ZGESVD
  105. * ..
  106. * .. External Functions ..
  107. LOGICAL LSAMEN, ZSLECT
  108. EXTERNAL LSAMEN, ZSLECT
  109. * ..
  110. * .. Intrinsic Functions ..
  111. INTRINSIC LEN_TRIM
  112. * ..
  113. * .. Arrays in Common ..
  114. LOGICAL SELVAL( 20 )
  115. DOUBLE PRECISION SELWI( 20 ), SELWR( 20 )
  116. * ..
  117. * .. Scalars in Common ..
  118. LOGICAL LERR, OK
  119. CHARACTER*32 SRNAMT
  120. INTEGER INFOT, NOUT, SELDIM, SELOPT
  121. * ..
  122. * .. Common blocks ..
  123. COMMON / INFOC / INFOT, NOUT, OK, LERR
  124. COMMON / SRNAMC / SRNAMT
  125. COMMON / SSLCT / SELOPT, SELDIM, SELVAL, SELWR, SELWI
  126. * ..
  127. * .. Executable Statements ..
  128. *
  129. NOUT = NUNIT
  130. WRITE( NOUT, FMT = * )
  131. C2 = PATH( 2: 3 )
  132. *
  133. * Initialize A
  134. *
  135. DO 20 J = 1, NMAX
  136. DO 10 I = 1, NMAX
  137. A( I, J ) = ZERO
  138. 10 CONTINUE
  139. 20 CONTINUE
  140. DO 30 I = 1, NMAX
  141. A( I, I ) = ONE
  142. 30 CONTINUE
  143. OK = .TRUE.
  144. NT = 0
  145. *
  146. IF( LSAMEN( 2, C2, 'EV' ) ) THEN
  147. *
  148. * Test ZGEEV
  149. *
  150. SRNAMT = 'ZGEEV '
  151. INFOT = 1
  152. CALL ZGEEV( 'X', 'N', 0, A, 1, X, VL, 1, VR, 1, W, 1, RW,
  153. $ INFO )
  154. CALL CHKXER( 'ZGEEV ', INFOT, NOUT, LERR, OK )
  155. INFOT = 2
  156. CALL ZGEEV( 'N', 'X', 0, A, 1, X, VL, 1, VR, 1, W, 1, RW,
  157. $ INFO )
  158. CALL CHKXER( 'ZGEEV ', INFOT, NOUT, LERR, OK )
  159. INFOT = 3
  160. CALL ZGEEV( 'N', 'N', -1, A, 1, X, VL, 1, VR, 1, W, 1, RW,
  161. $ INFO )
  162. CALL CHKXER( 'ZGEEV ', INFOT, NOUT, LERR, OK )
  163. INFOT = 5
  164. CALL ZGEEV( 'N', 'N', 2, A, 1, X, VL, 1, VR, 1, W, 4, RW,
  165. $ INFO )
  166. CALL CHKXER( 'ZGEEV ', INFOT, NOUT, LERR, OK )
  167. INFOT = 8
  168. CALL ZGEEV( 'V', 'N', 2, A, 2, X, VL, 1, VR, 1, W, 4, RW,
  169. $ INFO )
  170. CALL CHKXER( 'ZGEEV ', INFOT, NOUT, LERR, OK )
  171. INFOT = 10
  172. CALL ZGEEV( 'N', 'V', 2, A, 2, X, VL, 1, VR, 1, W, 4, RW,
  173. $ INFO )
  174. CALL CHKXER( 'ZGEEV ', INFOT, NOUT, LERR, OK )
  175. INFOT = 12
  176. CALL ZGEEV( 'V', 'V', 1, A, 1, X, VL, 1, VR, 1, W, 1, RW,
  177. $ INFO )
  178. CALL CHKXER( 'ZGEEV ', INFOT, NOUT, LERR, OK )
  179. NT = NT + 7
  180. *
  181. ELSE IF( LSAMEN( 2, C2, 'ES' ) ) THEN
  182. *
  183. * Test ZGEES
  184. *
  185. SRNAMT = 'ZGEES '
  186. INFOT = 1
  187. CALL ZGEES( 'X', 'N', ZSLECT, 0, A, 1, SDIM, X, VL, 1, W, 1,
  188. $ RW, B, INFO )
  189. CALL CHKXER( 'ZGEES ', INFOT, NOUT, LERR, OK )
  190. INFOT = 2
  191. CALL ZGEES( 'N', 'X', ZSLECT, 0, A, 1, SDIM, X, VL, 1, W, 1,
  192. $ RW, B, INFO )
  193. CALL CHKXER( 'ZGEES ', INFOT, NOUT, LERR, OK )
  194. INFOT = 4
  195. CALL ZGEES( 'N', 'S', ZSLECT, -1, A, 1, SDIM, X, VL, 1, W, 1,
  196. $ RW, B, INFO )
  197. CALL CHKXER( 'ZGEES ', INFOT, NOUT, LERR, OK )
  198. INFOT = 6
  199. CALL ZGEES( 'N', 'S', ZSLECT, 2, A, 1, SDIM, X, VL, 1, W, 4,
  200. $ RW, B, INFO )
  201. CALL CHKXER( 'ZGEES ', INFOT, NOUT, LERR, OK )
  202. INFOT = 10
  203. CALL ZGEES( 'V', 'S', ZSLECT, 2, A, 2, SDIM, X, VL, 1, W, 4,
  204. $ RW, B, INFO )
  205. CALL CHKXER( 'ZGEES ', INFOT, NOUT, LERR, OK )
  206. INFOT = 12
  207. CALL ZGEES( 'N', 'S', ZSLECT, 1, A, 1, SDIM, X, VL, 1, W, 1,
  208. $ RW, B, INFO )
  209. CALL CHKXER( 'ZGEES ', INFOT, NOUT, LERR, OK )
  210. NT = NT + 6
  211. *
  212. ELSE IF( LSAMEN( 2, C2, 'VX' ) ) THEN
  213. *
  214. * Test ZGEEVX
  215. *
  216. SRNAMT = 'ZGEEVX'
  217. INFOT = 1
  218. CALL ZGEEVX( 'X', 'N', 'N', 'N', 0, A, 1, X, VL, 1, VR, 1, ILO,
  219. $ IHI, S, ABNRM, R1, R2, W, 1, RW, INFO )
  220. CALL CHKXER( 'ZGEEVX', INFOT, NOUT, LERR, OK )
  221. INFOT = 2
  222. CALL ZGEEVX( 'N', 'X', 'N', 'N', 0, A, 1, X, VL, 1, VR, 1, ILO,
  223. $ IHI, S, ABNRM, R1, R2, W, 1, RW, INFO )
  224. CALL CHKXER( 'ZGEEVX', INFOT, NOUT, LERR, OK )
  225. INFOT = 3
  226. CALL ZGEEVX( 'N', 'N', 'X', 'N', 0, A, 1, X, VL, 1, VR, 1, ILO,
  227. $ IHI, S, ABNRM, R1, R2, W, 1, RW, INFO )
  228. CALL CHKXER( 'ZGEEVX', INFOT, NOUT, LERR, OK )
  229. INFOT = 4
  230. CALL ZGEEVX( 'N', 'N', 'N', 'X', 0, A, 1, X, VL, 1, VR, 1, ILO,
  231. $ IHI, S, ABNRM, R1, R2, W, 1, RW, INFO )
  232. CALL CHKXER( 'ZGEEVX', INFOT, NOUT, LERR, OK )
  233. INFOT = 5
  234. CALL ZGEEVX( 'N', 'N', 'N', 'N', -1, A, 1, X, VL, 1, VR, 1,
  235. $ ILO, IHI, S, ABNRM, R1, R2, W, 1, RW, INFO )
  236. CALL CHKXER( 'ZGEEVX', INFOT, NOUT, LERR, OK )
  237. INFOT = 7
  238. CALL ZGEEVX( 'N', 'N', 'N', 'N', 2, A, 1, X, VL, 1, VR, 1, ILO,
  239. $ IHI, S, ABNRM, R1, R2, W, 4, RW, INFO )
  240. CALL CHKXER( 'ZGEEVX', INFOT, NOUT, LERR, OK )
  241. INFOT = 10
  242. CALL ZGEEVX( 'N', 'V', 'N', 'N', 2, A, 2, X, VL, 1, VR, 1, ILO,
  243. $ IHI, S, ABNRM, R1, R2, W, 4, RW, INFO )
  244. CALL CHKXER( 'ZGEEVX', INFOT, NOUT, LERR, OK )
  245. INFOT = 12
  246. CALL ZGEEVX( 'N', 'N', 'V', 'N', 2, A, 2, X, VL, 1, VR, 1, ILO,
  247. $ IHI, S, ABNRM, R1, R2, W, 4, RW, INFO )
  248. CALL CHKXER( 'ZGEEVX', INFOT, NOUT, LERR, OK )
  249. INFOT = 20
  250. CALL ZGEEVX( 'N', 'N', 'N', 'N', 1, A, 1, X, VL, 1, VR, 1, ILO,
  251. $ IHI, S, ABNRM, R1, R2, W, 1, RW, INFO )
  252. CALL CHKXER( 'ZGEEVX', INFOT, NOUT, LERR, OK )
  253. INFOT = 20
  254. CALL ZGEEVX( 'N', 'N', 'V', 'V', 1, A, 1, X, VL, 1, VR, 1, ILO,
  255. $ IHI, S, ABNRM, R1, R2, W, 2, RW, INFO )
  256. CALL CHKXER( 'ZGEEVX', INFOT, NOUT, LERR, OK )
  257. NT = NT + 10
  258. *
  259. ELSE IF( LSAMEN( 2, C2, 'SX' ) ) THEN
  260. *
  261. * Test ZGEESX
  262. *
  263. SRNAMT = 'ZGEESX'
  264. INFOT = 1
  265. CALL ZGEESX( 'X', 'N', ZSLECT, 'N', 0, A, 1, SDIM, X, VL, 1,
  266. $ R1( 1 ), R2( 1 ), W, 1, RW, B, INFO )
  267. CALL CHKXER( 'ZGEESX', INFOT, NOUT, LERR, OK )
  268. INFOT = 2
  269. CALL ZGEESX( 'N', 'X', ZSLECT, 'N', 0, A, 1, SDIM, X, VL, 1,
  270. $ R1( 1 ), R2( 1 ), W, 1, RW, B, INFO )
  271. CALL CHKXER( 'ZGEESX', INFOT, NOUT, LERR, OK )
  272. INFOT = 4
  273. CALL ZGEESX( 'N', 'N', ZSLECT, 'X', 0, A, 1, SDIM, X, VL, 1,
  274. $ R1( 1 ), R2( 1 ), W, 1, RW, B, INFO )
  275. CALL CHKXER( 'ZGEESX', INFOT, NOUT, LERR, OK )
  276. INFOT = 5
  277. CALL ZGEESX( 'N', 'N', ZSLECT, 'N', -1, A, 1, SDIM, X, VL, 1,
  278. $ R1( 1 ), R2( 1 ), W, 1, RW, B, INFO )
  279. CALL CHKXER( 'ZGEESX', INFOT, NOUT, LERR, OK )
  280. INFOT = 7
  281. CALL ZGEESX( 'N', 'N', ZSLECT, 'N', 2, A, 1, SDIM, X, VL, 1,
  282. $ R1( 1 ), R2( 1 ), W, 4, RW, B, INFO )
  283. CALL CHKXER( 'ZGEESX', INFOT, NOUT, LERR, OK )
  284. INFOT = 11
  285. CALL ZGEESX( 'V', 'N', ZSLECT, 'N', 2, A, 2, SDIM, X, VL, 1,
  286. $ R1( 1 ), R2( 1 ), W, 4, RW, B, INFO )
  287. CALL CHKXER( 'ZGEESX', INFOT, NOUT, LERR, OK )
  288. INFOT = 15
  289. CALL ZGEESX( 'N', 'N', ZSLECT, 'N', 1, A, 1, SDIM, X, VL, 1,
  290. $ R1( 1 ), R2( 1 ), W, 1, RW, B, INFO )
  291. CALL CHKXER( 'ZGEESX', INFOT, NOUT, LERR, OK )
  292. NT = NT + 7
  293. *
  294. ELSE IF( LSAMEN( 2, C2, 'BD' ) ) THEN
  295. *
  296. * Test ZGESVD
  297. *
  298. SRNAMT = 'ZGESVD'
  299. INFOT = 1
  300. CALL ZGESVD( 'X', 'N', 0, 0, A, 1, S, U, 1, VT, 1, W, 1, RW,
  301. $ INFO )
  302. CALL CHKXER( 'ZGESVD', INFOT, NOUT, LERR, OK )
  303. INFOT = 2
  304. CALL ZGESVD( 'N', 'X', 0, 0, A, 1, S, U, 1, VT, 1, W, 1, RW,
  305. $ INFO )
  306. CALL CHKXER( 'ZGESVD', INFOT, NOUT, LERR, OK )
  307. INFOT = 2
  308. CALL ZGESVD( 'O', 'O', 0, 0, A, 1, S, U, 1, VT, 1, W, 1, RW,
  309. $ INFO )
  310. CALL CHKXER( 'ZGESVD', INFOT, NOUT, LERR, OK )
  311. INFOT = 3
  312. CALL ZGESVD( 'N', 'N', -1, 0, A, 1, S, U, 1, VT, 1, W, 1, RW,
  313. $ INFO )
  314. CALL CHKXER( 'ZGESVD', INFOT, NOUT, LERR, OK )
  315. INFOT = 4
  316. CALL ZGESVD( 'N', 'N', 0, -1, A, 1, S, U, 1, VT, 1, W, 1, RW,
  317. $ INFO )
  318. CALL CHKXER( 'ZGESVD', INFOT, NOUT, LERR, OK )
  319. INFOT = 6
  320. CALL ZGESVD( 'N', 'N', 2, 1, A, 1, S, U, 1, VT, 1, W, 5, RW,
  321. $ INFO )
  322. CALL CHKXER( 'ZGESVD', INFOT, NOUT, LERR, OK )
  323. INFOT = 9
  324. CALL ZGESVD( 'A', 'N', 2, 1, A, 2, S, U, 1, VT, 1, W, 5, RW,
  325. $ INFO )
  326. CALL CHKXER( 'ZGESVD', INFOT, NOUT, LERR, OK )
  327. INFOT = 11
  328. CALL ZGESVD( 'N', 'A', 1, 2, A, 1, S, U, 1, VT, 1, W, 5, RW,
  329. $ INFO )
  330. CALL CHKXER( 'ZGESVD', INFOT, NOUT, LERR, OK )
  331. NT = NT + 8
  332. IF( OK ) THEN
  333. WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
  334. $ NT
  335. ELSE
  336. WRITE( NOUT, FMT = 9998 )
  337. END IF
  338. *
  339. * Test ZGESDD
  340. *
  341. SRNAMT = 'ZGESDD'
  342. INFOT = 1
  343. CALL ZGESDD( 'X', 0, 0, A, 1, S, U, 1, VT, 1, W, 1, RW, IW,
  344. $ INFO )
  345. CALL CHKXER( 'ZGESDD', INFOT, NOUT, LERR, OK )
  346. INFOT = 2
  347. CALL ZGESDD( 'N', -1, 0, A, 1, S, U, 1, VT, 1, W, 1, RW, IW,
  348. $ INFO )
  349. CALL CHKXER( 'ZGESDD', INFOT, NOUT, LERR, OK )
  350. INFOT = 3
  351. CALL ZGESDD( 'N', 0, -1, A, 1, S, U, 1, VT, 1, W, 1, RW, IW,
  352. $ INFO )
  353. CALL CHKXER( 'ZGESDD', INFOT, NOUT, LERR, OK )
  354. INFOT = 5
  355. CALL ZGESDD( 'N', 2, 1, A, 1, S, U, 1, VT, 1, W, 5, RW, IW,
  356. $ INFO )
  357. CALL CHKXER( 'ZGESDD', INFOT, NOUT, LERR, OK )
  358. INFOT = 8
  359. CALL ZGESDD( 'A', 2, 1, A, 2, S, U, 1, VT, 1, W, 5, RW, IW,
  360. $ INFO )
  361. CALL CHKXER( 'ZGESDD', INFOT, NOUT, LERR, OK )
  362. INFOT = 10
  363. CALL ZGESDD( 'A', 1, 2, A, 1, S, U, 1, VT, 1, W, 5, RW, IW,
  364. $ INFO )
  365. CALL CHKXER( 'ZGESDD', INFOT, NOUT, LERR, OK )
  366. NT = NT - 2
  367. IF( OK ) THEN
  368. WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
  369. $ NT
  370. ELSE
  371. WRITE( NOUT, FMT = 9998 )
  372. END IF
  373. *
  374. * Test ZGEJSV
  375. *
  376. SRNAMT = 'ZGEJSV'
  377. INFOT = 1
  378. CALL ZGEJSV( 'X', 'U', 'V', 'R', 'N', 'N',
  379. $ 0, 0, A, 1, S, U, 1, VT, 1,
  380. $ W, 1, RW, 1, IW, INFO)
  381. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  382. INFOT = 2
  383. CALL ZGEJSV( 'G', 'X', 'V', 'R', 'N', 'N',
  384. $ 0, 0, A, 1, S, U, 1, VT, 1,
  385. $ W, 1, RW, 1, IW, INFO)
  386. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  387. INFOT = 3
  388. CALL ZGEJSV( 'G', 'U', 'X', 'R', 'N', 'N',
  389. $ 0, 0, A, 1, S, U, 1, VT, 1,
  390. $ W, 1, RW, 1, IW, INFO)
  391. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  392. INFOT = 4
  393. CALL ZGEJSV( 'G', 'U', 'V', 'X', 'N', 'N',
  394. $ 0, 0, A, 1, S, U, 1, VT, 1,
  395. $ W, 1, RW, 1, IW, INFO)
  396. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  397. INFOT = 5
  398. CALL ZGEJSV( 'G', 'U', 'V', 'R', 'X', 'N',
  399. $ 0, 0, A, 1, S, U, 1, VT, 1,
  400. $ W, 1, RW, 1, IW, INFO)
  401. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  402. INFOT = 6
  403. CALL ZGEJSV( 'G', 'U', 'V', 'R', 'N', 'X',
  404. $ 0, 0, A, 1, S, U, 1, VT, 1,
  405. $ W, 1, RW, 1, IW, INFO)
  406. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  407. INFOT = 7
  408. CALL ZGEJSV( 'G', 'U', 'V', 'R', 'N', 'N',
  409. $ -1, 0, A, 1, S, U, 1, VT, 1,
  410. $ W, 1, RW, 1, IW, INFO)
  411. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  412. INFOT = 8
  413. CALL ZGEJSV( 'G', 'U', 'V', 'R', 'N', 'N',
  414. $ 0, -1, A, 1, S, U, 1, VT, 1,
  415. $ W, 1, RW, 1, IW, INFO)
  416. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  417. INFOT = 10
  418. CALL ZGEJSV( 'G', 'U', 'V', 'R', 'N', 'N',
  419. $ 2, 1, A, 1, S, U, 1, VT, 1,
  420. $ W, 1, RW, 1, IW, INFO)
  421. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  422. INFOT = 13
  423. CALL ZGEJSV( 'G', 'U', 'V', 'R', 'N', 'N',
  424. $ 2, 2, A, 2, S, U, 1, VT, 2,
  425. $ W, 1, RW, 1, IW, INFO)
  426. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  427. INFOT = 15
  428. CALL ZGEJSV( 'G', 'U', 'V', 'R', 'N', 'N',
  429. $ 2, 2, A, 2, S, U, 2, VT, 1,
  430. $ W, 1, RW, 1, IW, INFO)
  431. CALL CHKXER( 'ZGEJSV', INFOT, NOUT, LERR, OK )
  432. NT = 11
  433. IF( OK ) THEN
  434. WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
  435. $ NT
  436. ELSE
  437. WRITE( NOUT, FMT = 9998 )
  438. END IF
  439. *
  440. * Test ZGESVDX
  441. *
  442. SRNAMT = 'ZGESVDX'
  443. INFOT = 1
  444. CALL ZGESVDX( 'X', 'N', 'A', 0, 0, A, 1, ZERO, ZERO,
  445. $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  446. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  447. INFOT = 2
  448. CALL ZGESVDX( 'N', 'X', 'A', 0, 0, A, 1, ZERO, ZERO,
  449. $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  450. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  451. INFOT = 3
  452. CALL ZGESVDX( 'N', 'N', 'X', 0, 0, A, 1, ZERO, ZERO,
  453. $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  454. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  455. INFOT = 4
  456. CALL ZGESVDX( 'N', 'N', 'A', -1, 0, A, 1, ZERO, ZERO,
  457. $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  458. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  459. INFOT = 5
  460. CALL ZGESVDX( 'N', 'N', 'A', 0, -1, A, 1, ZERO, ZERO,
  461. $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  462. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  463. INFOT = 7
  464. CALL ZGESVDX( 'N', 'N', 'A', 2, 1, A, 1, ZERO, ZERO,
  465. $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  466. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  467. INFOT = 8
  468. CALL ZGESVDX( 'N', 'N', 'V', 2, 1, A, 2, -ONE, ZERO,
  469. $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  470. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  471. INFOT = 9
  472. CALL ZGESVDX( 'N', 'N', 'V', 2, 1, A, 2, ONE, ZERO,
  473. $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  474. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  475. INFOT = 10
  476. CALL ZGESVDX( 'N', 'N', 'I', 2, 2, A, 2, ZERO, ZERO,
  477. $ 0, 1, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  478. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  479. INFOT = 11
  480. CALL ZGESVDX( 'V', 'N', 'I', 2, 2, A, 2, ZERO, ZERO,
  481. $ 1, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  482. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  483. INFOT = 15
  484. CALL ZGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
  485. $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  486. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  487. INFOT = 17
  488. CALL ZGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
  489. $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
  490. CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
  491. NT = 12
  492. IF( OK ) THEN
  493. WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
  494. $ NT
  495. ELSE
  496. WRITE( NOUT, FMT = 9998 )
  497. END IF
  498. END IF
  499. *
  500. * Print a summary line.
  501. *
  502. IF( .NOT.LSAMEN( 2, C2, 'BD' ) ) THEN
  503. IF( OK ) THEN
  504. WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
  505. $ NT
  506. ELSE
  507. WRITE( NOUT, FMT = 9998 )
  508. END IF
  509. END IF
  510. *
  511. 9999 FORMAT( 1X, A, ' passed the tests of the error exits (', I3,
  512. $ ' tests done)' )
  513. 9998 FORMAT( ' *** ', A, ' failed the tests of the error exits ***' )
  514. RETURN
  515. *
  516. * End of ZERRED
  517. *
  518. END