|
LAPACK 3.3.0
|
|
LAPACK 3.3.0
|
00001 SUBROUTINE CERRPO( PATH, NUNIT ) 00002 * 00003 * -- LAPACK test routine (version 3.2.1) -- 00004 * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. 00005 * April 2009 00006 * 00007 * .. Scalar Arguments .. 00008 CHARACTER*3 PATH 00009 INTEGER NUNIT 00010 * .. 00011 * 00012 * Purpose 00013 * ======= 00014 * 00015 * CERRPO tests the error exits for the COMPLEX routines 00016 * for Hermitian positive definite matrices. 00017 * 00018 * Note that this file is used only when the XBLAS are available, 00019 * otherwise cerrpo.f defines this subroutine. 00020 * 00021 * Arguments 00022 * ========= 00023 * 00024 * PATH (input) CHARACTER*3 00025 * The LAPACK path name for the routines to be tested. 00026 * 00027 * NUNIT (input) INTEGER 00028 * The unit number for output. 00029 * 00030 * ===================================================================== 00031 * 00032 * .. Parameters .. 00033 INTEGER NMAX 00034 PARAMETER ( NMAX = 4 ) 00035 * .. 00036 * .. Local Scalars .. 00037 CHARACTER EQ 00038 CHARACTER*2 C2 00039 INTEGER I, INFO, J, N_ERR_BNDS, NPARAMS 00040 REAL ANRM, RCOND, BERR 00041 * .. 00042 * .. Local Arrays .. 00043 REAL S( NMAX ), R( NMAX ), R1( NMAX ), R2( NMAX ), 00044 $ ERR_BNDS_N( NMAX, 3 ), ERR_BNDS_C( NMAX, 3 ), 00045 $ PARAMS( 1 ) 00046 COMPLEX A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ), 00047 $ W( 2*NMAX ), X( NMAX ) 00048 * .. 00049 * .. External Functions .. 00050 LOGICAL LSAMEN 00051 EXTERNAL LSAMEN 00052 * .. 00053 * .. External Subroutines .. 00054 EXTERNAL ALAESM, CHKXER, CPBCON, CPBEQU, CPBRFS, CPBTF2, 00055 $ CPBTRF, CPBTRS, CPOCON, CPOEQU, CPORFS, CPOTF2, 00056 $ CPOTRF, CPOTRI, CPOTRS, CPPCON, CPPEQU, CPPRFS, 00057 $ CPPTRF, CPPTRI, CPPTRS, CPOEQUB, CPORFSX 00058 * .. 00059 * .. Scalars in Common .. 00060 LOGICAL LERR, OK 00061 CHARACTER*32 SRNAMT 00062 INTEGER INFOT, NOUT 00063 * .. 00064 * .. Common blocks .. 00065 COMMON / INFOC / INFOT, NOUT, OK, LERR 00066 COMMON / SRNAMC / SRNAMT 00067 * .. 00068 * .. Intrinsic Functions .. 00069 INTRINSIC CMPLX, REAL 00070 * .. 00071 * .. Executable Statements .. 00072 * 00073 NOUT = NUNIT 00074 WRITE( NOUT, FMT = * ) 00075 C2 = PATH( 2: 3 ) 00076 * 00077 * Set the variables to innocuous values. 00078 * 00079 DO 20 J = 1, NMAX 00080 DO 10 I = 1, NMAX 00081 A( I, J ) = CMPLX( 1. / REAL( I+J ), -1. / REAL( I+J ) ) 00082 AF( I, J ) = CMPLX( 1. / REAL( I+J ), -1. / REAL( I+J ) ) 00083 10 CONTINUE 00084 B( J ) = 0. 00085 R1( J ) = 0. 00086 R2( J ) = 0. 00087 W( J ) = 0. 00088 X( J ) = 0. 00089 S( J ) = 0. 00090 20 CONTINUE 00091 ANRM = 1. 00092 OK = .TRUE. 00093 * 00094 * Test error exits of the routines that use the Cholesky 00095 * decomposition of a Hermitian positive definite matrix. 00096 * 00097 IF( LSAMEN( 2, C2, 'PO' ) ) THEN 00098 * 00099 * CPOTRF 00100 * 00101 SRNAMT = 'CPOTRF' 00102 INFOT = 1 00103 CALL CPOTRF( '/', 0, A, 1, INFO ) 00104 CALL CHKXER( 'CPOTRF', INFOT, NOUT, LERR, OK ) 00105 INFOT = 2 00106 CALL CPOTRF( 'U', -1, A, 1, INFO ) 00107 CALL CHKXER( 'CPOTRF', INFOT, NOUT, LERR, OK ) 00108 INFOT = 4 00109 CALL CPOTRF( 'U', 2, A, 1, INFO ) 00110 CALL CHKXER( 'CPOTRF', INFOT, NOUT, LERR, OK ) 00111 * 00112 * CPOTF2 00113 * 00114 SRNAMT = 'CPOTF2' 00115 INFOT = 1 00116 CALL CPOTF2( '/', 0, A, 1, INFO ) 00117 CALL CHKXER( 'CPOTF2', INFOT, NOUT, LERR, OK ) 00118 INFOT = 2 00119 CALL CPOTF2( 'U', -1, A, 1, INFO ) 00120 CALL CHKXER( 'CPOTF2', INFOT, NOUT, LERR, OK ) 00121 INFOT = 4 00122 CALL CPOTF2( 'U', 2, A, 1, INFO ) 00123 CALL CHKXER( 'CPOTF2', INFOT, NOUT, LERR, OK ) 00124 * 00125 * CPOTRI 00126 * 00127 SRNAMT = 'CPOTRI' 00128 INFOT = 1 00129 CALL CPOTRI( '/', 0, A, 1, INFO ) 00130 CALL CHKXER( 'CPOTRI', INFOT, NOUT, LERR, OK ) 00131 INFOT = 2 00132 CALL CPOTRI( 'U', -1, A, 1, INFO ) 00133 CALL CHKXER( 'CPOTRI', INFOT, NOUT, LERR, OK ) 00134 INFOT = 4 00135 CALL CPOTRI( 'U', 2, A, 1, INFO ) 00136 CALL CHKXER( 'CPOTRI', INFOT, NOUT, LERR, OK ) 00137 * 00138 * CPOTRS 00139 * 00140 SRNAMT = 'CPOTRS' 00141 INFOT = 1 00142 CALL CPOTRS( '/', 0, 0, A, 1, B, 1, INFO ) 00143 CALL CHKXER( 'CPOTRS', INFOT, NOUT, LERR, OK ) 00144 INFOT = 2 00145 CALL CPOTRS( 'U', -1, 0, A, 1, B, 1, INFO ) 00146 CALL CHKXER( 'CPOTRS', INFOT, NOUT, LERR, OK ) 00147 INFOT = 3 00148 CALL CPOTRS( 'U', 0, -1, A, 1, B, 1, INFO ) 00149 CALL CHKXER( 'CPOTRS', INFOT, NOUT, LERR, OK ) 00150 INFOT = 5 00151 CALL CPOTRS( 'U', 2, 1, A, 1, B, 2, INFO ) 00152 CALL CHKXER( 'CPOTRS', INFOT, NOUT, LERR, OK ) 00153 INFOT = 7 00154 CALL CPOTRS( 'U', 2, 1, A, 2, B, 1, INFO ) 00155 CALL CHKXER( 'CPOTRS', INFOT, NOUT, LERR, OK ) 00156 * 00157 * CPORFS 00158 * 00159 SRNAMT = 'CPORFS' 00160 INFOT = 1 00161 CALL CPORFS( '/', 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, R, 00162 $ INFO ) 00163 CALL CHKXER( 'CPORFS', INFOT, NOUT, LERR, OK ) 00164 INFOT = 2 00165 CALL CPORFS( 'U', -1, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, R, 00166 $ INFO ) 00167 CALL CHKXER( 'CPORFS', INFOT, NOUT, LERR, OK ) 00168 INFOT = 3 00169 CALL CPORFS( 'U', 0, -1, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, R, 00170 $ INFO ) 00171 CALL CHKXER( 'CPORFS', INFOT, NOUT, LERR, OK ) 00172 INFOT = 5 00173 CALL CPORFS( 'U', 2, 1, A, 1, AF, 2, B, 2, X, 2, R1, R2, W, R, 00174 $ INFO ) 00175 CALL CHKXER( 'CPORFS', INFOT, NOUT, LERR, OK ) 00176 INFOT = 7 00177 CALL CPORFS( 'U', 2, 1, A, 2, AF, 1, B, 2, X, 2, R1, R2, W, R, 00178 $ INFO ) 00179 CALL CHKXER( 'CPORFS', INFOT, NOUT, LERR, OK ) 00180 INFOT = 9 00181 CALL CPORFS( 'U', 2, 1, A, 2, AF, 2, B, 1, X, 2, R1, R2, W, R, 00182 $ INFO ) 00183 CALL CHKXER( 'CPORFS', INFOT, NOUT, LERR, OK ) 00184 INFOT = 11 00185 CALL CPORFS( 'U', 2, 1, A, 2, AF, 2, B, 2, X, 1, R1, R2, W, R, 00186 $ INFO ) 00187 CALL CHKXER( 'CPORFS', INFOT, NOUT, LERR, OK ) 00188 * 00189 * CPORFSX 00190 * 00191 N_ERR_BNDS = 3 00192 NPARAMS = 0 00193 SRNAMT = 'CPORFSX' 00194 INFOT = 1 00195 CALL CPORFSX( '/', EQ, 0, 0, A, 1, AF, 1, S, B, 1, X, 1, 00196 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS, 00197 $ PARAMS, W, R, INFO ) 00198 CALL CHKXER( 'CPORFSX', INFOT, NOUT, LERR, OK ) 00199 INFOT = 2 00200 CALL CPORFSX( 'U', EQ, -1, 0, A, 1, AF, 1, S, B, 1, X, 1, 00201 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS, 00202 $ PARAMS, W, R, INFO ) 00203 CALL CHKXER( 'CPORFSX', INFOT, NOUT, LERR, OK ) 00204 EQ = 'N' 00205 INFOT = 3 00206 CALL CPORFSX( 'U', EQ, -1, 0, A, 1, AF, 1, S, B, 1, X, 1, 00207 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS, 00208 $ PARAMS, W, R, INFO ) 00209 CALL CHKXER( 'CPORFSX', INFOT, NOUT, LERR, OK ) 00210 INFOT = 4 00211 CALL CPORFSX( 'U', EQ, 0, -1, A, 1, AF, 1, S, B, 1, X, 1, 00212 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS, 00213 $ PARAMS, W, R, INFO ) 00214 CALL CHKXER( 'CPORFSX', INFOT, NOUT, LERR, OK ) 00215 INFOT = 6 00216 CALL CPORFSX( 'U', EQ, 2, 1, A, 1, AF, 2, S, B, 2, X, 2, 00217 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS, 00218 $ PARAMS, W, R, INFO ) 00219 CALL CHKXER( 'CPORFSX', INFOT, NOUT, LERR, OK ) 00220 INFOT = 8 00221 CALL CPORFSX( 'U', EQ, 2, 1, A, 2, AF, 1, S, B, 2, X, 2, 00222 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS, 00223 $ PARAMS, W, R, INFO ) 00224 CALL CHKXER( 'CPORFSX', INFOT, NOUT, LERR, OK ) 00225 INFOT = 11 00226 CALL CPORFSX( 'U', EQ, 2, 1, A, 2, AF, 2, S, B, 1, X, 2, 00227 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS, 00228 $ PARAMS, W, R, INFO ) 00229 CALL CHKXER( 'CPORFSX', INFOT, NOUT, LERR, OK ) 00230 INFOT = 13 00231 CALL CPORFSX( 'U', EQ, 2, 1, A, 2, AF, 2, S, B, 2, X, 1, 00232 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS, 00233 $ PARAMS, W, R, INFO ) 00234 CALL CHKXER( 'CPORFSX', INFOT, NOUT, LERR, OK ) 00235 * 00236 * CPOCON 00237 * 00238 SRNAMT = 'CPOCON' 00239 INFOT = 1 00240 CALL CPOCON( '/', 0, A, 1, ANRM, RCOND, W, R, INFO ) 00241 CALL CHKXER( 'CPOCON', INFOT, NOUT, LERR, OK ) 00242 INFOT = 2 00243 CALL CPOCON( 'U', -1, A, 1, ANRM, RCOND, W, R, INFO ) 00244 CALL CHKXER( 'CPOCON', INFOT, NOUT, LERR, OK ) 00245 INFOT = 4 00246 CALL CPOCON( 'U', 2, A, 1, ANRM, RCOND, W, R, INFO ) 00247 CALL CHKXER( 'CPOCON', INFOT, NOUT, LERR, OK ) 00248 INFOT = 5 00249 CALL CPOCON( 'U', 1, A, 1, -ANRM, RCOND, W, R, INFO ) 00250 CALL CHKXER( 'CPOCON', INFOT, NOUT, LERR, OK ) 00251 * 00252 * CPOEQU 00253 * 00254 SRNAMT = 'CPOEQU' 00255 INFOT = 1 00256 CALL CPOEQU( -1, A, 1, R1, RCOND, ANRM, INFO ) 00257 CALL CHKXER( 'CPOEQU', INFOT, NOUT, LERR, OK ) 00258 INFOT = 3 00259 CALL CPOEQU( 2, A, 1, R1, RCOND, ANRM, INFO ) 00260 CALL CHKXER( 'CPOEQU', INFOT, NOUT, LERR, OK ) 00261 * 00262 * CPOEQUB 00263 * 00264 SRNAMT = 'CPOEQUB' 00265 INFOT = 1 00266 CALL CPOEQUB( -1, A, 1, R1, RCOND, ANRM, INFO ) 00267 CALL CHKXER( 'CPOEQUB', INFOT, NOUT, LERR, OK ) 00268 INFOT = 3 00269 CALL CPOEQUB( 2, A, 1, R1, RCOND, ANRM, INFO ) 00270 CALL CHKXER( 'CPOEQUB', INFOT, NOUT, LERR, OK ) 00271 * 00272 * Test error exits of the routines that use the Cholesky 00273 * decomposition of a Hermitian positive definite packed matrix. 00274 * 00275 ELSE IF( LSAMEN( 2, C2, 'PP' ) ) THEN 00276 * 00277 * CPPTRF 00278 * 00279 SRNAMT = 'CPPTRF' 00280 INFOT = 1 00281 CALL CPPTRF( '/', 0, A, INFO ) 00282 CALL CHKXER( 'CPPTRF', INFOT, NOUT, LERR, OK ) 00283 INFOT = 2 00284 CALL CPPTRF( 'U', -1, A, INFO ) 00285 CALL CHKXER( 'CPPTRF', INFOT, NOUT, LERR, OK ) 00286 * 00287 * CPPTRI 00288 * 00289 SRNAMT = 'CPPTRI' 00290 INFOT = 1 00291 CALL CPPTRI( '/', 0, A, INFO ) 00292 CALL CHKXER( 'CPPTRI', INFOT, NOUT, LERR, OK ) 00293 INFOT = 2 00294 CALL CPPTRI( 'U', -1, A, INFO ) 00295 CALL CHKXER( 'CPPTRI', INFOT, NOUT, LERR, OK ) 00296 * 00297 * CPPTRS 00298 * 00299 SRNAMT = 'CPPTRS' 00300 INFOT = 1 00301 CALL CPPTRS( '/', 0, 0, A, B, 1, INFO ) 00302 CALL CHKXER( 'CPPTRS', INFOT, NOUT, LERR, OK ) 00303 INFOT = 2 00304 CALL CPPTRS( 'U', -1, 0, A, B, 1, INFO ) 00305 CALL CHKXER( 'CPPTRS', INFOT, NOUT, LERR, OK ) 00306 INFOT = 3 00307 CALL CPPTRS( 'U', 0, -1, A, B, 1, INFO ) 00308 CALL CHKXER( 'CPPTRS', INFOT, NOUT, LERR, OK ) 00309 INFOT = 6 00310 CALL CPPTRS( 'U', 2, 1, A, B, 1, INFO ) 00311 CALL CHKXER( 'CPPTRS', INFOT, NOUT, LERR, OK ) 00312 * 00313 * CPPRFS 00314 * 00315 SRNAMT = 'CPPRFS' 00316 INFOT = 1 00317 CALL CPPRFS( '/', 0, 0, A, AF, B, 1, X, 1, R1, R2, W, R, INFO ) 00318 CALL CHKXER( 'CPPRFS', INFOT, NOUT, LERR, OK ) 00319 INFOT = 2 00320 CALL CPPRFS( 'U', -1, 0, A, AF, B, 1, X, 1, R1, R2, W, R, 00321 $ INFO ) 00322 CALL CHKXER( 'CPPRFS', INFOT, NOUT, LERR, OK ) 00323 INFOT = 3 00324 CALL CPPRFS( 'U', 0, -1, A, AF, B, 1, X, 1, R1, R2, W, R, 00325 $ INFO ) 00326 CALL CHKXER( 'CPPRFS', INFOT, NOUT, LERR, OK ) 00327 INFOT = 7 00328 CALL CPPRFS( 'U', 2, 1, A, AF, B, 1, X, 2, R1, R2, W, R, INFO ) 00329 CALL CHKXER( 'CPPRFS', INFOT, NOUT, LERR, OK ) 00330 INFOT = 9 00331 CALL CPPRFS( 'U', 2, 1, A, AF, B, 2, X, 1, R1, R2, W, R, INFO ) 00332 CALL CHKXER( 'CPPRFS', INFOT, NOUT, LERR, OK ) 00333 * 00334 * CPPCON 00335 * 00336 SRNAMT = 'CPPCON' 00337 INFOT = 1 00338 CALL CPPCON( '/', 0, A, ANRM, RCOND, W, R, INFO ) 00339 CALL CHKXER( 'CPPCON', INFOT, NOUT, LERR, OK ) 00340 INFOT = 2 00341 CALL CPPCON( 'U', -1, A, ANRM, RCOND, W, R, INFO ) 00342 CALL CHKXER( 'CPPCON', INFOT, NOUT, LERR, OK ) 00343 INFOT = 4 00344 CALL CPPCON( 'U', 1, A, -ANRM, RCOND, W, R, INFO ) 00345 CALL CHKXER( 'CPPCON', INFOT, NOUT, LERR, OK ) 00346 * 00347 * CPPEQU 00348 * 00349 SRNAMT = 'CPPEQU' 00350 INFOT = 1 00351 CALL CPPEQU( '/', 0, A, R1, RCOND, ANRM, INFO ) 00352 CALL CHKXER( 'CPPEQU', INFOT, NOUT, LERR, OK ) 00353 INFOT = 2 00354 CALL CPPEQU( 'U', -1, A, R1, RCOND, ANRM, INFO ) 00355 CALL CHKXER( 'CPPEQU', INFOT, NOUT, LERR, OK ) 00356 * 00357 * Test error exits of the routines that use the Cholesky 00358 * decomposition of a Hermitian positive definite band matrix. 00359 * 00360 ELSE IF( LSAMEN( 2, C2, 'PB' ) ) THEN 00361 * 00362 * CPBTRF 00363 * 00364 SRNAMT = 'CPBTRF' 00365 INFOT = 1 00366 CALL CPBTRF( '/', 0, 0, A, 1, INFO ) 00367 CALL CHKXER( 'CPBTRF', INFOT, NOUT, LERR, OK ) 00368 INFOT = 2 00369 CALL CPBTRF( 'U', -1, 0, A, 1, INFO ) 00370 CALL CHKXER( 'CPBTRF', INFOT, NOUT, LERR, OK ) 00371 INFOT = 3 00372 CALL CPBTRF( 'U', 1, -1, A, 1, INFO ) 00373 CALL CHKXER( 'CPBTRF', INFOT, NOUT, LERR, OK ) 00374 INFOT = 5 00375 CALL CPBTRF( 'U', 2, 1, A, 1, INFO ) 00376 CALL CHKXER( 'CPBTRF', INFOT, NOUT, LERR, OK ) 00377 * 00378 * CPBTF2 00379 * 00380 SRNAMT = 'CPBTF2' 00381 INFOT = 1 00382 CALL CPBTF2( '/', 0, 0, A, 1, INFO ) 00383 CALL CHKXER( 'CPBTF2', INFOT, NOUT, LERR, OK ) 00384 INFOT = 2 00385 CALL CPBTF2( 'U', -1, 0, A, 1, INFO ) 00386 CALL CHKXER( 'CPBTF2', INFOT, NOUT, LERR, OK ) 00387 INFOT = 3 00388 CALL CPBTF2( 'U', 1, -1, A, 1, INFO ) 00389 CALL CHKXER( 'CPBTF2', INFOT, NOUT, LERR, OK ) 00390 INFOT = 5 00391 CALL CPBTF2( 'U', 2, 1, A, 1, INFO ) 00392 CALL CHKXER( 'CPBTF2', INFOT, NOUT, LERR, OK ) 00393 * 00394 * CPBTRS 00395 * 00396 SRNAMT = 'CPBTRS' 00397 INFOT = 1 00398 CALL CPBTRS( '/', 0, 0, 0, A, 1, B, 1, INFO ) 00399 CALL CHKXER( 'CPBTRS', INFOT, NOUT, LERR, OK ) 00400 INFOT = 2 00401 CALL CPBTRS( 'U', -1, 0, 0, A, 1, B, 1, INFO ) 00402 CALL CHKXER( 'CPBTRS', INFOT, NOUT, LERR, OK ) 00403 INFOT = 3 00404 CALL CPBTRS( 'U', 1, -1, 0, A, 1, B, 1, INFO ) 00405 CALL CHKXER( 'CPBTRS', INFOT, NOUT, LERR, OK ) 00406 INFOT = 4 00407 CALL CPBTRS( 'U', 0, 0, -1, A, 1, B, 1, INFO ) 00408 CALL CHKXER( 'CPBTRS', INFOT, NOUT, LERR, OK ) 00409 INFOT = 6 00410 CALL CPBTRS( 'U', 2, 1, 1, A, 1, B, 1, INFO ) 00411 CALL CHKXER( 'CPBTRS', INFOT, NOUT, LERR, OK ) 00412 INFOT = 8 00413 CALL CPBTRS( 'U', 2, 0, 1, A, 1, B, 1, INFO ) 00414 CALL CHKXER( 'CPBTRS', INFOT, NOUT, LERR, OK ) 00415 * 00416 * CPBRFS 00417 * 00418 SRNAMT = 'CPBRFS' 00419 INFOT = 1 00420 CALL CPBRFS( '/', 0, 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, 00421 $ R, INFO ) 00422 CALL CHKXER( 'CPBRFS', INFOT, NOUT, LERR, OK ) 00423 INFOT = 2 00424 CALL CPBRFS( 'U', -1, 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, 00425 $ R, INFO ) 00426 CALL CHKXER( 'CPBRFS', INFOT, NOUT, LERR, OK ) 00427 INFOT = 3 00428 CALL CPBRFS( 'U', 1, -1, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, 00429 $ R, INFO ) 00430 CALL CHKXER( 'CPBRFS', INFOT, NOUT, LERR, OK ) 00431 INFOT = 4 00432 CALL CPBRFS( 'U', 0, 0, -1, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, 00433 $ R, INFO ) 00434 CALL CHKXER( 'CPBRFS', INFOT, NOUT, LERR, OK ) 00435 INFOT = 6 00436 CALL CPBRFS( 'U', 2, 1, 1, A, 1, AF, 2, B, 2, X, 2, R1, R2, W, 00437 $ R, INFO ) 00438 CALL CHKXER( 'CPBRFS', INFOT, NOUT, LERR, OK ) 00439 INFOT = 8 00440 CALL CPBRFS( 'U', 2, 1, 1, A, 2, AF, 1, B, 2, X, 2, R1, R2, W, 00441 $ R, INFO ) 00442 CALL CHKXER( 'CPBRFS', INFOT, NOUT, LERR, OK ) 00443 INFOT = 10 00444 CALL CPBRFS( 'U', 2, 0, 1, A, 1, AF, 1, B, 1, X, 2, R1, R2, W, 00445 $ R, INFO ) 00446 CALL CHKXER( 'CPBRFS', INFOT, NOUT, LERR, OK ) 00447 INFOT = 12 00448 CALL CPBRFS( 'U', 2, 0, 1, A, 1, AF, 1, B, 2, X, 1, R1, R2, W, 00449 $ R, INFO ) 00450 CALL CHKXER( 'CPBRFS', INFOT, NOUT, LERR, OK ) 00451 * 00452 * CPBCON 00453 * 00454 SRNAMT = 'CPBCON' 00455 INFOT = 1 00456 CALL CPBCON( '/', 0, 0, A, 1, ANRM, RCOND, W, R, INFO ) 00457 CALL CHKXER( 'CPBCON', INFOT, NOUT, LERR, OK ) 00458 INFOT = 2 00459 CALL CPBCON( 'U', -1, 0, A, 1, ANRM, RCOND, W, R, INFO ) 00460 CALL CHKXER( 'CPBCON', INFOT, NOUT, LERR, OK ) 00461 INFOT = 3 00462 CALL CPBCON( 'U', 1, -1, A, 1, ANRM, RCOND, W, R, INFO ) 00463 CALL CHKXER( 'CPBCON', INFOT, NOUT, LERR, OK ) 00464 INFOT = 5 00465 CALL CPBCON( 'U', 2, 1, A, 1, ANRM, RCOND, W, R, INFO ) 00466 CALL CHKXER( 'CPBCON', INFOT, NOUT, LERR, OK ) 00467 INFOT = 6 00468 CALL CPBCON( 'U', 1, 0, A, 1, -ANRM, RCOND, W, R, INFO ) 00469 CALL CHKXER( 'CPBCON', INFOT, NOUT, LERR, OK ) 00470 * 00471 * CPBEQU 00472 * 00473 SRNAMT = 'CPBEQU' 00474 INFOT = 1 00475 CALL CPBEQU( '/', 0, 0, A, 1, R1, RCOND, ANRM, INFO ) 00476 CALL CHKXER( 'CPBEQU', INFOT, NOUT, LERR, OK ) 00477 INFOT = 2 00478 CALL CPBEQU( 'U', -1, 0, A, 1, R1, RCOND, ANRM, INFO ) 00479 CALL CHKXER( 'CPBEQU', INFOT, NOUT, LERR, OK ) 00480 INFOT = 3 00481 CALL CPBEQU( 'U', 1, -1, A, 1, R1, RCOND, ANRM, INFO ) 00482 CALL CHKXER( 'CPBEQU', INFOT, NOUT, LERR, OK ) 00483 INFOT = 5 00484 CALL CPBEQU( 'U', 2, 1, A, 1, R1, RCOND, ANRM, INFO ) 00485 CALL CHKXER( 'CPBEQU', INFOT, NOUT, LERR, OK ) 00486 END IF 00487 * 00488 * Print a summary line. 00489 * 00490 CALL ALAESM( PATH, OK, NOUT ) 00491 * 00492 RETURN 00493 * 00494 * End of CERRPO 00495 * 00496 END
1.7.3 
