OpenXM_contrib2/asir2000/io/spexpr.c - diff

Return to spexpr.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / io

Diff for /OpenXM_contrib2/asir2000/io/spexpr.c between version 1.6 and 1.26

version 1.6, 2000/12/13 10:54:09

version 1.26, 2004/02/04 07:42:07

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/spexpr.c,v 1.5 2000/11/08 08:02:51 noro Exp $

#include "ca.h"

#include "al.h"

Line 52

Line 51

#include "comp.h"

#include "base.h"

#if defined(PARI)

#include "genpari.h"

#endif

#ifndef SPRINT

#define SPRINT

#endif

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

extern int outputstyle;

#ifdef FPRINT

FILE *asir_out;

Line 64 FILE *asir_out;

Line 68 FILE *asir_out;

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 76 int fortran_output;

Line 84 int fortran_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 88 int fortran_output;

Line 96 int fortran_output;

#define PRINTGF2MAT printgf2mat

#define PRINTGFMMAT printgfmmat

#define PRINTBYTEARRAY printbytearray

#define PRINTQUOTE printquote

#define PRINTERR printerr

#define PRINTLF printlf

#define PRINTLOP printlop

Line 95 int fortran_output;

Line 104 int fortran_output;

#define PRINTEOP printeop

#define PRINTQOP printqop

#define PRINTUP printup

#define PRINTUM printum

#define PRINTSF printsf

#define PRINTSYMBOL printsymbol

#endif

#ifdef SPRINT

Line 104 extern char DFORMAT[BUFSIZ];

Line 116 extern char DFORMAT[BUFSIZ];

extern int hex_output;

extern int fortran_output;

extern int double_output;

extern int real_digit;

extern int real_binary;

extern int print_quote;

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

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

#define PRINTF sprintf

Line 116 extern int double_output;

Line 132 extern int double_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 double_output;

Line 144 extern int double_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 double_output;

Line 152 extern int double_output;

#define PRINTEOP sprinteop

#define PRINTQOP sprintqop

#define PRINTUP sprintup

#define PRINTUM sprintum

#define PRINTSF sprintsf

#define PRINTSYMBOL sprintsymbol

#endif

void PRINTEXPR();

Line 153 void PRINTUI();

Line 173 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();

void PRINTSYMBOL();

#ifdef FPRINT

void output_init() {

Line 169 int mmono(p)

Line 198 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 213 P p;

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

}

#if PARI

#if defined(PARI)

void printbf(a)

BF a;

{

sor(a->body,'g',-1,0);

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

}

#endif

Line 195 char *s;

Line 229 char *s;

s[0] = 0; buf = s;

}

#if PARI

#if defined(PARI)

void myoutbrute(g)

GEN g;

{

bruteall(g,'f',-1,1);

}

void sprintbf(a)

BF a;

{

char *str;

char *GENtostr();

char *GENtostr0();

str = GENtostr(a->body);

if ( double_output ) {

str = GENtostr0(a->body,myoutbrute);

} else {

str = GENtostr(a->body);

}

TAIL PRINTF(OUT,"%s",str);

free(str);

}

Line 251 Obj p;

Line 295 Obj p;

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

case O_BYTEARRAY:

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

case O_QUOTE:

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

case O_SYMBOL:

PRINTSYMBOL((Symbol)p); break;

default:

break;

}

Line 287 N n;

Line 335 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 304 Num q;

Line 388 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__) || defined(__x86_64)

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 328 Num q;

Line 588 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 338 C a;

Line 604 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 415 V v;

Line 687 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 479 VECT vect;

Line 751 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 493 MAT mat;

Line 778 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 605 VL vl;

Line 908 VL vl;

GFMMAT mat;

{

int row,col,i,j;

unsigned int t;

unsigned int **b;

row = mat->row;

Line 625 VL vl;

Line 927 VL vl;

BYTEARRAY array;

{

int len,i;

unsigned int t;

unsigned char *b;

len = array->len;

b =array->body;

b = array->body;

PUTS("|");

for ( i = 0; i < len-1; i++ ) {

TAIL PRINTF(OUT,"%02x ",b[i]);

TAIL PRINTF(OUT,"%02x ",(unsigned int)b[i]);

}

TAIL PRINTF(OUT,"%02x",b[i]);

TAIL PRINTF(OUT,"%02x",(unsigned int)b[i]);

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 704 F f;

Line 1019 F f;

}

PRINTFOP(vl,f)

void PRINTFOP(vl,f)

VL vl;

F f;

{

Line 719 F f;

Line 1034 F f;

}

PRINTEOP(vl,f)

void PRINTEOP(vl,f)

VL vl;

F f;

{

Line 746 F f;

Line 1061 F f;

PUTS(")");

}

PRINTLOP(vl,f)

void PRINTLOP(vl,f)

VL vl;

F f;

{

Line 772 F f;

Line 1087 F f;

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

}

PRINTQOP(vl,f)

void PRINTQOP(vl,f)

VL vl;

F f;

{

Line 783 F f;

Line 1098 F f;

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

}

PRINTUP(n)

void PRINTUP(n)

UP n;

{

int i,d;

Line 815 UP n;

Line 1130 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));

}

void PRINTSYMBOL(Symbol p)

{

PUTS(p->name);

}

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