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.3 and 1.23

version 1.3, 2000/08/21 08:31:39

version 1.23, 2003/06/07 16:40:25

Line 23

* shall be made on your publication or presentation in any form of the

* results obtained by use of the SOFTWARE.

* (4) In the event that you modify the SOFTWARE, you shall notify FLL by

* e-mail at risa-admin@flab.fujitsu.co.jp of the detailed specification

* e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification

* for such modification or the source code of the modified part of the

* SOFTWARE.

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.2 2000/02/07 03:14:39 noro Exp $

* $OpenXM: OpenXM_contrib2/asir2000/io/pexpr.c,v 1.22 2003/05/29 16:45:01 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 64 FILE *asir_out;

Line 65 FILE *asir_out;

char DFORMAT[BUFSIZ];

int hex_output;

int fortran_output;

int double_output;

int real_digit;

int print_quote;

#define TAIL

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

Line 76 int fortran_output;

Line 80 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 87 int fortran_output;

Line 91 int fortran_output;

#define PRINTUI printui

#define PRINTGF2MAT printgf2mat

#define PRINTGFMMAT printgfmmat

#define PRINTBYTEARRAY printbytearray

#define PRINTQUOTE printquote

#define PRINTERR printerr

#define PRINTLF printlf

#define PRINTLOP printlop

Line 94 int fortran_output;

Line 100 int fortran_output;

#define PRINTEOP printeop

#define PRINTQOP printqop

#define PRINTUP printup

#define PRINTUM printum

#define PRINTSF printsf

#endif

#ifdef SPRINT

Line 102 static char *buf;

Line 110 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 114 extern int fortran_output;

Line 126 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 125 extern int fortran_output;

Line 137 extern int fortran_output;

#define PRINTUI sprintui

#define PRINTGF2MAT sprintgf2mat

#define PRINTGFMMAT sprintgfmmat

#define PRINTBYTEARRAY sprintbytearray

#define PRINTQUOTE sprintquote

#define PRINTERR sprinterr

#define PRINTLF sprintlf

#define PRINTLOP sprintlop

Line 132 extern int fortran_output;

Line 146 extern int fortran_output;

#define PRINTEOP sprinteop

#define PRINTQOP sprintqop

#define PRINTUP sprintup

#define PRINTUM sprintum

#define PRINTSF sprintsf

#endif

void PRINTEXPR();

Line 149 void PRINTDP();

Line 165 void PRINTDP();

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 165 int mmono(p)

Line 190 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 175 P p;

Line 205 P p;

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

}

#if PARI

#if defined(PARI)

void printbf(a)

BF a;

{

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

void sor();

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

}

#endif

Line 191 char *s;

Line 223 char *s;

s[0] = 0; buf = s;

}

#if PARI

#if defined(PARI)

#include "genpari.h"

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 245 Obj p;

Line 290 Obj p;

PRINTLF(vl,(F)p); break;

case O_GFMMAT:

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

case O_BYTEARRAY:

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

case O_QUOTE:

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

default:

break;

}

Line 281 N n;

Line 330 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 298 Num q;

Line 383 Num q;

}

break;

case N_R:

TAIL PRINTF(OUT,"%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_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 0:

PUTS("[");

PRINTNUM(INF((Itv)q));

PUTS(",");

PRINTNUM(SUP((Itv)q));

PUTS("]");

break;

case 1:

PUTS("intval(");

PRINTNUM(INF((Itv)q));

PUTS(",");

PRINTNUM(SUP((Itv)q));

PUTS(")");

break;

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 0:

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

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

break;

case 1:

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

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

break;

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 0:

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

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

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

break;

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;

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 0:

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

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

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

break;

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;

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 0:

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

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

break;

case 1:

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

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

break;

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

Line 534 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 332 C a;

Line 550 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 409 V v;

Line 633 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 473 VECT vect;

Line 697 VECT vect;

int i;

pointer *ptr;

PUTS("[ ");

switch ( outputstyle ) {

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

case 0:

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

PUTS("[ ");

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

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

}

PUTS("]");

break;

case 1:

PUTS("vect(");

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

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

PRINTEXPR(vl,ptr[i]);

}

PUTS(")");

break;

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 487 MAT mat;

Line 730 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 0:

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

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

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

PUTS("[ ");

}

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

PUTS("]");

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

if ( i < r - 1 )

}

PUTS("\n");

PUTS("]");

if ( i < r - 1 )

PUTS("\n");

}

break;

case 1:

PUTS("mat(\n");

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

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

PUTS("[ ");

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

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

PRINTEXPR(vl,ptr[j]);

}

PUTS("]");

}

PUTS(")");

break;

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 599 VL vl;

Line 870 VL vl;

GFMMAT mat;

{

int row,col,i,j;

unsigned int t;

unsigned int **b;

row = mat->row;

Line 614 GFMMAT mat;

Line 884 GFMMAT mat;

}

void PRINTBYTEARRAY(vl,array)

VL vl;

BYTEARRAY array;

{

int len,i;

unsigned char *b;

len = array->len;

b = array->body;

PUTS("|");

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

TAIL PRINTF(OUT,"%02x ",(unsigned int)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 680 F f;

Line 981 F f;

}

PRINTFOP(vl,f)

void PRINTFOP(vl,f)

VL vl;

F f;

{

Line 695 F f;

Line 996 F f;

}

PRINTEOP(vl,f)

void PRINTEOP(vl,f)

VL vl;

F f;

{

Line 722 F f;

Line 1023 F f;

PUTS(")");

}

PRINTLOP(vl,f)

void PRINTLOP(vl,f)

VL vl;

F f;

{

Line 748 F f;

Line 1049 F f;

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

}

PRINTQOP(vl,f)

void PRINTQOP(vl,f)

VL vl;

F f;

{

Line 759 F f;

Line 1060 F f;

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

}

PRINTUP(n)

void PRINTUP(n)

UP n;

{

int i,d;

Line 790 UP n;

Line 1091 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>