[BACK]Return to pexpr.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / io

Annotation of OpenXM_contrib2/asir2000/io/pexpr.c, Revision 1.33

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

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>