Annotation of OpenXM_contrib2/asir2000/io/pexpr.c, Revision 1.16
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.16 ! noro 47: * $OpenXM: OpenXM_contrib2/asir2000/io/pexpr.c,v 1.15 2001/08/28 05:46:26 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.14 noro 69: int print_quote;
1.1 noro 70:
71: #define TAIL
72: #define PUTS(s) fputs(s,OUT)
73: #define PRINTF fprintf
74: #define PRINTN printn
75: #define PRINTBF printbf
76: #define PRINTCPLX printcplx
77: #define PRINTLM printlm
78: #define PRINTUP2 printup2
79: #define PRINTV printv
80: #define PRINTEXPR printexpr
81: #define PRINTNUM printnum
82: #define PRINTP printp
83: #define PRINTR printr
84: #define PRINTLIST printlist
85: #define PRINTVECT printvect
86: #define PRINTMAT printmat
87: #define PRINTSTR printstr
88: #define PRINTCOMP printcomp
89: #define PRINTDP printdp
90: #define PRINTUI printui
91: #define PRINTGF2MAT printgf2mat
92: #define PRINTGFMMAT printgfmmat
1.5 noro 93: #define PRINTBYTEARRAY printbytearray
1.13 noro 94: #define PRINTQUOTE printquote
1.1 noro 95: #define PRINTERR printerr
96: #define PRINTLF printlf
97: #define PRINTLOP printlop
98: #define PRINTFOP printfop
99: #define PRINTEOP printeop
100: #define PRINTQOP printqop
101: #define PRINTUP printup
1.16 ! noro 102: #define PRINTUM printum
! 103: #define PRINTSF printsf
1.1 noro 104: #endif
105:
106: #ifdef SPRINT
107: static char *buf;
108: #define OUT buf
109: extern char DFORMAT[BUFSIZ];
1.2 noro 110: extern int hex_output;
111: extern int fortran_output;
1.14 noro 112: extern int double_output;
113: extern int real_digit;
114: extern int print_quote;
1.1 noro 115:
116: #define TAIL while ( *OUT ) OUT++;
117: #define PUTS(s) strcat(OUT,s)
118: #define PRINTF sprintf
119: #define PRINTN sprintn
120: #define PRINTBF sprintbf
121: #define PRINTCPLX sprintcplx
122: #define PRINTLM sprintlm
123: #define PRINTUP2 sprintup2
124: #define PRINTV sprintv
125: #define PRINTEXPR sprintexpr
126: #define PRINTNUM sprintnum
127: #define PRINTP sprintp
128: #define PRINTR sprintr
129: #define PRINTLIST sprintlist
130: #define PRINTVECT sprintvect
131: #define PRINTMAT sprintmat
132: #define PRINTSTR sprintstr
133: #define PRINTCOMP sprintcomp
134: #define PRINTDP sprintdp
135: #define PRINTUI sprintui
136: #define PRINTGF2MAT sprintgf2mat
137: #define PRINTGFMMAT sprintgfmmat
1.5 noro 138: #define PRINTBYTEARRAY sprintbytearray
1.13 noro 139: #define PRINTQUOTE sprintquote
1.1 noro 140: #define PRINTERR sprinterr
141: #define PRINTLF sprintlf
142: #define PRINTLOP sprintlop
143: #define PRINTFOP sprintfop
144: #define PRINTEOP sprinteop
145: #define PRINTQOP sprintqop
1.2 noro 146: #define PRINTUP sprintup
1.16 ! noro 147: #define PRINTUM sprintum
! 148: #define PRINTSF sprintsf
1.1 noro 149: #endif
150:
151: void PRINTEXPR();
152: void PRINTNUM();
153: void PRINTN();
154: void PRINTV();
155: void PRINTP();
156: void PRINTR();
157: void PRINTLIST();
158: void PRINTVECT();
159: void PRINTMAT();
160: void PRINTSTR();
161: void PRINTCOMP();
162: void PRINTDP();
163: void PRINTUI();
164: void PRINTGF2MAT();
165: void PRINTGFMMAT();
1.5 noro 166: void PRINTBYTEARRAY();
1.13 noro 167: void PRINTQUOTE();
1.1 noro 168: void PRINTERR();
169: void PRINTCPLX();
170: void PRINTLM();
171: void PRINTLF();
1.2 noro 172: void PRINTUP2();
1.1 noro 173:
174: #ifdef FPRINT
175: void output_init() {
176: OUT = stdout;
177: sprintf(DFORMAT,"%%0%dd",DLENGTH);
178: }
179:
180: int mmono(p)
181: P p;
182: {
183: if ( NUM(p) )
1.9 saito 184: #if defined(INTERVAL)
185: if ( NID(p) != N_IP && NID(p) != N_ID && NID(p) != N_IT && NID(p) != N_IF
186: && compnum(CO,(Num)p,0) < 0 )
187: #else
1.1 noro 188: if ( compnum(CO,(Num)p,0) < 0 )
1.9 saito 189: #endif
1.1 noro 190: return ( 1 );
191: else
192: return ( 0 );
193: else if ( NEXT(DC(p)) )
194: return ( 0 );
195: else
196: return (mmono(COEF(DC(p))));
197: }
198:
199: #if PARI
200: void printbf(a)
201: BF a;
202: {
1.8 noro 203: sor(a->body,double_output ? 'f' : 'g',-1,0);
1.1 noro 204: }
205: #endif
206: #endif
207:
208: #ifdef SPRINT
209: void soutput_init(s)
210: char *s;
211: {
212: s[0] = 0; buf = s;
213: }
214:
215: #if PARI
1.8 noro 216: #include "genpari.h"
217:
218: void myoutbrute(g)
219: GEN g;
220: {
221: bruteall(g,'f',-1,1);
222: }
223:
1.1 noro 224: void sprintbf(a)
225: BF a;
226: {
1.2 noro 227: char *str;
228: char *GENtostr();
1.8 noro 229: char *GENtostr0();
1.2 noro 230:
1.8 noro 231: if ( double_output ) {
232: str = GENtostr0(a->body,myoutbrute);
233: } else {
234: str = GENtostr(a->body);
235: }
1.2 noro 236: TAIL PRINTF(OUT,"%s",str);
237: free(str);
1.1 noro 238: }
239: #endif
240: #endif
241:
242: void PRINTEXPR(vl,p)
243: VL vl;
244: Obj p;
245: {
246: if ( !p ) {
247: PRINTR(vl,(R)p);
248: return;
249: }
250:
251: switch ( OID(p) ) {
252: case O_N:
253: PRINTNUM((Num)p); break;
254: case O_P:
255: PRINTP(vl,(P)p); break;
256: case O_R:
257: PRINTR(vl,(R)p); break;
258: case O_LIST:
259: PRINTLIST(vl,(LIST)p); break;
260: case O_VECT:
261: PRINTVECT(vl,(VECT)p); break;
262: case O_MAT:
263: PRINTMAT(vl,(MAT)p); break;
264: case O_STR:
265: PRINTSTR((STRING)p); break;
266: case O_COMP:
267: PRINTCOMP(vl,(COMP)p); break;
268: case O_DP:
269: PRINTDP(vl,(DP)p); break;
270: case O_USINT:
271: PRINTUI(vl,(USINT)p); break;
272: case O_GF2MAT:
273: PRINTGF2MAT(vl,(GF2MAT)p); break;
274: case O_ERR:
275: PRINTERR(vl,(ERR)p); break;
276: case O_MATHCAP:
277: PRINTLIST(vl,((MATHCAP)p)->body); break;
278: case O_F:
279: PRINTLF(vl,(F)p); break;
280: case O_GFMMAT:
281: PRINTGFMMAT(vl,(GFMMAT)p); break;
1.5 noro 282: case O_BYTEARRAY:
283: PRINTBYTEARRAY(vl,(BYTEARRAY)p); break;
1.13 noro 284: case O_QUOTE:
285: PRINTQUOTE(vl,(QUOTE)p); break;
1.1 noro 286: default:
287: break;
288: }
289: }
290:
291: void PRINTN(n)
292: N n;
293: {
294: register int i,*ptr;
295: N tn;
296:
297: if ( !n ) {
298: PUTS("0");
299: return;
300: }
301: if ( hex_output ) {
302: ptr = BD(n);
303: TAIL PRINTF(OUT,"0x%x",ptr[PL(n)-1]);
304: if ( hex_output < 0 )
305: for ( i = PL(n) - 2; i >= 0; i-- ) {
306: TAIL PRINTF(OUT,"|%08x",ptr[i]);
307: }
308: else
309: for ( i = PL(n) - 2; i >= 0; i-- ) {
310: TAIL PRINTF(OUT,"%08x",ptr[i]);
311: }
312: } else {
313: ntobn(DBASE,n,&tn);
314: ptr = BD(tn);
315: TAIL PRINTF(OUT,"%d",ptr[PL(tn) - 1]);
316: for ( i = PL(tn) - 2; i >= 0; i-- ) {
317: TAIL PRINTF(OUT,DFORMAT,ptr[i]);
318: }
319: }
320: }
321:
1.15 noro 322: #if defined(FPRINT)
1.9 saito 323: #if defined(ITVDEBUG)
324: void printbin(double z)
325: {
326: int i, j, mask;
327: union {
328: double x;
329: char c[8];
330: } a;
331:
332: a.x = z;
333: for(i=7;i>=0;i--) {
334: mask = 0x80;
335: for(j=0;j<8;j++) {
336: if (a.c[i] & mask) fprintf(stderr,"1");
337: else fprintf(stderr,"0");
338: mask >>= 1;
339: }
340: }
341: fprintf(stderr,"\n");
342: }
343: #endif
1.15 noro 344: #endif
1.9 saito 345:
1.15 noro 346: #if defined(FPRINT)
1.9 saito 347: #if 0
348: int printmode = PRINTF_E;
349: #else
350: int printmode = PRINTF_G;
1.15 noro 351: #endif
352: #else
353: extern int printmode;
1.9 saito 354: #endif
355:
1.1 noro 356: void PRINTNUM(q)
357: Num q;
358: {
1.10 noro 359: char real_format[20];
360:
1.1 noro 361: if ( !q ) {
362: PUTS("0");
363: return;
364: }
365: switch ( NID(q) ) {
366: case N_Q:
367: if ( SGN((Q)q) == -1 )
368: PUTS("-");
369: PRINTN(NM((Q)q));
370: if ( !INT((Q)q) ) {
371: PUTS("/"); PRINTN(DN((Q)q));
372: }
373: break;
374: case N_R:
1.9 saito 375: switch (printmode) {
376: case PRINTF_E:
377: #if defined(INTERVAL)
378: case MID_PRINTF_E:
379: #endif
380: TAIL PRINTF(OUT,"%.16e",BDY((Real)q));
381: break;
382: case PRINTF_G:
383: #if defined(INTERVAL)
384: case MID_PRINTF_G:
385: #endif
386: default:
1.10 noro 387: if ( real_digit ) {
388: sprintf(real_format,
389: double_output?"%%.%df":"%%.%dg",real_digit);
390: TAIL PRINTF(OUT,real_format,BDY((Real)q));
391: } else {
392: TAIL PRINTF(OUT,double_output?"%f":"%g",BDY((Real)q));
393: }
1.9 saito 394: break;
395: }
1.1 noro 396: break;
397: case N_A:
398: PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")");
399: break;
400: #if PARI
401: case N_B:
402: PRINTBF((BF)q); break;
403: #endif
1.9 saito 404: #if defined(INTERVAL)
405: case N_IP:
406: case N_IF:
407: PUTS("[");
408: PRINTNUM(INF((Itv)q));
409: PUTS(",");
410: PRINTNUM(SUP((Itv)q));
411: PUTS("]");
412: break;
413: case N_ID:
414: switch (printmode) {
415: case PRINTF_E:
416: TAIL PRINTF(OUT, "[%.16e,%.16e]",INF((ItvD)q),SUP((ItvD)q));
417: #if defined(ITVDEBUG)
418: printbin(INF((ItvD)q));
419: printbin(SUP((ItvD)q));
420: #endif
421: break;
422: case MID_PRINTF_G:
423: TAIL PRINTF(OUT, "<%g,%g>", (SUP((ItvD)q)+INF((ItvD)q))*0.5,(SUP((ItvD)q)-INF((ItvD)q))*0.5);
424: break;
425: case MID_PRINTF_E:
426: TAIL PRINTF(OUT, "<%.16e,%.16e>", (SUP((ItvD)q)+INF((ItvD)q))*0.5,(SUP((ItvD)q)-INF((ItvD)q))*0.5);
427: break;
428: case PRINTF_G:
429: default:
430: TAIL PRINTF(OUT, "[%g,%g]",INF((ItvD)q),SUP((ItvD)q));
431: break;
432: }
433: break;
434: #endif
1.1 noro 435: case N_C:
436: PRINTCPLX((C)q); break;
437: case N_M:
438: TAIL PRINTF(OUT,"%d",CONT((MQ)q)); break;
439: case N_LM:
440: PRINTN(((LM)q)->body); break;
441: case N_GF2N:
442: if ( hex_output )
443: PRINTN((N)(((GF2N)q)->body));
444: else
445: PRINTUP2(((GF2N)q)->body);
446: break;
447: case N_GFPN:
448: PRINTUP((UP)(((GFPN)q)->body));
1.11 noro 449: break;
450: case N_GFS:
1.14 noro 451: TAIL PRINTF(OUT,"@_%d",CONT((GFS)q));
1.1 noro 452: break;
1.16 ! noro 453: case N_GFSN:
! 454: PRINTUM(BDY((GFSN)q));
! 455: break;
1.1 noro 456: }
457: }
458:
459: void PRINTCPLX(a)
460: C a;
461: {
462: PUTS("(");
463: if ( a->r )
464: PRINTNUM(a->r);
465: if ( a->i ) {
1.9 saito 466: #if defined(INTERVAL)
467: if ( a->r && ((compnum(0,a->i,0) > 0)
468: || NID(a->i) == N_IP || NID(a->i) == N_ID
469: || NID(a->i) == N_IT || NID(a->i) == N_IF) )
470: #else
1.1 noro 471: if ( a->r && (compnum(0,a->i,0) > 0) )
1.9 saito 472: #endif
1.1 noro 473: PUTS("+");
474: PRINTNUM(a->i); PUTS("*@i");
475: }
476: PUTS(")");
477: }
478:
479: void PRINTP(vl,p)
480: VL vl;
481: P p;
482: {
483: V v;
484: DCP dc;
485:
486: if ( !p )
487: PUTS("0");
488: else if ( NUM(p) )
489: PRINTNUM((Num)p);
490: else
491: for ( dc = DC(p), v = VR(p); dc; dc = NEXT(dc) ) {
492: if ( !DEG(dc) )
493: PRINTP(vl,COEF(dc));
494: else {
495: if ( NUM(COEF(dc)) && UNIQ((Q)COEF(dc)) ) {
496: ;
497: } else if ( NUM(COEF(dc)) && MUNIQ((Q)COEF(dc)) )
498: PUTS("-");
499: else if ( NUM(COEF(dc)) || !NEXT(DC(COEF(dc)))) {
500: PRINTP(vl,COEF(dc)); PUTS("*");
501: } else {
502: PUTS("("); PRINTP(vl,COEF(dc)); PUTS(")*");
503: }
504: PRINTV(vl,v);
505: if ( cmpq(DEG(dc),ONE) ) {
1.2 noro 506: PRINTHAT;
1.1 noro 507: if ( INT(DEG(dc)) && SGN(DEG(dc))>0 )
508: PRINTNUM((Num)DEG(dc));
509: else {
510: PUTS("("); PRINTNUM((Num)DEG(dc)); PUTS(")");
511: }
512: }
513: }
514: if ( NEXT(dc) ) {
515: P t;
516:
517: t = COEF(NEXT(dc));
518: if (!DEG(NEXT(dc))) {
519: if ( NUM(t) ) {
520: if ( !mmono(t) )
521: PUTS("+");
522: } else {
523: if (!mmono(COEF(DC(t))))
524: PUTS("+");
525: }
526: } else {
527: if ( !mmono(t) )
528: PUTS("+");
529: }
530: }
531: }
532: }
533:
534: int hideargs;
535:
536: void PRINTV(vl,v)
537: VL vl;
538: V v;
539: {
540: PF pf;
541: PFAD ad;
542: int i;
543:
544: if ( NAME(v) )
545: PUTS(NAME(v));
546: else if ( (vid)v->attr == V_PF ) {
547: pf = ((PFINS)v->priv)->pf; ad = ((PFINS)v->priv)->ad;
548: if ( !strcmp(NAME(pf),"pow") ) {
1.2 noro 549: PUTS("("); PRINTR(vl,(R)ad[0].arg); PUTS(")"); PRINTHAT; PUTS("(");
1.1 noro 550: PRINTR(vl,(R)ad[1].arg); PUTS(")");
551: } else if ( !pf->argc ) {
552: TAIL PRINTF(OUT,"%s",NAME(pf));
553: } else {
554: if ( hideargs ) {
555: for ( i = 0; i < pf->argc; i++ )
556: if ( ad[i].d )
557: break;
558: if ( i < pf->argc ) {
559: TAIL PRINTF(OUT,"%s{%d",NAME(pf),ad[0].d);
560: for ( i = 1; i < pf->argc; i++ ) {
561: TAIL PRINTF(OUT,",%d",ad[i].d);
562: }
563: PUTS("}");
564: } else {
565: TAIL PRINTF(OUT,"%s",NAME(pf));
566: }
567: } else {
568: for ( i = 0; i < pf->argc; i++ )
569: if ( ad[i].d )
570: break;
571: if ( i < pf->argc ) {
572: TAIL PRINTF(OUT,"%s{%d",NAME(pf),ad[0].d);
573: for ( i = 1; i < pf->argc; i++ ) {
574: TAIL PRINTF(OUT,",%d",ad[i].d);
575: }
576: PUTS("}(");
577: } else {
578: TAIL PRINTF(OUT,"%s(",NAME(pf));
579: }
580: PRINTR(vl,(R)ad[0].arg);
581: for ( i = 1; i < pf->argc; i++ ) {
582: PUTS(","); PRINTR(vl,(R)ad[i].arg);
583: }
584: PUTS(")");
585: }
586: }
587: }
588: }
589:
590: void PRINTR(vl,a)
591: VL vl;
592: R a;
593: {
594: if ( !a )
595: PUTS("0");
596: else
597: switch (OID(a)) {
598: case O_N: case O_P:
599: PRINTP(vl,(P)a); break;
600: default:
601: PUTS("("); PRINTP(vl,NM((R)a)); PUTS(")/("); PRINTP(vl,DN((R)a)); PUTS(")");
602: break;
603: }
604: }
605:
606: void PRINTVECT(vl,vect)
607: VL vl;
608: VECT vect;
609: {
610: int i;
611: pointer *ptr;
612:
613: PUTS("[ ");
614: for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {
615: PRINTEXPR(vl,ptr[i]); PUTS(" ");
616: }
617: PUTS("]");
618: }
619:
620: void PRINTMAT(vl,mat)
621: VL vl;
622: MAT mat;
623: {
624: int i,j,r,c;
625: pointer *ptr;
626:
627: for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) {
628: PUTS("[ ");
629: for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) {
630: PRINTEXPR(vl,ptr[j]); PUTS(" ");
631: }
632: PUTS("]");
633: if ( i < r - 1 )
634: PUTS("\n");
635: }
636: }
637:
638: void PRINTLIST(vl,list)
639: VL vl;
640: LIST list;
641: {
642: NODE tnode;
643:
644: PUTS("[");
645: for ( tnode = (NODE)list->body; tnode; tnode = NEXT(tnode) ) {
646: PRINTEXPR(vl,tnode->body);
647: if ( NEXT(tnode) )
648: PUTS(",");
649: }
650: PUTS("]");
651: }
652:
653: void PRINTSTR(str)
654: STRING str;
655: {
656: char *p;
657:
658: for ( p = BDY(str); *p; p++ )
659: if ( *p == '"' )
660: PUTS("\"");
661: else {
662: TAIL PRINTF(OUT,"%c",*p);
663: }
664: }
665:
666: void PRINTCOMP(vl,c)
667: VL vl;
668: COMP c;
669: {
670: int n,i;
671:
672: n = getcompsize((int)c->type);
673: PUTS("{");
674: for ( i = 0; i < n; i++ ) {
675: PRINTEXPR(vl,(pointer)c->member[i]);
676: if ( i < n-1 )
677: PUTS(",");
678: }
679: PUTS("}");
680: }
681:
682: void PRINTDP(vl,d)
683: VL vl;
684: DP d;
685: {
686: int n,i;
687: MP m;
688: DL dl;
689:
690: for ( n = d->nv, m = BDY(d); m; m = NEXT(m) ) {
691: PUTS("("); PRINTEXPR(vl,(pointer)m->c); PUTS(")*<<");
692: for ( i = 0, dl = m->dl; i < n-1; i++ ) {
693: TAIL PRINTF(OUT,"%d,",dl->d[i]);
694: }
695: TAIL PRINTF(OUT,"%d",dl->d[i]);
696: PUTS(">>");
697: if ( NEXT(m) )
698: PUTS("+");
699: }
700: }
701:
702: void PRINTUI(vl,u)
703: VL vl;
704: USINT u;
705: {
706: TAIL PRINTF(OUT,"%u",BDY(u));
707: }
708:
709: void PRINTGF2MAT(vl,mat)
710: VL vl;
711: GF2MAT mat;
712: {
713: int row,col,w,i,j,k,m;
714: unsigned int t;
715: unsigned int **b;
716:
717: row = mat->row;
718: col = mat->col;
719: w = (col+BSH-1)/BSH;
720: b = mat->body;
721: for ( i = 0; i < row; i++ ) {
722: for ( j = 0, m = 0; j < w; j++ ) {
723: t = b[i][j];
724: for ( k = 0; m < col && k < BSH; k++, m++ )
725: if ( t & (1<<k) )
726: PUTS("1");
727: else
728: PUTS("0");
729: }
730: PUTS("\n");
731: }
732: }
733:
734: void PRINTGFMMAT(vl,mat)
735: VL vl;
736: GFMMAT mat;
737: {
738: int row,col,i,j;
739: unsigned int t;
740: unsigned int **b;
741:
742: row = mat->row;
743: col = mat->col;
744: b = mat->body;
745: for ( i = 0; i < row; i++ ) {
746: PUTS("[");
747: for ( j = 0; j < col; j++ ) {
748: TAIL PRINTF(OUT,"%8d",b[i][j]);
749: }
750: PUTS("]\n");
751: }
1.5 noro 752: }
753:
754: void PRINTBYTEARRAY(vl,array)
755: VL vl;
756: BYTEARRAY array;
757: {
758: int len,i;
759: unsigned int t;
760: unsigned char *b;
761:
762: len = array->len;
763: b = array->body;
764: PUTS("|");
765: for ( i = 0; i < len-1; i++ ) {
766: TAIL PRINTF(OUT,"%02x ",(unsigned int)b[i]);
767: }
768: TAIL PRINTF(OUT,"%02x",(unsigned int)b[i]);
769: PUTS("|");
1.13 noro 770: }
771:
772: void PRINTQUOTE(vl,quote)
773: VL vl;
774: QUOTE quote;
775: {
1.14 noro 776: LIST list;
777:
778: if ( print_quote ) {
779: fnodetotree(BDY(quote),&list);
780: PRINTEXPR(vl,(Obj)list);
781: } else {
782: PUTS("<...quoted...>");
783: }
1.1 noro 784: }
785:
786: void PRINTERR(vl,e)
787: VL vl;
788: ERR e;
789: {
790: PUTS("error("); PRINTEXPR(vl,e->body); PUTS(")");
791: }
792:
793: void PRINTUP2(p)
794: UP2 p;
795: {
796: int d,i;
797:
1.2 noro 798: if ( !p ) {
799: TAIL PRINTF(OUT,"0");
800: } else {
1.1 noro 801: d = degup2(p);
1.2 noro 802: TAIL PRINTF(OUT,"(");
803: if ( !d ) {
804: TAIL PRINTF(OUT,"1");
805: } else if ( d == 1 ) {
806: TAIL PRINTF(OUT,"@");
807: } else {
808: PUTS("@"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
809: }
1.1 noro 810: for ( i = d-1; i >= 0; i-- ) {
811: if ( p->b[i/BSH] & (1<<(i%BSH)) )
1.2 noro 812: if ( !i ) {
813: TAIL PRINTF(OUT,"+1");
814: } else if ( i == 1 ) {
815: TAIL PRINTF(OUT,"+@");
816: } else {
817: PUTS("+@"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
818: }
1.1 noro 819: }
1.2 noro 820: TAIL PRINTF(OUT,")");
1.1 noro 821: }
822: }
823:
824: void PRINTLF(vl,f)
825: VL vl;
826: F f;
827: {
828: switch ( FOP(f) ) {
829: case AL_TRUE:
830: TAIL PRINTF(OUT,"@true");
831: break;
832: case AL_FALSE:
833: TAIL PRINTF(OUT,"@false");
834: break;
835:
836: case AL_OR: case AL_AND:
837: PRINTFOP(vl,f); break;
838: case AL_NOT: case AL_IMPL: case AL_REPL: case AL_EQUIV:
839: PRINTEOP(vl,f); break;
840:
841: case AL_EQUAL: case AL_NEQ: case AL_LESSP:
842: case AL_GREATERP: case AL_LEQ: case AL_GEQ:
843: PRINTLOP(vl,f); break;
844:
845: case AL_EX: case AL_ALL:
846: PRINTQOP(vl,f); break;
847: default:
848: break;
849: }
850: }
851:
852: PRINTFOP(vl,f)
853: VL vl;
854: F f;
855: {
856: char *op;
857: NODE n;
858:
859: op = FOP(f)==AL_OR?" @|| ":" @&& ";
860: n = FJARG(f);
861: PUTS("("); PRINTEXPR(vl,BDY(n)); PUTS(")");
862: for ( n = NEXT(n); n; n = NEXT(n) ) {
863: PUTS(op); PUTS("("); PRINTEXPR(vl,BDY(n)); PUTS(")");
864: }
865: }
866:
867: PRINTEOP(vl,f)
868: VL vl;
869: F f;
870: {
871: oFOP op;
872: char *sop;
873:
874: if ( (op = FOP(f)) == AL_NOT ) {
875: PUTS("(@! "); PRINTEXPR(vl,(Obj)FARG(f)); PUTS(")"); return;
876: }
877: switch ( op ) {
878: case AL_IMPL:
879: sop = " @impl "; break;
880: case AL_REPL:
881: sop = " @repl "; break;
882: case AL_EQUIV:
883: sop = " @equiv "; break;
884: default:
885: break;
886: }
887: PUTS("(");
888: PRINTEXPR(vl,(Obj)FLHS(f));
889: PUTS(sop);
890: PRINTEXPR(vl,(Obj)FRHS(f));
891: PUTS(")");
892: }
893:
894: PRINTLOP(vl,f)
895: VL vl;
896: F f;
897: {
898: char *op;
899:
900: switch ( FOP(f) ) {
901: case AL_EQUAL:
902: op = " @== "; break;
903: case AL_NEQ:
904: op = " @!= "; break;
905: case AL_LESSP:
906: op = " @< "; break;
907: case AL_GREATERP:
908: op = " @> "; break;
909: case AL_LEQ:
910: op = " @<= "; break;
911: case AL_GEQ:
912: op = " @>= "; break;
913: default:
914: error("PRINTLOP : invalid operator");
915: break;
916: }
917: PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0");
918: }
919:
920: PRINTQOP(vl,f)
921: VL vl;
922: F f;
923: {
924: char *op;
925:
926: op = FOP(f)==AL_EX?"ex":"all";
927: TAIL PRINTF(OUT,"%s(%s,",op,NAME(FQVR(f)));
928: PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")");
929: }
930:
931: PRINTUP(n)
932: UP n;
933: {
934: int i,d;
935:
936: if ( !n )
937: PUTS("0");
938: else if ( !n->d )
939: PRINTNUM(n->c[0]);
940: else {
941: d = n->d;
942: PUTS("(");
943: if ( !d ) {
944: PRINTNUM(n->c[d]);
945: } else if ( d == 1 ) {
946: PRINTNUM(n->c[d]);
947: PUTS("*@p");
948: } else {
949: PRINTNUM(n->c[d]);
1.2 noro 950: PUTS("*@p"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
1.1 noro 951: }
952: for ( i = d-1; i >= 0; i-- ) {
953: if ( n->c[i] ) {
954: PUTS("+("); PRINTNUM(n->c[i]); PUTS(")");
955: if ( i >= 2 ) {
1.2 noro 956: PUTS("*@p"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
1.1 noro 957: } else if ( i == 1 )
958: PUTS("*@p");
959: }
960: }
961: PUTS(")");
1.16 ! noro 962: }
! 963: }
! 964:
! 965: PRINTUM(n)
! 966: UM n;
! 967: {
! 968: int i,d;
! 969:
! 970: if ( !n )
! 971: PUTS("0");
! 972: else if ( !n->d )
! 973: PRINTSF(n->c[0]);
! 974: else {
! 975: d = n->d;
! 976: PUTS("(");
! 977: if ( !d ) {
! 978: PRINTSF(n->c[d]);
! 979: } else if ( d == 1 ) {
! 980: PRINTSF(n->c[d]);
! 981: PUTS("*@s");
! 982: } else {
! 983: PRINTSF(n->c[d]);
! 984: PUTS("*@s"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
! 985: }
! 986: for ( i = d-1; i >= 0; i-- ) {
! 987: if ( n->c[i] ) {
! 988: PUTS("+("); PRINTSF(n->c[i]); PUTS(")");
! 989: if ( i >= 2 ) {
! 990: PUTS("*@s"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
! 991: } else if ( i == 1 )
! 992: PUTS("*@s");
! 993: }
! 994: }
! 995: PUTS(")");
! 996: }
! 997: }
! 998:
! 999: PRINTSF(i)
! 1000: unsigned int i;
! 1001: {
! 1002: if ( !i ) {
! 1003: PUTS("0");
! 1004: } else {
! 1005: TAIL PRINTF(OUT,"@_%d",IFTOF(i));
1.1 noro 1006: }
1007: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>