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