derrhs.c
Go to the documentation of this file.
00001 /* derrhs.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 static integer c__3 = 3;
00038 static integer c__4 = 4;
00039 
00040 /* Subroutine */ int derrhs_(char *path, integer *nunit)
00041 {
00042     /* Format strings */
00043     static char fmt_9999[] = "(1x,a3,\002 routines passed the tests of the e"
00044             "rror exits\002,\002 (\002,i3,\002 tests done)\002)";
00045     static char fmt_9998[] = "(\002 *** \002,a3,\002 routines failed the tes"
00046             "ts of the error \002,\002exits ***\002)";
00047 
00048     /* Builtin functions */
00049     integer s_wsle(cilist *), e_wsle(void);
00050     /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
00051     integer s_wsfe(cilist *), do_fio(integer *, char *, ftnlen), e_wsfe(void);
00052 
00053     /* Local variables */
00054     doublereal a[9]     /* was [3][3] */, c__[9]        /* was [3][3] */;
00055     integer i__, j, m;
00056     doublereal s[3], w[28];
00057     char c2[2];
00058     doublereal wi[3];
00059     integer nt;
00060     doublereal vl[9]    /* was [3][3] */, vr[9] /* was [3][3] */, wr[3];
00061     integer ihi, ilo;
00062     logical sel[3];
00063     doublereal tau[3];
00064     integer info;
00065     extern /* Subroutine */ int dgebak_(char *, char *, integer *, integer *, 
00066             integer *, doublereal *, integer *, doublereal *, integer *, 
00067             integer *), dgebal_(char *, integer *, doublereal 
00068             *, integer *, integer *, integer *, doublereal *, integer *), dgehrd_(integer *, integer *, integer *, doublereal *, 
00069             integer *, doublereal *, doublereal *, integer *, integer *);
00070     integer ifaill[3], ifailr[3];
00071     extern /* Subroutine */ int dhsein_(char *, char *, char *, logical *, 
00072             integer *, doublereal *, integer *, doublereal *, doublereal *, 
00073             doublereal *, integer *, doublereal *, integer *, integer *, 
00074             integer *, doublereal *, integer *, integer *, integer *);
00075     extern logical lsamen_(integer *, char *, char *);
00076     extern /* Subroutine */ int chkxer_(char *, integer *, integer *, logical 
00077             *, logical *), dorghr_(integer *, integer *, integer *, 
00078             doublereal *, integer *, doublereal *, doublereal *, integer *, 
00079             integer *), dhseqr_(char *, char *, integer *, integer *, integer 
00080             *, doublereal *, integer *, doublereal *, doublereal *, 
00081             doublereal *, integer *, doublereal *, integer *, integer *), dtrevc_(char *, char *, logical *, integer *, 
00082             doublereal *, integer *, doublereal *, integer *, doublereal *, 
00083             integer *, integer *, integer *, doublereal *, integer *), dormhr_(char *, char *, integer *, integer *, integer *, 
00084             integer *, doublereal *, integer *, doublereal *, doublereal *, 
00085             integer *, doublereal *, integer *, integer *);
00086 
00087     /* Fortran I/O blocks */
00088     static cilist io___1 = { 0, 0, 0, 0, 0 };
00089     static cilist io___22 = { 0, 0, 0, fmt_9999, 0 };
00090     static cilist io___23 = { 0, 0, 0, fmt_9998, 0 };
00091 
00092 
00093 
00094 /*  -- LAPACK test routine (version 3.1) -- */
00095 /*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
00096 /*     November 2006 */
00097 
00098 /*     .. Scalar Arguments .. */
00099 /*     .. */
00100 
00101 /*  Purpose */
00102 /*  ======= */
00103 
00104 /*  DERRHS tests the error exits for DGEBAK, SGEBAL, SGEHRD, DORGHR, */
00105 /*  DORMHR, DHSEQR, SHSEIN, and DTREVC. */
00106 
00107 /*  Arguments */
00108 /*  ========= */
00109 
00110 /*  PATH    (input) CHARACTER*3 */
00111 /*          The LAPACK path name for the routines to be tested. */
00112 
00113 /*  NUNIT   (input) INTEGER */
00114 /*          The unit number for output. */
00115 
00116 /*  ===================================================================== */
00117 
00118 /*     .. Parameters .. */
00119 /*     .. */
00120 /*     .. Local Scalars .. */
00121 /*     .. */
00122 /*     .. Local Arrays .. */
00123 /*     .. */
00124 /*     .. External Functions .. */
00125 /*     .. */
00126 /*     .. External Subroutines .. */
00127 /*     .. */
00128 /*     .. Intrinsic Functions .. */
00129 /*     .. */
00130 /*     .. Scalars in Common .. */
00131 /*     .. */
00132 /*     .. Common blocks .. */
00133 /*     .. */
00134 /*     .. Executable Statements .. */
00135 
00136     infoc_1.nout = *nunit;
00137     io___1.ciunit = infoc_1.nout;
00138     s_wsle(&io___1);
00139     e_wsle();
00140     s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00141 
00142 /*     Set the variables to innocuous values. */
00143 
00144     for (j = 1; j <= 3; ++j) {
00145         for (i__ = 1; i__ <= 3; ++i__) {
00146             a[i__ + j * 3 - 4] = 1. / (doublereal) (i__ + j);
00147 /* L10: */
00148         }
00149         wi[j - 1] = (doublereal) j;
00150         sel[j - 1] = TRUE_;
00151 /* L20: */
00152     }
00153     infoc_1.ok = TRUE_;
00154     nt = 0;
00155 
00156 /*     Test error exits of the nonsymmetric eigenvalue routines. */
00157 
00158     if (lsamen_(&c__2, c2, "HS")) {
00159 
00160 /*        DGEBAL */
00161 
00162         s_copy(srnamc_1.srnamt, "DGEBAL", (ftnlen)32, (ftnlen)6);
00163         infoc_1.infot = 1;
00164         dgebal_("/", &c__0, a, &c__1, &ilo, &ihi, s, &info);
00165         chkxer_("DGEBAL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00166                 infoc_1.ok);
00167         infoc_1.infot = 2;
00168         dgebal_("N", &c_n1, a, &c__1, &ilo, &ihi, s, &info);
00169         chkxer_("DGEBAL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00170                 infoc_1.ok);
00171         infoc_1.infot = 4;
00172         dgebal_("N", &c__2, a, &c__1, &ilo, &ihi, s, &info);
00173         chkxer_("DGEBAL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00174                 infoc_1.ok);
00175         nt += 3;
00176 
00177 /*        DGEBAK */
00178 
00179         s_copy(srnamc_1.srnamt, "DGEBAK", (ftnlen)32, (ftnlen)6);
00180         infoc_1.infot = 1;
00181         dgebak_("/", "R", &c__0, &c__1, &c__0, s, &c__0, a, &c__1, &info);
00182         chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00183                 infoc_1.ok);
00184         infoc_1.infot = 2;
00185         dgebak_("N", "/", &c__0, &c__1, &c__0, s, &c__0, a, &c__1, &info);
00186         chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00187                 infoc_1.ok);
00188         infoc_1.infot = 3;
00189         dgebak_("N", "R", &c_n1, &c__1, &c__0, s, &c__0, a, &c__1, &info);
00190         chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00191                 infoc_1.ok);
00192         infoc_1.infot = 4;
00193         dgebak_("N", "R", &c__0, &c__0, &c__0, s, &c__0, a, &c__1, &info);
00194         chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00195                 infoc_1.ok);
00196         infoc_1.infot = 4;
00197         dgebak_("N", "R", &c__0, &c__2, &c__0, s, &c__0, a, &c__1, &info);
00198         chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00199                 infoc_1.ok);
00200         infoc_1.infot = 5;
00201         dgebak_("N", "R", &c__2, &c__2, &c__1, s, &c__0, a, &c__2, &info);
00202         chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00203                 infoc_1.ok);
00204         infoc_1.infot = 5;
00205         dgebak_("N", "R", &c__0, &c__1, &c__1, s, &c__0, a, &c__1, &info);
00206         chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00207                 infoc_1.ok);
00208         infoc_1.infot = 7;
00209         dgebak_("N", "R", &c__0, &c__1, &c__0, s, &c_n1, a, &c__1, &info);
00210         chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00211                 infoc_1.ok);
00212         infoc_1.infot = 9;
00213         dgebak_("N", "R", &c__2, &c__1, &c__2, s, &c__0, a, &c__1, &info);
00214         chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00215                 infoc_1.ok);
00216         nt += 9;
00217 
00218 /*        DGEHRD */
00219 
00220         s_copy(srnamc_1.srnamt, "DGEHRD", (ftnlen)32, (ftnlen)6);
00221         infoc_1.infot = 1;
00222         dgehrd_(&c_n1, &c__1, &c__1, a, &c__1, tau, w, &c__1, &info);
00223         chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00224                 infoc_1.ok);
00225         infoc_1.infot = 2;
00226         dgehrd_(&c__0, &c__0, &c__0, a, &c__1, tau, w, &c__1, &info);
00227         chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00228                 infoc_1.ok);
00229         infoc_1.infot = 2;
00230         dgehrd_(&c__0, &c__2, &c__0, a, &c__1, tau, w, &c__1, &info);
00231         chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00232                 infoc_1.ok);
00233         infoc_1.infot = 3;
00234         dgehrd_(&c__1, &c__1, &c__0, a, &c__1, tau, w, &c__1, &info);
00235         chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00236                 infoc_1.ok);
00237         infoc_1.infot = 3;
00238         dgehrd_(&c__0, &c__1, &c__1, a, &c__1, tau, w, &c__1, &info);
00239         chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00240                 infoc_1.ok);
00241         infoc_1.infot = 5;
00242         dgehrd_(&c__2, &c__1, &c__1, a, &c__1, tau, w, &c__2, &info);
00243         chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00244                 infoc_1.ok);
00245         infoc_1.infot = 8;
00246         dgehrd_(&c__2, &c__1, &c__2, a, &c__2, tau, w, &c__1, &info);
00247         chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00248                 infoc_1.ok);
00249         nt += 7;
00250 
00251 /*        DORGHR */
00252 
00253         s_copy(srnamc_1.srnamt, "DORGHR", (ftnlen)32, (ftnlen)6);
00254         infoc_1.infot = 1;
00255         dorghr_(&c_n1, &c__1, &c__1, a, &c__1, tau, w, &c__1, &info);
00256         chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00257                 infoc_1.ok);
00258         infoc_1.infot = 2;
00259         dorghr_(&c__0, &c__0, &c__0, a, &c__1, tau, w, &c__1, &info);
00260         chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00261                 infoc_1.ok);
00262         infoc_1.infot = 2;
00263         dorghr_(&c__0, &c__2, &c__0, a, &c__1, tau, w, &c__1, &info);
00264         chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00265                 infoc_1.ok);
00266         infoc_1.infot = 3;
00267         dorghr_(&c__1, &c__1, &c__0, a, &c__1, tau, w, &c__1, &info);
00268         chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00269                 infoc_1.ok);
00270         infoc_1.infot = 3;
00271         dorghr_(&c__0, &c__1, &c__1, a, &c__1, tau, w, &c__1, &info);
00272         chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00273                 infoc_1.ok);
00274         infoc_1.infot = 5;
00275         dorghr_(&c__2, &c__1, &c__1, a, &c__1, tau, w, &c__1, &info);
00276         chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00277                 infoc_1.ok);
00278         infoc_1.infot = 8;
00279         dorghr_(&c__3, &c__1, &c__3, a, &c__3, tau, w, &c__1, &info);
00280         chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00281                 infoc_1.ok);
00282         nt += 7;
00283 
00284 /*        DORMHR */
00285 
00286         s_copy(srnamc_1.srnamt, "DORMHR", (ftnlen)32, (ftnlen)6);
00287         infoc_1.infot = 1;
00288         dormhr_("/", "N", &c__0, &c__0, &c__1, &c__0, a, &c__1, tau, c__, &
00289                 c__1, w, &c__1, &info);
00290         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00291                 infoc_1.ok);
00292         infoc_1.infot = 2;
00293         dormhr_("L", "/", &c__0, &c__0, &c__1, &c__0, a, &c__1, tau, c__, &
00294                 c__1, w, &c__1, &info);
00295         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00296                 infoc_1.ok);
00297         infoc_1.infot = 3;
00298         dormhr_("L", "N", &c_n1, &c__0, &c__1, &c__0, a, &c__1, tau, c__, &
00299                 c__1, w, &c__1, &info);
00300         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00301                 infoc_1.ok);
00302         infoc_1.infot = 4;
00303         dormhr_("L", "N", &c__0, &c_n1, &c__1, &c__0, a, &c__1, tau, c__, &
00304                 c__1, w, &c__1, &info);
00305         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00306                 infoc_1.ok);
00307         infoc_1.infot = 5;
00308         dormhr_("L", "N", &c__0, &c__0, &c__0, &c__0, a, &c__1, tau, c__, &
00309                 c__1, w, &c__1, &info);
00310         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00311                 infoc_1.ok);
00312         infoc_1.infot = 5;
00313         dormhr_("L", "N", &c__0, &c__0, &c__2, &c__0, a, &c__1, tau, c__, &
00314                 c__1, w, &c__1, &info);
00315         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00316                 infoc_1.ok);
00317         infoc_1.infot = 5;
00318         dormhr_("L", "N", &c__1, &c__2, &c__2, &c__1, a, &c__1, tau, c__, &
00319                 c__1, w, &c__2, &info);
00320         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00321                 infoc_1.ok);
00322         infoc_1.infot = 5;
00323         dormhr_("R", "N", &c__2, &c__1, &c__2, &c__1, a, &c__1, tau, c__, &
00324                 c__2, w, &c__2, &info);
00325         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00326                 infoc_1.ok);
00327         infoc_1.infot = 6;
00328         dormhr_("L", "N", &c__1, &c__1, &c__1, &c__0, a, &c__1, tau, c__, &
00329                 c__1, w, &c__1, &info);
00330         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00331                 infoc_1.ok);
00332         infoc_1.infot = 6;
00333         dormhr_("L", "N", &c__0, &c__1, &c__1, &c__1, a, &c__1, tau, c__, &
00334                 c__1, w, &c__1, &info);
00335         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00336                 infoc_1.ok);
00337         infoc_1.infot = 6;
00338         dormhr_("R", "N", &c__1, &c__0, &c__1, &c__1, a, &c__1, tau, c__, &
00339                 c__1, w, &c__1, &info);
00340         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00341                 infoc_1.ok);
00342         infoc_1.infot = 8;
00343         dormhr_("L", "N", &c__2, &c__1, &c__1, &c__1, a, &c__1, tau, c__, &
00344                 c__2, w, &c__1, &info);
00345         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00346                 infoc_1.ok);
00347         infoc_1.infot = 8;
00348         dormhr_("R", "N", &c__1, &c__2, &c__1, &c__1, a, &c__1, tau, c__, &
00349                 c__1, w, &c__1, &info);
00350         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00351                 infoc_1.ok);
00352         infoc_1.infot = 11;
00353         dormhr_("L", "N", &c__2, &c__1, &c__1, &c__1, a, &c__2, tau, c__, &
00354                 c__1, w, &c__1, &info);
00355         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00356                 infoc_1.ok);
00357         infoc_1.infot = 13;
00358         dormhr_("L", "N", &c__1, &c__2, &c__1, &c__1, a, &c__1, tau, c__, &
00359                 c__1, w, &c__1, &info);
00360         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00361                 infoc_1.ok);
00362         infoc_1.infot = 13;
00363         dormhr_("R", "N", &c__2, &c__1, &c__1, &c__1, a, &c__1, tau, c__, &
00364                 c__2, w, &c__1, &info);
00365         chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00366                 infoc_1.ok);
00367         nt += 16;
00368 
00369 /*        DHSEQR */
00370 
00371         s_copy(srnamc_1.srnamt, "DHSEQR", (ftnlen)32, (ftnlen)6);
00372         infoc_1.infot = 1;
00373         dhseqr_("/", "N", &c__0, &c__1, &c__0, a, &c__1, wr, wi, c__, &c__1, 
00374                 w, &c__1, &info);
00375         chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00376                 infoc_1.ok);
00377         infoc_1.infot = 2;
00378         dhseqr_("E", "/", &c__0, &c__1, &c__0, a, &c__1, wr, wi, c__, &c__1, 
00379                 w, &c__1, &info);
00380         chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00381                 infoc_1.ok);
00382         infoc_1.infot = 3;
00383         dhseqr_("E", "N", &c_n1, &c__1, &c__0, a, &c__1, wr, wi, c__, &c__1, 
00384                 w, &c__1, &info);
00385         chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00386                 infoc_1.ok);
00387         infoc_1.infot = 4;
00388         dhseqr_("E", "N", &c__0, &c__0, &c__0, a, &c__1, wr, wi, c__, &c__1, 
00389                 w, &c__1, &info);
00390         chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00391                 infoc_1.ok);
00392         infoc_1.infot = 4;
00393         dhseqr_("E", "N", &c__0, &c__2, &c__0, a, &c__1, wr, wi, c__, &c__1, 
00394                 w, &c__1, &info);
00395         chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00396                 infoc_1.ok);
00397         infoc_1.infot = 5;
00398         dhseqr_("E", "N", &c__1, &c__1, &c__0, a, &c__1, wr, wi, c__, &c__1, 
00399                 w, &c__1, &info);
00400         chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00401                 infoc_1.ok);
00402         infoc_1.infot = 5;
00403         dhseqr_("E", "N", &c__1, &c__1, &c__2, a, &c__1, wr, wi, c__, &c__1, 
00404                 w, &c__1, &info);
00405         chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00406                 infoc_1.ok);
00407         infoc_1.infot = 7;
00408         dhseqr_("E", "N", &c__2, &c__1, &c__2, a, &c__1, wr, wi, c__, &c__2, 
00409                 w, &c__1, &info);
00410         chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00411                 infoc_1.ok);
00412         infoc_1.infot = 11;
00413         dhseqr_("E", "V", &c__2, &c__1, &c__2, a, &c__2, wr, wi, c__, &c__1, 
00414                 w, &c__1, &info);
00415         chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00416                 infoc_1.ok);
00417         nt += 9;
00418 
00419 /*        DHSEIN */
00420 
00421         s_copy(srnamc_1.srnamt, "DHSEIN", (ftnlen)32, (ftnlen)6);
00422         infoc_1.infot = 1;
00423         dhsein_("/", "N", "N", sel, &c__0, a, &c__1, wr, wi, vl, &c__1, vr, &
00424                 c__1, &c__0, &m, w, ifaill, ifailr, &info);
00425         chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00426                 infoc_1.ok);
00427         infoc_1.infot = 2;
00428         dhsein_("R", "/", "N", sel, &c__0, a, &c__1, wr, wi, vl, &c__1, vr, &
00429                 c__1, &c__0, &m, w, ifaill, ifailr, &info);
00430         chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00431                 infoc_1.ok);
00432         infoc_1.infot = 3;
00433         dhsein_("R", "N", "/", sel, &c__0, a, &c__1, wr, wi, vl, &c__1, vr, &
00434                 c__1, &c__0, &m, w, ifaill, ifailr, &info);
00435         chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00436                 infoc_1.ok);
00437         infoc_1.infot = 5;
00438         dhsein_("R", "N", "N", sel, &c_n1, a, &c__1, wr, wi, vl, &c__1, vr, &
00439                 c__1, &c__0, &m, w, ifaill, ifailr, &info);
00440         chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00441                 infoc_1.ok);
00442         infoc_1.infot = 7;
00443         dhsein_("R", "N", "N", sel, &c__2, a, &c__1, wr, wi, vl, &c__1, vr, &
00444                 c__2, &c__4, &m, w, ifaill, ifailr, &info);
00445         chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00446                 infoc_1.ok);
00447         infoc_1.infot = 11;
00448         dhsein_("L", "N", "N", sel, &c__2, a, &c__2, wr, wi, vl, &c__1, vr, &
00449                 c__1, &c__4, &m, w, ifaill, ifailr, &info);
00450         chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00451                 infoc_1.ok);
00452         infoc_1.infot = 13;
00453         dhsein_("R", "N", "N", sel, &c__2, a, &c__2, wr, wi, vl, &c__1, vr, &
00454                 c__1, &c__4, &m, w, ifaill, ifailr, &info);
00455         chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00456                 infoc_1.ok);
00457         infoc_1.infot = 14;
00458         dhsein_("R", "N", "N", sel, &c__2, a, &c__2, wr, wi, vl, &c__1, vr, &
00459                 c__2, &c__1, &m, w, ifaill, ifailr, &info);
00460         chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00461                 infoc_1.ok);
00462         nt += 8;
00463 
00464 /*        DTREVC */
00465 
00466         s_copy(srnamc_1.srnamt, "DTREVC", (ftnlen)32, (ftnlen)6);
00467         infoc_1.infot = 1;
00468         dtrevc_("/", "A", sel, &c__0, a, &c__1, vl, &c__1, vr, &c__1, &c__0, &
00469                 m, w, &info);
00470         chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00471                 infoc_1.ok);
00472         infoc_1.infot = 2;
00473         dtrevc_("L", "/", sel, &c__0, a, &c__1, vl, &c__1, vr, &c__1, &c__0, &
00474                 m, w, &info);
00475         chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00476                 infoc_1.ok);
00477         infoc_1.infot = 4;
00478         dtrevc_("L", "A", sel, &c_n1, a, &c__1, vl, &c__1, vr, &c__1, &c__0, &
00479                 m, w, &info);
00480         chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00481                 infoc_1.ok);
00482         infoc_1.infot = 6;
00483         dtrevc_("L", "A", sel, &c__2, a, &c__1, vl, &c__2, vr, &c__1, &c__4, &
00484                 m, w, &info);
00485         chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00486                 infoc_1.ok);
00487         infoc_1.infot = 8;
00488         dtrevc_("L", "A", sel, &c__2, a, &c__2, vl, &c__1, vr, &c__1, &c__4, &
00489                 m, w, &info);
00490         chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00491                 infoc_1.ok);
00492         infoc_1.infot = 10;
00493         dtrevc_("R", "A", sel, &c__2, a, &c__2, vl, &c__1, vr, &c__1, &c__4, &
00494                 m, w, &info);
00495         chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00496                 infoc_1.ok);
00497         infoc_1.infot = 11;
00498         dtrevc_("L", "A", sel, &c__2, a, &c__2, vl, &c__2, vr, &c__1, &c__1, &
00499                 m, w, &info);
00500         chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00501                 infoc_1.ok);
00502         nt += 7;
00503     }
00504 
00505 /*     Print a summary line. */
00506 
00507     if (infoc_1.ok) {
00508         io___22.ciunit = infoc_1.nout;
00509         s_wsfe(&io___22);
00510         do_fio(&c__1, path, (ftnlen)3);
00511         do_fio(&c__1, (char *)&nt, (ftnlen)sizeof(integer));
00512         e_wsfe();
00513     } else {
00514         io___23.ciunit = infoc_1.nout;
00515         s_wsfe(&io___23);
00516         do_fio(&c__1, path, (ftnlen)3);
00517         e_wsfe();
00518     }
00519 
00520 
00521     return 0;
00522 
00523 /*     End of DERRHS */
00524 
00525 } /* derrhs_ */


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