5 #include "sphinxbase/f2c.h"
10 extern doublereal slamch_(
char *);
11 #define EPSILON slamch_("Epsilon")
12 #define SAFEMINIMUM slamch_("Safe minimum")
13 #define PRECISION slamch_("Precision")
14 #define BASE slamch_("Base")
18 extern doublereal slapy2_(real *, real *);
24 static integer c__1 = 1;
26 logical lsame_(
char *ca,
char *cb)
32 static integer inta, intb, zcode;
61 ret_val = *(
unsigned char *)ca == *(
unsigned char *)cb;
77 inta = *(
unsigned char *)ca;
78 intb = *(
unsigned char *)cb;
80 if (zcode == 90 || zcode == 122) {
87 if (inta >= 97 && inta <= 122) {
90 if (intb >= 97 && intb <= 122) {
94 }
else if (zcode == 233 || zcode == 169) {
101 if (inta >= 129 && inta <= 137 || inta >= 145 && inta <= 153 || inta
102 >= 162 && inta <= 169) {
105 if (intb >= 129 && intb <= 137 || intb >= 145 && intb <= 153 || intb
106 >= 162 && intb <= 169) {
110 }
else if (zcode == 218 || zcode == 250) {
117 if (inta >= 225 && inta <= 250) {
120 if (intb >= 225 && intb <= 250) {
124 ret_val = inta == intb;
135 doublereal sdot_(integer *n, real *sx, integer *incx, real *sy, integer *incy)
142 static integer i__, m, ix, iy, mp1;
164 if (*incx == 1 && *incy == 1) {
176 ix = (-(*n) + 1) * *incx + 1;
179 iy = (-(*n) + 1) * *incy + 1;
182 for (i__ = 1; i__ <= i__1; ++i__) {
183 stemp += sx[ix] * sy[iy];
204 for (i__ = 1; i__ <= i__1; ++i__) {
205 stemp += sx[i__] * sy[i__];
214 for (i__ = mp1; i__ <= i__1; i__ += 5) {
215 stemp = stemp + sx[i__] * sy[i__] + sx[i__ + 1] * sy[i__ + 1] + sx[
216 i__ + 2] * sy[i__ + 2] + sx[i__ + 3] * sy[i__ + 3] + sx[i__ +
225 int sgemm_(
char *transa,
char *transb, integer *m, integer *
226 n, integer *k, real *alpha, real *a, integer *lda, real *b, integer *
227 ldb, real *beta, real *c__, integer *ldc)
230 integer a_dim1, a_offset, b_dim1, b_offset, c_dim1, c_offset, i__1, i__2,
234 static integer i__, j, l, info;
235 static logical nota, notb;
237 static integer ncola;
238 extern logical lsame_(
char *,
char *);
239 static integer nrowa, nrowb;
240 extern int xerbla_(
char *, integer *);
371 a_offset = 1 + a_dim1;
374 b_offset = 1 + b_dim1;
377 c_offset = 1 + c_dim1;
381 nota = lsame_(transa,
"N");
382 notb = lsame_(transb,
"N");
399 if (! nota && ! lsame_(transa,
"C") && ! lsame_(
402 }
else if (! notb && ! lsame_(transb,
"C") && !
403 lsame_(transb,
"T")) {
411 }
else if (*lda < max(1,nrowa)) {
413 }
else if (*ldb < max(1,nrowb)) {
415 }
else if (*ldc < max(1,*m)) {
419 xerbla_(
"SGEMM ", &info);
425 if (*m == 0 || *n == 0 || (*alpha == 0.f || *k == 0) && *beta == 1.f) {
434 for (j = 1; j <= i__1; ++j) {
436 for (i__ = 1; i__ <= i__2; ++i__) {
437 c__[i__ + j * c_dim1] = 0.f;
444 for (j = 1; j <= i__1; ++j) {
446 for (i__ = 1; i__ <= i__2; ++i__) {
447 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1];
464 for (j = 1; j <= i__1; ++j) {
467 for (i__ = 1; i__ <= i__2; ++i__) {
468 c__[i__ + j * c_dim1] = 0.f;
471 }
else if (*beta != 1.f) {
473 for (i__ = 1; i__ <= i__2; ++i__) {
474 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1];
479 for (l = 1; l <= i__2; ++l) {
480 if (b[l + j * b_dim1] != 0.f) {
481 temp = *alpha * b[l + j * b_dim1];
483 for (i__ = 1; i__ <= i__3; ++i__) {
484 c__[i__ + j * c_dim1] += temp * a[i__ + l *
498 for (j = 1; j <= i__1; ++j) {
500 for (i__ = 1; i__ <= i__2; ++i__) {
503 for (l = 1; l <= i__3; ++l) {
504 temp += a[l + i__ * a_dim1] * b[l + j * b_dim1];
508 c__[i__ + j * c_dim1] = *alpha * temp;
510 c__[i__ + j * c_dim1] = *alpha * temp + *beta * c__[
524 for (j = 1; j <= i__1; ++j) {
527 for (i__ = 1; i__ <= i__2; ++i__) {
528 c__[i__ + j * c_dim1] = 0.f;
531 }
else if (*beta != 1.f) {
533 for (i__ = 1; i__ <= i__2; ++i__) {
534 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1];
539 for (l = 1; l <= i__2; ++l) {
540 if (b[j + l * b_dim1] != 0.f) {
541 temp = *alpha * b[j + l * b_dim1];
543 for (i__ = 1; i__ <= i__3; ++i__) {
544 c__[i__ + j * c_dim1] += temp * a[i__ + l *
558 for (j = 1; j <= i__1; ++j) {
560 for (i__ = 1; i__ <= i__2; ++i__) {
563 for (l = 1; l <= i__3; ++l) {
564 temp += a[l + i__ * a_dim1] * b[j + l * b_dim1];
568 c__[i__ + j * c_dim1] = *alpha * temp;
570 c__[i__ + j * c_dim1] = *alpha * temp + *beta * c__[
586 int sgemv_(
char *trans, integer *m, integer *n, real *alpha,
587 real *a, integer *lda, real *x, integer *incx, real *beta, real *y,
591 integer a_dim1, a_offset, i__1, i__2;
594 static integer i__, j, ix, iy, jx, jy, kx, ky, info;
596 static integer lenx, leny;
597 extern logical lsame_(
char *,
char *);
598 extern int xerbla_(
char *, integer *);
698 a_offset = 1 + a_dim1;
705 if (! lsame_(trans,
"N") && ! lsame_(trans,
"T") && ! lsame_(trans,
"C")
712 }
else if (*lda < max(1,*m)) {
714 }
else if (*incx == 0) {
716 }
else if (*incy == 0) {
720 xerbla_(
"SGEMV ", &info);
726 if (*m == 0 || *n == 0 || *alpha == 0.f && *beta == 1.f) {
735 if (lsame_(trans,
"N")) {
745 kx = 1 - (lenx - 1) * *incx;
750 ky = 1 - (leny - 1) * *incy;
764 for (i__ = 1; i__ <= i__1; ++i__) {
770 for (i__ = 1; i__ <= i__1; ++i__) {
771 y[i__] = *beta * y[i__];
779 for (i__ = 1; i__ <= i__1; ++i__) {
786 for (i__ = 1; i__ <= i__1; ++i__) {
787 y[iy] = *beta * y[iy];
797 if (lsame_(trans,
"N")) {
804 for (j = 1; j <= i__1; ++j) {
806 temp = *alpha * x[jx];
808 for (i__ = 1; i__ <= i__2; ++i__) {
809 y[i__] += temp * a[i__ + j * a_dim1];
818 for (j = 1; j <= i__1; ++j) {
820 temp = *alpha * x[jx];
823 for (i__ = 1; i__ <= i__2; ++i__) {
824 y[iy] += temp * a[i__ + j * a_dim1];
840 for (j = 1; j <= i__1; ++j) {
843 for (i__ = 1; i__ <= i__2; ++i__) {
844 temp += a[i__ + j * a_dim1] * x[i__];
847 y[jy] += *alpha * temp;
853 for (j = 1; j <= i__1; ++j) {
857 for (i__ = 1; i__ <= i__2; ++i__) {
858 temp += a[i__ + j * a_dim1] * x[ix];
862 y[jy] += *alpha * temp;
875 int sscal_(integer *n, real *sa, real *sx, integer *incx)
881 static integer i__, m, mp1, nincx;
897 if (*n <= 0 || *incx <= 0) {
909 for (i__ = 1; i__2 < 0 ? i__ >= i__1 : i__ <= i__1; i__ += i__2) {
910 sx[i__] = *sa * sx[i__];
928 for (i__ = 1; i__ <= i__2; ++i__) {
929 sx[i__] = *sa * sx[i__];
938 for (i__ = mp1; i__ <= i__2; i__ += 5) {
939 sx[i__] = *sa * sx[i__];
940 sx[i__ + 1] = *sa * sx[i__ + 1];
941 sx[i__ + 2] = *sa * sx[i__ + 2];
942 sx[i__ + 3] = *sa * sx[i__ + 3];
943 sx[i__ + 4] = *sa * sx[i__ + 4];
949 int ssymm_(
char *side,
char *uplo, integer *m, integer *n,
950 real *alpha, real *a, integer *lda, real *b, integer *ldb, real *beta,
951 real *c__, integer *ldc)
954 integer a_dim1, a_offset, b_dim1, b_offset, c_dim1, c_offset, i__1, i__2,
958 static integer i__, j, k, info;
959 static real temp1, temp2;
960 extern logical lsame_(
char *,
char *);
961 static integer nrowa;
962 static logical upper;
963 extern int xerbla_(
char *, integer *);
1096 a_offset = 1 + a_dim1;
1099 b_offset = 1 + b_dim1;
1102 c_offset = 1 + c_dim1;
1106 if (lsame_(side,
"L")) {
1111 upper = lsame_(uplo,
"U");
1116 if (! lsame_(side,
"L") && ! lsame_(side,
"R")) {
1118 }
else if (! upper && ! lsame_(uplo,
"L")) {
1120 }
else if (*m < 0) {
1122 }
else if (*n < 0) {
1124 }
else if (*lda < max(1,nrowa)) {
1126 }
else if (*ldb < max(1,*m)) {
1128 }
else if (*ldc < max(1,*m)) {
1132 xerbla_(
"SSYMM ", &info);
1138 if (*m == 0 || *n == 0 || *alpha == 0.f && *beta == 1.f) {
1144 if (*alpha == 0.f) {
1147 for (j = 1; j <= i__1; ++j) {
1149 for (i__ = 1; i__ <= i__2; ++i__) {
1150 c__[i__ + j * c_dim1] = 0.f;
1157 for (j = 1; j <= i__1; ++j) {
1159 for (i__ = 1; i__ <= i__2; ++i__) {
1160 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1];
1171 if (lsame_(side,
"L")) {
1177 for (j = 1; j <= i__1; ++j) {
1179 for (i__ = 1; i__ <= i__2; ++i__) {
1180 temp1 = *alpha * b[i__ + j * b_dim1];
1183 for (k = 1; k <= i__3; ++k) {
1184 c__[k + j * c_dim1] += temp1 * a[k + i__ * a_dim1];
1185 temp2 += b[k + j * b_dim1] * a[k + i__ * a_dim1];
1189 c__[i__ + j * c_dim1] = temp1 * a[i__ + i__ * a_dim1]
1192 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1]
1193 + temp1 * a[i__ + i__ * a_dim1] + *alpha *
1202 for (j = 1; j <= i__1; ++j) {
1203 for (i__ = *m; i__ >= 1; --i__) {
1204 temp1 = *alpha * b[i__ + j * b_dim1];
1207 for (k = i__ + 1; k <= i__2; ++k) {
1208 c__[k + j * c_dim1] += temp1 * a[k + i__ * a_dim1];
1209 temp2 += b[k + j * b_dim1] * a[k + i__ * a_dim1];
1213 c__[i__ + j * c_dim1] = temp1 * a[i__ + i__ * a_dim1]
1216 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1]
1217 + temp1 * a[i__ + i__ * a_dim1] + *alpha *
1230 for (j = 1; j <= i__1; ++j) {
1231 temp1 = *alpha * a[j + j * a_dim1];
1234 for (i__ = 1; i__ <= i__2; ++i__) {
1235 c__[i__ + j * c_dim1] = temp1 * b[i__ + j * b_dim1];
1240 for (i__ = 1; i__ <= i__2; ++i__) {
1241 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1] +
1242 temp1 * b[i__ + j * b_dim1];
1247 for (k = 1; k <= i__2; ++k) {
1249 temp1 = *alpha * a[k + j * a_dim1];
1251 temp1 = *alpha * a[j + k * a_dim1];
1254 for (i__ = 1; i__ <= i__3; ++i__) {
1255 c__[i__ + j * c_dim1] += temp1 * b[i__ + k * b_dim1];
1261 for (k = j + 1; k <= i__2; ++k) {
1263 temp1 = *alpha * a[j + k * a_dim1];
1265 temp1 = *alpha * a[k + j * a_dim1];
1268 for (i__ = 1; i__ <= i__3; ++i__) {
1269 c__[i__ + j * c_dim1] += temp1 * b[i__ + k * b_dim1];
1284 int ssyrk_(
char *uplo,
char *trans, integer *n, integer *k,
1285 real *alpha, real *a, integer *lda, real *beta, real *c__, integer *
1289 integer a_dim1, a_offset, c_dim1, c_offset, i__1, i__2, i__3;
1292 static integer i__, j, l, info;
1294 extern logical lsame_(
char *,
char *);
1295 static integer nrowa;
1296 static logical upper;
1297 extern int xerbla_(
char *, integer *);
1414 a_offset = 1 + a_dim1;
1417 c_offset = 1 + c_dim1;
1421 if (lsame_(trans,
"N")) {
1426 upper = lsame_(uplo,
"U");
1429 if (! upper && ! lsame_(uplo,
"L")) {
1431 }
else if (! lsame_(trans,
"N") && ! lsame_(trans,
1432 "T") && ! lsame_(trans,
"C")) {
1434 }
else if (*n < 0) {
1436 }
else if (*k < 0) {
1438 }
else if (*lda < max(1,nrowa)) {
1440 }
else if (*ldc < max(1,*n)) {
1444 xerbla_(
"SSYRK ", &info);
1450 if (*n == 0 || (*alpha == 0.f || *k == 0) && *beta == 1.f) {
1456 if (*alpha == 0.f) {
1460 for (j = 1; j <= i__1; ++j) {
1462 for (i__ = 1; i__ <= i__2; ++i__) {
1463 c__[i__ + j * c_dim1] = 0.f;
1470 for (j = 1; j <= i__1; ++j) {
1472 for (i__ = 1; i__ <= i__2; ++i__) {
1473 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1];
1482 for (j = 1; j <= i__1; ++j) {
1484 for (i__ = j; i__ <= i__2; ++i__) {
1485 c__[i__ + j * c_dim1] = 0.f;
1492 for (j = 1; j <= i__1; ++j) {
1494 for (i__ = j; i__ <= i__2; ++i__) {
1495 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1];
1507 if (lsame_(trans,
"N")) {
1513 for (j = 1; j <= i__1; ++j) {
1516 for (i__ = 1; i__ <= i__2; ++i__) {
1517 c__[i__ + j * c_dim1] = 0.f;
1520 }
else if (*beta != 1.f) {
1522 for (i__ = 1; i__ <= i__2; ++i__) {
1523 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1];
1528 for (l = 1; l <= i__2; ++l) {
1529 if (a[j + l * a_dim1] != 0.f) {
1530 temp = *alpha * a[j + l * a_dim1];
1532 for (i__ = 1; i__ <= i__3; ++i__) {
1533 c__[i__ + j * c_dim1] += temp * a[i__ + l *
1544 for (j = 1; j <= i__1; ++j) {
1547 for (i__ = j; i__ <= i__2; ++i__) {
1548 c__[i__ + j * c_dim1] = 0.f;
1551 }
else if (*beta != 1.f) {
1553 for (i__ = j; i__ <= i__2; ++i__) {
1554 c__[i__ + j * c_dim1] = *beta * c__[i__ + j * c_dim1];
1559 for (l = 1; l <= i__2; ++l) {
1560 if (a[j + l * a_dim1] != 0.f) {
1561 temp = *alpha * a[j + l * a_dim1];
1563 for (i__ = j; i__ <= i__3; ++i__) {
1564 c__[i__ + j * c_dim1] += temp * a[i__ + l *
1580 for (j = 1; j <= i__1; ++j) {
1582 for (i__ = 1; i__ <= i__2; ++i__) {
1585 for (l = 1; l <= i__3; ++l) {
1586 temp += a[l + i__ * a_dim1] * a[l + j * a_dim1];
1590 c__[i__ + j * c_dim1] = *alpha * temp;
1592 c__[i__ + j * c_dim1] = *alpha * temp + *beta * c__[
1601 for (j = 1; j <= i__1; ++j) {
1603 for (i__ = j; i__ <= i__2; ++i__) {
1606 for (l = 1; l <= i__3; ++l) {
1607 temp += a[l + i__ * a_dim1] * a[l + j * a_dim1];
1611 c__[i__ + j * c_dim1] = *alpha * temp;
1613 c__[i__ + j * c_dim1] = *alpha * temp + *beta * c__[
1629 int strsm_(
char *side,
char *uplo,
char *transa,
char *diag,
1630 integer *m, integer *n, real *alpha, real *a, integer *lda, real *b,
1634 integer a_dim1, a_offset, b_dim1, b_offset, i__1, i__2, i__3;
1637 static integer i__, j, k, info;
1639 static logical lside;
1640 extern logical lsame_(
char *,
char *);
1641 static integer nrowa;
1642 static logical upper;
1643 extern int xerbla_(
char *, integer *);
1644 static logical nounit;
1772 a_offset = 1 + a_dim1;
1775 b_offset = 1 + b_dim1;
1779 lside = lsame_(side,
"L");
1785 nounit = lsame_(diag,
"N");
1786 upper = lsame_(uplo,
"U");
1789 if (! lside && ! lsame_(side,
"R")) {
1791 }
else if (! upper && ! lsame_(uplo,
"L")) {
1793 }
else if (! lsame_(transa,
"N") && ! lsame_(transa,
1794 "T") && ! lsame_(transa,
"C")) {
1796 }
else if (! lsame_(diag,
"U") && ! lsame_(diag,
1799 }
else if (*m < 0) {
1801 }
else if (*n < 0) {
1803 }
else if (*lda < max(1,nrowa)) {
1805 }
else if (*ldb < max(1,*m)) {
1809 xerbla_(
"STRSM ", &info);
1821 if (*alpha == 0.f) {
1823 for (j = 1; j <= i__1; ++j) {
1825 for (i__ = 1; i__ <= i__2; ++i__) {
1826 b[i__ + j * b_dim1] = 0.f;
1837 if (lsame_(transa,
"N")) {
1843 for (j = 1; j <= i__1; ++j) {
1844 if (*alpha != 1.f) {
1846 for (i__ = 1; i__ <= i__2; ++i__) {
1847 b[i__ + j * b_dim1] = *alpha * b[i__ + j * b_dim1]
1852 for (k = *m; k >= 1; --k) {
1853 if (b[k + j * b_dim1] != 0.f) {
1855 b[k + j * b_dim1] /= a[k + k * a_dim1];
1858 for (i__ = 1; i__ <= i__2; ++i__) {
1859 b[i__ + j * b_dim1] -= b[k + j * b_dim1] * a[
1870 for (j = 1; j <= i__1; ++j) {
1871 if (*alpha != 1.f) {
1873 for (i__ = 1; i__ <= i__2; ++i__) {
1874 b[i__ + j * b_dim1] = *alpha * b[i__ + j * b_dim1]
1880 for (k = 1; k <= i__2; ++k) {
1881 if (b[k + j * b_dim1] != 0.f) {
1883 b[k + j * b_dim1] /= a[k + k * a_dim1];
1886 for (i__ = k + 1; i__ <= i__3; ++i__) {
1887 b[i__ + j * b_dim1] -= b[k + j * b_dim1] * a[
1903 for (j = 1; j <= i__1; ++j) {
1905 for (i__ = 1; i__ <= i__2; ++i__) {
1906 temp = *alpha * b[i__ + j * b_dim1];
1908 for (k = 1; k <= i__3; ++k) {
1909 temp -= a[k + i__ * a_dim1] * b[k + j * b_dim1];
1913 temp /= a[i__ + i__ * a_dim1];
1915 b[i__ + j * b_dim1] = temp;
1922 for (j = 1; j <= i__1; ++j) {
1923 for (i__ = *m; i__ >= 1; --i__) {
1924 temp = *alpha * b[i__ + j * b_dim1];
1926 for (k = i__ + 1; k <= i__2; ++k) {
1927 temp -= a[k + i__ * a_dim1] * b[k + j * b_dim1];
1931 temp /= a[i__ + i__ * a_dim1];
1933 b[i__ + j * b_dim1] = temp;
1941 if (lsame_(transa,
"N")) {
1947 for (j = 1; j <= i__1; ++j) {
1948 if (*alpha != 1.f) {
1950 for (i__ = 1; i__ <= i__2; ++i__) {
1951 b[i__ + j * b_dim1] = *alpha * b[i__ + j * b_dim1]
1957 for (k = 1; k <= i__2; ++k) {
1958 if (a[k + j * a_dim1] != 0.f) {
1960 for (i__ = 1; i__ <= i__3; ++i__) {
1961 b[i__ + j * b_dim1] -= a[k + j * a_dim1] * b[
1969 temp = 1.f / a[j + j * a_dim1];
1971 for (i__ = 1; i__ <= i__2; ++i__) {
1972 b[i__ + j * b_dim1] = temp * b[i__ + j * b_dim1];
1979 for (j = *n; j >= 1; --j) {
1980 if (*alpha != 1.f) {
1982 for (i__ = 1; i__ <= i__1; ++i__) {
1983 b[i__ + j * b_dim1] = *alpha * b[i__ + j * b_dim1]
1989 for (k = j + 1; k <= i__1; ++k) {
1990 if (a[k + j * a_dim1] != 0.f) {
1992 for (i__ = 1; i__ <= i__2; ++i__) {
1993 b[i__ + j * b_dim1] -= a[k + j * a_dim1] * b[
2001 temp = 1.f / a[j + j * a_dim1];
2003 for (i__ = 1; i__ <= i__1; ++i__) {
2004 b[i__ + j * b_dim1] = temp * b[i__ + j * b_dim1];
2016 for (k = *n; k >= 1; --k) {
2018 temp = 1.f / a[k + k * a_dim1];
2020 for (i__ = 1; i__ <= i__1; ++i__) {
2021 b[i__ + k * b_dim1] = temp * b[i__ + k * b_dim1];
2026 for (j = 1; j <= i__1; ++j) {
2027 if (a[j + k * a_dim1] != 0.f) {
2028 temp = a[j + k * a_dim1];
2030 for (i__ = 1; i__ <= i__2; ++i__) {
2031 b[i__ + j * b_dim1] -= temp * b[i__ + k *
2038 if (*alpha != 1.f) {
2040 for (i__ = 1; i__ <= i__1; ++i__) {
2041 b[i__ + k * b_dim1] = *alpha * b[i__ + k * b_dim1]
2050 for (k = 1; k <= i__1; ++k) {
2052 temp = 1.f / a[k + k * a_dim1];
2054 for (i__ = 1; i__ <= i__2; ++i__) {
2055 b[i__ + k * b_dim1] = temp * b[i__ + k * b_dim1];
2060 for (j = k + 1; j <= i__2; ++j) {
2061 if (a[j + k * a_dim1] != 0.f) {
2062 temp = a[j + k * a_dim1];
2064 for (i__ = 1; i__ <= i__3; ++i__) {
2065 b[i__ + j * b_dim1] -= temp * b[i__ + k *
2072 if (*alpha != 1.f) {
2074 for (i__ = 1; i__ <= i__2; ++i__) {
2075 b[i__ + k * b_dim1] = *alpha * b[i__ + k * b_dim1]
2092 int xerbla_(
char *srname, integer *info)
2095 static char fmt_9999[] =
"(\002 ** On entry to \002,a6,\002 parameter nu"
2096 "mber \002,i2,\002 had \002,\002an illegal value\002)";
2099 integer s_wsfe(
cilist *), do_fio(integer *,
char *, ftnlen), e_wsfe(
void);
2100 int s_stop(
char *, ftnlen);
2103 static cilist io___60 = { 0, 6, 0, fmt_9999, 0 };
2136 do_fio(&c__1, srname, (ftnlen)6);
2137 do_fio(&c__1, (
char *)&(*info), (ftnlen)
sizeof(integer));
2140 s_stop(
"", (ftnlen)0);