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.21 and 1.29

version 1.21, 2003/02/14 22:29:15

version 1.29, 2004/02/09 08:23:30

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.20 2002/09/11 07:23:27 noro Exp $

* $OpenXM: OpenXM_contrib2/asir2000/io/pexpr.c,v 1.28 2004/02/04 07:42:07 noro Exp $

#include "ca.h"

#include "al.h"

Line 52

#include "comp.h"

#include "base.h"

#if defined(PARI)

#include "genpari.h"

#endif

#ifndef FPRINT

#define FPRINT

#endif

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

extern int outputstyle;

#ifdef FPRINT

FILE *asir_out;

Line 66 int hex_output;

Line 71 int hex_output;

int fortran_output;

int double_output;

int real_digit;

int real_binary;

int print_quote;

#define TAIL

Line 101 int print_quote;

Line 107 int print_quote;

#define PRINTUP printup

#define PRINTUM printum

#define PRINTSF printsf

#define PRINTSYMBOL printsymbol

#define PRINTRANGE printrange

#endif

#ifdef SPRINT

Line 111 extern int hex_output;

Line 119 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 146 extern int print_quote;

Line 156 extern int print_quote;

#define PRINTUP sprintup

#define PRINTUM sprintum

#define PRINTSF sprintsf

#define PRINTSYMBOL sprintsymbol

#define PRINTRANGE sprintrange

#endif

void PRINTEXPR();

Line 177 void PRINTEOP();

Line 189 void PRINTEOP();

void PRINTLOP();

void PRINTQOP();

void PRINTSF();

void PRINTSYMBOL();

void PRINTRANGE();

#ifdef FPRINT

void output_init() {

Line 207 P p;

Line 221 P p;

void printbf(a)

BF a;

{

void sor();

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

}

#endif

Line 222 char *s;

Line 234 char *s;

}

#if defined(PARI)

#include "genpari.h"

void myoutbrute(g)

GEN g;

{

Line 234 void sprintbf(a)

Line 244 void sprintbf(a)

BF a;

{

char *str;

char *GENtostr();

char *GENtostr0();

if ( double_output ) {

Line 292 Obj p;

Line 301 Obj p;

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

case O_QUOTE:

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

case O_SYMBOL:

PRINTSYMBOL((SYMBOL)p); break;

case O_RANGE:

PRINTRANGE(vl,(RANGE)p); break;

default:

break;

}

Line 393 Num q;

Line 406 Num q;

case MID_PRINTF_G:

#endif

default:

if ( real_digit ) {

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

Line 413 Num q;

Line 495 Num q;

#if defined(INTERVAL)

case N_IP:

case N_IntervalBigFloat:

PUTS("[");

switch ( outputstyle ) {

PRINTNUM(INF((Itv)q));

case 1:

PUTS(",");

PUTS("intval(");

PRINTNUM(SUP((Itv)q));

PRINTNUM(INF((Itv)q));

PUTS("]");

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:

TAIL PRINTF(OUT, "[%.16e,%.16e]",INF((IntervalDouble)q),SUP((IntervalDouble)q));

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:

TAIL PRINTF(OUT, "<%g,%g>", (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);

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:

TAIL PRINTF(OUT, "<%.16e,%.16e>", (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);

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:

TAIL PRINTF(OUT, "[%g,%g]",INF((IntervalDouble)q),SUP((IntervalDouble)q));

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;

Line 619 VECT vect;

Line 757 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 633 MAT mat;

Line 784 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 1011 unsigned int i;

Line 1180 unsigned int i;

} else {

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

}

void PRINTSYMBOL(SYMBOL sym)

{

PUTS(sym->name);

}

void PRINTRANGE(VL vl,RANGE p)

{

PUTS("range(");

PRINTEXPR(vl,p->start);

PUTS(",");

PRINTEXPR(vl,p->end);

PUTS(")");

}

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