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