| 1 | #include "f2c.h"
|
|---|
| 2 | #include "fio.h"
|
|---|
| 3 | #include "fmt.h"
|
|---|
| 4 | #include "lio.h"
|
|---|
| 5 |
|
|---|
| 6 | ftnint L_len;
|
|---|
| 7 | int f__Aquote;
|
|---|
| 8 |
|
|---|
| 9 | static void
|
|---|
| 10 | donewrec (void)
|
|---|
| 11 | {
|
|---|
| 12 | if (f__recpos)
|
|---|
| 13 | (*f__donewrec) ();
|
|---|
| 14 | }
|
|---|
| 15 |
|
|---|
| 16 | static void
|
|---|
| 17 | lwrt_I (longint n)
|
|---|
| 18 | {
|
|---|
| 19 | char *p;
|
|---|
| 20 | int ndigit, sign;
|
|---|
| 21 |
|
|---|
| 22 | p = f__icvt (n, &ndigit, &sign, 10);
|
|---|
| 23 | if (f__recpos + ndigit >= L_len)
|
|---|
| 24 | donewrec ();
|
|---|
| 25 | PUT (' ');
|
|---|
| 26 | if (sign)
|
|---|
| 27 | PUT ('-');
|
|---|
| 28 | while (*p)
|
|---|
| 29 | PUT (*p++);
|
|---|
| 30 | }
|
|---|
| 31 | static void
|
|---|
| 32 | lwrt_L (ftnint n, ftnlen len)
|
|---|
| 33 | {
|
|---|
| 34 | if (f__recpos + LLOGW >= L_len)
|
|---|
| 35 | donewrec ();
|
|---|
| 36 | wrt_L ((Uint *) & n, LLOGW, len);
|
|---|
| 37 | }
|
|---|
| 38 | static void
|
|---|
| 39 | lwrt_A (char *p, ftnlen len)
|
|---|
| 40 | {
|
|---|
| 41 | int a;
|
|---|
| 42 | char *p1, *pe;
|
|---|
| 43 |
|
|---|
| 44 | a = 0;
|
|---|
| 45 | pe = p + len;
|
|---|
| 46 | if (f__Aquote)
|
|---|
| 47 | {
|
|---|
| 48 | a = 3;
|
|---|
| 49 | if (len > 1 && p[len - 1] == ' ')
|
|---|
| 50 | {
|
|---|
| 51 | while (--len > 1 && p[len - 1] == ' ');
|
|---|
| 52 | pe = p + len;
|
|---|
| 53 | }
|
|---|
| 54 | p1 = p;
|
|---|
| 55 | while (p1 < pe)
|
|---|
| 56 | if (*p1++ == '\'')
|
|---|
| 57 | a++;
|
|---|
| 58 | }
|
|---|
| 59 | if (f__recpos + len + a >= L_len)
|
|---|
| 60 | donewrec ();
|
|---|
| 61 | if (a
|
|---|
| 62 | #ifndef OMIT_BLANK_CC
|
|---|
| 63 | || !f__recpos
|
|---|
| 64 | #endif
|
|---|
| 65 | )
|
|---|
| 66 | PUT (' ');
|
|---|
| 67 | if (a)
|
|---|
| 68 | {
|
|---|
| 69 | PUT ('\'');
|
|---|
| 70 | while (p < pe)
|
|---|
| 71 | {
|
|---|
| 72 | if (*p == '\'')
|
|---|
| 73 | PUT ('\'');
|
|---|
| 74 | PUT (*p++);
|
|---|
| 75 | }
|
|---|
| 76 | PUT ('\'');
|
|---|
| 77 | }
|
|---|
| 78 | else
|
|---|
| 79 | while (p < pe)
|
|---|
| 80 | PUT (*p++);
|
|---|
| 81 | }
|
|---|
| 82 |
|
|---|
| 83 | static int
|
|---|
| 84 | l_g (char *buf, double n)
|
|---|
| 85 | {
|
|---|
| 86 | #ifdef Old_list_output
|
|---|
| 87 | doublereal absn;
|
|---|
| 88 | char *fmt;
|
|---|
| 89 |
|
|---|
| 90 | absn = n;
|
|---|
| 91 | if (absn < 0)
|
|---|
| 92 | absn = -absn;
|
|---|
| 93 | fmt = LLOW <= absn && absn < LHIGH ? LFFMT : LEFMT;
|
|---|
| 94 | #ifdef USE_STRLEN
|
|---|
| 95 | sprintf (buf, fmt, n);
|
|---|
| 96 | return strlen (buf);
|
|---|
| 97 | #else
|
|---|
| 98 | return sprintf (buf, fmt, n);
|
|---|
| 99 | #endif
|
|---|
| 100 |
|
|---|
| 101 | #else
|
|---|
| 102 | register char *b, c, c1;
|
|---|
| 103 |
|
|---|
| 104 | b = buf;
|
|---|
| 105 | *b++ = ' ';
|
|---|
| 106 | if (n < 0)
|
|---|
| 107 | {
|
|---|
| 108 | *b++ = '-';
|
|---|
| 109 | n = -n;
|
|---|
| 110 | }
|
|---|
| 111 | else
|
|---|
| 112 | *b++ = ' ';
|
|---|
| 113 | if (n == 0)
|
|---|
| 114 | {
|
|---|
| 115 | *b++ = '0';
|
|---|
| 116 | *b++ = '.';
|
|---|
| 117 | *b = 0;
|
|---|
| 118 | goto f__ret;
|
|---|
| 119 | }
|
|---|
| 120 | sprintf (b, LGFMT, n);
|
|---|
| 121 | switch (*b)
|
|---|
| 122 | {
|
|---|
| 123 | #ifndef WANT_LEAD_0
|
|---|
| 124 | case '0':
|
|---|
| 125 | while (b[0] = b[1])
|
|---|
| 126 | b++;
|
|---|
| 127 | break;
|
|---|
| 128 | #endif
|
|---|
| 129 | case 'i':
|
|---|
| 130 | case 'I':
|
|---|
| 131 | /* Infinity */
|
|---|
| 132 | case 'n':
|
|---|
| 133 | case 'N':
|
|---|
| 134 | /* NaN */
|
|---|
| 135 | while (*++b);
|
|---|
| 136 | break;
|
|---|
| 137 |
|
|---|
| 138 | default:
|
|---|
| 139 | /* Fortran 77 insists on having a decimal point... */
|
|---|
| 140 | for (;; b++)
|
|---|
| 141 | switch (*b)
|
|---|
| 142 | {
|
|---|
| 143 | case 0:
|
|---|
| 144 | *b++ = '.';
|
|---|
| 145 | *b = 0;
|
|---|
| 146 | goto f__ret;
|
|---|
| 147 | case '.':
|
|---|
| 148 | while (*++b);
|
|---|
| 149 | goto f__ret;
|
|---|
| 150 | case 'E':
|
|---|
| 151 | for (c1 = '.', c = 'E'; (*b = c1); c1 = c, c = *++b);
|
|---|
| 152 | goto f__ret;
|
|---|
| 153 | }
|
|---|
| 154 | }
|
|---|
| 155 | f__ret:
|
|---|
| 156 | return b - buf;
|
|---|
| 157 | #endif
|
|---|
| 158 | }
|
|---|
| 159 |
|
|---|
| 160 | static void
|
|---|
| 161 | l_put (register char *s)
|
|---|
| 162 | {
|
|---|
| 163 | register void (*pn) (int) = f__putn;
|
|---|
| 164 | register int c;
|
|---|
| 165 |
|
|---|
| 166 | while ((c = *s++))
|
|---|
| 167 | (*pn) (c);
|
|---|
| 168 | }
|
|---|
| 169 |
|
|---|
| 170 | static void
|
|---|
| 171 | lwrt_F (double n)
|
|---|
| 172 | {
|
|---|
| 173 | char buf[LEFBL];
|
|---|
| 174 |
|
|---|
| 175 | if (f__recpos + l_g (buf, n) >= L_len)
|
|---|
| 176 | donewrec ();
|
|---|
| 177 | l_put (buf);
|
|---|
| 178 | }
|
|---|
| 179 | static void
|
|---|
| 180 | lwrt_C (double a, double b)
|
|---|
| 181 | {
|
|---|
| 182 | char *ba, *bb, bufa[LEFBL], bufb[LEFBL];
|
|---|
| 183 | int al, bl;
|
|---|
| 184 |
|
|---|
| 185 | al = l_g (bufa, a);
|
|---|
| 186 | for (ba = bufa; *ba == ' '; ba++)
|
|---|
| 187 | --al;
|
|---|
| 188 | bl = l_g (bufb, b) + 1; /* intentionally high by 1 */
|
|---|
| 189 | for (bb = bufb; *bb == ' '; bb++)
|
|---|
| 190 | --bl;
|
|---|
| 191 | if (f__recpos + al + bl + 3 >= L_len)
|
|---|
| 192 | donewrec ();
|
|---|
| 193 | #ifdef OMIT_BLANK_CC
|
|---|
| 194 | else
|
|---|
| 195 | #endif
|
|---|
| 196 | PUT (' ');
|
|---|
| 197 | PUT ('(');
|
|---|
| 198 | l_put (ba);
|
|---|
| 199 | PUT (',');
|
|---|
| 200 | if (f__recpos + bl >= L_len)
|
|---|
| 201 | {
|
|---|
| 202 | (*f__donewrec) ();
|
|---|
| 203 | #ifndef OMIT_BLANK_CC
|
|---|
| 204 | PUT (' ');
|
|---|
| 205 | #endif
|
|---|
| 206 | }
|
|---|
| 207 | l_put (bb);
|
|---|
| 208 | PUT (')');
|
|---|
| 209 | }
|
|---|
| 210 |
|
|---|
| 211 | int
|
|---|
| 212 | l_write (ftnint * number, char *ptr, ftnlen len, ftnint type)
|
|---|
| 213 | {
|
|---|
| 214 | #define Ptr ((flex *)ptr)
|
|---|
| 215 | int i;
|
|---|
| 216 | longint x;
|
|---|
| 217 | double y, z;
|
|---|
| 218 | real *xx;
|
|---|
| 219 | doublereal *yy;
|
|---|
| 220 | for (i = 0; i < *number; i++)
|
|---|
| 221 | {
|
|---|
| 222 | switch ((int) type)
|
|---|
| 223 | {
|
|---|
| 224 | default:
|
|---|
| 225 | f__fatal (204, "unknown type in lio");
|
|---|
| 226 | case TYINT1:
|
|---|
| 227 | x = Ptr->flchar;
|
|---|
| 228 | goto xint;
|
|---|
| 229 | case TYSHORT:
|
|---|
| 230 | x = Ptr->flshort;
|
|---|
| 231 | goto xint;
|
|---|
| 232 | #ifdef Allow_TYQUAD
|
|---|
| 233 | case TYQUAD:
|
|---|
| 234 | x = Ptr->fllongint;
|
|---|
| 235 | goto xint;
|
|---|
| 236 | #endif
|
|---|
| 237 | case TYLONG:
|
|---|
| 238 | x = Ptr->flint;
|
|---|
| 239 | xint:lwrt_I (x);
|
|---|
| 240 | break;
|
|---|
| 241 | case TYREAL:
|
|---|
| 242 | y = Ptr->flreal;
|
|---|
| 243 | goto xfloat;
|
|---|
| 244 | case TYDREAL:
|
|---|
| 245 | y = Ptr->fldouble;
|
|---|
| 246 | xfloat:lwrt_F (y);
|
|---|
| 247 | break;
|
|---|
| 248 | case TYCOMPLEX:
|
|---|
| 249 | xx = &Ptr->flreal;
|
|---|
| 250 | y = *xx++;
|
|---|
| 251 | z = *xx;
|
|---|
| 252 | goto xcomplex;
|
|---|
| 253 | case TYDCOMPLEX:
|
|---|
| 254 | yy = &Ptr->fldouble;
|
|---|
| 255 | y = *yy++;
|
|---|
| 256 | z = *yy;
|
|---|
| 257 | xcomplex:
|
|---|
| 258 | lwrt_C (y, z);
|
|---|
| 259 | break;
|
|---|
| 260 | case TYLOGICAL1:
|
|---|
| 261 | x = Ptr->flchar;
|
|---|
| 262 | goto xlog;
|
|---|
| 263 | case TYLOGICAL2:
|
|---|
| 264 | x = Ptr->flshort;
|
|---|
| 265 | goto xlog;
|
|---|
| 266 | case TYLOGICAL:
|
|---|
| 267 | x = Ptr->flint;
|
|---|
| 268 | xlog:lwrt_L (Ptr->flint, len);
|
|---|
| 269 | break;
|
|---|
| 270 | case TYCHAR:
|
|---|
| 271 | lwrt_A (ptr, len);
|
|---|
| 272 | break;
|
|---|
| 273 | }
|
|---|
| 274 | ptr += len;
|
|---|
| 275 | }
|
|---|
| 276 | return (0);
|
|---|
| 277 | }
|
|---|