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