Annotation of OpenXM_contrib2/asir2000/io/pexpr_body.c, Revision 1.6
1.6 ! noro 1: /* $OpenXM: OpenXM_contrib2/asir2000/io/pexpr_body.c,v 1.5 2004/12/02 13:48:43 noro Exp $ */
1.1 noro 2:
3: #define PRINTHAT (fortran_output?PUTS("**"):PUTS("^"))
4:
5: void PRINTEXPR();
6: void PRINTNUM();
7: void PRINTN();
8: void PRINTV();
9: void PRINTP();
10: void PRINTR();
11: void PRINTLIST();
12: void PRINTVECT();
13: void PRINTMAT();
14: void PRINTSTR();
15: void PRINTCOMP();
16: void PRINTDP();
17: void PRINTUI();
18: void PRINTGF2MAT();
19: void PRINTGFMMAT();
20: void PRINTBYTEARRAY();
21: void PRINTQUOTE();
1.4 noro 22: void PRINTQUOTEARG();
1.1 noro 23: void PRINTERR();
24: void PRINTCPLX();
25: void PRINTLM();
26: void PRINTLF();
27: void PRINTUP2();
28: void PRINTUP();
29: void PRINTUM();
30: void PRINTFOP();
31: void PRINTEOP();
32: void PRINTLOP();
33: void PRINTQOP();
34: void PRINTSF();
35: void PRINTSYMBOL();
36: void PRINTRANGE();
37: void PRINTTB();
1.3 noro 38: void PRINTDPV();
1.1 noro 39: void PRINTFNODE();
40: void PRINTBF();
1.5 noro 41: void PRINTDAlg();
1.1 noro 42:
43: void PRINTEXPR(vl,p)
44: VL vl;
45: Obj p;
46: {
47: if ( !p ) {
48: PRINTR(vl,(R)p);
49: } else
50: switch ( OID(p) ) {
51: case O_N:
52: PRINTNUM((Num)p); break;
53: case O_P:
54: PRINTP(vl,(P)p); break;
55: case O_R:
56: PRINTR(vl,(R)p); break;
57: case O_LIST:
58: PRINTLIST(vl,(LIST)p); break;
59: case O_VECT:
60: PRINTVECT(vl,(VECT)p); break;
61: case O_MAT:
62: PRINTMAT(vl,(MAT)p); break;
63: case O_STR:
64: PRINTSTR((STRING)p); break;
65: case O_COMP:
66: PRINTCOMP(vl,(COMP)p); break;
67: case O_DP:
68: PRINTDP(vl,(DP)p); break;
69: case O_USINT:
70: PRINTUI(vl,(USINT)p); break;
71: case O_GF2MAT:
72: PRINTGF2MAT(vl,(GF2MAT)p); break;
73: case O_ERR:
74: PRINTERR(vl,(ERR)p); break;
75: case O_MATHCAP:
76: PRINTLIST(vl,((MATHCAP)p)->body); break;
77: case O_F:
78: PRINTLF(vl,(F)p); break;
79: case O_GFMMAT:
80: PRINTGFMMAT(vl,(GFMMAT)p); break;
81: case O_BYTEARRAY:
82: PRINTBYTEARRAY(vl,(BYTEARRAY)p); break;
83: case O_QUOTE:
84: PRINTQUOTE(vl,(QUOTE)p); break;
1.4 noro 85: case O_QUOTEARG:
86: PRINTQUOTEARG(vl,(QUOTEARG)p); break;
1.1 noro 87: case O_SYMBOL:
88: PRINTSYMBOL((SYMBOL)p); break;
89: case O_RANGE:
90: PRINTRANGE(vl,(RANGE)p); break;
91: case O_TB:
92: PRINTTB(vl,(TB)p); break;
1.3 noro 93: case O_DPV:
94: PRINTDPV(vl,(DPV)p); break;
1.1 noro 95: default:
96: break;
97: }
98: }
99:
100: #if defined(FPRINT)
101: #if defined(ITVDEBUG)
102: void printbin(double z)
103: {
104: int i, j, mask;
105: union {
106: double x;
107: char c[8];
108: } a;
109:
110: a.x = z;
111: for(i=7;i>=0;i--) {
112: mask = 0x80;
113: for(j=0;j<8;j++) {
114: if (a.c[i] & mask) fprintf(stderr,"1");
115: else fprintf(stderr,"0");
116: mask >>= 1;
117: }
118: }
119: fprintf(stderr,"\n");
120: }
121: #endif
122: #endif
123:
124: void PRINTCPLX(a)
125: C a;
126: {
127: PUTS("(");
128: if ( a->r )
129: PRINTNUM(a->r);
130: if ( a->i ) {
131: #if defined(INTERVAL)
132: if ( a->r && ((compnum(0,a->i,0) > 0)
133: || NID(a->i) == N_IP || NID(a->i) == N_IntervalDouble
134: || NID(a->i) == N_IntervalQuad || NID(a->i) == N_IntervalBigFloat) )
135: #else
136: if ( a->r && (compnum(0,a->i,0) > 0) )
137: #endif
138: PUTS("+");
139: PRINTNUM(a->i); PUTS("*@i");
140: }
141: PUTS(")");
142: }
143:
144: void PRINTP(vl,p)
145: VL vl;
146: P p;
147: {
148: V v;
149: DCP dc;
150:
151: if ( !p )
152: PUTS("0");
153: else if ( NUM(p) )
154: PRINTNUM((Num)p);
155: else
156: for ( dc = DC(p), v = VR(p); dc; dc = NEXT(dc) ) {
157: if ( !DEG(dc) )
158: PRINTP(vl,COEF(dc));
159: else {
160: if ( NUM(COEF(dc)) && UNIQ((Q)COEF(dc)) ) {
161: ;
162: } else if ( NUM(COEF(dc)) && MUNIQ((Q)COEF(dc)) )
163: PUTS("-");
164: else if ( NUM(COEF(dc)) || !NEXT(DC(COEF(dc)))) {
165: PRINTP(vl,COEF(dc)); PUTS("*");
166: } else {
167: PUTS("("); PRINTP(vl,COEF(dc)); PUTS(")*");
168: }
169: PRINTV(vl,v);
170: if ( cmpq(DEG(dc),ONE) ) {
171: PRINTHAT;
172: if ( INT(DEG(dc)) && SGN(DEG(dc))>0 )
173: PRINTNUM((Num)DEG(dc));
174: else {
175: PUTS("("); PRINTNUM((Num)DEG(dc)); PUTS(")");
176: }
177: }
178: }
179: if ( NEXT(dc) ) {
180: P t;
181:
182: t = COEF(NEXT(dc));
183: if (!DEG(NEXT(dc))) {
184: if ( NUM(t) ) {
185: if ( !mmono(t) )
186: PUTS("+");
187: } else {
188: if (!mmono(COEF(DC(t))))
189: PUTS("+");
190: }
191: } else {
192: if ( !mmono(t) )
193: PUTS("+");
194: }
195: }
196: }
197: }
198:
199: void PRINTR(vl,a)
200: VL vl;
201: R a;
202: {
203: if ( !a )
204: PUTS("0");
205: else
206: switch (OID(a)) {
207: case O_N: case O_P:
208: PRINTP(vl,(P)a); break;
209: default:
210: PUTS("("); PRINTP(vl,NM((R)a)); PUTS(")/("); PRINTP(vl,DN((R)a)); PUTS(")");
211: break;
212: }
213: }
214:
215: void PRINTVECT(vl,vect)
216: VL vl;
217: VECT vect;
218: {
219: int i;
220: pointer *ptr;
221:
222: switch ( outputstyle ) {
223: case 1:
224: PUTS("vect(");
225: for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {
226: if ( i != 0 ) PUTS(",");
227: PRINTEXPR(vl,ptr[i]);
228: }
229: PUTS(")");
230: break;
231: case 0:
232: default:
233: PUTS("[ ");
234: for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {
235: PRINTEXPR(vl,ptr[i]); PUTS(" ");
236: }
237: PUTS("]");
238: break;
239: }
1.3 noro 240: }
241:
242: void PRINTDPV(vl,d)
243: VL vl;
244: DPV d;
245: {
246: int i;
247: DP *ptr;
248:
249: PUTS("modvect(");
250: for ( i = 0, ptr = BDY(d); i < d->len; i++ ) {
251: if ( i != 0 ) PUTS(",");
252: PRINTEXPR(vl,(Obj)ptr[i]);
253: }
254: PUTS(")");
1.1 noro 255: }
256:
257: void PRINTMAT(vl,mat)
258: VL vl;
259: MAT mat;
260: {
261: int i,j,r,c;
262: pointer *ptr;
263:
264: switch ( outputstyle ) {
265: case 1:
266: PUTS("mat(\n");
267: for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) {
268: if ( i != 0 ) PUTS(",\n");
269: PUTS("[ ");
270: for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) {
271: if ( j != 0 ) PUTS(",");
272: PRINTEXPR(vl,ptr[j]);
273: }
274: PUTS(" ]");
275: }
276: PUTS(")");
277: break;
278: case 0:
279: default:
280: for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) {
281: PUTS("[ ");
282: for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) {
283: PRINTEXPR(vl,ptr[j]); PUTS(" ");
284: }
285: PUTS("]");
286: if ( i < r - 1 )
287: PUTS("\n");
288: }
289: break;
290: }
291: }
292:
293: void PRINTLIST(vl,list)
294: VL vl;
295: LIST list;
296: {
297: NODE tnode;
298:
299: PUTS("[");
300: for ( tnode = (NODE)list->body; tnode; tnode = NEXT(tnode) ) {
301: PRINTEXPR(vl,tnode->body);
302: if ( NEXT(tnode) )
303: PUTS(",");
304: }
305: PUTS("]");
306: }
307:
308: void PRINTCOMP(vl,c)
309: VL vl;
310: COMP c;
311: {
312: int n,i;
313:
314: n = getcompsize((int)c->type);
315: PUTS("{");
316: for ( i = 0; i < n; i++ ) {
317: PRINTEXPR(vl,(pointer)c->member[i]);
318: if ( i < n-1 )
319: PUTS(",");
320: }
321: PUTS("}");
322: }
323:
324: void PRINTGF2MAT(vl,mat)
325: VL vl;
326: GF2MAT mat;
327: {
328: int row,col,w,i,j,k,m;
329: unsigned int t;
330: unsigned int **b;
331:
332: row = mat->row;
333: col = mat->col;
334: w = (col+BSH-1)/BSH;
335: b = mat->body;
336: for ( i = 0; i < row; i++ ) {
337: for ( j = 0, m = 0; j < w; j++ ) {
338: t = b[i][j];
339: for ( k = 0; m < col && k < BSH; k++, m++ )
340: if ( t & (1<<k) )
341: PUTS("1");
342: else
343: PUTS("0");
344: }
345: PUTS("\n");
346: }
347: }
348: void PRINTERR(vl,e)
349: VL vl;
350: ERR e;
351: {
352: PUTS("error("); PRINTEXPR(vl,e->body); PUTS(")");
353: }
354:
355: void PRINTLF(vl,f)
356: VL vl;
357: F f;
358: {
359: switch ( FOP(f) ) {
360: case AL_TRUE:
361: PUTS("@true");
362: break;
363: case AL_FALSE:
364: PUTS("@false");
365: break;
366:
367: case AL_OR: case AL_AND:
368: PRINTFOP(vl,f); break;
369: case AL_NOT: case AL_IMPL: case AL_REPL: case AL_EQUIV:
370: PRINTEOP(vl,f); break;
371:
372: case AL_EQUAL: case AL_NEQ: case AL_LESSP:
373: case AL_GREATERP: case AL_LEQ: case AL_GEQ:
374: PRINTLOP(vl,f); break;
375:
376: case AL_EX: case AL_ALL:
377: PRINTQOP(vl,f); break;
378: default:
379: break;
380: }
381: }
382:
383: void PRINTFOP(vl,f)
384: VL vl;
385: F f;
386: {
387: char *op;
388: NODE n;
389:
390: op = FOP(f)==AL_OR?" @|| ":" @&& ";
391: n = FJARG(f);
392: PUTS("("); PRINTEXPR(vl,BDY(n)); PUTS(")");
393: for ( n = NEXT(n); n; n = NEXT(n) ) {
394: PUTS(op); PUTS("("); PRINTEXPR(vl,BDY(n)); PUTS(")");
395: }
396: }
397:
398: void PRINTEOP(vl,f)
399: VL vl;
400: F f;
401: {
402: oFOP op;
403: char *sop;
404:
405: if ( (op = FOP(f)) == AL_NOT ) {
406: PUTS("(@! "); PRINTEXPR(vl,(Obj)FARG(f)); PUTS(")"); return;
407: }
408: switch ( op ) {
409: case AL_IMPL:
410: sop = " @impl "; break;
411: case AL_REPL:
412: sop = " @repl "; break;
413: case AL_EQUIV:
414: sop = " @equiv "; break;
415: default:
416: break;
417: }
418: PUTS("(");
419: PRINTEXPR(vl,(Obj)FLHS(f));
420: PUTS(sop);
421: PRINTEXPR(vl,(Obj)FRHS(f));
422: PUTS(")");
423: }
424:
425: void PRINTLOP(vl,f)
426: VL vl;
427: F f;
428: {
429: char *op;
430:
431: switch ( FOP(f) ) {
432: case AL_EQUAL:
433: op = " @== "; break;
434: case AL_NEQ:
435: op = " @!= "; break;
436: case AL_LESSP:
437: op = " @< "; break;
438: case AL_GREATERP:
439: op = " @> "; break;
440: case AL_LEQ:
441: op = " @<= "; break;
442: case AL_GEQ:
443: op = " @>= "; break;
444: default:
445: error("PRINTLOP : invalid operator");
446: break;
447: }
448: PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0");
449: }
450:
451: void PRINTRANGE(VL vl,RANGE p)
452: {
453: PUTS("range(");
454: PRINTEXPR(vl,p->start);
455: PUTS(",");
456: PRINTEXPR(vl,p->end);
457: PUTS(")");
458: }
459:
460: void PRINTFNODENODE(NODE n)
461: {
462: for ( ; n; n = NEXT(n) ) {
463: PRINTFNODE((FNODE)BDY(n),0);
464: if ( NEXT(n) ) PUTS(",");
465: }
466: }
467:
468: void PRINTFARGS(FNODE f)
469: {
470: NODE n;
471:
472: if ( f->id == I_LIST ) {
473: n = (NODE)FA0(f);
474: PRINTFNODENODE(n);
475: } else
476: PRINTFNODE(f,0);
477: }
478:
479: void PRINTFNODE(FNODE f,int paren)
480: {
481: NODE n,t,t0;
482: char vname[BUFSIZ],prefix[BUFSIZ];
483: char *opname,*vname_conv,*prefix_conv;
484: Obj obj;
485: int i,len,allzero,elen,elen2;
486: C cplx;
487: char *r;
488: FNODE fi,f2;
489:
490: if ( !f ) {
491: PUTS("(0)");
492: return;
493: }
494: if ( paren ) PUTS("(");
495: switch ( f->id ) {
496: /* unary operators */
497: case I_NOT: PUTS("!"); PRINTFNODE((FNODE)FA0(f),1); break;
498: case I_PAREN: PRINTFNODE((FNODE)FA0(f),0); break;
499: case I_MINUS: PUTS("-"); PRINTFNODE((FNODE)FA0(f),1); break;
500: /* binary operators */
501: /* arg list */
502: /* I_AND, I_OR => FA0(f), FA1(f) */
503: /* otherwise => FA1(f), FA2(f) */
504: case I_BOP:
505: PRINTFNODE((FNODE)FA1(f),1);
506: PUTS(((ARF)FA0(f))->name);
507: PRINTFNODE((FNODE)FA2(f),1);
508: break;
509: case I_COP:
510: switch( (cid)FA0(f) ) {
511: case C_EQ: opname = ("=="); break;
512: case C_NE: opname = ("!="); break;
513: case C_GT: opname = (">"); break;
514: case C_LT: opname = ("<"); break;
515: case C_GE: opname = (">="); break;
516: case C_LE: opname = ("<="); break;
517: }
518: PRINTFNODE((FNODE)FA1(f),1);
519: PUTS(opname);
520: PRINTFNODE((FNODE)FA2(f),1);
521: break;
522: case I_LOP:
523: switch( (lid)FA0(f) ) {
524: case L_EQ: opname = ("@=="); break;
525: case L_NE: opname = ("@!="); break;
526: case L_GT: opname = ("@>"); break;
527: case L_LT: opname = ("@<"); break;
528: case L_GE: opname = ("@>="); break;
529: case L_LE: opname = ("@<="); break;
530: case L_AND: opname = ("@&&"); break;
531: case L_OR: opname = ("@||"); break;
532: case L_NOT: opname = ("@!"); break;
533: }
534: if ( (lid)FA0(f)==L_NOT ) {
535: PUTS(opname); PRINTFNODE((FNODE)FA1(f),1);
536: } else {
537: PRINTFNODE((FNODE)FA1(f),1);
538: PUTS(opname);
539: PRINTFNODE((FNODE)FA2(f),1);
540: }
541: break;
542: case I_AND:
543: PRINTFNODE((FNODE)FA0(f),1);
544: PUTS("&&");
545: PRINTFNODE((FNODE)FA1(f),1);
546: break;
547: case I_OR:
548: PRINTFNODE((FNODE)FA0(f),1);
549: PUTS("!!");
550: PRINTFNODE((FNODE)FA1(f),1);
551: break;
552: /* ternary operators */
553: case I_CE:
554: PRINTFNODE((FNODE)FA0(f),1); PUTS("?"); PRINTFNODE((FNODE)FA1(f),1);
555: PUTS(":"); PRINTFNODE((FNODE)FA2(f),1);
556: break;
557: /* lists */
558: case I_LIST: PUTS("["); PRINTFNODENODE((NODE)FA0(f)); PUTS("]"); break;
559: /* function */
560: case I_FUNC:
561: if ( !strcmp(((FUNC)FA0(f))->name,"@pi") ) PUTS("@pi");
562: else if ( !strcmp(((FUNC)FA0(f))->name,"@e") ) PUTS("@e");
563: else {
564: PUTS(((FUNC)FA0(f))->name);
565: PUTS("("); PRINTFARGS(FA1(f)); PUTS(")");
566: }
567: break;
568: /* XXX */
569: case I_CAR: PUTS("car("); PRINTFNODE(FA0(f),0); PUTS(")"); break;
570: case I_CDR: PUTS("cdr("); PRINTFNODE(FA0(f),0); PUTS(")"); break;
571: /* exponent vector */
572: case I_EV: PUTS("<<"); PRINTFNODENODE((NODE)FA0(f)); PUTS(">>"); break;
573: /* string */
574: case I_STR: PUTS((char *)FA0(f)); break;
575: /* internal object */
576: case I_FORMULA: obj = (Obj)FA0(f); PRINTEXPR(CO,obj); break;
577: /* program variable */
578: case I_PVAR:
579: if ( FA1(f) )
580: error("printfnode : not implemented yet");
581: GETPVNAME(FA0(f),opname);
582: PUTS(opname);
583: break;
584: default: error("printfnode : not implemented yet");
585: }
586: if ( paren ) PUTS(")");
587: }
588:
589: /* functions not used in cpexpr.c */
590:
591: #if !defined(CPRINT)
592:
593: void PRINTNUM(q)
594: Num q;
595: {
596: char real_format[20];
1.5 noro 597: DAlg d;
598: DP nm;
599: Q dn;
1.1 noro 600:
601: if ( !q ) {
602: PUTS("0");
603: return;
604: }
605: switch ( NID(q) ) {
606: case N_Q:
607: if ( SGN((Q)q) == -1 )
608: PUTS("-");
609: PRINTN(NM((Q)q));
610: if ( !INT((Q)q) ) {
611: PUTS("/"); PRINTN(DN((Q)q));
612: }
613: break;
614: case N_R:
615: switch (printmode) {
616: case PRINTF_E:
617: #if defined(INTERVAL)
618: case MID_PRINTF_E:
619: #endif
620: TAIL PRINTF(OUT,"%.16e",BDY((Real)q));
621: break;
622: case PRINTF_G:
623: #if defined(INTERVAL)
624: case MID_PRINTF_G:
625: #endif
626: default:
627: if ( real_binary ) {
628: unsigned int *m;
629: unsigned int u,l,mask;
630: int i,expo;
631:
632: m = (unsigned int *)&BDY((Real)q);
633: #if defined(__i386__) || defined(MIPSEL) || defined(VISUAL) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) || defined(__x86_64)
634: u = m[1]; l = m[0];
635: #else
636: u = m[0]; l = m[1];
637: #endif
638: if (u&0x80000000) {
639: TAIL PRINTF(OUT,"-");
640: }
641: u &= 0x7fffffff;
642: expo = ((int)(u>>20));
643: u &= 0xfffff;
644: if ( expo == 2047 ) {
645: if ( u || l ) {
646: TAIL PRINTF(OUT,"NaN");
647: } else {
648: TAIL PRINTF(OUT,"Inf");
649: }
650: } else if ( expo == 0 ) {
651: if ( u || l ) {
652: TAIL PRINTF(OUT,"0b0.");
653: for ( i = 0, mask = 0x80000; i < 20;
654: i++, mask >>= 1) {
655: TAIL
656: if ( u&mask )
657: PRINTF(OUT,"1");
658: else
659: PRINTF(OUT,"0");
660: }
661: for ( i = 0, mask = 0x80000000; i < 32;
662: i++, mask >>= 1) {
663: TAIL
664: if ( l&mask )
665: PRINTF(OUT,"1");
666: else
667: PRINTF(OUT,"0");
668: }
669: TAIL PRINTF(OUT,"*2^%d",-1022);
670: } else {
671: TAIL PRINTF(OUT,"0");
672: }
673: } else {
674: expo -= 1023;
675: TAIL PRINTF(OUT,"0b1.");
676: for ( i = 0, mask = 0x80000; i < 20;
677: i++, mask >>= 1) {
678: TAIL
679: if ( u&mask )
680: PRINTF(OUT,"1");
681: else
682: PRINTF(OUT,"0");
683: }
684: for ( i = 0, mask = 0x80000000; i < 32;
685: i++, mask >>= 1) {
686: TAIL
687: if ( l&mask )
688: PRINTF(OUT,"1");
689: else
690: PRINTF(OUT,"0");
691: }
692: if ( expo ) {
693: TAIL PRINTF(OUT,"*2^%d",expo);
694: }
695: }
696: } else if ( real_digit ) {
697: sprintf(real_format,
698: double_output?"%%.%df":"%%.%dg",real_digit);
699: TAIL PRINTF(OUT,real_format,BDY((Real)q));
700: } else {
701: TAIL PRINTF(OUT,double_output?"%f":"%g",BDY((Real)q));
702: }
703: break;
704: }
705: break;
706: case N_A:
707: PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")");
708: break;
709: #if defined(PARI)
710: case N_B:
711: PRINTBF((BF)q); break;
712: #endif
713: #if defined(INTERVAL)
714: case N_IP:
715: case N_IntervalBigFloat:
716: switch ( outputstyle ) {
717: case 1:
718: PUTS("intval(");
719: PRINTNUM(INF((Itv)q));
720: PUTS(",");
721: PRINTNUM(SUP((Itv)q));
722: PUTS(")");
723: break;
724: case 0:
725: default:
726: PUTS("[");
727: PRINTNUM(INF((Itv)q));
728: PUTS(",");
729: PRINTNUM(SUP((Itv)q));
730: PUTS("]");
731: break;
732: }
733: break;
734: case N_IntervalDouble:
735: switch (printmode) {
736: case PRINTF_E:
737: switch ( outputstyle ) {
738: case 1:
739: TAIL PRINTF(OUT, "intval(%.16e,%.16e)",
740: INF((IntervalDouble)q),SUP((IntervalDouble)q));
741: break;
742: case 0:
743: default:
744: TAIL PRINTF(OUT, "[%.16e,%.16e]",
745: INF((IntervalDouble)q),SUP((IntervalDouble)q));
746: break;
747: }
748: #if defined(ITVDEBUG)
749: printbin(INF((IntervalDouble)q));
750: printbin(SUP((IntervalDouble)q));
751: #endif
752: break;
753: case MID_PRINTF_G:
754: switch ( outputstyle ) {
755: case 1:
756: TAIL PRINTF(OUT, "intvalm(%g,%g)",
757: (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,
758: (SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);
759: break;
760: case 0:
761: default:
762: TAIL PRINTF(OUT, "<%g,%g>",
763: (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,
764: (SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);
765: break;
766: }
767: break;
768: case MID_PRINTF_E:
769: switch ( outputstyle ) {
770: case 1:
771: TAIL PRINTF(OUT, "intvalm(%.16e,%.16e)",
772: (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,
773: (SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);
774: break;
775: case 0:
776: default:
777: TAIL PRINTF(OUT, "<%.16e,%.16e>",
778: (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,
779: (SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);
780: break;
781: }
782: break;
783: case PRINTF_G:
784: default:
785: switch ( outputstyle ) {
786: case 1:
787: TAIL PRINTF(OUT, "intval(%g,%g)",
788: INF((IntervalDouble)q),SUP((IntervalDouble)q));
789: break;
790: case 0:
791: default:
792: TAIL PRINTF(OUT, "[%g,%g]",
793: INF((IntervalDouble)q),SUP((IntervalDouble)q));
794: break;
795: }
796: break;
797: }
798: break;
799: #endif
800: case N_C:
801: PRINTCPLX((C)q); break;
802: case N_M:
803: TAIL PRINTF(OUT,"%d",CONT((MQ)q)); break;
804: case N_LM:
805: PRINTN(((LM)q)->body); break;
806: case N_GF2N:
807: if ( hex_output )
808: PRINTN((N)(((GF2N)q)->body));
809: else
810: PRINTUP2(((GF2N)q)->body);
811: break;
812: case N_GFPN:
813: PRINTUP((UP)(((GFPN)q)->body));
814: break;
815: case N_GFS:
816: TAIL PRINTF(OUT,"@_%d",CONT((GFS)q));
817: break;
818: case N_GFSN:
819: PRINTUM(BDY((GFSN)q));
1.5 noro 820: break;
821: case N_DA:
822: d = (DAlg)q;
823: nm = d->nm;
824: dn = d->dn;
825: if ( SGN((Q)dn) == -1 ) PUTS("-");
826: PUTS("(");
827: PRINTDP(CO,((DAlg)q)->nm);
828: PUTS(")");
829: if ( !UNIN(NM(dn)) ) {
830: PUTS("/");
831: PRINTN(NM(dn));
832: }
1.1 noro 833: break;
834: }
835: }
836:
837: void PRINTN(n)
838: N n;
839: {
840: register int i,*ptr;
841: N tn;
842:
843: if ( !n ) {
844: PUTS("0");
845: return;
846: }
847: if ( hex_output ) {
848: ptr = BD(n);
849: TAIL PRINTF(OUT,"0x%x",ptr[PL(n)-1]);
850: if ( hex_output < 0 )
851: for ( i = PL(n) - 2; i >= 0; i-- ) {
852: TAIL PRINTF(OUT,"|%08x",ptr[i]);
853: }
854: else
855: for ( i = PL(n) - 2; i >= 0; i-- ) {
856: TAIL PRINTF(OUT,"%08x",ptr[i]);
857: }
858: } else {
859: ntobn(DBASE,n,&tn);
860: ptr = BD(tn);
861: TAIL PRINTF(OUT,"%d",ptr[PL(tn) - 1]);
862: for ( i = PL(tn) - 2; i >= 0; i-- ) {
863: TAIL PRINTF(OUT,DFORMAT,ptr[i]);
864: }
865: }
866: }
867:
868: void PRINTV(vl,v)
869: VL vl;
870: V v;
871: {
872: PF pf;
873: PFAD ad;
874: int i;
875:
876: if ( NAME(v) )
877: PUTS(NAME(v));
878: else if ( (vid)v->attr == V_PF ) {
879: pf = ((PFINS)v->priv)->pf; ad = ((PFINS)v->priv)->ad;
880: if ( !strcmp(NAME(pf),"pow") ) {
881: PUTS("(("); PRINTR(vl,(R)ad[0].arg); PUTS(")"); PRINTHAT; PUTS("(");
882: PRINTR(vl,(R)ad[1].arg); PUTS("))");
883: } else if ( !pf->argc ) {
884: TAIL PRINTF(OUT,"%s",NAME(pf));
885: } else {
886: if ( hideargs ) {
887: for ( i = 0; i < pf->argc; i++ )
888: if ( ad[i].d )
889: break;
890: if ( i < pf->argc ) {
891: TAIL PRINTF(OUT,"%s{%d",NAME(pf),ad[0].d);
892: for ( i = 1; i < pf->argc; i++ ) {
893: TAIL PRINTF(OUT,",%d",ad[i].d);
894: }
895: PUTS("}");
896: } else {
897: TAIL PRINTF(OUT,"%s",NAME(pf));
898: }
899: } else {
900: for ( i = 0; i < pf->argc; i++ )
901: if ( ad[i].d )
902: break;
903: if ( i < pf->argc ) {
904: TAIL PRINTF(OUT,"%s{%d",NAME(pf),ad[0].d);
905: for ( i = 1; i < pf->argc; i++ ) {
906: TAIL PRINTF(OUT,",%d",ad[i].d);
907: }
908: PUTS("}(");
909: } else {
910: TAIL PRINTF(OUT,"%s(",NAME(pf));
911: }
912: PRINTR(vl,(R)ad[0].arg);
913: for ( i = 1; i < pf->argc; i++ ) {
914: PUTS(","); PRINTR(vl,(R)ad[i].arg);
915: }
916: PUTS(")");
917: }
918: }
919: }
920: }
921:
922: void PRINTSTR(str)
923: STRING str;
924: {
925: char *p;
926:
927: for ( p = BDY(str); *p; p++ )
928: if ( *p == '"' )
929: PUTS("\"");
930: else {
931: TAIL PRINTF(OUT,"%c",*p);
932: }
933: }
934:
935: void PRINTDP(vl,d)
936: VL vl;
937: DP d;
938: {
939: int n,i;
940: MP m;
941: DL dl;
942:
943: for ( n = d->nv, m = BDY(d); m; m = NEXT(m) ) {
944: PUTS("("); PRINTEXPR(vl,(pointer)m->c); PUTS(")*<<");
945: for ( i = 0, dl = m->dl; i < n-1; i++ ) {
946: TAIL PRINTF(OUT,"%d,",dl->d[i]);
947: }
948: TAIL PRINTF(OUT,"%d",dl->d[i]);
949: PUTS(">>");
950: if ( NEXT(m) )
951: PUTS("+");
952: }
953: }
954:
955: void PRINTUI(vl,u)
956: VL vl;
957: USINT u;
958: {
959: TAIL PRINTF(OUT,"%u",BDY(u));
960: }
961:
962: void PRINTGFMMAT(vl,mat)
963: VL vl;
964: GFMMAT mat;
965: {
966: int row,col,i,j;
967: unsigned int **b;
968:
969: row = mat->row;
970: col = mat->col;
971: b = mat->body;
972: for ( i = 0; i < row; i++ ) {
973: PUTS("[");
974: for ( j = 0; j < col; j++ ) {
975: TAIL PRINTF(OUT,"%8d",b[i][j]);
976: }
977: PUTS("]\n");
978: }
979: }
980:
981: void PRINTBYTEARRAY(vl,array)
982: VL vl;
983: BYTEARRAY array;
984: {
985: int len,i;
986: unsigned char *b;
987:
988: len = array->len;
989: b = array->body;
990: PUTS("|");
991: for ( i = 0; i < len-1; i++ ) {
992: TAIL PRINTF(OUT,"%02x ",(unsigned int)b[i]);
993: }
994: TAIL PRINTF(OUT,"%02x",(unsigned int)b[i]);
995: PUTS("|");
996: }
997:
998: void PRINTQUOTE(vl,quote)
999: VL vl;
1000: QUOTE quote;
1001: {
1002: LIST list;
1003:
1004: if ( print_quote == 1 ) {
1005: fnodetotree(BDY(quote),&list);
1006: PRINTEXPR(vl,(Obj)list);
1007: } else if ( print_quote == 2 ) {
1008: PRINTFNODE(BDY(quote),0);
1009: } else {
1010: PUTS("<...quoted...>");
1011: }
1.4 noro 1012: }
1013:
1014: void PRINTQUOTEARG(vl,quote)
1015: VL vl;
1016: QUOTEARG quote;
1017: {
1018: PUTS("<...quoted...>");
1.1 noro 1019: }
1020:
1021: void PRINTTB(VL vl,TB p)
1022: {
1023: int i;
1.2 noro 1024: #if defined(SPRINT)
1025: char *ptr,*e;
1.1 noro 1026:
1.2 noro 1027: for ( i = 0, ptr = buf+strlen(OUT); i < p->next; i++ )
1028: for ( e = p->body[i]; *e; *ptr++ = *e++ );
1029: *ptr = 0;
1030: #else
1.1 noro 1031: for ( i = 0; i < p->next; i++ ) {
1032: PUTS(p->body[i]);
1033: }
1.2 noro 1034: #endif
1.1 noro 1035: }
1036:
1037: void PRINTUP2(p)
1038: UP2 p;
1039: {
1040: int d,i;
1041:
1042: if ( !p ) {
1043: TAIL PRINTF(OUT,"0");
1044: } else {
1045: d = degup2(p);
1046: TAIL PRINTF(OUT,"(");
1047: if ( !d ) {
1048: TAIL PRINTF(OUT,"1");
1049: } else if ( d == 1 ) {
1050: TAIL PRINTF(OUT,"@");
1051: } else {
1052: PUTS("@"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
1053: }
1054: for ( i = d-1; i >= 0; i-- ) {
1055: if ( p->b[i/BSH] & (1<<(i%BSH)) )
1056: if ( !i ) {
1057: TAIL PRINTF(OUT,"+1");
1058: } else if ( i == 1 ) {
1059: TAIL PRINTF(OUT,"+@");
1060: } else {
1061: PUTS("+@"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
1062: }
1063: }
1064: TAIL PRINTF(OUT,")");
1065: }
1066: }
1067:
1068: void PRINTQOP(vl,f)
1069: VL vl;
1070: F f;
1071: {
1072: char *op;
1073:
1074: op = FOP(f)==AL_EX?"ex":"all";
1075: TAIL PRINTF(OUT,"%s(%s,",op,NAME(FQVR(f)));
1076: PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")");
1077: }
1078:
1079: void PRINTUP(n)
1080: UP n;
1081: {
1082: int i,d;
1083:
1084: if ( !n )
1085: PUTS("0");
1086: else if ( !n->d )
1087: PRINTNUM(n->c[0]);
1088: else {
1089: d = n->d;
1090: PUTS("(");
1091: if ( !d ) {
1092: PRINTNUM(n->c[d]);
1093: } else if ( d == 1 ) {
1094: PRINTNUM(n->c[d]);
1095: PUTS("*@p");
1096: } else {
1097: PRINTNUM(n->c[d]);
1098: PUTS("*@p"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
1099: }
1100: for ( i = d-1; i >= 0; i-- ) {
1101: if ( n->c[i] ) {
1102: PUTS("+("); PRINTNUM(n->c[i]); PUTS(")");
1103: if ( i >= 2 ) {
1104: PUTS("*@p"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
1105: } else if ( i == 1 )
1106: PUTS("*@p");
1107: }
1108: }
1109: PUTS(")");
1110: }
1111: }
1112:
1113: void PRINTUM(n)
1114: UM n;
1115: {
1116: int i,d;
1117:
1118: if ( !n )
1119: PUTS("0");
1120: else if ( !n->d )
1121: PRINTSF(n->c[0]);
1122: else {
1123: d = n->d;
1124: PUTS("(");
1125: if ( !d ) {
1126: PRINTSF(n->c[d]);
1127: } else if ( d == 1 ) {
1128: PRINTSF(n->c[d]);
1129: PUTS("*@s");
1130: } else {
1131: PRINTSF(n->c[d]);
1132: PUTS("*@s"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
1133: }
1134: for ( i = d-1; i >= 0; i-- ) {
1135: if ( n->c[i] ) {
1136: PUTS("+("); PRINTSF(n->c[i]); PUTS(")");
1137: if ( i >= 2 ) {
1138: PUTS("*@s"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
1139: } else if ( i == 1 )
1140: PUTS("*@s");
1141: }
1142: }
1143: PUTS(")");
1144: }
1145: }
1146:
1147: void PRINTSF(i)
1148: unsigned int i;
1149: {
1150: if ( !i ) {
1151: PUTS("0");
1152: } else {
1153: TAIL PRINTF(OUT,"@_%d",IFTOF(i));
1154: }
1155: }
1156:
1157: void PRINTSYMBOL(SYMBOL sym)
1158: {
1159: PUTS(sym->name);
1160: }
1161:
1162: #endif /* CPRINT */
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>