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