cerrpo.c
Go to the documentation of this file.
00001 /* cerrpo.f -- translated by f2c (version 20061008).
00002    You must link the resulting object file with libf2c:
00003         on Microsoft Windows system, link with libf2c.lib;
00004         on Linux or Unix systems, link with .../path/to/libf2c.a -lm
00005         or, if you install libf2c.a in a standard place, with -lf2c -lm
00006         -- in that order, at the end of the command line, as in
00007                 cc *.o -lf2c -lm
00008         Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
00009 
00010                 http://www.netlib.org/f2c/libf2c.zip
00011 */
00012 
00013 #include "f2c.h"
00014 #include "blaswrap.h"
00015 
00016 /* Common Block Declarations */
00017 
00018 struct {
00019     integer infot, nout;
00020     logical ok, lerr;
00021 } infoc_;
00022 
00023 #define infoc_1 infoc_
00024 
00025 struct {
00026     char srnamt[32];
00027 } srnamc_;
00028 
00029 #define srnamc_1 srnamc_
00030 
00031 /* Table of constant values */
00032 
00033 static integer c__2 = 2;
00034 static integer c__0 = 0;
00035 static integer c__1 = 1;
00036 static integer c_n1 = -1;
00037 
00038 /* Subroutine */ int cerrpo_(char *path, integer *nunit)
00039 {
00040     /* System generated locals */
00041     integer i__1;
00042     real r__1, r__2;
00043     complex q__1;
00044 
00045     /* Builtin functions */
00046     integer s_wsle(cilist *), e_wsle(void);
00047     /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
00048 
00049     /* Local variables */
00050     complex a[16]       /* was [4][4] */, b[4];
00051     integer i__, j;
00052     real r__[4];
00053     complex w[8], x[4];
00054     char c2[2];
00055     real r1[4], r2[4];
00056     complex af[16]      /* was [4][4] */;
00057     integer info;
00058     real anrm, rcond;
00059     extern /* Subroutine */ int cpbtf2_(char *, integer *, integer *, complex 
00060             *, integer *, integer *), cpotf2_(char *, integer *, 
00061             complex *, integer *, integer *), alaesm_(char *, logical 
00062             *, integer *), cpbcon_(char *, integer *, integer *, 
00063             complex *, integer *, real *, real *, complex *, real *, integer *
00064 );
00065     extern logical lsamen_(integer *, char *, char *);
00066     extern /* Subroutine */ int cpbequ_(char *, integer *, integer *, complex 
00067             *, integer *, real *, real *, real *, integer *), cpbrfs_(
00068             char *, integer *, integer *, integer *, complex *, integer *, 
00069             complex *, integer *, complex *, integer *, complex *, integer *, 
00070             real *, real *, complex *, real *, integer *), cpbtrf_(
00071             char *, integer *, integer *, complex *, integer *, integer *), cpocon_(char *, integer *, complex *, integer *, real *, 
00072             real *, complex *, real *, integer *), chkxer_(char *, 
00073             integer *, integer *, logical *, logical *), cppcon_(char 
00074             *, integer *, complex *, real *, real *, complex *, real *, 
00075             integer *), cpoequ_(integer *, complex *, integer *, real 
00076             *, real *, real *, integer *), cpbtrs_(char *, integer *, integer 
00077             *, integer *, complex *, integer *, complex *, integer *, integer 
00078             *), cporfs_(char *, integer *, integer *, complex *, 
00079             integer *, complex *, integer *, complex *, integer *, complex *, 
00080             integer *, real *, real *, complex *, real *, integer *), 
00081             cpotrf_(char *, integer *, complex *, integer *, integer *), cpotri_(char *, integer *, complex *, integer *, integer 
00082             *), cppequ_(char *, integer *, complex *, real *, real *, 
00083             real *, integer *), cpprfs_(char *, integer *, integer *, 
00084             complex *, complex *, complex *, integer *, complex *, integer *, 
00085             real *, real *, complex *, real *, integer *), cpptrf_(
00086             char *, integer *, complex *, integer *), cpptri_(char *, 
00087             integer *, complex *, integer *), cpotrs_(char *, integer 
00088             *, integer *, complex *, integer *, complex *, integer *, integer 
00089             *), cpptrs_(char *, integer *, integer *, complex *, 
00090             complex *, integer *, integer *);
00091 
00092     /* Fortran I/O blocks */
00093     static cilist io___1 = { 0, 0, 0, 0, 0 };
00094 
00095 
00096 
00097 /*  -- LAPACK test routine (version 3.1) -- */
00098 /*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
00099 /*     November 2006 */
00100 
00101 /*     .. Scalar Arguments .. */
00102 /*     .. */
00103 
00104 /*  Purpose */
00105 /*  ======= */
00106 
00107 /*  CERRPO tests the error exits for the COMPLEX routines */
00108 /*  for Hermitian positive definite matrices. */
00109 
00110 /*  Arguments */
00111 /*  ========= */
00112 
00113 /*  PATH    (input) CHARACTER*3 */
00114 /*          The LAPACK path name for the routines to be tested. */
00115 
00116 /*  NUNIT   (input) INTEGER */
00117 /*          The unit number for output. */
00118 
00119 /*  ===================================================================== */
00120 
00121 /*     .. Parameters .. */
00122 /*     .. */
00123 /*     .. Local Scalars .. */
00124 /*     .. */
00125 /*     .. Local Arrays .. */
00126 /*     .. */
00127 /*     .. External Functions .. */
00128 /*     .. */
00129 /*     .. External Subroutines .. */
00130 /*     .. */
00131 /*     .. Scalars in Common .. */
00132 /*     .. */
00133 /*     .. Common blocks .. */
00134 /*     .. */
00135 /*     .. Intrinsic Functions .. */
00136 /*     .. */
00137 /*     .. Executable Statements .. */
00138 
00139     infoc_1.nout = *nunit;
00140     io___1.ciunit = infoc_1.nout;
00141     s_wsle(&io___1);
00142     e_wsle();
00143     s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00144 
00145 /*     Set the variables to innocuous values. */
00146 
00147     for (j = 1; j <= 4; ++j) {
00148         for (i__ = 1; i__ <= 4; ++i__) {
00149             i__1 = i__ + (j << 2) - 5;
00150             r__1 = 1.f / (real) (i__ + j);
00151             r__2 = -1.f / (real) (i__ + j);
00152             q__1.r = r__1, q__1.i = r__2;
00153             a[i__1].r = q__1.r, a[i__1].i = q__1.i;
00154             i__1 = i__ + (j << 2) - 5;
00155             r__1 = 1.f / (real) (i__ + j);
00156             r__2 = -1.f / (real) (i__ + j);
00157             q__1.r = r__1, q__1.i = r__2;
00158             af[i__1].r = q__1.r, af[i__1].i = q__1.i;
00159 /* L10: */
00160         }
00161         i__1 = j - 1;
00162         b[i__1].r = 0.f, b[i__1].i = 0.f;
00163         r1[j - 1] = 0.f;
00164         r2[j - 1] = 0.f;
00165         i__1 = j - 1;
00166         w[i__1].r = 0.f, w[i__1].i = 0.f;
00167         i__1 = j - 1;
00168         x[i__1].r = 0.f, x[i__1].i = 0.f;
00169 /* L20: */
00170     }
00171     anrm = 1.f;
00172     infoc_1.ok = TRUE_;
00173 
00174 /*     Test error exits of the routines that use the Cholesky */
00175 /*     decomposition of a Hermitian positive definite matrix. */
00176 
00177     if (lsamen_(&c__2, c2, "PO")) {
00178 
00179 /*        CPOTRF */
00180 
00181         s_copy(srnamc_1.srnamt, "CPOTRF", (ftnlen)32, (ftnlen)6);
00182         infoc_1.infot = 1;
00183         cpotrf_("/", &c__0, a, &c__1, &info);
00184         chkxer_("CPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00185                 infoc_1.ok);
00186         infoc_1.infot = 2;
00187         cpotrf_("U", &c_n1, a, &c__1, &info);
00188         chkxer_("CPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00189                 infoc_1.ok);
00190         infoc_1.infot = 4;
00191         cpotrf_("U", &c__2, a, &c__1, &info);
00192         chkxer_("CPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00193                 infoc_1.ok);
00194 
00195 /*        CPOTF2 */
00196 
00197         s_copy(srnamc_1.srnamt, "CPOTF2", (ftnlen)32, (ftnlen)6);
00198         infoc_1.infot = 1;
00199         cpotf2_("/", &c__0, a, &c__1, &info);
00200         chkxer_("CPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00201                 infoc_1.ok);
00202         infoc_1.infot = 2;
00203         cpotf2_("U", &c_n1, a, &c__1, &info);
00204         chkxer_("CPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00205                 infoc_1.ok);
00206         infoc_1.infot = 4;
00207         cpotf2_("U", &c__2, a, &c__1, &info);
00208         chkxer_("CPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00209                 infoc_1.ok);
00210 
00211 /*        CPOTRI */
00212 
00213         s_copy(srnamc_1.srnamt, "CPOTRI", (ftnlen)32, (ftnlen)6);
00214         infoc_1.infot = 1;
00215         cpotri_("/", &c__0, a, &c__1, &info);
00216         chkxer_("CPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00217                 infoc_1.ok);
00218         infoc_1.infot = 2;
00219         cpotri_("U", &c_n1, a, &c__1, &info);
00220         chkxer_("CPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00221                 infoc_1.ok);
00222         infoc_1.infot = 4;
00223         cpotri_("U", &c__2, a, &c__1, &info);
00224         chkxer_("CPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00225                 infoc_1.ok);
00226 
00227 /*        CPOTRS */
00228 
00229         s_copy(srnamc_1.srnamt, "CPOTRS", (ftnlen)32, (ftnlen)6);
00230         infoc_1.infot = 1;
00231         cpotrs_("/", &c__0, &c__0, a, &c__1, b, &c__1, &info);
00232         chkxer_("CPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00233                 infoc_1.ok);
00234         infoc_1.infot = 2;
00235         cpotrs_("U", &c_n1, &c__0, a, &c__1, b, &c__1, &info);
00236         chkxer_("CPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00237                 infoc_1.ok);
00238         infoc_1.infot = 3;
00239         cpotrs_("U", &c__0, &c_n1, a, &c__1, b, &c__1, &info);
00240         chkxer_("CPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00241                 infoc_1.ok);
00242         infoc_1.infot = 5;
00243         cpotrs_("U", &c__2, &c__1, a, &c__1, b, &c__2, &info);
00244         chkxer_("CPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00245                 infoc_1.ok);
00246         infoc_1.infot = 7;
00247         cpotrs_("U", &c__2, &c__1, a, &c__2, b, &c__1, &info);
00248         chkxer_("CPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00249                 infoc_1.ok);
00250 
00251 /*        CPORFS */
00252 
00253         s_copy(srnamc_1.srnamt, "CPORFS", (ftnlen)32, (ftnlen)6);
00254         infoc_1.infot = 1;
00255         cporfs_("/", &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00256                 r1, r2, w, r__, &info);
00257         chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00258                 infoc_1.ok);
00259         infoc_1.infot = 2;
00260         cporfs_("U", &c_n1, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00261                 r1, r2, w, r__, &info);
00262         chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00263                 infoc_1.ok);
00264         infoc_1.infot = 3;
00265         cporfs_("U", &c__0, &c_n1, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00266                 r1, r2, w, r__, &info);
00267         chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00268                 infoc_1.ok);
00269         infoc_1.infot = 5;
00270         cporfs_("U", &c__2, &c__1, a, &c__1, af, &c__2, b, &c__2, x, &c__2, 
00271                 r1, r2, w, r__, &info);
00272         chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00273                 infoc_1.ok);
00274         infoc_1.infot = 7;
00275         cporfs_("U", &c__2, &c__1, a, &c__2, af, &c__1, b, &c__2, x, &c__2, 
00276                 r1, r2, w, r__, &info);
00277         chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00278                 infoc_1.ok);
00279         infoc_1.infot = 9;
00280         cporfs_("U", &c__2, &c__1, a, &c__2, af, &c__2, b, &c__1, x, &c__2, 
00281                 r1, r2, w, r__, &info);
00282         chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00283                 infoc_1.ok);
00284         infoc_1.infot = 11;
00285         cporfs_("U", &c__2, &c__1, a, &c__2, af, &c__2, b, &c__2, x, &c__1, 
00286                 r1, r2, w, r__, &info);
00287         chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00288                 infoc_1.ok);
00289 
00290 /*        CPOCON */
00291 
00292         s_copy(srnamc_1.srnamt, "CPOCON", (ftnlen)32, (ftnlen)6);
00293         infoc_1.infot = 1;
00294         cpocon_("/", &c__0, a, &c__1, &anrm, &rcond, w, r__, &info)
00295                 ;
00296         chkxer_("CPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00297                 infoc_1.ok);
00298         infoc_1.infot = 2;
00299         cpocon_("U", &c_n1, a, &c__1, &anrm, &rcond, w, r__, &info)
00300                 ;
00301         chkxer_("CPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00302                 infoc_1.ok);
00303         infoc_1.infot = 4;
00304         cpocon_("U", &c__2, a, &c__1, &anrm, &rcond, w, r__, &info)
00305                 ;
00306         chkxer_("CPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00307                 infoc_1.ok);
00308         infoc_1.infot = 5;
00309         r__1 = -anrm;
00310         cpocon_("U", &c__1, a, &c__1, &r__1, &rcond, w, r__, &info)
00311                 ;
00312         chkxer_("CPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00313                 infoc_1.ok);
00314 
00315 /*        CPOEQU */
00316 
00317         s_copy(srnamc_1.srnamt, "CPOEQU", (ftnlen)32, (ftnlen)6);
00318         infoc_1.infot = 1;
00319         cpoequ_(&c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00320         chkxer_("CPOEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00321                 infoc_1.ok);
00322         infoc_1.infot = 3;
00323         cpoequ_(&c__2, a, &c__1, r1, &rcond, &anrm, &info);
00324         chkxer_("CPOEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00325                 infoc_1.ok);
00326 
00327 /*     Test error exits of the routines that use the Cholesky */
00328 /*     decomposition of a Hermitian positive definite packed matrix. */
00329 
00330     } else if (lsamen_(&c__2, c2, "PP")) {
00331 
00332 /*        CPPTRF */
00333 
00334         s_copy(srnamc_1.srnamt, "CPPTRF", (ftnlen)32, (ftnlen)6);
00335         infoc_1.infot = 1;
00336         cpptrf_("/", &c__0, a, &info);
00337         chkxer_("CPPTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00338                 infoc_1.ok);
00339         infoc_1.infot = 2;
00340         cpptrf_("U", &c_n1, a, &info);
00341         chkxer_("CPPTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00342                 infoc_1.ok);
00343 
00344 /*        CPPTRI */
00345 
00346         s_copy(srnamc_1.srnamt, "CPPTRI", (ftnlen)32, (ftnlen)6);
00347         infoc_1.infot = 1;
00348         cpptri_("/", &c__0, a, &info);
00349         chkxer_("CPPTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00350                 infoc_1.ok);
00351         infoc_1.infot = 2;
00352         cpptri_("U", &c_n1, a, &info);
00353         chkxer_("CPPTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00354                 infoc_1.ok);
00355 
00356 /*        CPPTRS */
00357 
00358         s_copy(srnamc_1.srnamt, "CPPTRS", (ftnlen)32, (ftnlen)6);
00359         infoc_1.infot = 1;
00360         cpptrs_("/", &c__0, &c__0, a, b, &c__1, &info);
00361         chkxer_("CPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00362                 infoc_1.ok);
00363         infoc_1.infot = 2;
00364         cpptrs_("U", &c_n1, &c__0, a, b, &c__1, &info);
00365         chkxer_("CPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00366                 infoc_1.ok);
00367         infoc_1.infot = 3;
00368         cpptrs_("U", &c__0, &c_n1, a, b, &c__1, &info);
00369         chkxer_("CPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00370                 infoc_1.ok);
00371         infoc_1.infot = 6;
00372         cpptrs_("U", &c__2, &c__1, a, b, &c__1, &info);
00373         chkxer_("CPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00374                 infoc_1.ok);
00375 
00376 /*        CPPRFS */
00377 
00378         s_copy(srnamc_1.srnamt, "CPPRFS", (ftnlen)32, (ftnlen)6);
00379         infoc_1.infot = 1;
00380         cpprfs_("/", &c__0, &c__0, a, af, b, &c__1, x, &c__1, r1, r2, w, r__, 
00381                 &info);
00382         chkxer_("CPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00383                 infoc_1.ok);
00384         infoc_1.infot = 2;
00385         cpprfs_("U", &c_n1, &c__0, a, af, b, &c__1, x, &c__1, r1, r2, w, r__, 
00386                 &info);
00387         chkxer_("CPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00388                 infoc_1.ok);
00389         infoc_1.infot = 3;
00390         cpprfs_("U", &c__0, &c_n1, a, af, b, &c__1, x, &c__1, r1, r2, w, r__, 
00391                 &info);
00392         chkxer_("CPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00393                 infoc_1.ok);
00394         infoc_1.infot = 7;
00395         cpprfs_("U", &c__2, &c__1, a, af, b, &c__1, x, &c__2, r1, r2, w, r__, 
00396                 &info);
00397         chkxer_("CPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00398                 infoc_1.ok);
00399         infoc_1.infot = 9;
00400         cpprfs_("U", &c__2, &c__1, a, af, b, &c__2, x, &c__1, r1, r2, w, r__, 
00401                 &info);
00402         chkxer_("CPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00403                 infoc_1.ok);
00404 
00405 /*        CPPCON */
00406 
00407         s_copy(srnamc_1.srnamt, "CPPCON", (ftnlen)32, (ftnlen)6);
00408         infoc_1.infot = 1;
00409         cppcon_("/", &c__0, a, &anrm, &rcond, w, r__, &info);
00410         chkxer_("CPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00411                 infoc_1.ok);
00412         infoc_1.infot = 2;
00413         cppcon_("U", &c_n1, a, &anrm, &rcond, w, r__, &info);
00414         chkxer_("CPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00415                 infoc_1.ok);
00416         infoc_1.infot = 4;
00417         r__1 = -anrm;
00418         cppcon_("U", &c__1, a, &r__1, &rcond, w, r__, &info);
00419         chkxer_("CPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00420                 infoc_1.ok);
00421 
00422 /*        CPPEQU */
00423 
00424         s_copy(srnamc_1.srnamt, "CPPEQU", (ftnlen)32, (ftnlen)6);
00425         infoc_1.infot = 1;
00426         cppequ_("/", &c__0, a, r1, &rcond, &anrm, &info);
00427         chkxer_("CPPEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00428                 infoc_1.ok);
00429         infoc_1.infot = 2;
00430         cppequ_("U", &c_n1, a, r1, &rcond, &anrm, &info);
00431         chkxer_("CPPEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00432                 infoc_1.ok);
00433 
00434 /*     Test error exits of the routines that use the Cholesky */
00435 /*     decomposition of a Hermitian positive definite band matrix. */
00436 
00437     } else if (lsamen_(&c__2, c2, "PB")) {
00438 
00439 /*        CPBTRF */
00440 
00441         s_copy(srnamc_1.srnamt, "CPBTRF", (ftnlen)32, (ftnlen)6);
00442         infoc_1.infot = 1;
00443         cpbtrf_("/", &c__0, &c__0, a, &c__1, &info);
00444         chkxer_("CPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00445                 infoc_1.ok);
00446         infoc_1.infot = 2;
00447         cpbtrf_("U", &c_n1, &c__0, a, &c__1, &info);
00448         chkxer_("CPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00449                 infoc_1.ok);
00450         infoc_1.infot = 3;
00451         cpbtrf_("U", &c__1, &c_n1, a, &c__1, &info);
00452         chkxer_("CPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00453                 infoc_1.ok);
00454         infoc_1.infot = 5;
00455         cpbtrf_("U", &c__2, &c__1, a, &c__1, &info);
00456         chkxer_("CPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00457                 infoc_1.ok);
00458 
00459 /*        CPBTF2 */
00460 
00461         s_copy(srnamc_1.srnamt, "CPBTF2", (ftnlen)32, (ftnlen)6);
00462         infoc_1.infot = 1;
00463         cpbtf2_("/", &c__0, &c__0, a, &c__1, &info);
00464         chkxer_("CPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00465                 infoc_1.ok);
00466         infoc_1.infot = 2;
00467         cpbtf2_("U", &c_n1, &c__0, a, &c__1, &info);
00468         chkxer_("CPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00469                 infoc_1.ok);
00470         infoc_1.infot = 3;
00471         cpbtf2_("U", &c__1, &c_n1, a, &c__1, &info);
00472         chkxer_("CPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00473                 infoc_1.ok);
00474         infoc_1.infot = 5;
00475         cpbtf2_("U", &c__2, &c__1, a, &c__1, &info);
00476         chkxer_("CPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00477                 infoc_1.ok);
00478 
00479 /*        CPBTRS */
00480 
00481         s_copy(srnamc_1.srnamt, "CPBTRS", (ftnlen)32, (ftnlen)6);
00482         infoc_1.infot = 1;
00483         cpbtrs_("/", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &info);
00484         chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00485                 infoc_1.ok);
00486         infoc_1.infot = 2;
00487         cpbtrs_("U", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, &info);
00488         chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00489                 infoc_1.ok);
00490         infoc_1.infot = 3;
00491         cpbtrs_("U", &c__1, &c_n1, &c__0, a, &c__1, b, &c__1, &info);
00492         chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00493                 infoc_1.ok);
00494         infoc_1.infot = 4;
00495         cpbtrs_("U", &c__0, &c__0, &c_n1, a, &c__1, b, &c__1, &info);
00496         chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00497                 infoc_1.ok);
00498         infoc_1.infot = 6;
00499         cpbtrs_("U", &c__2, &c__1, &c__1, a, &c__1, b, &c__1, &info);
00500         chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00501                 infoc_1.ok);
00502         infoc_1.infot = 8;
00503         cpbtrs_("U", &c__2, &c__0, &c__1, a, &c__1, b, &c__1, &info);
00504         chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00505                 infoc_1.ok);
00506 
00507 /*        CPBRFS */
00508 
00509         s_copy(srnamc_1.srnamt, "CPBRFS", (ftnlen)32, (ftnlen)6);
00510         infoc_1.infot = 1;
00511         cpbrfs_("/", &c__0, &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00512                 c__1, r1, r2, w, r__, &info);
00513         chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00514                 infoc_1.ok);
00515         infoc_1.infot = 2;
00516         cpbrfs_("U", &c_n1, &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00517                 c__1, r1, r2, w, r__, &info);
00518         chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00519                 infoc_1.ok);
00520         infoc_1.infot = 3;
00521         cpbrfs_("U", &c__1, &c_n1, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00522                 c__1, r1, r2, w, r__, &info);
00523         chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00524                 infoc_1.ok);
00525         infoc_1.infot = 4;
00526         cpbrfs_("U", &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, b, &c__1, x, &
00527                 c__1, r1, r2, w, r__, &info);
00528         chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00529                 infoc_1.ok);
00530         infoc_1.infot = 6;
00531         cpbrfs_("U", &c__2, &c__1, &c__1, a, &c__1, af, &c__2, b, &c__2, x, &
00532                 c__2, r1, r2, w, r__, &info);
00533         chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00534                 infoc_1.ok);
00535         infoc_1.infot = 8;
00536         cpbrfs_("U", &c__2, &c__1, &c__1, a, &c__2, af, &c__1, b, &c__2, x, &
00537                 c__2, r1, r2, w, r__, &info);
00538         chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00539                 infoc_1.ok);
00540         infoc_1.infot = 10;
00541         cpbrfs_("U", &c__2, &c__0, &c__1, a, &c__1, af, &c__1, b, &c__1, x, &
00542                 c__2, r1, r2, w, r__, &info);
00543         chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00544                 infoc_1.ok);
00545         infoc_1.infot = 12;
00546         cpbrfs_("U", &c__2, &c__0, &c__1, a, &c__1, af, &c__1, b, &c__2, x, &
00547                 c__1, r1, r2, w, r__, &info);
00548         chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00549                 infoc_1.ok);
00550 
00551 /*        CPBCON */
00552 
00553         s_copy(srnamc_1.srnamt, "CPBCON", (ftnlen)32, (ftnlen)6);
00554         infoc_1.infot = 1;
00555         cpbcon_("/", &c__0, &c__0, a, &c__1, &anrm, &rcond, w, r__, &info);
00556         chkxer_("CPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00557                 infoc_1.ok);
00558         infoc_1.infot = 2;
00559         cpbcon_("U", &c_n1, &c__0, a, &c__1, &anrm, &rcond, w, r__, &info);
00560         chkxer_("CPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00561                 infoc_1.ok);
00562         infoc_1.infot = 3;
00563         cpbcon_("U", &c__1, &c_n1, a, &c__1, &anrm, &rcond, w, r__, &info);
00564         chkxer_("CPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00565                 infoc_1.ok);
00566         infoc_1.infot = 5;
00567         cpbcon_("U", &c__2, &c__1, a, &c__1, &anrm, &rcond, w, r__, &info);
00568         chkxer_("CPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00569                 infoc_1.ok);
00570         infoc_1.infot = 6;
00571         r__1 = -anrm;
00572         cpbcon_("U", &c__1, &c__0, a, &c__1, &r__1, &rcond, w, r__, &info);
00573         chkxer_("CPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00574                 infoc_1.ok);
00575 
00576 /*        CPBEQU */
00577 
00578         s_copy(srnamc_1.srnamt, "CPBEQU", (ftnlen)32, (ftnlen)6);
00579         infoc_1.infot = 1;
00580         cpbequ_("/", &c__0, &c__0, a, &c__1, r1, &rcond, &anrm, &info);
00581         chkxer_("CPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00582                 infoc_1.ok);
00583         infoc_1.infot = 2;
00584         cpbequ_("U", &c_n1, &c__0, a, &c__1, r1, &rcond, &anrm, &info);
00585         chkxer_("CPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00586                 infoc_1.ok);
00587         infoc_1.infot = 3;
00588         cpbequ_("U", &c__1, &c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00589         chkxer_("CPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00590                 infoc_1.ok);
00591         infoc_1.infot = 5;
00592         cpbequ_("U", &c__2, &c__1, a, &c__1, r1, &rcond, &anrm, &info);
00593         chkxer_("CPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00594                 infoc_1.ok);
00595     }
00596 
00597 /*     Print a summary line. */
00598 
00599     alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00600 
00601     return 0;
00602 
00603 /*     End of CERRPO */
00604 
00605 } /* cerrpo_ */


swiftnav
Author(s):
autogenerated on Sat Jun 8 2019 18:55:24