vif.h
1 //! @file vif.h
2 //! @author J. Marcel van der Veer
3 //
4 //! @section Copyright
5 //
6 // This file is part of VIF - vintage FORTRAN compiler.
7 // Copyright 2020-2025 J. Marcel van der Veer <algol68g@xs4all.nl>.
8 //
9 //! @section License
10 //
11 // This program is free software; you can redistribute it and/or modify it
12 // under the terms of the GNU General Public License as published by the
13 // Free Software Foundation; either version 3 of the License, or
14 // (at your option) any later version.
15 //
16 // This program is distributed in the hope that it will be useful, but
17 // WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
18 // or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for
19 // more details. You should have received a copy of the GNU General Public
20 // License along with this program. If not, see <http://www.gnu.org/licenses/>.
21
22 //! @section Synopsis
23 //!
24 //! VIF include file.
25
26 // This is a monolithic include file to avoid cluttering the installation directory.
27
28 #if ! defined (__VIF_H__)
29 #define __VIF_H__
30
31 #define _GNU_SOURCE
32
33 #include <complex.h>
34 #include <ctype.h>
35 #include <endian.h>
36 #include <errno.h>
37 #include <float.h>
38 #include <inttypes.h>
39 #include <libgen.h>
40 #include <limits.h>
41 #include <math.h>
42 #include <quadmath.h>
43 #include <signal.h>
44 #include <stdarg.h>
45 #include <stddef.h>
46 #include <stdint.h>
47 #include <stdio.h>
48 #include <stdlib.h>
49 #include <string.h>
50 #include <sys/stat.h>
51 #include <sys/time.h>
52 #include <sys/utsname.h>
53 #include <time.h>
54 #include <unistd.h>
55
56 #define prototype
57
58 #define PACKAGE "vif"
59
60 #define NO_IDENT ((IDENT *) NULL)
61 #define NO_INTRINS ((INTRINS *) NULL)
62 #define NO_LABEL ((LBL *) NULL)
63 #define NO_MODE ((MODE *) NULL)
64 #define NO_TEXT ((char *) NULL)
65 #define NO_FUN NULL
66 #define NO_FTNFILE ((FTNFILE *) NULL)
67 #define NO_FTN_LINE ((FTN_LINE *) NULL)
68 #define NO_FILE ((FILE *) NULL)
69 #define NO_REF_INTEGER ((int_4 *) NULL)
70 #define NO_EXPR ((EXPR *) NULL)
71 #define NO_REF_TEXT ((char **) NULL)
72
73 #define STDF_IN 5
74 #define STDF_OUT 6
75 #define STDF_PUN 6
76 #define STDF_ERR 7
77
78 typedef char *FORMAT;
79
80 // Flags for gcc when compiling generated code from FORTRAN.
81
82 #if defined (BOOTSTRAP)
83 #define CFLAGS "-pipe -ggdb -I. -Isrc/include -L. -I../include -Wno-implicit-function-declaration -Wno-incompatible-pointer-types -fdiagnostics-plain-output -fdiagnostics-show-location=once"
84 #define LD_FLAGS "-rdynamic -Lsrc/lib -L../lib -lvif -lquadmath -lm"
85 #else
86 #define CFLAGS "-pipe -ggdb -I. -Isrc/include -L. -I../include -Wno-implicit-function-declaration -Wno-incompatible-pointer-types -fdiagnostics-plain-output -fdiagnostics-show-location=once -z execstack"
87 #define LD_FLAGS "-rdynamic -L./src/lib -lvif -lquadmath -lm"
88 #endif
89 #define OFLAGS "-O2"
90
91 enum {NOTYPE, ETYPE, INTEGER, LOGICAL, REAL, COMPLEX, CHARACTER};
92 enum {UNFORMATTED = 1, STDFORMAT, FORMATTED };
93 enum {SOURCE = 0, TEMP, MACRO };
94 enum {HEADER = 0, BANNER, LIST, MESSAGES, JCL, SYMBOLS, CONSTANTS, STRINGS, PARAMETERS, TYPEDEF, FREQ, COMMON, PROTOTYPE, STAT, TITLE, PRE, DECL, REFDECL, EQUIV, FMT, DATA, NESTED, BODY, POST, MAXPHA };
95 enum {EXPR_OTHER = 0, EXPR_VAR, EXPR_CONST, EXPR_SLICE, EXPR_SUBSTR };
96 enum {LEXEME = 0, DECLAR, WORD, LABEL, INT_NUMBER, FLT_NUMBER, TEXT };
97 enum {STATIC = 0, AUTOMATIC };
98
99 // Shorthands to improve code legibility.
100
101 #define int_2 int16_t
102 #define int_4 int32_t
103 #define int_8 int64_t
104
105 #define unt_2 uint16_t
106 #define unt_4 uint32_t
107 #define unt_8 uint64_t
108 #define unt_16 unsigned __int128
109
110 #define logical_4 unsigned
111
112 #define real_4 float
113 #define real_8 double
114 #define real_16 __float128
115
116 #define complex_8 float complex
117 #define complex_16 double complex
118 #define complex_32 __complex128
119
120 // MAX_STRLEN = 2 ** MAX_STRLENS - 1
121 #define MAX_STRLEN 4095
122 #define MAX_STRLENS 12
123 extern int_4 strlens[MAX_STRLENS];
124
125 #define ERR (-32)
126 #define SKIP //
127 #define TRUE 1
128 #define FALSE 0
129 #define MAX_LRECL 32760
130 #define MAX_NEST 9
131 #define INCREMENT 500
132 #define INDENT 2
133 #define LINES_PER_PAGE 50
134 #define LINE_WIDTH 122
135 #define MAX_ARGS 16
136 #define MAX_COMMONS 100
137 #define MAX_DIMS 6
138 #define MAX_IDENTS 1000
139 #define MAX_LABELS 1000
140 #define MAX_MODULES 5000
141 #define MAX_PRIO 9
142 #define MAX_ERROR 5
143 #define MAX_WARNS MAX_ERROR
144 #define RETURN "goto _l0;\n"
145 #define M_LOG10_2 0.30102999566398119521373889472449q
146
147 #define NOT_LOCAL(idf) ((idf)->common != LOCAL)
148 #define LOCAL 0
149 #define EXTERN 1
150
151 #define WITHIN (curret != END_OF_LINE && curret != END_OF_MODULE)
152
153 #define END_OF_MODULE (-3)
154 #define END_OF_LINE (-2)
155 #define EXPECT_NONE (NULL)
156 #define START_OF_LINE (-1)
157 #define ARG TRUE
158 #define NOARG FALSE
159 #define CONST TRUE
160 #define NOCONST FALSE
161 #define FUN TRUE
162 #define NOFUN FALSE
163 #define UNIQ TRUE
164 #define NOUNIQ FALSE
165 #define NOPATCH 0
166 #define NOTYPE 0
167 #define NOLEN 0
168 #define FORMAL TRUE
169 #define ACTUAL FALSE
170 #define CAST TRUE
171 #define NOCAST FALSE
172 #define PROTEST TRUE
173 #define QUIET FALSE
174
175 #define EXPECT_LABEL "label"
176 #define INTERNAL_CONSISTENCY "internal consistency"
177
178 #define ABS(n) ((n) > 0 ? (n) : -(n))
179 #define MAXIMISE(u, v) ((u) = _max (u, v))
180 #define MINIMISE(u, v) ((u) = _min (u, v))
181
182 #define IS_JCL(c) (c == '/' || CUR_LIN.jcl)
183 #define IS_COMMENT(c) (strchr ("cd*!#", tolower (c)) != NULL || IS_JCL (c))
184 #define IS_VAR(n) (isalpha(n[0]) || n[0] == '_' || n[0] == '$')
185 #define UNSCAN {curlin = prelin; curcol = precol; strcpy (curlex, prelex); curret = preret;}
186 #define SAVE_POS int_4 _l_ = curlin, _c_ = curcol;
187 #define SAVE_PRE int_4 _l_ = prelin, _c_ = precol;
188 #define RESTORE_POS {curlin = _l_; curcol = _c_;}
189 #define _EXPCHAR(n) (tolower (n) != '\0' && strchr ("edqx", tolower (n)) != NULL)
190 #define _IDFCHAR(n) (tolower (n) != '\0' && strchr ("abcdefghijklmnopqrstuvwxyz0123456789_ $", tolower (n)) != NULL)
191 #define CUR_LIN (source[curlin])
192 #define POS(n) (tolower (CUR_LIN.text[n]))
193 #define EXPONENT(n) _EXPCHAR (POS (n))
194 #define CUR_COL (POS (curcol))
195
196 #define MSG(sev, num, text, info)\
197 if (prelin < 1) {\
198 message (&source[curlin], curcol, (sev), (num), (text), (info));\
199 } else {\
200 message (&source[prelin], precol, (sev), (num), (text), (info));\
201 }
202
203 extern int_4 _srecordf (char *, const char *, ...);
204 extern void _vif_backtr(int_4);
205
206 #define ERROR(num, s, t) {MSG ("error", (num), (s), (t));}
207 #define MODE_ERROR(num, s, t) {\
208 RECORD _txt_;\
209 _srecordf(_txt_, "%s to %s", (s), (t));\
210 ERROR (num, "cannot convert", _txt_);\
211 }
212 #define PRECISION_LOSS(num, s, t) {\
213 RECORD _txt_;\
214 _srecordf(_txt_, "%s to %s", (s), (t));\
215 WARNING (num, "possible precision loss", _txt_);\
216 }
217 #define SYNTAX(num, s) {ERROR ((num), "syntax", (s));}
218 #define EXPECT(num, s) {ERROR ((num), "expected", (s));}
219 #define ADJACENT(num, s) {ERROR ((num), "adjacent", (s));}
220 #define WARNING(num, s, t) {MSG ("warning", (num), (s), (t));}
221 #define ECHO(num, s, t) {MSG ("info", (num), (s), (t));}
222 #define SCANER(num, s, t) {MSG ("fatal", (num), (s), (t)); exit (EXIT_FAILURE);}
223 #define FATAL(num, s, t) {message (NULL, ERR, "fatal", (num), (s), (t)); exit (EXIT_FAILURE);}
224 #define OVERFLOW(num, t) {fprintf (stderr, "\n** exception ** f%d: fatal: overflow, (%s)\n", (num), (t)); _vif_backtr (SIGTERM);}
225
226 #define TOKEN(s) EQUAL (curlex, (s))
227 #define CHECKPOINT(num, s) {if (!EQUAL (curlex, (s))) {EXPECT ((num), (s));}}
228 #define CHECKDIGIT(num, s) {if (!isdigit (s)) {EXPECT ((num), "digit");}}
229 #define IS_NUMBER(u) ((u) == INT_NUMBER || (u) == FLT_NUMBER)
230
231 #define FMT_TEXT ((char *) 1)
232 #define FMT_INT ((char *) 2)
233 #define FMT_REAL ((char *) 3)
234 #define FMT_LOGICAL ((char *) 4)
235 #define FMT_CHAR ((char *) 5)
236 #define FMT_TERM ((char *) 6)
237 extern int_4 __scale__;
238
239 // REAL*32
240
241 #define FLT128_LEN 7 // Do NOT change this!
242 #define FLT256_LEN 15 // Do NOT change this!
243 #define FLT256_GUARD 2 // Guard digits.
244 #define FLT256_DIG (72 - FLT256_GUARD) // 240 * log10 (2) minus guard digits.
245
246 typedef unt_2 REAL16[8];
247 typedef unt_2 REAL32[FLT256_LEN + 1]; // One for exponent.
248
249 struct __real_32__
250 {
251 REAL32 value;
252 };
253 typedef struct __real_32__ real_32;
254
255 struct __complex_64__
256 {
257 struct __real_32__ re, im;
258 };
259 typedef struct __complex_64__ complex_64;
260
261 #define xsub(a, b) xadd (a, b, 1)
262 #define xsum(a, b) xadd (a, b, 0)
263
264 // In Fortran, only EQUIVALENCE can alias names.
265
266 #define _p_ * restrict
267
268 //
269
270 #define RECLN 1280
271 typedef char RECORD[RECLN + 1];
272 #define RECCLR(z) memset ((z), '\0', RECLN + 1)
273 #define RECCPY(dst, src) bufcpy (dst, src, RECLN);
274 #define NEW_RECORD(z) RECORD z; RECCLR(z)
275 #define NEW_STATIC_RECORD(z) static RECORD z; RECCLR(z)
276
277 typedef struct FTNFILE FTNFILE;
278 struct FTNFILE {
279 FILE *unit;
280 char *name, *form, *action, *disp, *buff, *rewind;
281 int_4 lrecl, vers, in_stream, record, records, redirect, buff_init, buff_pos, buff_len, memopen_len;
282 };
283
284 #define NEW_FTN_FILE(fp, fmt, act, recl)\
285 ((FTNFILE) {.unit = fp, .form = fmt, .action = act, .lrecl = recl, .in_stream = FALSE})
286
287 #define MAX_SOURCE_FILES 1500
288
289 // MAX_FILES should be 100.
290 #define MAX_FILES 100
291
292 extern FTNFILE _ffile[MAX_FILES];
293
294 typedef struct CALLS CALLS;
295 struct CALLS {
296 char *name;
297 int_8 calls;
298 };
299
300 typedef struct C_SRC C_SRC;
301 struct C_SRC
302 {
303 int_4 proc, phase;
304 char *text;
305 };
306 extern C_SRC *object;
307
308 typedef struct FTN_LINE FTN_LINE;
309 struct FTN_LINE
310 {
311 FTN_LINE *file;
312 int_4 num, isn, len, label, jcl, diag, proc, cpp;
313 char *text;
314 };
315 extern FTN_LINE *source, *files;
316 extern int_4 MAX_FTN_LINES;
317
318 typedef struct MODE MODE;
319 struct MODE
320 {
321 int_4 type, len, dim, save, fun;
322 };
323
324 #define PLAIN_MODE(t, l) ((MODE) {.type = type, .len = len, .dim = 0, .fun = FALSE})
325 #define TYPE(z, t, l) ((z)->mode.type == t && (z)->mode.len == l)
326 #define IS_ROW(m) ((m).dim > 0)
327 #define IS_SCALAR(m) ((m).dim == 0)
328
329 typedef struct IDENT IDENT;
330 struct IDENT
331 {
332 char *name, *fname, *parm;
333 int_4 arg, call, common, const_ref, external, intrinsic, line, macro, nest, save, source, used, variable;
334 IDENT *alias, *equiv;
335 MODE mode;
336 int_4 patch1, patch2;
337 char *lwb[MAX_DIMS], *upb[MAX_DIMS], *len[MAX_DIMS];
338 char *size;
339 };
340 #define C_NAME(z) ((z)->name)
341 #define FTN_NAME(z) ((z)->fname)
342
343 typedef struct IMPLICIT IMPLICIT;
344 struct IMPLICIT
345 {
346 MODE mode;
347 };
348
349 typedef struct LBL LBL;
350 struct LBL
351 {
352 int_4 index, num, line, jumped, patch, nonexe, data, format;
353 };
354 extern LBL *lbl;
355
356 typedef struct EXPR EXPR;
357 struct EXPR
358 {
359 RECORD str, elem;
360 int_4 variant, value, opt;
361 IDENT *idf;
362 MODE mode;
363 };
364
365 typedef struct INTRINS INTRINS;
366 struct INTRINS
367 {
368 char *ffun, *bfun, *cfun;
369 int_4 spec, rtyp, rlen, alen;
370 int_4 anum, atyp;
371 real_32 (*f1) (real_32);
372 real_32 (*f2) (real_32, real_32);
373 complex_64 (*f3) (complex_64);
374 };
375
376 extern char *action_default;
377 extern char *action_read;
378 extern char *action_readwrite;
379 extern char *action_write;
380 extern char *commons[], *modules[];
381 extern char *disp_delete;
382 extern char *disp_keep;
383 extern char *disp_new;
384 extern char *disp_old;
385 extern char *form_formatted;
386 extern char *form_unformatted;
387
388 extern IDENT globals[], locals[];
389
390 extern IMPLICIT implic[];
391
392 extern int_4 curlin, curcol;
393 extern int_4 curret, preret;
394 extern int_4 end_statements;
395 extern int_4 indent;
396 extern int_4 jcllin;
397 extern int_4 lhs_factor;
398 extern int_4 macro_nest;
399 extern int_4 MAX_C_SRC;
400 extern int_4 MAX_C_SRC;
401 extern int_4 MAX_FTN_LINES;
402 extern int_4 n_c_src;
403 extern int_4 nerrors, merrors, nwarns;
404 extern int_4 nloctmps, nglobtmps, func;
405 extern int_4 optimise;
406 extern int_4 page, line, ncommons, nmodules;
407 extern int_4 pnprocs, nprocs, nglobals, nlocals, nlabels, nftnlines, nfiles;
408 extern int_4 prelin, precol;
409
410 extern LBL labels[];
411
412 extern logical_4 abend;
413 extern logical_4 brief;
414 extern logical_4 use_strcasecmp;
415 extern logical_4 compile_only;
416 extern logical_4 f4_do_loops;
417 extern logical_4 gcc_ftn_lines;
418 extern logical_4 hollerith;
419 extern logical_4 implicit_r8;
420 extern logical_4 keep;
421 extern logical_4 no_warnings;
422 extern logical_4 pretty;
423 extern logical_4 quiet_mode;
424 extern logical_4 syntax_only;
425 extern logical_4 trace;
426
427 extern RECORD curlex, prelex;
428 extern RECORD hdate;
429 extern RECORD hmodule, hsection;
430 extern RECORD libnam, modnam, procnam;
431 extern RECORD program, block;
432 extern RECORD retnam;
433 extern RECORD stat_start;
434
435 extern char *bufcat (char *, char *, int_4);
436 extern char *bufcpy (char *, char *, int_4);
437 extern char *bufrep (char *, char *);
438 extern char *_bufsub (char *, char *, int_4, int_4);
439 extern char *c_name (char *);
440 extern char *concat (char *, char *, char *);
441 extern char *edit_f (char *);
442 extern char *edit_fmt (int_4);
443 extern char *edit_i (char *);
444 extern char *edit_tmp (int_4);
445 extern char *edit_unit (int_4);
446 extern char *edit_v (char *);
447 extern char *edit_vn (char *, int_4);
448 extern char *encode (char *, char *);
449 extern char *f2c_type (char *, MODE *, int_4, int_4);
450 extern char *f_strallocat (char *, char *);
451 extern char *f_stralloc (char *);
452 extern char *get_uniq_str (char *, char *);
453 extern char *idf_full_c_name (RECORD, IDENT *);
454 extern char *intnot (char *, int_8, int_4);
455 extern char *newpage (char *, char *);
456 extern char *pretty_float (char *);
457 extern char *ptr_to_array (IDENT *, int_4, int_4, int_4);
458 extern char *qtype (MODE *);
459 extern char *_strlower (char *);
460 extern char *wtype (MODE *, int_4, int_4);
461 extern char *xfixed (char *, real_32, int_4, int_4, int_4);
462 extern char *xfloat (char *, real_32, int_4, int_4, int_4, int_4, int_4, char);
463 extern char *xsubfixed (char *, real_32, logical_4, int_4);
464 extern char *xtoa (char *, real_32, int_4);
465
466 extern IDENT *add_local (char *, int_4, int_4, int_4, int_4, int_4, int_4, int_4);
467 extern IDENT *add_nest (char *, int_4, MODE *);
468 extern IDENT *extf_decl (char *, MODE *);
469 extern IDENT *find_local (char *, MODE *);
470 extern IDENT *impl_decl (char *, MODE *);
471 extern IDENT *void_decl (char *, MODE *);
472
473 extern int_4 add_block (char *);
474 extern int_4 code (int_4, int_4, char *);
475 extern int_4 code_real_32_const (char *);
476 extern int_4 code_uniq_str (char *);
477 extern int_4 format_str (char *str);
478 extern int_4 impl_do (void);
479 extern int_4 mix_len (EXPR *, EXPR *);
480 extern int_4 new_charlen (int_4);
481 extern int_4 ord (char);
482 extern int_4 patch (int_4, char *);
483 extern int_4 scan (char *);
484 extern int_4 scan_fmt (void);
485 extern int_4 scan_fmt (void);
486
487 extern LBL *find_label (char *);
488
489 extern logical_4 accept_mode (int_4, int_4, int_4, int_4);
490 extern logical_4 coerce (EXPR *, EXPR *);
491 extern logical_4 express (EXPR *, int_4, int_4);
492 extern logical_4 find_module (char *);
493 extern logical_4 fold_expr (EXPR *, int_4);
494 extern logical_4 fold_intrinsic (INTRINS *, EXPR *, EXPR *);
495 extern logical_4 intrinsic_call (char *, EXPR *);
496 extern logical_4 is_intrins (char *, MODE *);
497 extern logical_4 is_int4 (char *, int_4 *);
498 extern logical_4 is_macro_decl (char *);
499 extern logical_4 is_specific (char *);
500 extern logical_4 lookahead (char *);
501 extern logical_4 reserved (char *);
502 extern logical_4 same_name (char *, char *);
503 extern logical_4 valid_expr (EXPR *);
504
505 extern void assign (EXPR *);
506 extern void banner (int_4, int_4, char *);
507 extern void call (void);
508 extern void code_comment (void);
509 extern void code_common (void);
510 extern void code_exts (IDENT *, int_4, int_4, int_4, int_4);
511 extern void code_idfs (IDENT *, int_4, int_4, int_4, int_4);
512 extern void code_index (RECORD, IDENT *, int_4);
513 extern void code_parms (RECORD);
514 extern void code_row_len (IDENT *);
515 extern void common (void);
516 extern void common (void);
517 extern void compute_row_size (RECORD, IDENT *);
518 extern void cpp_direct (int_4, int_4, int_4);
519 extern void decl_autosave (void);
520 extern void decl_data (void);
521 extern void decl_equiv (void);
522 extern void decl_macros (void);
523 extern void default_impl (void);
524 extern void diagnostic (int_4, char *);
525 extern void dimension (void);
526 extern void do_data (int_4 *);
527 extern void do_io (char *, int_4 *);
528 extern void equivalence (void);
529 extern void executable (void);
530 extern void exprio (EXPR *, int_4, logical_4);
531 extern void externals (void);
532 extern void factor (EXPR *);
533 extern void factor_function_call (EXPR *, RECORD);
534 extern void factor_integer_number (EXPR *, char *);
535 extern void factor_slice_char (EXPR *, IDENT *);
536 extern void factor_slice (EXPR *, IDENT *);
537 extern void factor_variable (EXPR *, IDENT *, MODE *, RECORD);
538 extern void *f_malloc (size_t);
539 extern void fold_int_4 (char *, char *);
540 extern void format (LBL *);
541 extern void *f_realloc (void *, size_t);
542 extern void compile_nested_intrinsic (char *);
543 extern void gen_statements (LBL *, int_4);
544 extern void get_decls (void);
545 extern void get_dims (IDENT *, int_4);
546 extern void get_impl (void);
547 extern void get_source (char *, int_4);
548 extern void idfs_impl (void);
549 extern void idfs_unused (void);
550 extern void implicit (void);
551 extern void impl_type (char *, MODE *);
552 extern void inline_args (RECORD *, int_4 *);
553 extern void intrinsics (void);
554 extern void jump (void);
555 extern void macro (EXPR *, IDENT *);
556 extern void merge_commons (void);
557 extern void message (FTN_LINE *, int_4, char *, int_4, char *, char *);
558 extern void norm_mode (MODE *);
559 extern void option (char *);
560 extern void parameter (void);
561 extern void patch_args (void);
562 extern void patch_args (void);
563 extern void proc_listing (int_4);
564 extern void recursion (EXPR *, RECORD, IDENT *);
565 extern void relabel (char *);
566 extern void RTE (const char *, const char *);
567 extern void RTW (const char *, const char *);
568 extern void scan_modules (void);
569 extern void skip_card (int_4);
570 extern void subprograms (void);
571 extern void vif_backspace (void);
572 extern void vif_close (void);
573 extern void vif_endfile (void);
574 extern void vif_jcl (void);
575 extern void vif_open (void);
576 extern void vif_rewind (void);
577 extern void write_object (char *);
578
579 // ++++ MACROS
580
581 #define ln(x) log(x)
582
583 static inline logical_4 EQUAL (char *s, char *t)
584 {
585 if (s == NULL || t == NULL) {
586 return s == t;
587 } else {
588 return strcasecmp (s, t) == 0;
589 }
590 }
591
592 static inline logical_4 MATCH (char *t)
593 {
594 if (EQUAL (curlex, t)) {
595 return TRUE;
596 } else {
597 NEW_RECORD (str);
598 _srecordf (str, "\"%s\"", t); // Stringize t
599 return EQUAL (curlex, str);
600 }
601 }
602
603 static logical_4 LEQUAL (char *s, char *t)
604 {
605 if (s == NULL || t == NULL) {
606 return s == t;
607 } else {
608 return strncasecmp (s, t, strlen (s)) == 0;
609 }
610 }
611
612 static inline void BUG (char *s)
613 {
614 fprintf (stderr, "%d %s\n", curlin, source[curlin].text);
615 message (NULL, ERR, "fatal", 4001, "compiler bug", (s));
616 exit (EXIT_FAILURE);
617 }
618
619 #define _write_err(rc, funit, action)\
620 fflush (_ffile[(funit)].unit);\
621 if ((rc) < 1 || ferror (_ffile[(funit)].unit)) {\
622 action;\
623 }
624
625 #define _read_err(rc, funit, action_end, action_err)\
626 if (feof (_ffile[(funit)].unit)) {\
627 action_end;\
628 }\
629 if ((rc) < 1 || ferror (_ffile[(funit)].unit)) {\
630 action_err;\
631 }
632
633 #define _write_eol(funit) {\
634 (void) _vif_printf ((funit), "\n", NULL, NOTYPE, 0);\
635 fflush (_ffile[(funit)].unit);\
636 }
637
638 #define _read_eol(funit) {\
639 (void) _vif_scanf ((funit), NULL, NULL, NOTYPE, 0);\
640 }
641
642 #define _abs(a) ({__typeof__ (a) _u = (a); _u >= 0 ? _u : -_u;})
643 #define _dim_(a, b) ({__typeof__ (a) _a = (a); __typeof__ (b) _b = (b); _a > _b ? _a - _b : 0;})
644 #define _max(a, b) ({__typeof__ (a) _a = (a); __typeof__ (b) _b = (b); _a > _b ? _a : _b;})
645 #define _min(a, b) ({__typeof__ (a) _a = (a); __typeof__ (b) _b = (b); _a < _b ? _a : _b;})
646 #define _imod(a, b) ({__typeof__ (a) _a = (a), _b = (b), _q = (__typeof__ (a)) (_a / _b); (_a - _q * _b);})
647 #define _sign(a, b) ({__typeof__ (a) _a = _abs (a); b >= 0 ? _a : -_a;})
648
649 #define _ichar(s) ((int_4) ((s)[0]))
650
651 // AINT
652 static inline real_4 _aintf (real_4 x) {return truncf (x);}
653 static inline real_8 _aintd (real_8 x) {return trunc (x);}
654 static inline real_16 _aintq (real_16 x) {return truncq (x);}
655
656 // NINT
657 static inline int_4 _nintf (real_4 x) {return (int_4) (x) >= 0.0 ? floorf (x + 0.5) : -floorf (0.5 - x);}
658 static inline int_4 _nintd (real_8 x) {return (int_4) (x >= 0.0 ? floor (x + 0.5) : -floor (0.5 - x));}
659 static inline int_4 _nintq (real_16 x) {return (int_4) (x >= 0.0q ? floorq (x + 0.5q) : -floorq (0.5q - x));}
660
661 // ANINT
662 static inline real_4 _anintf (real_4 x) {return (x) >= 0.0 ? floorf (x + 0.5) : -floorf (0.5 - x);}
663 static inline real_8 _anintd (real_8 x) {return (x >= 0.0 ? floor (x + 0.5) : -floor (0.5 - x));}
664 static inline real_16 _anintq (real_8 x) {return (x >= 0.0q ? floor (x + 0.5q) : -floor (0.5q - x));}
665
666 // COMPLEX*32
667
668 static inline complex_32 CMPLXQ (real_16 re, real_16 im) {complex_32 z; __real__ z = re; __imag__ z = im; return z;}
669 extern complex_32 qxcmplx (real_32, real_32);
670
671 // COMPLEX*64
672
673 #define CMPLXX(re, im) (complex_64){re, im}
674 #define CMPLXZ(re) (complex_64){re, X_0}
675
676 // RTS
677
678 extern char *_char (int);
679 extern char *__strtok_r (char *, const char *, char **);
680 extern char *_strupper (char *);
681 extern char **_vif_jit (char *, char *);
682
683 extern void _pi4 (real_4 *);
684 extern void _pi8 (real_8 *);
685 extern void _pi16 (real_16 *);
686 extern void _pi32 (real_32 *);
687
688 extern complex_16 _dcmplxq (complex_32);
689 extern complex_16 _up_complex (complex_16, int_4);
690
691 extern complex_32 _qcmplxd (complex_16);
692 extern complex_32 _up_complex_32 (complex_32, int_4);
693
694 extern complex_8 _cmplxd (complex_16);
695 extern complex_8 _up_complex_8 (complex_8, int_4);
696
697 extern int_4 _backspace (char *, int_4);
698 extern int_4 _i1mach (int_4 *);
699 extern int_4 _index (char *, char *);
700 extern int_4 _init_file_buffer (int_4);
701 extern int_4 _rewind (char *, int_4);
702 extern int_4 _set_record (char *, int_4, int_4);
703 extern int_4 _sys (char *, char *, char *, const char *, ...);
704 extern int_4 _str_to_int4 (char *);
705 extern int_4 _up_int_4 (int_4, int_4);
706 extern int_4 _vif_printf (int_4, char *, void *, int_4, int_4);
707 extern int_4 _vif_scanf (int_4, char *, void *, int_4, int_4);
708 extern int_4 _xerclr (void);
709 extern int_4 _xermsg (char *, char *, char *, int_4 *, int_4 *);
710 extern int_4 _xgetf (int_4 *);
711 extern int_4 _xint4 (real_32);
712 extern int_4 _xsetf (int_4 *);
713
714 extern int_8 _up_int_8 (int_8, int_4);
715 extern int_8 _xint8 (real_32);
716 extern int_8 _xnint8 (real_32);
717
718 extern real_16 acotanq (real_16);
719 extern real_16 cotanq (real_16);
720 extern real_16 _qext (real_8);
721 extern real_16 _up_real_16 (real_16, int_4);
722 extern real_16 _strtoquad (char *, char **);
723 extern real_16 xtoquad (real_32);
724 extern real_16 _zabs_32 (real_16 re, real_16 im);
725 extern real_16 cximagq (complex_64);
726 extern real_16 cxrealq (complex_64);
727
728 extern real_32 _quadtop (real_32 *, real_16);
729 extern real_32 quadtox (real_16);
730 extern real_32 _xerf (real_32);
731 extern real_32 _aintx (real_32);
732 extern real_32 _xI (real_32);
733 extern real_32 _xmod (real_32, real_32);
734 extern real_32 _xdimx (real_32, real_32);
735 extern real_32 _anintx (real_32);
736 extern real_32 _zabs_64 (real_32 re, real_32 im);
737
738 extern real_4 acotanf (real_4);
739 extern real_4 cotanf (real_4);
740 extern real_4 _up_real_4 (real_4, int_4);
741 extern real_4 _zabs_8 (real_4 re, real_4 im);
742
743 extern real_8 acotan (real_8);
744 extern real_8 cotan (real_8);
745 extern real_8 _drand48 (void);
746 extern real_8 _seconds (void);
747 extern real_8 _str_to_real8 (char *);
748 extern real_8 _up_real_8 (real_8, int_4);
749 extern real_8 _zabs_16 (real_8 re, real_8 im);
750
751 extern void _cputim (real_8 *);
752 extern void _cputyd (int_4 *);
753 extern void _fcheck (char *, int_4, char *, char *);
754 extern void _fclose (int_4);
755 extern void _fprintf_real_32 (char *, char *, real_32, int_4, int_4);
756 extern void _fregister (char *, int_4, int_4, char *, char *, char *, char *);
757 extern void _funregister (char *, int_4);
758 extern void _ioend_read (char *, int_4);
759 extern void _ioerr (char *, int_4);
760 extern void _ioerr_read (char *, int_4);
761 extern void _ioerr_write (char *, int_4);
762 extern void _merfi (real_8 *, real_8 *, int_4 *);
763 extern void _qhex (real_16 *);
764 extern void _skip_eol (FILE *);
765 extern void _srand48 (int_4 *);
766 extern void _vif_exit (void);
767 extern void _vif_freq (CALLS *);
768 extern void _vif_init (void);
769 extern void _xhex (real_32 *);
770
771 extern complex_64 _coctotop (complex_64 *, real_32);
772 extern complex_64 _cquadtop (complex_64 *, complex_32);
773 extern complex_64 cxcos (complex_64);
774 extern complex_64 cxflt(complex_8);
775 extern complex_64 cxdbl(complex_16);
776 extern complex_64 cxdiv (complex_64, complex_64);
777 extern complex_64 cxexp (complex_64);
778 extern complex_64 cxlog (complex_64);
779 extern complex_64 cxmul (complex_64, complex_64);
780 extern complex_64 cxneg (complex_64);
781 extern complex_64 cxquad(complex_32);
782 extern complex_64 cxreal32(real_32);
783 extern complex_64 cxsin (complex_64);
784 extern complex_64 cxsqrt (complex_64);
785 extern complex_64 cxsub (complex_64, complex_64);
786 extern complex_64 cxsum (complex_64, complex_64);
787 extern complex_64 cxtan (complex_64);
788
789 extern int_4 xsgn (const real_32 *);
790
791 extern logical_4 xeq (real_32, real_32);
792 extern logical_4 xge (real_32, real_32);
793 extern logical_4 xgt (real_32, real_32);
794 extern logical_4 xis0 (const real_32 *);
795 extern logical_4 xis_minf (const real_32 *);
796 extern logical_4 xis_nan (const real_32 *);
797 extern logical_4 xis_pinf (const real_32 *);
798 extern logical_4 xle (real_32, real_32);
799 extern logical_4 xlt (real_32, real_32);
800 extern logical_4 xneq (real_32, real_32);
801 extern logical_4 xnot0 (const real_32 *);
802
803 extern real_32 atox (char *);
804 extern real_32 cximag (complex_64);
805 extern real_32 cxreal (complex_64);
806 extern real_32 dbltox (real_8);
807 extern real_32 flttox (real_4);
808 extern real_32 inttox (int_4);
809 extern real_32 strtox (char *, char **);
810 extern real_32 xabs (real_32);
811 extern real_32 xacosh (real_32);
812 extern real_32 xacos (real_32);
813 extern real_32 xacotan (real_32);
814 extern real_32 xadd (real_32, real_32, int_4);
815 extern real_32 xasinh (real_32);
816 extern real_32 xasin (real_32);
817 extern real_32 xatan2 (real_32, real_32);
818 extern real_32 xatanh (real_32);
819 extern real_32 xatan (real_32);
820 extern real_32 xcosh (real_32);
821 extern real_32 xcos (real_32);
822 extern real_32 xcotan (real_32);
823 extern real_32 xdiv (real_32, real_32);
824 extern real_32 xexp (real_32);
825 extern real_32 xfrac (real_32);
826 extern real_32 xlog10 (real_32);
827 extern real_32 xlog (real_32);
828 extern real_32 xmul (real_32, real_32);
829 extern real_32 xneg (real_32);
830 extern real_32 xpwr (real_32, int_4);
831 extern real_32 xround (real_32);
832 extern real_32 xsfmod (real_32, int_4 *);
833 extern real_32 xsinh (real_32);
834 extern real_32 xsin (real_32);
835 extern real_32 xsqrt (real_32);
836 extern real_32 xtanh (real_32);
837 extern real_32 xtan (real_32);
838 extern real_32 xtenup (int_4);
839 extern real_32 xtrunc (real_32);
840 extern real_4 xtoflt (real_32);
841 extern real_8 xtodbl (real_32);
842
843 // SLATEC message handling routines.
844
845 extern int_4 _j4save (int_4 _p_, int_4 _p_, logical_4 _p_);
846 extern int_4 _xerabt (char _p_, int_4 *);
847 extern int_4 _xerbla (char _p_, int_4 _p_);
848 extern int_4 _xerclr (void);
849 extern int_4 _xerdmp (void);
850 extern int_4 _xermax (int_4 _p_);
851 extern int_4 _xermsg (char _p_, char _p_, char _p_, int_4 _p_, int_4 _p_);
852 extern int_4 _xerprn (char _p_, int_4 _p_, char _p_, int_4 _p_);
853 extern int_4 _xersve (char _p_, char _p_, char _p_, int_4 _p_, int_4 _p_, int_4 _p_, int_4 _p_);
854 extern int_4 _xgetf (int_4 _p_);
855 extern int_4 _xgetua (int_4 _p_, int_4 _p_);
856 extern int_4 _xgetun (int_4 _p_);
857 extern int_4 _xsetf (int_4 _p_);
858 extern int_4 _xsetua (int_4 _p_, int_4 _p_);
859 extern int_4 _xsetun (int_4 _p_);
860
861 #endif
|
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|
