00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013 #include "f2c.h"
00014 #include "blaswrap.h"
00015
00016
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
00032
00033 static integer c_n1 = -1;
00034 static integer c__0 = 0;
00035 static integer c__1 = 1;
00036 static integer c__2 = 2;
00037
00038 int serrlq_(char *path, integer *nunit)
00039 {
00040
00041 integer s_wsle(cilist *), e_wsle(void);
00042 int s_copy(char *, char *, ftnlen, ftnlen);
00043
00044
00045 real a[4] , b[2];
00046 integer i__, j;
00047 real w[2], x[2], af[4] ;
00048 integer info;
00049 extern int sgelq2_(integer *, integer *, real *, integer
00050 *, real *, real *, integer *), sorgl2_(integer *, integer *,
00051 integer *, real *, integer *, real *, real *, integer *), sorml2_(
00052 char *, char *, integer *, integer *, integer *, real *, integer *
00053 , real *, real *, integer *, real *, integer *),
00054 alaesm_(char *, logical *, integer *), sgelqf_(integer *,
00055 integer *, real *, integer *, real *, real *, integer *, integer *
00056 ), chkxer_(char *, integer *, integer *, logical *, logical *), sgelqs_(integer *, integer *, integer *, real *, integer
00057 *, real *, real *, integer *, real *, integer *, integer *),
00058 sorglq_(integer *, integer *, integer *, real *, integer *, real *
00059 , real *, integer *, integer *), sormlq_(char *, char *, integer *
00060 , integer *, integer *, real *, integer *, real *, real *,
00061 integer *, real *, integer *, integer *);
00062
00063
00064 static cilist io___1 = { 0, 0, 0, 0, 0 };
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102
00103
00104
00105
00106
00107
00108 infoc_1.nout = *nunit;
00109 io___1.ciunit = infoc_1.nout;
00110 s_wsle(&io___1);
00111 e_wsle();
00112
00113
00114
00115 for (j = 1; j <= 2; ++j) {
00116 for (i__ = 1; i__ <= 2; ++i__) {
00117 a[i__ + (j << 1) - 3] = 1.f / (real) (i__ + j);
00118 af[i__ + (j << 1) - 3] = 1.f / (real) (i__ + j);
00119
00120 }
00121 b[j - 1] = 0.f;
00122 w[j - 1] = 0.f;
00123 x[j - 1] = 0.f;
00124
00125 }
00126 infoc_1.ok = TRUE_;
00127
00128
00129
00130
00131
00132 s_copy(srnamc_1.srnamt, "SGELQF", (ftnlen)32, (ftnlen)6);
00133 infoc_1.infot = 1;
00134 sgelqf_(&c_n1, &c__0, a, &c__1, b, w, &c__1, &info);
00135 chkxer_("SGELQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00136 infoc_1.ok);
00137 infoc_1.infot = 2;
00138 sgelqf_(&c__0, &c_n1, a, &c__1, b, w, &c__1, &info);
00139 chkxer_("SGELQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00140 infoc_1.ok);
00141 infoc_1.infot = 4;
00142 sgelqf_(&c__2, &c__1, a, &c__1, b, w, &c__2, &info);
00143 chkxer_("SGELQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00144 infoc_1.ok);
00145 infoc_1.infot = 7;
00146 sgelqf_(&c__2, &c__1, a, &c__2, b, w, &c__1, &info);
00147 chkxer_("SGELQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00148 infoc_1.ok);
00149
00150
00151
00152 s_copy(srnamc_1.srnamt, "SGELQ2", (ftnlen)32, (ftnlen)6);
00153 infoc_1.infot = 1;
00154 sgelq2_(&c_n1, &c__0, a, &c__1, b, w, &info);
00155 chkxer_("SGELQ2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00156 infoc_1.ok);
00157 infoc_1.infot = 2;
00158 sgelq2_(&c__0, &c_n1, a, &c__1, b, w, &info);
00159 chkxer_("SGELQ2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00160 infoc_1.ok);
00161 infoc_1.infot = 4;
00162 sgelq2_(&c__2, &c__1, a, &c__1, b, w, &info);
00163 chkxer_("SGELQ2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00164 infoc_1.ok);
00165
00166
00167
00168 s_copy(srnamc_1.srnamt, "SGELQS", (ftnlen)32, (ftnlen)6);
00169 infoc_1.infot = 1;
00170 sgelqs_(&c_n1, &c__0, &c__0, a, &c__1, x, b, &c__1, w, &c__1, &info);
00171 chkxer_("SGELQS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00172 infoc_1.ok);
00173 infoc_1.infot = 2;
00174 sgelqs_(&c__0, &c_n1, &c__0, a, &c__1, x, b, &c__1, w, &c__1, &info);
00175 chkxer_("SGELQS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00176 infoc_1.ok);
00177 infoc_1.infot = 2;
00178 sgelqs_(&c__2, &c__1, &c__0, a, &c__2, x, b, &c__1, w, &c__1, &info);
00179 chkxer_("SGELQS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00180 infoc_1.ok);
00181 infoc_1.infot = 3;
00182 sgelqs_(&c__0, &c__0, &c_n1, a, &c__1, x, b, &c__1, w, &c__1, &info);
00183 chkxer_("SGELQS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00184 infoc_1.ok);
00185 infoc_1.infot = 5;
00186 sgelqs_(&c__2, &c__2, &c__0, a, &c__1, x, b, &c__2, w, &c__1, &info);
00187 chkxer_("SGELQS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00188 infoc_1.ok);
00189 infoc_1.infot = 8;
00190 sgelqs_(&c__1, &c__2, &c__0, a, &c__1, x, b, &c__1, w, &c__1, &info);
00191 chkxer_("SGELQS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00192 infoc_1.ok);
00193 infoc_1.infot = 10;
00194 sgelqs_(&c__1, &c__1, &c__2, a, &c__1, x, b, &c__1, w, &c__1, &info);
00195 chkxer_("SGELQS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00196 infoc_1.ok);
00197
00198
00199
00200 s_copy(srnamc_1.srnamt, "SORGLQ", (ftnlen)32, (ftnlen)6);
00201 infoc_1.infot = 1;
00202 sorglq_(&c_n1, &c__0, &c__0, a, &c__1, x, w, &c__1, &info);
00203 chkxer_("SORGLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00204 infoc_1.ok);
00205 infoc_1.infot = 2;
00206 sorglq_(&c__0, &c_n1, &c__0, a, &c__1, x, w, &c__1, &info);
00207 chkxer_("SORGLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00208 infoc_1.ok);
00209 infoc_1.infot = 2;
00210 sorglq_(&c__2, &c__1, &c__0, a, &c__2, x, w, &c__2, &info);
00211 chkxer_("SORGLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00212 infoc_1.ok);
00213 infoc_1.infot = 3;
00214 sorglq_(&c__0, &c__0, &c_n1, a, &c__1, x, w, &c__1, &info);
00215 chkxer_("SORGLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00216 infoc_1.ok);
00217 infoc_1.infot = 3;
00218 sorglq_(&c__1, &c__1, &c__2, a, &c__1, x, w, &c__1, &info);
00219 chkxer_("SORGLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00220 infoc_1.ok);
00221 infoc_1.infot = 5;
00222 sorglq_(&c__2, &c__2, &c__0, a, &c__1, x, w, &c__2, &info);
00223 chkxer_("SORGLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00224 infoc_1.ok);
00225 infoc_1.infot = 8;
00226 sorglq_(&c__2, &c__2, &c__0, a, &c__2, x, w, &c__1, &info);
00227 chkxer_("SORGLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00228 infoc_1.ok);
00229
00230
00231
00232 s_copy(srnamc_1.srnamt, "SORGL2", (ftnlen)32, (ftnlen)6);
00233 infoc_1.infot = 1;
00234 sorgl2_(&c_n1, &c__0, &c__0, a, &c__1, x, w, &info);
00235 chkxer_("SORGL2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00236 infoc_1.ok);
00237 infoc_1.infot = 2;
00238 sorgl2_(&c__0, &c_n1, &c__0, a, &c__1, x, w, &info);
00239 chkxer_("SORGL2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00240 infoc_1.ok);
00241 infoc_1.infot = 2;
00242 sorgl2_(&c__2, &c__1, &c__0, a, &c__2, x, w, &info);
00243 chkxer_("SORGL2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00244 infoc_1.ok);
00245 infoc_1.infot = 3;
00246 sorgl2_(&c__0, &c__0, &c_n1, a, &c__1, x, w, &info);
00247 chkxer_("SORGL2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00248 infoc_1.ok);
00249 infoc_1.infot = 3;
00250 sorgl2_(&c__1, &c__1, &c__2, a, &c__1, x, w, &info);
00251 chkxer_("SORGL2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00252 infoc_1.ok);
00253 infoc_1.infot = 5;
00254 sorgl2_(&c__2, &c__2, &c__0, a, &c__1, x, w, &info);
00255 chkxer_("SORGL2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00256 infoc_1.ok);
00257
00258
00259
00260 s_copy(srnamc_1.srnamt, "SORMLQ", (ftnlen)32, (ftnlen)6);
00261 infoc_1.infot = 1;
00262 sormlq_("/", "N", &c__0, &c__0, &c__0, a, &c__1, x, af, &c__1, w, &c__1, &
00263 info);
00264 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00265 infoc_1.ok);
00266 infoc_1.infot = 2;
00267 sormlq_("L", "/", &c__0, &c__0, &c__0, a, &c__1, x, af, &c__1, w, &c__1, &
00268 info);
00269 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00270 infoc_1.ok);
00271 infoc_1.infot = 3;
00272 sormlq_("L", "N", &c_n1, &c__0, &c__0, a, &c__1, x, af, &c__1, w, &c__1, &
00273 info);
00274 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00275 infoc_1.ok);
00276 infoc_1.infot = 4;
00277 sormlq_("L", "N", &c__0, &c_n1, &c__0, a, &c__1, x, af, &c__1, w, &c__1, &
00278 info);
00279 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00280 infoc_1.ok);
00281 infoc_1.infot = 5;
00282 sormlq_("L", "N", &c__0, &c__0, &c_n1, a, &c__1, x, af, &c__1, w, &c__1, &
00283 info);
00284 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00285 infoc_1.ok);
00286 infoc_1.infot = 5;
00287 sormlq_("L", "N", &c__0, &c__1, &c__1, a, &c__1, x, af, &c__1, w, &c__1, &
00288 info);
00289 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00290 infoc_1.ok);
00291 infoc_1.infot = 5;
00292 sormlq_("R", "N", &c__1, &c__0, &c__1, a, &c__1, x, af, &c__1, w, &c__1, &
00293 info);
00294 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00295 infoc_1.ok);
00296 infoc_1.infot = 7;
00297 sormlq_("L", "N", &c__2, &c__0, &c__2, a, &c__1, x, af, &c__2, w, &c__1, &
00298 info);
00299 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00300 infoc_1.ok);
00301 infoc_1.infot = 7;
00302 sormlq_("R", "N", &c__0, &c__2, &c__2, a, &c__1, x, af, &c__1, w, &c__1, &
00303 info);
00304 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00305 infoc_1.ok);
00306 infoc_1.infot = 10;
00307 sormlq_("L", "N", &c__2, &c__1, &c__0, a, &c__2, x, af, &c__1, w, &c__1, &
00308 info);
00309 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00310 infoc_1.ok);
00311 infoc_1.infot = 12;
00312 sormlq_("L", "N", &c__1, &c__2, &c__0, a, &c__1, x, af, &c__1, w, &c__1, &
00313 info);
00314 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00315 infoc_1.ok);
00316 infoc_1.infot = 12;
00317 sormlq_("R", "N", &c__2, &c__1, &c__0, a, &c__1, x, af, &c__2, w, &c__1, &
00318 info);
00319 chkxer_("SORMLQ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00320 infoc_1.ok);
00321
00322
00323
00324 s_copy(srnamc_1.srnamt, "SORML2", (ftnlen)32, (ftnlen)6);
00325 infoc_1.infot = 1;
00326 sorml2_("/", "N", &c__0, &c__0, &c__0, a, &c__1, x, af, &c__1, w, &info);
00327 chkxer_("SORML2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00328 infoc_1.ok);
00329 infoc_1.infot = 2;
00330 sorml2_("L", "/", &c__0, &c__0, &c__0, a, &c__1, x, af, &c__1, w, &info);
00331 chkxer_("SORML2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00332 infoc_1.ok);
00333 infoc_1.infot = 3;
00334 sorml2_("L", "N", &c_n1, &c__0, &c__0, a, &c__1, x, af, &c__1, w, &info);
00335 chkxer_("SORML2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00336 infoc_1.ok);
00337 infoc_1.infot = 4;
00338 sorml2_("L", "N", &c__0, &c_n1, &c__0, a, &c__1, x, af, &c__1, w, &info);
00339 chkxer_("SORML2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00340 infoc_1.ok);
00341 infoc_1.infot = 5;
00342 sorml2_("L", "N", &c__0, &c__0, &c_n1, a, &c__1, x, af, &c__1, w, &info);
00343 chkxer_("SORML2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00344 infoc_1.ok);
00345 infoc_1.infot = 5;
00346 sorml2_("L", "N", &c__0, &c__1, &c__1, a, &c__1, x, af, &c__1, w, &info);
00347 chkxer_("SORML2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00348 infoc_1.ok);
00349 infoc_1.infot = 5;
00350 sorml2_("R", "N", &c__1, &c__0, &c__1, a, &c__1, x, af, &c__1, w, &info);
00351 chkxer_("SORML2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00352 infoc_1.ok);
00353 infoc_1.infot = 7;
00354 sorml2_("L", "N", &c__2, &c__1, &c__2, a, &c__1, x, af, &c__2, w, &info);
00355 chkxer_("SORML2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00356 infoc_1.ok);
00357 infoc_1.infot = 7;
00358 sorml2_("R", "N", &c__1, &c__2, &c__2, a, &c__1, x, af, &c__1, w, &info);
00359 chkxer_("SORML2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00360 infoc_1.ok);
00361 infoc_1.infot = 10;
00362 sorml2_("L", "N", &c__2, &c__1, &c__0, a, &c__2, x, af, &c__1, w, &info);
00363 chkxer_("SORML2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00364 infoc_1.ok);
00365
00366
00367
00368 alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00369
00370 return 0;
00371
00372
00373
00374 }