Annotation of OpenXM_contrib2/asir2000/io/pexpr.c, Revision 1.13
1.3 noro 1: /*
2: * Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED
3: * All rights reserved.
4: *
5: * FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited,
6: * non-exclusive and royalty-free license to use, copy, modify and
7: * redistribute, solely for non-commercial and non-profit purposes, the
8: * computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and
9: * conditions of this Agreement. For the avoidance of doubt, you acquire
10: * only a limited right to use the SOFTWARE hereunder, and FLL or any
11: * third party developer retains all rights, including but not limited to
12: * copyrights, in and to the SOFTWARE.
13: *
14: * (1) FLL does not grant you a license in any way for commercial
15: * purposes. You may use the SOFTWARE only for non-commercial and
16: * non-profit purposes only, such as academic, research and internal
17: * business use.
18: * (2) The SOFTWARE is protected by the Copyright Law of Japan and
19: * international copyright treaties. If you make copies of the SOFTWARE,
20: * with or without modification, as permitted hereunder, you shall affix
21: * to all such copies of the SOFTWARE the above copyright notice.
22: * (3) An explicit reference to this SOFTWARE and its copyright owner
23: * shall be made on your publication or presentation in any form of the
24: * results obtained by use of the SOFTWARE.
25: * (4) In the event that you modify the SOFTWARE, you shall notify FLL by
1.4 noro 26: * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification
1.3 noro 27: * for such modification or the source code of the modified part of the
28: * SOFTWARE.
29: *
30: * THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL
31: * MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND
32: * EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS
33: * FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES'
34: * RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY
35: * MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY.
36: * UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT,
37: * OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY
38: * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL
39: * DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES
40: * ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES
41: * FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY
42: * DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF
43: * SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART
44: * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY
45: * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
46: * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
1.13 ! noro 47: * $OpenXM: OpenXM_contrib2/asir2000/io/pexpr.c,v 1.12 2001/03/15 05:52:12 noro Exp $
1.3 noro 48: */
1.1 noro 49: #include "ca.h"
50: #include "al.h"
51: #include "parse.h"
52: #include "comp.h"
53: #include "base.h"
54:
55: #ifndef FPRINT
56: #define FPRINT
57: #endif
58:
1.2 noro 59: #define PRINTHAT (fortran_output?PUTS("**"):PUTS("^"))
1.1 noro 60:
61: #ifdef FPRINT
62: FILE *asir_out;
63: #define OUT asir_out
64: char DFORMAT[BUFSIZ];
1.2 noro 65: int hex_output;
66: int fortran_output;
1.6 noro 67: int double_output;
1.10 noro 68: int real_digit;
1.1 noro 69:
70: #define TAIL
71: #define PUTS(s) fputs(s,OUT)
72: #define PRINTF fprintf
73: #define PRINTN printn
74: #define PRINTBF printbf
75: #define PRINTCPLX printcplx
76: #define PRINTLM printlm
77: #define PRINTUP2 printup2
78: #define PRINTV printv
79: #define PRINTEXPR printexpr
80: #define PRINTNUM printnum
81: #define PRINTP printp
82: #define PRINTR printr
83: #define PRINTLIST printlist
84: #define PRINTVECT printvect
85: #define PRINTMAT printmat
86: #define PRINTSTR printstr
87: #define PRINTCOMP printcomp
88: #define PRINTDP printdp
89: #define PRINTUI printui
90: #define PRINTGF2MAT printgf2mat
91: #define PRINTGFMMAT printgfmmat
1.5 noro 92: #define PRINTBYTEARRAY printbytearray
1.13 ! noro 93: #define PRINTQUOTE printquote
1.1 noro 94: #define PRINTERR printerr
95: #define PRINTLF printlf
96: #define PRINTLOP printlop
97: #define PRINTFOP printfop
98: #define PRINTEOP printeop
99: #define PRINTQOP printqop
100: #define PRINTUP printup
101: #endif
102:
103: #ifdef SPRINT
104: static char *buf;
105: #define OUT buf
106: extern char DFORMAT[BUFSIZ];
1.2 noro 107: extern int hex_output;
108: extern int fortran_output;
1.1 noro 109:
110: #define TAIL while ( *OUT ) OUT++;
111: #define PUTS(s) strcat(OUT,s)
112: #define PRINTF sprintf
113: #define PRINTN sprintn
114: #define PRINTBF sprintbf
115: #define PRINTCPLX sprintcplx
116: #define PRINTLM sprintlm
117: #define PRINTUP2 sprintup2
118: #define PRINTV sprintv
119: #define PRINTEXPR sprintexpr
120: #define PRINTNUM sprintnum
121: #define PRINTP sprintp
122: #define PRINTR sprintr
123: #define PRINTLIST sprintlist
124: #define PRINTVECT sprintvect
125: #define PRINTMAT sprintmat
126: #define PRINTSTR sprintstr
127: #define PRINTCOMP sprintcomp
128: #define PRINTDP sprintdp
129: #define PRINTUI sprintui
130: #define PRINTGF2MAT sprintgf2mat
131: #define PRINTGFMMAT sprintgfmmat
1.5 noro 132: #define PRINTBYTEARRAY sprintbytearray
1.13 ! noro 133: #define PRINTQUOTE sprintquote
1.1 noro 134: #define PRINTERR sprinterr
135: #define PRINTLF sprintlf
136: #define PRINTLOP sprintlop
137: #define PRINTFOP sprintfop
138: #define PRINTEOP sprinteop
139: #define PRINTQOP sprintqop
1.2 noro 140: #define PRINTUP sprintup
1.1 noro 141: #endif
142:
143: void PRINTEXPR();
144: void PRINTNUM();
145: void PRINTN();
146: void PRINTV();
147: void PRINTP();
148: void PRINTR();
149: void PRINTLIST();
150: void PRINTVECT();
151: void PRINTMAT();
152: void PRINTSTR();
153: void PRINTCOMP();
154: void PRINTDP();
155: void PRINTUI();
156: void PRINTGF2MAT();
157: void PRINTGFMMAT();
1.5 noro 158: void PRINTBYTEARRAY();
1.13 ! noro 159: void PRINTQUOTE();
1.1 noro 160: void PRINTERR();
161: void PRINTCPLX();
162: void PRINTLM();
163: void PRINTLF();
1.2 noro 164: void PRINTUP2();
1.1 noro 165:
166: #ifdef FPRINT
167: void output_init() {
168: OUT = stdout;
169: sprintf(DFORMAT,"%%0%dd",DLENGTH);
170: }
171:
172: int mmono(p)
173: P p;
174: {
175: if ( NUM(p) )
1.9 saito 176: #if defined(INTERVAL)
177: if ( NID(p) != N_IP && NID(p) != N_ID && NID(p) != N_IT && NID(p) != N_IF
178: && compnum(CO,(Num)p,0) < 0 )
179: #else
1.1 noro 180: if ( compnum(CO,(Num)p,0) < 0 )
1.9 saito 181: #endif
1.1 noro 182: return ( 1 );
183: else
184: return ( 0 );
185: else if ( NEXT(DC(p)) )
186: return ( 0 );
187: else
188: return (mmono(COEF(DC(p))));
189: }
190:
191: #if PARI
192: void printbf(a)
193: BF a;
194: {
1.8 noro 195: sor(a->body,double_output ? 'f' : 'g',-1,0);
1.1 noro 196: }
197: #endif
198: #endif
199:
200: #ifdef SPRINT
201: void soutput_init(s)
202: char *s;
203: {
204: s[0] = 0; buf = s;
205: }
206:
207: #if PARI
1.8 noro 208: #include "genpari.h"
209:
210: void myoutbrute(g)
211: GEN g;
212: {
213: bruteall(g,'f',-1,1);
214: }
215:
1.1 noro 216: void sprintbf(a)
217: BF a;
218: {
1.2 noro 219: char *str;
220: char *GENtostr();
1.8 noro 221: char *GENtostr0();
1.2 noro 222:
1.8 noro 223: if ( double_output ) {
224: str = GENtostr0(a->body,myoutbrute);
225: } else {
226: str = GENtostr(a->body);
227: }
1.2 noro 228: TAIL PRINTF(OUT,"%s",str);
229: free(str);
1.1 noro 230: }
231: #endif
232: #endif
233:
234: void PRINTEXPR(vl,p)
235: VL vl;
236: Obj p;
237: {
238: if ( !p ) {
239: PRINTR(vl,(R)p);
240: return;
241: }
242:
243: switch ( OID(p) ) {
244: case O_N:
245: PRINTNUM((Num)p); break;
246: case O_P:
247: PRINTP(vl,(P)p); break;
248: case O_R:
249: PRINTR(vl,(R)p); break;
250: case O_LIST:
251: PRINTLIST(vl,(LIST)p); break;
252: case O_VECT:
253: PRINTVECT(vl,(VECT)p); break;
254: case O_MAT:
255: PRINTMAT(vl,(MAT)p); break;
256: case O_STR:
257: PRINTSTR((STRING)p); break;
258: case O_COMP:
259: PRINTCOMP(vl,(COMP)p); break;
260: case O_DP:
261: PRINTDP(vl,(DP)p); break;
262: case O_USINT:
263: PRINTUI(vl,(USINT)p); break;
264: case O_GF2MAT:
265: PRINTGF2MAT(vl,(GF2MAT)p); break;
266: case O_ERR:
267: PRINTERR(vl,(ERR)p); break;
268: case O_MATHCAP:
269: PRINTLIST(vl,((MATHCAP)p)->body); break;
270: case O_F:
271: PRINTLF(vl,(F)p); break;
272: case O_GFMMAT:
273: PRINTGFMMAT(vl,(GFMMAT)p); break;
1.5 noro 274: case O_BYTEARRAY:
275: PRINTBYTEARRAY(vl,(BYTEARRAY)p); break;
1.13 ! noro 276: case O_QUOTE:
! 277: PRINTQUOTE(vl,(QUOTE)p); break;
1.1 noro 278: default:
279: break;
280: }
281: }
282:
283: void PRINTN(n)
284: N n;
285: {
286: register int i,*ptr;
287: N tn;
288:
289: if ( !n ) {
290: PUTS("0");
291: return;
292: }
293: if ( hex_output ) {
294: ptr = BD(n);
295: TAIL PRINTF(OUT,"0x%x",ptr[PL(n)-1]);
296: if ( hex_output < 0 )
297: for ( i = PL(n) - 2; i >= 0; i-- ) {
298: TAIL PRINTF(OUT,"|%08x",ptr[i]);
299: }
300: else
301: for ( i = PL(n) - 2; i >= 0; i-- ) {
302: TAIL PRINTF(OUT,"%08x",ptr[i]);
303: }
304: } else {
305: ntobn(DBASE,n,&tn);
306: ptr = BD(tn);
307: TAIL PRINTF(OUT,"%d",ptr[PL(tn) - 1]);
308: for ( i = PL(tn) - 2; i >= 0; i-- ) {
309: TAIL PRINTF(OUT,DFORMAT,ptr[i]);
310: }
311: }
312: }
313:
1.9 saito 314: #if defined(ITVDEBUG)
315: void printbin(double z)
316: {
317: int i, j, mask;
318: union {
319: double x;
320: char c[8];
321: } a;
322:
323: a.x = z;
324: for(i=7;i>=0;i--) {
325: mask = 0x80;
326: for(j=0;j<8;j++) {
327: if (a.c[i] & mask) fprintf(stderr,"1");
328: else fprintf(stderr,"0");
329: mask >>= 1;
330: }
331: }
332: fprintf(stderr,"\n");
333: }
334: #endif
335:
336: #if 0
337: int printmode = PRINTF_E;
338: #else
339: int printmode = PRINTF_G;
340: #endif
341:
1.1 noro 342: void PRINTNUM(q)
343: Num q;
344: {
1.10 noro 345: char real_format[20];
346:
1.1 noro 347: if ( !q ) {
348: PUTS("0");
349: return;
350: }
351: switch ( NID(q) ) {
352: case N_Q:
353: if ( SGN((Q)q) == -1 )
354: PUTS("-");
355: PRINTN(NM((Q)q));
356: if ( !INT((Q)q) ) {
357: PUTS("/"); PRINTN(DN((Q)q));
358: }
359: break;
360: case N_R:
1.9 saito 361: switch (printmode) {
362: case PRINTF_E:
363: #if defined(INTERVAL)
364: case MID_PRINTF_E:
365: #endif
366: TAIL PRINTF(OUT,"%.16e",BDY((Real)q));
367: break;
368: case PRINTF_G:
369: #if defined(INTERVAL)
370: case MID_PRINTF_G:
371: #endif
372: default:
1.10 noro 373: if ( real_digit ) {
374: sprintf(real_format,
375: double_output?"%%.%df":"%%.%dg",real_digit);
376: TAIL PRINTF(OUT,real_format,BDY((Real)q));
377: } else {
378: TAIL PRINTF(OUT,double_output?"%f":"%g",BDY((Real)q));
379: }
1.9 saito 380: break;
381: }
1.1 noro 382: break;
383: case N_A:
384: PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")");
385: break;
386: #if PARI
387: case N_B:
388: PRINTBF((BF)q); break;
389: #endif
1.9 saito 390: #if defined(INTERVAL)
391: case N_IP:
392: case N_IF:
393: PUTS("[");
394: PRINTNUM(INF((Itv)q));
395: PUTS(",");
396: PRINTNUM(SUP((Itv)q));
397: PUTS("]");
398: break;
399: case N_ID:
400: switch (printmode) {
401: case PRINTF_E:
402: TAIL PRINTF(OUT, "[%.16e,%.16e]",INF((ItvD)q),SUP((ItvD)q));
403: #if defined(ITVDEBUG)
404: printbin(INF((ItvD)q));
405: printbin(SUP((ItvD)q));
406: #endif
407: break;
408: case MID_PRINTF_G:
409: TAIL PRINTF(OUT, "<%g,%g>", (SUP((ItvD)q)+INF((ItvD)q))*0.5,(SUP((ItvD)q)-INF((ItvD)q))*0.5);
410: break;
411: case MID_PRINTF_E:
412: TAIL PRINTF(OUT, "<%.16e,%.16e>", (SUP((ItvD)q)+INF((ItvD)q))*0.5,(SUP((ItvD)q)-INF((ItvD)q))*0.5);
413: break;
414: case PRINTF_G:
415: default:
416: TAIL PRINTF(OUT, "[%g,%g]",INF((ItvD)q),SUP((ItvD)q));
417: break;
418: }
419: break;
420: #endif
1.1 noro 421: case N_C:
422: PRINTCPLX((C)q); break;
423: case N_M:
424: TAIL PRINTF(OUT,"%d",CONT((MQ)q)); break;
425: case N_LM:
426: PRINTN(((LM)q)->body); break;
427: case N_GF2N:
428: if ( hex_output )
429: PRINTN((N)(((GF2N)q)->body));
430: else
431: PRINTUP2(((GF2N)q)->body);
432: break;
433: case N_GFPN:
434: PRINTUP((UP)(((GFPN)q)->body));
1.11 noro 435: break;
436: case N_GFS:
1.12 noro 437: TAIL PRINTF(OUT,"@_%d",CONT((GFS)q)); break;
1.1 noro 438: break;
439: }
440: }
441:
442: void PRINTCPLX(a)
443: C a;
444: {
445: PUTS("(");
446: if ( a->r )
447: PRINTNUM(a->r);
448: if ( a->i ) {
1.9 saito 449: #if defined(INTERVAL)
450: if ( a->r && ((compnum(0,a->i,0) > 0)
451: || NID(a->i) == N_IP || NID(a->i) == N_ID
452: || NID(a->i) == N_IT || NID(a->i) == N_IF) )
453: #else
1.1 noro 454: if ( a->r && (compnum(0,a->i,0) > 0) )
1.9 saito 455: #endif
1.1 noro 456: PUTS("+");
457: PRINTNUM(a->i); PUTS("*@i");
458: }
459: PUTS(")");
460: }
461:
462: void PRINTP(vl,p)
463: VL vl;
464: P p;
465: {
466: V v;
467: DCP dc;
468:
469: if ( !p )
470: PUTS("0");
471: else if ( NUM(p) )
472: PRINTNUM((Num)p);
473: else
474: for ( dc = DC(p), v = VR(p); dc; dc = NEXT(dc) ) {
475: if ( !DEG(dc) )
476: PRINTP(vl,COEF(dc));
477: else {
478: if ( NUM(COEF(dc)) && UNIQ((Q)COEF(dc)) ) {
479: ;
480: } else if ( NUM(COEF(dc)) && MUNIQ((Q)COEF(dc)) )
481: PUTS("-");
482: else if ( NUM(COEF(dc)) || !NEXT(DC(COEF(dc)))) {
483: PRINTP(vl,COEF(dc)); PUTS("*");
484: } else {
485: PUTS("("); PRINTP(vl,COEF(dc)); PUTS(")*");
486: }
487: PRINTV(vl,v);
488: if ( cmpq(DEG(dc),ONE) ) {
1.2 noro 489: PRINTHAT;
1.1 noro 490: if ( INT(DEG(dc)) && SGN(DEG(dc))>0 )
491: PRINTNUM((Num)DEG(dc));
492: else {
493: PUTS("("); PRINTNUM((Num)DEG(dc)); PUTS(")");
494: }
495: }
496: }
497: if ( NEXT(dc) ) {
498: P t;
499:
500: t = COEF(NEXT(dc));
501: if (!DEG(NEXT(dc))) {
502: if ( NUM(t) ) {
503: if ( !mmono(t) )
504: PUTS("+");
505: } else {
506: if (!mmono(COEF(DC(t))))
507: PUTS("+");
508: }
509: } else {
510: if ( !mmono(t) )
511: PUTS("+");
512: }
513: }
514: }
515: }
516:
517: int hideargs;
518:
519: void PRINTV(vl,v)
520: VL vl;
521: V v;
522: {
523: PF pf;
524: PFAD ad;
525: int i;
526:
527: if ( NAME(v) )
528: PUTS(NAME(v));
529: else if ( (vid)v->attr == V_PF ) {
530: pf = ((PFINS)v->priv)->pf; ad = ((PFINS)v->priv)->ad;
531: if ( !strcmp(NAME(pf),"pow") ) {
1.2 noro 532: PUTS("("); PRINTR(vl,(R)ad[0].arg); PUTS(")"); PRINTHAT; PUTS("(");
1.1 noro 533: PRINTR(vl,(R)ad[1].arg); PUTS(")");
534: } else if ( !pf->argc ) {
535: TAIL PRINTF(OUT,"%s",NAME(pf));
536: } else {
537: if ( hideargs ) {
538: for ( i = 0; i < pf->argc; i++ )
539: if ( ad[i].d )
540: break;
541: if ( i < pf->argc ) {
542: TAIL PRINTF(OUT,"%s{%d",NAME(pf),ad[0].d);
543: for ( i = 1; i < pf->argc; i++ ) {
544: TAIL PRINTF(OUT,",%d",ad[i].d);
545: }
546: PUTS("}");
547: } else {
548: TAIL PRINTF(OUT,"%s",NAME(pf));
549: }
550: } else {
551: for ( i = 0; i < pf->argc; i++ )
552: if ( ad[i].d )
553: break;
554: if ( i < pf->argc ) {
555: TAIL PRINTF(OUT,"%s{%d",NAME(pf),ad[0].d);
556: for ( i = 1; i < pf->argc; i++ ) {
557: TAIL PRINTF(OUT,",%d",ad[i].d);
558: }
559: PUTS("}(");
560: } else {
561: TAIL PRINTF(OUT,"%s(",NAME(pf));
562: }
563: PRINTR(vl,(R)ad[0].arg);
564: for ( i = 1; i < pf->argc; i++ ) {
565: PUTS(","); PRINTR(vl,(R)ad[i].arg);
566: }
567: PUTS(")");
568: }
569: }
570: }
571: }
572:
573: void PRINTR(vl,a)
574: VL vl;
575: R a;
576: {
577: if ( !a )
578: PUTS("0");
579: else
580: switch (OID(a)) {
581: case O_N: case O_P:
582: PRINTP(vl,(P)a); break;
583: default:
584: PUTS("("); PRINTP(vl,NM((R)a)); PUTS(")/("); PRINTP(vl,DN((R)a)); PUTS(")");
585: break;
586: }
587: }
588:
589: void PRINTVECT(vl,vect)
590: VL vl;
591: VECT vect;
592: {
593: int i;
594: pointer *ptr;
595:
596: PUTS("[ ");
597: for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {
598: PRINTEXPR(vl,ptr[i]); PUTS(" ");
599: }
600: PUTS("]");
601: }
602:
603: void PRINTMAT(vl,mat)
604: VL vl;
605: MAT mat;
606: {
607: int i,j,r,c;
608: pointer *ptr;
609:
610: for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) {
611: PUTS("[ ");
612: for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) {
613: PRINTEXPR(vl,ptr[j]); PUTS(" ");
614: }
615: PUTS("]");
616: if ( i < r - 1 )
617: PUTS("\n");
618: }
619: }
620:
621: void PRINTLIST(vl,list)
622: VL vl;
623: LIST list;
624: {
625: NODE tnode;
626:
627: PUTS("[");
628: for ( tnode = (NODE)list->body; tnode; tnode = NEXT(tnode) ) {
629: PRINTEXPR(vl,tnode->body);
630: if ( NEXT(tnode) )
631: PUTS(",");
632: }
633: PUTS("]");
634: }
635:
636: void PRINTSTR(str)
637: STRING str;
638: {
639: char *p;
640:
641: for ( p = BDY(str); *p; p++ )
642: if ( *p == '"' )
643: PUTS("\"");
644: else {
645: TAIL PRINTF(OUT,"%c",*p);
646: }
647: }
648:
649: void PRINTCOMP(vl,c)
650: VL vl;
651: COMP c;
652: {
653: int n,i;
654:
655: n = getcompsize((int)c->type);
656: PUTS("{");
657: for ( i = 0; i < n; i++ ) {
658: PRINTEXPR(vl,(pointer)c->member[i]);
659: if ( i < n-1 )
660: PUTS(",");
661: }
662: PUTS("}");
663: }
664:
665: void PRINTDP(vl,d)
666: VL vl;
667: DP d;
668: {
669: int n,i;
670: MP m;
671: DL dl;
672:
673: for ( n = d->nv, m = BDY(d); m; m = NEXT(m) ) {
674: PUTS("("); PRINTEXPR(vl,(pointer)m->c); PUTS(")*<<");
675: for ( i = 0, dl = m->dl; i < n-1; i++ ) {
676: TAIL PRINTF(OUT,"%d,",dl->d[i]);
677: }
678: TAIL PRINTF(OUT,"%d",dl->d[i]);
679: PUTS(">>");
680: if ( NEXT(m) )
681: PUTS("+");
682: }
683: }
684:
685: void PRINTUI(vl,u)
686: VL vl;
687: USINT u;
688: {
689: TAIL PRINTF(OUT,"%u",BDY(u));
690: }
691:
692: void PRINTGF2MAT(vl,mat)
693: VL vl;
694: GF2MAT mat;
695: {
696: int row,col,w,i,j,k,m;
697: unsigned int t;
698: unsigned int **b;
699:
700: row = mat->row;
701: col = mat->col;
702: w = (col+BSH-1)/BSH;
703: b = mat->body;
704: for ( i = 0; i < row; i++ ) {
705: for ( j = 0, m = 0; j < w; j++ ) {
706: t = b[i][j];
707: for ( k = 0; m < col && k < BSH; k++, m++ )
708: if ( t & (1<<k) )
709: PUTS("1");
710: else
711: PUTS("0");
712: }
713: PUTS("\n");
714: }
715: }
716:
717: void PRINTGFMMAT(vl,mat)
718: VL vl;
719: GFMMAT mat;
720: {
721: int row,col,i,j;
722: unsigned int t;
723: unsigned int **b;
724:
725: row = mat->row;
726: col = mat->col;
727: b = mat->body;
728: for ( i = 0; i < row; i++ ) {
729: PUTS("[");
730: for ( j = 0; j < col; j++ ) {
731: TAIL PRINTF(OUT,"%8d",b[i][j]);
732: }
733: PUTS("]\n");
734: }
1.5 noro 735: }
736:
737: void PRINTBYTEARRAY(vl,array)
738: VL vl;
739: BYTEARRAY array;
740: {
741: int len,i;
742: unsigned int t;
743: unsigned char *b;
744:
745: len = array->len;
746: b = array->body;
747: PUTS("|");
748: for ( i = 0; i < len-1; i++ ) {
749: TAIL PRINTF(OUT,"%02x ",(unsigned int)b[i]);
750: }
751: TAIL PRINTF(OUT,"%02x",(unsigned int)b[i]);
752: PUTS("|");
1.13 ! noro 753: }
! 754:
! 755: void PRINTQUOTE(vl,quote)
! 756: VL vl;
! 757: QUOTE quote;
! 758: {
! 759: PUTS("<...quoted...>");
1.1 noro 760: }
761:
762: void PRINTERR(vl,e)
763: VL vl;
764: ERR e;
765: {
766: PUTS("error("); PRINTEXPR(vl,e->body); PUTS(")");
767: }
768:
769: void PRINTUP2(p)
770: UP2 p;
771: {
772: int d,i;
773:
1.2 noro 774: if ( !p ) {
775: TAIL PRINTF(OUT,"0");
776: } else {
1.1 noro 777: d = degup2(p);
1.2 noro 778: TAIL PRINTF(OUT,"(");
779: if ( !d ) {
780: TAIL PRINTF(OUT,"1");
781: } else if ( d == 1 ) {
782: TAIL PRINTF(OUT,"@");
783: } else {
784: PUTS("@"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
785: }
1.1 noro 786: for ( i = d-1; i >= 0; i-- ) {
787: if ( p->b[i/BSH] & (1<<(i%BSH)) )
1.2 noro 788: if ( !i ) {
789: TAIL PRINTF(OUT,"+1");
790: } else if ( i == 1 ) {
791: TAIL PRINTF(OUT,"+@");
792: } else {
793: PUTS("+@"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
794: }
1.1 noro 795: }
1.2 noro 796: TAIL PRINTF(OUT,")");
1.1 noro 797: }
798: }
799:
800: void PRINTLF(vl,f)
801: VL vl;
802: F f;
803: {
804: switch ( FOP(f) ) {
805: case AL_TRUE:
806: TAIL PRINTF(OUT,"@true");
807: break;
808: case AL_FALSE:
809: TAIL PRINTF(OUT,"@false");
810: break;
811:
812: case AL_OR: case AL_AND:
813: PRINTFOP(vl,f); break;
814: case AL_NOT: case AL_IMPL: case AL_REPL: case AL_EQUIV:
815: PRINTEOP(vl,f); break;
816:
817: case AL_EQUAL: case AL_NEQ: case AL_LESSP:
818: case AL_GREATERP: case AL_LEQ: case AL_GEQ:
819: PRINTLOP(vl,f); break;
820:
821: case AL_EX: case AL_ALL:
822: PRINTQOP(vl,f); break;
823: default:
824: break;
825: }
826: }
827:
828: PRINTFOP(vl,f)
829: VL vl;
830: F f;
831: {
832: char *op;
833: NODE n;
834:
835: op = FOP(f)==AL_OR?" @|| ":" @&& ";
836: n = FJARG(f);
837: PUTS("("); PRINTEXPR(vl,BDY(n)); PUTS(")");
838: for ( n = NEXT(n); n; n = NEXT(n) ) {
839: PUTS(op); PUTS("("); PRINTEXPR(vl,BDY(n)); PUTS(")");
840: }
841: }
842:
843: PRINTEOP(vl,f)
844: VL vl;
845: F f;
846: {
847: oFOP op;
848: char *sop;
849:
850: if ( (op = FOP(f)) == AL_NOT ) {
851: PUTS("(@! "); PRINTEXPR(vl,(Obj)FARG(f)); PUTS(")"); return;
852: }
853: switch ( op ) {
854: case AL_IMPL:
855: sop = " @impl "; break;
856: case AL_REPL:
857: sop = " @repl "; break;
858: case AL_EQUIV:
859: sop = " @equiv "; break;
860: default:
861: break;
862: }
863: PUTS("(");
864: PRINTEXPR(vl,(Obj)FLHS(f));
865: PUTS(sop);
866: PRINTEXPR(vl,(Obj)FRHS(f));
867: PUTS(")");
868: }
869:
870: PRINTLOP(vl,f)
871: VL vl;
872: F f;
873: {
874: char *op;
875:
876: switch ( FOP(f) ) {
877: case AL_EQUAL:
878: op = " @== "; break;
879: case AL_NEQ:
880: op = " @!= "; break;
881: case AL_LESSP:
882: op = " @< "; break;
883: case AL_GREATERP:
884: op = " @> "; break;
885: case AL_LEQ:
886: op = " @<= "; break;
887: case AL_GEQ:
888: op = " @>= "; break;
889: default:
890: error("PRINTLOP : invalid operator");
891: break;
892: }
893: PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0");
894: }
895:
896: PRINTQOP(vl,f)
897: VL vl;
898: F f;
899: {
900: char *op;
901:
902: op = FOP(f)==AL_EX?"ex":"all";
903: TAIL PRINTF(OUT,"%s(%s,",op,NAME(FQVR(f)));
904: PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")");
905: }
906:
907: PRINTUP(n)
908: UP n;
909: {
910: int i,d;
911:
912: if ( !n )
913: PUTS("0");
914: else if ( !n->d )
915: PRINTNUM(n->c[0]);
916: else {
917: d = n->d;
918: PUTS("(");
919: if ( !d ) {
920: PRINTNUM(n->c[d]);
921: } else if ( d == 1 ) {
922: PRINTNUM(n->c[d]);
923: PUTS("*@p");
924: } else {
925: PRINTNUM(n->c[d]);
1.2 noro 926: PUTS("*@p"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
1.1 noro 927: }
928: for ( i = d-1; i >= 0; i-- ) {
929: if ( n->c[i] ) {
930: PUTS("+("); PRINTNUM(n->c[i]); PUTS(")");
931: if ( i >= 2 ) {
1.2 noro 932: PUTS("*@p"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
1.1 noro 933: } else if ( i == 1 )
934: PUTS("*@p");
935: }
936: }
937: PUTS(")");
938: }
939: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>