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