00001 SUBROUTINE DERRPO( PATH, NUNIT )
00002
00003
00004
00005
00006
00007
00008 CHARACTER*3 PATH
00009 INTEGER NUNIT
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033 INTEGER NMAX
00034 PARAMETER ( NMAX = 4 )
00035
00036
00037 CHARACTER EQ
00038 CHARACTER*2 C2
00039 INTEGER I, INFO, J, N_ERR_BNDS, NPARAMS
00040 DOUBLE PRECISION ANRM, RCOND, BERR
00041
00042
00043 INTEGER IW( NMAX )
00044 DOUBLE PRECISION A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
00045 $ R1( NMAX ), R2( NMAX ), W( 3*NMAX ), X( NMAX ),
00046 $ S( NMAX ), ERR_BNDS_N( NMAX, 3 ),
00047 $ ERR_BNDS_C( NMAX, 3), PARAMS( 1 )
00048
00049
00050 LOGICAL LSAMEN
00051 EXTERNAL LSAMEN
00052
00053
00054 EXTERNAL ALAESM, CHKXER, DPBCON, DPBEQU, DPBRFS, DPBTF2,
00055 $ DPBTRF, DPBTRS, DPOCON, DPOEQU, DPORFS, DPOTF2,
00056 $ DPOTRF, DPOTRI, DPOTRS, DPPCON, DPPEQU, DPPRFS,
00057 $ DPPTRF, DPPTRI, DPPTRS, DPOEQUB, DPORFSX
00058
00059
00060 LOGICAL LERR, OK
00061 CHARACTER*32 SRNAMT
00062 INTEGER INFOT, NOUT
00063
00064
00065 COMMON / INFOC / INFOT, NOUT, OK, LERR
00066 COMMON / SRNAMC / SRNAMT
00067
00068
00069 INTRINSIC DBLE
00070
00071
00072
00073 NOUT = NUNIT
00074 WRITE( NOUT, FMT = * )
00075 C2 = PATH( 2: 3 )
00076
00077
00078
00079 DO 20 J = 1, NMAX
00080 DO 10 I = 1, NMAX
00081 A( I, J ) = 1.D0 / DBLE( I+J )
00082 AF( I, J ) = 1.D0 / DBLE( I+J )
00083 10 CONTINUE
00084 B( J ) = 0.D0
00085 R1( J ) = 0.D0
00086 R2( J ) = 0.D0
00087 W( J ) = 0.D0
00088 X( J ) = 0.D0
00089 S( J ) = 0.D0
00090 IW( J ) = J
00091 20 CONTINUE
00092 OK = .TRUE.
00093
00094 IF( LSAMEN( 2, C2, 'PO' ) ) THEN
00095
00096
00097
00098
00099
00100
00101 SRNAMT = 'DPOTRF'
00102 INFOT = 1
00103 CALL DPOTRF( '/', 0, A, 1, INFO )
00104 CALL CHKXER( 'DPOTRF', INFOT, NOUT, LERR, OK )
00105 INFOT = 2
00106 CALL DPOTRF( 'U', -1, A, 1, INFO )
00107 CALL CHKXER( 'DPOTRF', INFOT, NOUT, LERR, OK )
00108 INFOT = 4
00109 CALL DPOTRF( 'U', 2, A, 1, INFO )
00110 CALL CHKXER( 'DPOTRF', INFOT, NOUT, LERR, OK )
00111
00112
00113
00114 SRNAMT = 'DPOTF2'
00115 INFOT = 1
00116 CALL DPOTF2( '/', 0, A, 1, INFO )
00117 CALL CHKXER( 'DPOTF2', INFOT, NOUT, LERR, OK )
00118 INFOT = 2
00119 CALL DPOTF2( 'U', -1, A, 1, INFO )
00120 CALL CHKXER( 'DPOTF2', INFOT, NOUT, LERR, OK )
00121 INFOT = 4
00122 CALL DPOTF2( 'U', 2, A, 1, INFO )
00123 CALL CHKXER( 'DPOTF2', INFOT, NOUT, LERR, OK )
00124
00125
00126
00127 SRNAMT = 'DPOTRI'
00128 INFOT = 1
00129 CALL DPOTRI( '/', 0, A, 1, INFO )
00130 CALL CHKXER( 'DPOTRI', INFOT, NOUT, LERR, OK )
00131 INFOT = 2
00132 CALL DPOTRI( 'U', -1, A, 1, INFO )
00133 CALL CHKXER( 'DPOTRI', INFOT, NOUT, LERR, OK )
00134 INFOT = 4
00135 CALL DPOTRI( 'U', 2, A, 1, INFO )
00136 CALL CHKXER( 'DPOTRI', INFOT, NOUT, LERR, OK )
00137
00138
00139
00140 SRNAMT = 'DPOTRS'
00141 INFOT = 1
00142 CALL DPOTRS( '/', 0, 0, A, 1, B, 1, INFO )
00143 CALL CHKXER( 'DPOTRS', INFOT, NOUT, LERR, OK )
00144 INFOT = 2
00145 CALL DPOTRS( 'U', -1, 0, A, 1, B, 1, INFO )
00146 CALL CHKXER( 'DPOTRS', INFOT, NOUT, LERR, OK )
00147 INFOT = 3
00148 CALL DPOTRS( 'U', 0, -1, A, 1, B, 1, INFO )
00149 CALL CHKXER( 'DPOTRS', INFOT, NOUT, LERR, OK )
00150 INFOT = 5
00151 CALL DPOTRS( 'U', 2, 1, A, 1, B, 2, INFO )
00152 CALL CHKXER( 'DPOTRS', INFOT, NOUT, LERR, OK )
00153 INFOT = 7
00154 CALL DPOTRS( 'U', 2, 1, A, 2, B, 1, INFO )
00155 CALL CHKXER( 'DPOTRS', INFOT, NOUT, LERR, OK )
00156
00157
00158
00159 SRNAMT = 'DPORFS'
00160 INFOT = 1
00161 CALL DPORFS( '/', 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, IW,
00162 $ INFO )
00163 CALL CHKXER( 'DPORFS', INFOT, NOUT, LERR, OK )
00164 INFOT = 2
00165 CALL DPORFS( 'U', -1, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
00166 $ IW, INFO )
00167 CALL CHKXER( 'DPORFS', INFOT, NOUT, LERR, OK )
00168 INFOT = 3
00169 CALL DPORFS( 'U', 0, -1, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
00170 $ IW, INFO )
00171 CALL CHKXER( 'DPORFS', INFOT, NOUT, LERR, OK )
00172 INFOT = 5
00173 CALL DPORFS( 'U', 2, 1, A, 1, AF, 2, B, 2, X, 2, R1, R2, W, IW,
00174 $ INFO )
00175 CALL CHKXER( 'DPORFS', INFOT, NOUT, LERR, OK )
00176 INFOT = 7
00177 CALL DPORFS( 'U', 2, 1, A, 2, AF, 1, B, 2, X, 2, R1, R2, W, IW,
00178 $ INFO )
00179 CALL CHKXER( 'DPORFS', INFOT, NOUT, LERR, OK )
00180 INFOT = 9
00181 CALL DPORFS( 'U', 2, 1, A, 2, AF, 2, B, 1, X, 2, R1, R2, W, IW,
00182 $ INFO )
00183 CALL CHKXER( 'DPORFS', INFOT, NOUT, LERR, OK )
00184 INFOT = 11
00185 CALL DPORFS( 'U', 2, 1, A, 2, AF, 2, B, 2, X, 1, R1, R2, W, IW,
00186 $ INFO )
00187 CALL CHKXER( 'DPORFS', INFOT, NOUT, LERR, OK )
00188
00189
00190
00191 N_ERR_BNDS = 3
00192 NPARAMS = 0
00193 SRNAMT = 'DPORFSX'
00194 INFOT = 1
00195 CALL DPORFSX( '/', 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, IW, INFO )
00198 CALL CHKXER( 'DPORFSX', INFOT, NOUT, LERR, OK )
00199 INFOT = 2
00200 CALL DPORFSX( '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, IW, INFO )
00203 CALL CHKXER( 'DPORFSX', INFOT, NOUT, LERR, OK )
00204 EQ = 'N'
00205 INFOT = 3
00206 CALL DPORFSX( '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, IW, INFO )
00209 CALL CHKXER( 'DPORFSX', INFOT, NOUT, LERR, OK )
00210 INFOT = 4
00211 CALL DPORFSX( '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, IW, INFO )
00214 CALL CHKXER( 'DPORFSX', INFOT, NOUT, LERR, OK )
00215 INFOT = 6
00216 CALL DPORFSX( '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, IW, INFO )
00219 CALL CHKXER( 'DPORFSX', INFOT, NOUT, LERR, OK )
00220 INFOT = 8
00221 CALL DPORFSX( '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, IW, INFO )
00224 CALL CHKXER( 'DPORFSX', INFOT, NOUT, LERR, OK )
00225 INFOT = 11
00226 CALL DPORFSX( '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, IW, INFO )
00229 CALL CHKXER( 'DPORFSX', INFOT, NOUT, LERR, OK )
00230 INFOT = 13
00231 CALL DPORFSX( '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, IW, INFO )
00234 CALL CHKXER( 'DPORFSX', INFOT, NOUT, LERR, OK )
00235
00236
00237
00238 SRNAMT = 'DPOCON'
00239 INFOT = 1
00240 CALL DPOCON( '/', 0, A, 1, ANRM, RCOND, W, IW, INFO )
00241 CALL CHKXER( 'DPOCON', INFOT, NOUT, LERR, OK )
00242 INFOT = 2
00243 CALL DPOCON( 'U', -1, A, 1, ANRM, RCOND, W, IW, INFO )
00244 CALL CHKXER( 'DPOCON', INFOT, NOUT, LERR, OK )
00245 INFOT = 4
00246 CALL DPOCON( 'U', 2, A, 1, ANRM, RCOND, W, IW, INFO )
00247 CALL CHKXER( 'DPOCON', INFOT, NOUT, LERR, OK )
00248
00249
00250
00251 SRNAMT = 'DPOEQU'
00252 INFOT = 1
00253 CALL DPOEQU( -1, A, 1, R1, RCOND, ANRM, INFO )
00254 CALL CHKXER( 'DPOEQU', INFOT, NOUT, LERR, OK )
00255 INFOT = 3
00256 CALL DPOEQU( 2, A, 1, R1, RCOND, ANRM, INFO )
00257 CALL CHKXER( 'DPOEQU', INFOT, NOUT, LERR, OK )
00258
00259
00260
00261 SRNAMT = 'DPOEQUB'
00262 INFOT = 1
00263 CALL DPOEQUB( -1, A, 1, R1, RCOND, ANRM, INFO )
00264 CALL CHKXER( 'DPOEQUB', INFOT, NOUT, LERR, OK )
00265 INFOT = 3
00266 CALL DPOEQUB( 2, A, 1, R1, RCOND, ANRM, INFO )
00267 CALL CHKXER( 'DPOEQUB', INFOT, NOUT, LERR, OK )
00268
00269 ELSE IF( LSAMEN( 2, C2, 'PP' ) ) THEN
00270
00271
00272
00273
00274
00275
00276 SRNAMT = 'DPPTRF'
00277 INFOT = 1
00278 CALL DPPTRF( '/', 0, A, INFO )
00279 CALL CHKXER( 'DPPTRF', INFOT, NOUT, LERR, OK )
00280 INFOT = 2
00281 CALL DPPTRF( 'U', -1, A, INFO )
00282 CALL CHKXER( 'DPPTRF', INFOT, NOUT, LERR, OK )
00283
00284
00285
00286 SRNAMT = 'DPPTRI'
00287 INFOT = 1
00288 CALL DPPTRI( '/', 0, A, INFO )
00289 CALL CHKXER( 'DPPTRI', INFOT, NOUT, LERR, OK )
00290 INFOT = 2
00291 CALL DPPTRI( 'U', -1, A, INFO )
00292 CALL CHKXER( 'DPPTRI', INFOT, NOUT, LERR, OK )
00293
00294
00295
00296 SRNAMT = 'DPPTRS'
00297 INFOT = 1
00298 CALL DPPTRS( '/', 0, 0, A, B, 1, INFO )
00299 CALL CHKXER( 'DPPTRS', INFOT, NOUT, LERR, OK )
00300 INFOT = 2
00301 CALL DPPTRS( 'U', -1, 0, A, B, 1, INFO )
00302 CALL CHKXER( 'DPPTRS', INFOT, NOUT, LERR, OK )
00303 INFOT = 3
00304 CALL DPPTRS( 'U', 0, -1, A, B, 1, INFO )
00305 CALL CHKXER( 'DPPTRS', INFOT, NOUT, LERR, OK )
00306 INFOT = 6
00307 CALL DPPTRS( 'U', 2, 1, A, B, 1, INFO )
00308 CALL CHKXER( 'DPPTRS', INFOT, NOUT, LERR, OK )
00309
00310
00311
00312 SRNAMT = 'DPPRFS'
00313 INFOT = 1
00314 CALL DPPRFS( '/', 0, 0, A, AF, B, 1, X, 1, R1, R2, W, IW,
00315 $ INFO )
00316 CALL CHKXER( 'DPPRFS', INFOT, NOUT, LERR, OK )
00317 INFOT = 2
00318 CALL DPPRFS( 'U', -1, 0, A, AF, B, 1, X, 1, R1, R2, W, IW,
00319 $ INFO )
00320 CALL CHKXER( 'DPPRFS', INFOT, NOUT, LERR, OK )
00321 INFOT = 3
00322 CALL DPPRFS( 'U', 0, -1, A, AF, B, 1, X, 1, R1, R2, W, IW,
00323 $ INFO )
00324 CALL CHKXER( 'DPPRFS', INFOT, NOUT, LERR, OK )
00325 INFOT = 7
00326 CALL DPPRFS( 'U', 2, 1, A, AF, B, 1, X, 2, R1, R2, W, IW,
00327 $ INFO )
00328 CALL CHKXER( 'DPPRFS', INFOT, NOUT, LERR, OK )
00329 INFOT = 9
00330 CALL DPPRFS( 'U', 2, 1, A, AF, B, 2, X, 1, R1, R2, W, IW,
00331 $ INFO )
00332 CALL CHKXER( 'DPPRFS', INFOT, NOUT, LERR, OK )
00333
00334
00335
00336 SRNAMT = 'DPPCON'
00337 INFOT = 1
00338 CALL DPPCON( '/', 0, A, ANRM, RCOND, W, IW, INFO )
00339 CALL CHKXER( 'DPPCON', INFOT, NOUT, LERR, OK )
00340 INFOT = 2
00341 CALL DPPCON( 'U', -1, A, ANRM, RCOND, W, IW, INFO )
00342 CALL CHKXER( 'DPPCON', INFOT, NOUT, LERR, OK )
00343
00344
00345
00346 SRNAMT = 'DPPEQU'
00347 INFOT = 1
00348 CALL DPPEQU( '/', 0, A, R1, RCOND, ANRM, INFO )
00349 CALL CHKXER( 'DPPEQU', INFOT, NOUT, LERR, OK )
00350 INFOT = 2
00351 CALL DPPEQU( 'U', -1, A, R1, RCOND, ANRM, INFO )
00352 CALL CHKXER( 'DPPEQU', INFOT, NOUT, LERR, OK )
00353
00354 ELSE IF( LSAMEN( 2, C2, 'PB' ) ) THEN
00355
00356
00357
00358
00359
00360
00361 SRNAMT = 'DPBTRF'
00362 INFOT = 1
00363 CALL DPBTRF( '/', 0, 0, A, 1, INFO )
00364 CALL CHKXER( 'DPBTRF', INFOT, NOUT, LERR, OK )
00365 INFOT = 2
00366 CALL DPBTRF( 'U', -1, 0, A, 1, INFO )
00367 CALL CHKXER( 'DPBTRF', INFOT, NOUT, LERR, OK )
00368 INFOT = 3
00369 CALL DPBTRF( 'U', 1, -1, A, 1, INFO )
00370 CALL CHKXER( 'DPBTRF', INFOT, NOUT, LERR, OK )
00371 INFOT = 5
00372 CALL DPBTRF( 'U', 2, 1, A, 1, INFO )
00373 CALL CHKXER( 'DPBTRF', INFOT, NOUT, LERR, OK )
00374
00375
00376
00377 SRNAMT = 'DPBTF2'
00378 INFOT = 1
00379 CALL DPBTF2( '/', 0, 0, A, 1, INFO )
00380 CALL CHKXER( 'DPBTF2', INFOT, NOUT, LERR, OK )
00381 INFOT = 2
00382 CALL DPBTF2( 'U', -1, 0, A, 1, INFO )
00383 CALL CHKXER( 'DPBTF2', INFOT, NOUT, LERR, OK )
00384 INFOT = 3
00385 CALL DPBTF2( 'U', 1, -1, A, 1, INFO )
00386 CALL CHKXER( 'DPBTF2', INFOT, NOUT, LERR, OK )
00387 INFOT = 5
00388 CALL DPBTF2( 'U', 2, 1, A, 1, INFO )
00389 CALL CHKXER( 'DPBTF2', INFOT, NOUT, LERR, OK )
00390
00391
00392
00393 SRNAMT = 'DPBTRS'
00394 INFOT = 1
00395 CALL DPBTRS( '/', 0, 0, 0, A, 1, B, 1, INFO )
00396 CALL CHKXER( 'DPBTRS', INFOT, NOUT, LERR, OK )
00397 INFOT = 2
00398 CALL DPBTRS( 'U', -1, 0, 0, A, 1, B, 1, INFO )
00399 CALL CHKXER( 'DPBTRS', INFOT, NOUT, LERR, OK )
00400 INFOT = 3
00401 CALL DPBTRS( 'U', 1, -1, 0, A, 1, B, 1, INFO )
00402 CALL CHKXER( 'DPBTRS', INFOT, NOUT, LERR, OK )
00403 INFOT = 4
00404 CALL DPBTRS( 'U', 0, 0, -1, A, 1, B, 1, INFO )
00405 CALL CHKXER( 'DPBTRS', INFOT, NOUT, LERR, OK )
00406 INFOT = 6
00407 CALL DPBTRS( 'U', 2, 1, 1, A, 1, B, 1, INFO )
00408 CALL CHKXER( 'DPBTRS', INFOT, NOUT, LERR, OK )
00409 INFOT = 8
00410 CALL DPBTRS( 'U', 2, 0, 1, A, 1, B, 1, INFO )
00411 CALL CHKXER( 'DPBTRS', INFOT, NOUT, LERR, OK )
00412
00413
00414
00415 SRNAMT = 'DPBRFS'
00416 INFOT = 1
00417 CALL DPBRFS( '/', 0, 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
00418 $ IW, INFO )
00419 CALL CHKXER( 'DPBRFS', INFOT, NOUT, LERR, OK )
00420 INFOT = 2
00421 CALL DPBRFS( 'U', -1, 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
00422 $ IW, INFO )
00423 CALL CHKXER( 'DPBRFS', INFOT, NOUT, LERR, OK )
00424 INFOT = 3
00425 CALL DPBRFS( 'U', 1, -1, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
00426 $ IW, INFO )
00427 CALL CHKXER( 'DPBRFS', INFOT, NOUT, LERR, OK )
00428 INFOT = 4
00429 CALL DPBRFS( 'U', 0, 0, -1, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
00430 $ IW, INFO )
00431 CALL CHKXER( 'DPBRFS', INFOT, NOUT, LERR, OK )
00432 INFOT = 6
00433 CALL DPBRFS( 'U', 2, 1, 1, A, 1, AF, 2, B, 2, X, 2, R1, R2, W,
00434 $ IW, INFO )
00435 CALL CHKXER( 'DPBRFS', INFOT, NOUT, LERR, OK )
00436 INFOT = 8
00437 CALL DPBRFS( 'U', 2, 1, 1, A, 2, AF, 1, B, 2, X, 2, R1, R2, W,
00438 $ IW, INFO )
00439 CALL CHKXER( 'DPBRFS', INFOT, NOUT, LERR, OK )
00440 INFOT = 10
00441 CALL DPBRFS( 'U', 2, 0, 1, A, 1, AF, 1, B, 1, X, 2, R1, R2, W,
00442 $ IW, INFO )
00443 CALL CHKXER( 'DPBRFS', INFOT, NOUT, LERR, OK )
00444 INFOT = 12
00445 CALL DPBRFS( 'U', 2, 0, 1, A, 1, AF, 1, B, 2, X, 1, R1, R2, W,
00446 $ IW, INFO )
00447 CALL CHKXER( 'DPBRFS', INFOT, NOUT, LERR, OK )
00448
00449
00450
00451 SRNAMT = 'DPBCON'
00452 INFOT = 1
00453 CALL DPBCON( '/', 0, 0, A, 1, ANRM, RCOND, W, IW, INFO )
00454 CALL CHKXER( 'DPBCON', INFOT, NOUT, LERR, OK )
00455 INFOT = 2
00456 CALL DPBCON( 'U', -1, 0, A, 1, ANRM, RCOND, W, IW, INFO )
00457 CALL CHKXER( 'DPBCON', INFOT, NOUT, LERR, OK )
00458 INFOT = 3
00459 CALL DPBCON( 'U', 1, -1, A, 1, ANRM, RCOND, W, IW, INFO )
00460 CALL CHKXER( 'DPBCON', INFOT, NOUT, LERR, OK )
00461 INFOT = 5
00462 CALL DPBCON( 'U', 2, 1, A, 1, ANRM, RCOND, W, IW, INFO )
00463 CALL CHKXER( 'DPBCON', INFOT, NOUT, LERR, OK )
00464
00465
00466
00467 SRNAMT = 'DPBEQU'
00468 INFOT = 1
00469 CALL DPBEQU( '/', 0, 0, A, 1, R1, RCOND, ANRM, INFO )
00470 CALL CHKXER( 'DPBEQU', INFOT, NOUT, LERR, OK )
00471 INFOT = 2
00472 CALL DPBEQU( 'U', -1, 0, A, 1, R1, RCOND, ANRM, INFO )
00473 CALL CHKXER( 'DPBEQU', INFOT, NOUT, LERR, OK )
00474 INFOT = 3
00475 CALL DPBEQU( 'U', 1, -1, A, 1, R1, RCOND, ANRM, INFO )
00476 CALL CHKXER( 'DPBEQU', INFOT, NOUT, LERR, OK )
00477 INFOT = 5
00478 CALL DPBEQU( 'U', 2, 1, A, 1, R1, RCOND, ANRM, INFO )
00479 CALL CHKXER( 'DPBEQU', INFOT, NOUT, LERR, OK )
00480 END IF
00481
00482
00483
00484 CALL ALAESM( PATH, OK, NOUT )
00485
00486 RETURN
00487
00488
00489
00490 END