version 1.7, 2000/12/15 05:30:08 |
version 1.14, 2001/08/06 01:48:33 |
|
|
* OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY |
* OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* $OpenXM: OpenXM_contrib2/asir2000/io/spexpr.c,v 1.6 2000/12/13 10:54:09 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/io/spexpr.c,v 1.13 2001/04/20 02:34:23 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "al.h" |
#include "al.h" |
|
|
char DFORMAT[BUFSIZ]; |
char DFORMAT[BUFSIZ]; |
int hex_output; |
int hex_output; |
int fortran_output; |
int fortran_output; |
|
int double_output; |
|
int real_digit; |
|
int print_quote; |
|
|
#define TAIL |
#define TAIL |
#define PUTS(s) fputs(s,OUT) |
#define PUTS(s) fputs(s,OUT) |
Line 88 int fortran_output; |
|
Line 91 int fortran_output; |
|
#define PRINTGF2MAT printgf2mat |
#define PRINTGF2MAT printgf2mat |
#define PRINTGFMMAT printgfmmat |
#define PRINTGFMMAT printgfmmat |
#define PRINTBYTEARRAY printbytearray |
#define PRINTBYTEARRAY printbytearray |
|
#define PRINTQUOTE printquote |
#define PRINTERR printerr |
#define PRINTERR printerr |
#define PRINTLF printlf |
#define PRINTLF printlf |
#define PRINTLOP printlop |
#define PRINTLOP printlop |
Line 104 extern char DFORMAT[BUFSIZ]; |
|
Line 108 extern char DFORMAT[BUFSIZ]; |
|
extern int hex_output; |
extern int hex_output; |
extern int fortran_output; |
extern int fortran_output; |
extern int double_output; |
extern int double_output; |
|
extern int real_digit; |
|
extern int print_quote; |
|
|
#define TAIL while ( *OUT ) OUT++; |
#define TAIL while ( *OUT ) OUT++; |
#define PUTS(s) strcat(OUT,s) |
#define PUTS(s) strcat(OUT,s) |
Line 128 extern int double_output; |
|
Line 134 extern int double_output; |
|
#define PRINTGF2MAT sprintgf2mat |
#define PRINTGF2MAT sprintgf2mat |
#define PRINTGFMMAT sprintgfmmat |
#define PRINTGFMMAT sprintgfmmat |
#define PRINTBYTEARRAY sprintbytearray |
#define PRINTBYTEARRAY sprintbytearray |
|
#define PRINTQUOTE sprintquote |
#define PRINTERR sprinterr |
#define PRINTERR sprinterr |
#define PRINTLF sprintlf |
#define PRINTLF sprintlf |
#define PRINTLOP sprintlop |
#define PRINTLOP sprintlop |
|
|
void PRINTGF2MAT(); |
void PRINTGF2MAT(); |
void PRINTGFMMAT(); |
void PRINTGFMMAT(); |
void PRINTBYTEARRAY(); |
void PRINTBYTEARRAY(); |
|
void PRINTQUOTE(); |
void PRINTERR(); |
void PRINTERR(); |
void PRINTCPLX(); |
void PRINTCPLX(); |
void PRINTLM(); |
void PRINTLM(); |
|
|
P p; |
P p; |
{ |
{ |
if ( NUM(p) ) |
if ( NUM(p) ) |
|
#if defined(INTERVAL) |
|
if ( NID(p) != N_IP && NID(p) != N_ID && NID(p) != N_IT && NID(p) != N_IF |
|
&& compnum(CO,(Num)p,0) < 0 ) |
|
#else |
if ( compnum(CO,(Num)p,0) < 0 ) |
if ( compnum(CO,(Num)p,0) < 0 ) |
|
#endif |
return ( 1 ); |
return ( 1 ); |
else |
else |
return ( 0 ); |
return ( 0 ); |
|
|
} |
} |
|
|
#if PARI |
#if PARI |
|
#include "genpari.h" |
|
|
|
void myoutbrute(g) |
|
GEN g; |
|
{ |
|
bruteall(g,'f',-1,1); |
|
} |
|
|
void sprintbf(a) |
void sprintbf(a) |
BF a; |
BF a; |
{ |
{ |
char *str; |
char *str; |
char *GENtostr(); |
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); |
TAIL PRINTF(OUT,"%s",str); |
free(str); |
free(str); |
} |
} |
|
|
PRINTGFMMAT(vl,(GFMMAT)p); break; |
PRINTGFMMAT(vl,(GFMMAT)p); break; |
case O_BYTEARRAY: |
case O_BYTEARRAY: |
PRINTBYTEARRAY(vl,(BYTEARRAY)p); break; |
PRINTBYTEARRAY(vl,(BYTEARRAY)p); break; |
|
case O_QUOTE: |
|
PRINTQUOTE(vl,(QUOTE)p); break; |
default: |
default: |
break; |
break; |
} |
} |
|
|
} |
} |
} |
} |
|
|
|
extern int printmode; |
|
|
void PRINTNUM(q) |
void PRINTNUM(q) |
Num q; |
Num q; |
{ |
{ |
|
char real_format[20]; |
|
|
if ( !q ) { |
if ( !q ) { |
PUTS("0"); |
PUTS("0"); |
return; |
return; |
|
|
} |
} |
break; |
break; |
case N_R: |
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_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; |
break; |
case N_A: |
case N_A: |
PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")"); |
PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")"); |
|
|
case N_B: |
case N_B: |
PRINTBF((BF)q); break; |
PRINTBF((BF)q); break; |
#endif |
#endif |
|
#if defined(INTERVAL) |
|
case N_IP: |
|
case N_IF: |
|
PUTS("["); |
|
PRINTNUM(INF((Itv)q)); |
|
PUTS(","); |
|
PRINTNUM(SUP((Itv)q)); |
|
PUTS("]"); |
|
break; |
|
case N_ID: |
|
switch (printmode) { |
|
case PRINTF_E: |
|
TAIL PRINTF(OUT, "[%.16e,%.16e]",INF((ItvD)q),SUP((ItvD)q)); |
|
#if defined(ITVDEBUG) |
|
printbin(INF((ItvD)q)); |
|
printbin(SUP((ItvD)q)); |
|
#endif |
|
break; |
|
case MID_PRINTF_G: |
|
TAIL PRINTF(OUT, "<%g,%g>", (SUP((ItvD)q)+INF((ItvD)q))*0.5,(SUP((ItvD)q)-INF((ItvD)q))*0.5); |
|
break; |
|
case MID_PRINTF_E: |
|
TAIL PRINTF(OUT, "<%.16e,%.16e>", (SUP((ItvD)q)+INF((ItvD)q))*0.5,(SUP((ItvD)q)-INF((ItvD)q))*0.5); |
|
break; |
|
case PRINTF_G: |
|
default: |
|
TAIL PRINTF(OUT, "[%g,%g]",INF((ItvD)q),SUP((ItvD)q)); |
|
break; |
|
} |
|
break; |
|
#endif |
case N_C: |
case N_C: |
PRINTCPLX((C)q); break; |
PRINTCPLX((C)q); break; |
case N_M: |
case N_M: |
|
|
case N_GFPN: |
case N_GFPN: |
PRINTUP((UP)(((GFPN)q)->body)); |
PRINTUP((UP)(((GFPN)q)->body)); |
break; |
break; |
|
case N_GFS: |
|
TAIL PRINTF(OUT,"@_%d",CONT((GFS)q)); |
|
break; |
} |
} |
} |
} |
|
|
|
|
if ( a->r ) |
if ( a->r ) |
PRINTNUM(a->r); |
PRINTNUM(a->r); |
if ( a->i ) { |
if ( a->i ) { |
|
#if defined(INTERVAL) |
|
if ( a->r && ((compnum(0,a->i,0) > 0) |
|
|| NID(a->i) == N_IP || NID(a->i) == N_ID |
|
|| NID(a->i) == N_IT || NID(a->i) == N_IF) ) |
|
#else |
if ( a->r && (compnum(0,a->i,0) > 0) ) |
if ( a->r && (compnum(0,a->i,0) > 0) ) |
|
#endif |
PUTS("+"); |
PUTS("+"); |
PRINTNUM(a->i); PUTS("*@i"); |
PRINTNUM(a->i); PUTS("*@i"); |
} |
} |
Line 629 BYTEARRAY array; |
|
Line 721 BYTEARRAY array; |
|
unsigned char *b; |
unsigned char *b; |
|
|
len = array->len; |
len = array->len; |
b =array->body; |
b = array->body; |
PUTS("|"); |
PUTS("|"); |
for ( i = 0; i < len-1; i++ ) { |
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("|"); |
PUTS("|"); |
|
} |
|
|
|
void PRINTQUOTE(vl,quote) |
|
VL vl; |
|
QUOTE quote; |
|
{ |
|
PUTS("<...quoted...>"); |
} |
} |
|
|
void PRINTERR(vl,e) |
void PRINTERR(vl,e) |