OpenXM_contrib2/asir2000/io/pexpr.c - diff

Return to pexpr.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / io

Diff for /OpenXM_contrib2/asir2000/io/pexpr.c between version 1.8 and 1.25

version 1.8, 2000/12/16 06:16:10

version 1.25, 2003/12/02 06:56:48

Line 44

* OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY

* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,

* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.

* $OpenXM: OpenXM_contrib2/asir2000/io/pexpr.c,v 1.7 2000/12/15 05:30:08 noro Exp $

* $OpenXM: OpenXM_contrib2/asir2000/io/pexpr.c,v 1.24 2003/06/09 16:18:10 saito Exp $

#include "ca.h"

#include "al.h"

Line 57

#endif

#define PRINTHAT (fortran_output?PUTS("**"):PUTS("^"))

extern int outputstyle;

#ifdef FPRINT

FILE *asir_out;

Line 65 char DFORMAT[BUFSIZ];

Line 66 char DFORMAT[BUFSIZ];

int hex_output;

int fortran_output;

int double_output;

int real_digit;

int real_binary;

int print_quote;

#define TAIL

#define PUTS(s) fputs(s,OUT)

Line 77 int double_output;

Line 81 int double_output;

#define PRINTV printv

#define PRINTEXPR printexpr

#define PRINTNUM printnum

#define PRINTP printp

#define PRINTP asir_printp

#define PRINTR printr

#define PRINTLIST printlist

#define PRINTVECT printvect

Line 89 int double_output;

Line 93 int double_output;

#define PRINTGF2MAT printgf2mat

#define PRINTGFMMAT printgfmmat

#define PRINTBYTEARRAY printbytearray

#define PRINTQUOTE printquote

#define PRINTERR printerr

#define PRINTLF printlf

#define PRINTLOP printlop

Line 96 int double_output;

Line 101 int double_output;

#define PRINTEOP printeop

#define PRINTQOP printqop

#define PRINTUP printup

#define PRINTUM printum

#define PRINTSF printsf

#endif

#ifdef SPRINT

Line 104 static char *buf;

Line 111 static char *buf;

extern char DFORMAT[BUFSIZ];

extern int hex_output;

extern int fortran_output;

extern int double_output;

extern int real_digit;

extern int print_quote;

#define TAIL while ( *OUT ) OUT++;

#define PUTS(s) strcat(OUT,s)

#define PRINTF sprintf

Line 116 extern int fortran_output;

Line 127 extern int fortran_output;

#define PRINTV sprintv

#define PRINTEXPR sprintexpr

#define PRINTNUM sprintnum

#define PRINTP sprintp

#define PRINTP asir_sprintp

#define PRINTR sprintr

#define PRINTLIST sprintlist

#define PRINTVECT sprintvect

Line 128 extern int fortran_output;

Line 139 extern int fortran_output;

#define PRINTGF2MAT sprintgf2mat

#define PRINTGFMMAT sprintgfmmat

#define PRINTBYTEARRAY sprintbytearray

#define PRINTQUOTE sprintquote

#define PRINTERR sprinterr

#define PRINTLF sprintlf

#define PRINTLOP sprintlop

Line 135 extern int fortran_output;

Line 147 extern int fortran_output;

#define PRINTEOP sprinteop

#define PRINTQOP sprintqop

#define PRINTUP sprintup

#define PRINTUM sprintum

#define PRINTSF sprintsf

#endif

void PRINTEXPR();

Line 153 void PRINTUI();

Line 167 void PRINTUI();

void PRINTGF2MAT();

void PRINTGFMMAT();

void PRINTBYTEARRAY();

void PRINTQUOTE();

void PRINTERR();

void PRINTCPLX();

void PRINTLM();

void PRINTLF();

void PRINTUP2();

void PRINTUP();

void PRINTUM();

void PRINTFOP();

void PRINTEOP();

void PRINTLOP();

void PRINTQOP();

void PRINTSF();

#ifdef FPRINT

void output_init() {

Line 169 int mmono(p)

Line 191 int mmono(p)

P p;

{

if ( NUM(p) )

#if defined(INTERVAL)

if ( NID(p) != N_IP && NID(p) != N_IntervalDouble && NID(p) != N_IntervalQuad && NID(p) != N_IntervalBigFloat

&& compnum(CO,(Num)p,0) < 0 )

#else

if ( compnum(CO,(Num)p,0) < 0 )

#endif

return ( 1 );

else

return ( 0 );

Line 179 P p;

Line 206 P p;

return (mmono(COEF(DC(p))));

}

#if PARI

#if defined(PARI)

void printbf(a)

BF a;

{

void sor();

sor(a->body,double_output ? 'f' : 'g',-1,0);

}

#endif

Line 195 char *s;

Line 224 char *s;

s[0] = 0; buf = s;

}

#if PARI

#if defined(PARI)

#include "genpari.h"

void myoutbrute(g)

Line 264 Obj p;

Line 293 Obj p;

PRINTGFMMAT(vl,(GFMMAT)p); break;

case O_BYTEARRAY:

PRINTBYTEARRAY(vl,(BYTEARRAY)p); break;

case O_QUOTE:

PRINTQUOTE(vl,(QUOTE)p); break;

default:

break;

}

Line 300 N n;

Line 331 N n;

}

#if defined(FPRINT)

#if defined(ITVDEBUG)

void printbin(double z)

{

int i, j, mask;

union {

double x;

char c[8];

} a;

a.x = z;

for(i=7;i>=0;i--) {

mask = 0x80;

for(j=0;j<8;j++) {

if (a.c[i] & mask) fprintf(stderr,"1");

else fprintf(stderr,"0");

mask >>= 1;

}

fprintf(stderr,"\n");

}

#endif

#if defined(FPRINT)

#if 0

int printmode = PRINTF_E;

#else

int printmode = PRINTF_G;

#endif

#else

extern int printmode;

#endif

void PRINTNUM(q)

Num q;

{

char real_format[20];

if ( !q ) {

PUTS("0");

return;

Line 317 Num q;

Line 384 Num q;

}

break;

case N_R:

TAIL PRINTF(OUT,double_output?"%f":"%g",BDY((Real)q));

switch (printmode) {

case PRINTF_E:

#if defined(INTERVAL)

case MID_PRINTF_E:

#endif

TAIL PRINTF(OUT,"%.16e",BDY((Real)q));

break;

case PRINTF_G:

#if defined(INTERVAL)

case MID_PRINTF_G:

#endif

default:

if ( real_binary ) {

unsigned int *m;

unsigned int u,l,mask;

int i,expo;

m = (unsigned int *)&BDY((Real)q);

#if defined(__i386__) || defined(MIPSEL) || defined(VISUAL) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__)

u = m[1]; l = m[0];

#else

u = m[0]; l = m[1];

#endif

if (u&0x80000000) {

TAIL PRINTF(OUT,"-");

}

u &= 0x7fffffff;

expo = ((int)(u>>20));

u &= 0xfffff;

if ( expo == 2047 ) {

if ( u || l ) {

TAIL PRINTF(OUT,"NaN");

} else {

TAIL PRINTF(OUT,"Inf");

}

} else if ( expo == 0 ) {

if ( u || l ) {

TAIL PRINTF(OUT,"0b0.");

for ( i = 0, mask = 0x80000; i < 20;

i++, mask >>= 1) {

TAIL

if ( u&mask )

PRINTF(OUT,"1");

else

PRINTF(OUT,"0");

}

for ( i = 0, mask = 0x80000000; i < 32;

i++, mask >>= 1) {

TAIL

if ( l&mask )

PRINTF(OUT,"1");

else

PRINTF(OUT,"0");

}

TAIL PRINTF(OUT,"*2^%d",-1022);

} else {

TAIL PRINTF(OUT,"0");

}

} else {

expo -= 1023;

TAIL PRINTF(OUT,"0b1.");

for ( i = 0, mask = 0x80000; i < 20;

i++, mask >>= 1) {

TAIL

if ( u&mask )

PRINTF(OUT,"1");

else

PRINTF(OUT,"0");

}

for ( i = 0, mask = 0x80000000; i < 32;

i++, mask >>= 1) {

TAIL

if ( l&mask )

PRINTF(OUT,"1");

else

PRINTF(OUT,"0");

}

if ( expo ) {

TAIL PRINTF(OUT,"*2^%d",expo);

}

} else if ( real_digit ) {

sprintf(real_format,

double_output?"%%.%df":"%%.%dg",real_digit);

TAIL PRINTF(OUT,real_format,BDY((Real)q));

} else {

TAIL PRINTF(OUT,double_output?"%f":"%g",BDY((Real)q));

}

break;

}

break;

case N_A:

PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")");

break;

#if PARI

#if defined(PARI)

case N_B:

PRINTBF((BF)q); break;

#endif

#if defined(INTERVAL)

case N_IP:

case N_IntervalBigFloat:

switch ( outputstyle ) {

case 1:

PUTS("intval(");

PRINTNUM(INF((Itv)q));

PUTS(",");

PRINTNUM(SUP((Itv)q));

PUTS(")");

break;

case 0:

default:

PUTS("[");

PRINTNUM(INF((Itv)q));

PUTS(",");

PRINTNUM(SUP((Itv)q));

PUTS("]");

break;

}

break;

case N_IntervalDouble:

switch (printmode) {

case PRINTF_E:

switch ( outputstyle ) {

case 1:

TAIL PRINTF(OUT, "intval(%.16e,%.16e)",

INF((IntervalDouble)q),SUP((IntervalDouble)q));

break;

case 0:

default:

TAIL PRINTF(OUT, "[%.16e,%.16e]",

INF((IntervalDouble)q),SUP((IntervalDouble)q));

break;

}

#if defined(ITVDEBUG)

printbin(INF((IntervalDouble)q));

printbin(SUP((IntervalDouble)q));

#endif

break;

case MID_PRINTF_G:

switch ( outputstyle ) {

case 1:

TAIL PRINTF(OUT, "intvalm(%g,%g)",

(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,

(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);

break;

case 0:

default:

TAIL PRINTF(OUT, "<%g,%g>",

(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,

(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);

break;

}

break;

case MID_PRINTF_E:

switch ( outputstyle ) {

case 1:

TAIL PRINTF(OUT, "intvalm(%.16e,%.16e)",

(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,

(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);

break;

case 0:

default:

TAIL PRINTF(OUT, "<%.16e,%.16e>",

(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,

(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);

break;

}

break;

case PRINTF_G:

default:

switch ( outputstyle ) {

case 1:

TAIL PRINTF(OUT, "intval(%g,%g)",

INF((IntervalDouble)q),SUP((IntervalDouble)q));

break;

case 0:

default:

TAIL PRINTF(OUT, "[%g,%g]",

INF((IntervalDouble)q),SUP((IntervalDouble)q));

break;

}

break;

}

break;

#endif

case N_C:

PRINTCPLX((C)q); break;

case N_M:

Line 341 Num q;

Line 584 Num q;

case N_GFPN:

PRINTUP((UP)(((GFPN)q)->body));

break;

case N_GFS:

TAIL PRINTF(OUT,"@_%d",CONT((GFS)q));

break;

case N_GFSN:

PRINTUM(BDY((GFSN)q));

break;

}

Line 351 C a;

Line 600 C a;

if ( a->r )

PRINTNUM(a->r);

if ( a->i ) {

#if defined(INTERVAL)

if ( a->r && ((compnum(0,a->i,0) > 0)

|| NID(a->i) == N_IP || NID(a->i) == N_IntervalDouble

|| NID(a->i) == N_IntervalQuad || NID(a->i) == N_IntervalBigFloat) )

#else

if ( a->r && (compnum(0,a->i,0) > 0) )

#endif

PUTS("+");

PRINTNUM(a->i); PUTS("*@i");

}

Line 428 V v;

Line 683 V v;

else if ( (vid)v->attr == V_PF ) {

pf = ((PFINS)v->priv)->pf; ad = ((PFINS)v->priv)->ad;

if ( !strcmp(NAME(pf),"pow") ) {

PUTS("("); PRINTR(vl,(R)ad[0].arg); PUTS(")"); PRINTHAT; PUTS("(");

PUTS("(("); PRINTR(vl,(R)ad[0].arg); PUTS(")"); PRINTHAT; PUTS("(");

PRINTR(vl,(R)ad[1].arg); PUTS(")");

PRINTR(vl,(R)ad[1].arg); PUTS("))");

} else if ( !pf->argc ) {

TAIL PRINTF(OUT,"%s",NAME(pf));

} else {

Line 492 VECT vect;

Line 747 VECT vect;

int i;

pointer *ptr;

PUTS("[ ");

switch ( outputstyle ) {

for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {

case 1:

PRINTEXPR(vl,ptr[i]); PUTS(" ");

PUTS("vect(");

for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {

if ( i != 0 ) PUTS(",");

PRINTEXPR(vl,ptr[i]);

}

PUTS(")");

break;

case 0:

default:

PUTS("[ ");

for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {

PRINTEXPR(vl,ptr[i]); PUTS(" ");

}

PUTS("]");

break;

}

PUTS("]");

}

void PRINTMAT(vl,mat)

Line 506 MAT mat;

Line 774 MAT mat;

int i,j,r,c;

pointer *ptr;

for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) {

switch ( outputstyle ) {

PUTS("[ ");

case 1:

for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) {

PUTS("mat(\n");

PRINTEXPR(vl,ptr[j]); PUTS(" ");

for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) {

}

if ( i != 0 ) PUTS(",\n");

PUTS("]");

PUTS("[ ");

if ( i < r - 1 )

for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) {

PUTS("\n");

if ( j != 0 ) PUTS(",");

PRINTEXPR(vl,ptr[j]);

}

PUTS(" ]");

}

PUTS(")");

break;

case 0:

default:

for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) {

PUTS("[ ");

for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) {

PRINTEXPR(vl,ptr[j]); PUTS(" ");

}

PUTS("]");

if ( i < r - 1 )

PUTS("\n");

}

break;

}

Line 618 VL vl;

Line 904 VL vl;

GFMMAT mat;

{

int row,col,i,j;

unsigned int t;

unsigned int **b;

row = mat->row;

Line 638 VL vl;

Line 923 VL vl;

BYTEARRAY array;

{

int len,i;

unsigned int t;

unsigned char *b;

len = array->len;

Line 651 BYTEARRAY array;

Line 935 BYTEARRAY array;

PUTS("|");

}

void PRINTQUOTE(vl,quote)

VL vl;

QUOTE quote;

{

LIST list;

if ( print_quote ) {

fnodetotree(BDY(quote),&list);

PRINTEXPR(vl,(Obj)list);

} else {

PUTS("<...quoted...>");

}

void PRINTERR(vl,e)

VL vl;

ERR e;

Line 717 F f;

Line 1015 F f;

}

PRINTFOP(vl,f)

void PRINTFOP(vl,f)

VL vl;

F f;

{

Line 732 F f;

Line 1030 F f;

}

PRINTEOP(vl,f)

void PRINTEOP(vl,f)

VL vl;

F f;

{

Line 759 F f;

Line 1057 F f;

PUTS(")");

}

PRINTLOP(vl,f)

void PRINTLOP(vl,f)

VL vl;

F f;

{

Line 785 F f;

Line 1083 F f;

PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0");

}

PRINTQOP(vl,f)

void PRINTQOP(vl,f)

VL vl;

F f;

{

Line 796 F f;

Line 1094 F f;

PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")");

}

PRINTUP(n)

void PRINTUP(n)

UP n;

{

int i,d;

Line 827 UP n;

Line 1125 UP n;

}

PUTS(")");

}

void PRINTUM(n)

UM n;

{

int i,d;

if ( !n )

PUTS("0");

else if ( !n->d )

PRINTSF(n->c[0]);

else {

d = n->d;

PUTS("(");

if ( !d ) {

PRINTSF(n->c[d]);

} else if ( d == 1 ) {

PRINTSF(n->c[d]);

PUTS("*@s");

} else {

PRINTSF(n->c[d]);

PUTS("*@s"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);

}

for ( i = d-1; i >= 0; i-- ) {

if ( n->c[i] ) {

PUTS("+("); PRINTSF(n->c[i]); PUTS(")");

if ( i >= 2 ) {

PUTS("*@s"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);

} else if ( i == 1 )

PUTS("*@s");

}

PUTS(")");

}

void PRINTSF(i)

unsigned int i;

{

if ( !i ) {

PUTS("0");

} else {

TAIL PRINTF(OUT,"@_%d",IFTOF(i));

}

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