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