version 1.1, 1999/12/03 07:39:11 |
version 1.14, 2001/08/06 01:48:33 |
|
|
/* $OpenXM: OpenXM/src/asir99/io/pexpr.c,v 1.1.1.1 1999/11/10 08:12:30 noro Exp $ */ |
/* |
|
* Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED |
|
* All rights reserved. |
|
* |
|
* FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited, |
|
* non-exclusive and royalty-free license to use, copy, modify and |
|
* redistribute, solely for non-commercial and non-profit purposes, the |
|
* computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and |
|
* conditions of this Agreement. For the avoidance of doubt, you acquire |
|
* only a limited right to use the SOFTWARE hereunder, and FLL or any |
|
* third party developer retains all rights, including but not limited to |
|
* copyrights, in and to the SOFTWARE. |
|
* |
|
* (1) FLL does not grant you a license in any way for commercial |
|
* purposes. You may use the SOFTWARE only for non-commercial and |
|
* non-profit purposes only, such as academic, research and internal |
|
* business use. |
|
* (2) The SOFTWARE is protected by the Copyright Law of Japan and |
|
* international copyright treaties. If you make copies of the SOFTWARE, |
|
* with or without modification, as permitted hereunder, you shall affix |
|
* to all such copies of the SOFTWARE the above copyright notice. |
|
* (3) An explicit reference to this SOFTWARE and its copyright owner |
|
* 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@sec.flab.fujitsu.co.jp of the detailed specification |
|
* for such modification or the source code of the modified part of the |
|
* SOFTWARE. |
|
* |
|
* THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL |
|
* MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND |
|
* EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS |
|
* FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES' |
|
* RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY |
|
* MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY. |
|
* UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT, |
|
* OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY |
|
* DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL |
|
* DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES |
|
* ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES |
|
* FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY |
|
* DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF |
|
* SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART |
|
* 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.13 2001/04/20 02:34:23 noro Exp $ |
|
*/ |
#include "ca.h" |
#include "ca.h" |
#include "al.h" |
#include "al.h" |
#include "parse.h" |
#include "parse.h" |
|
|
#define FPRINT |
#define FPRINT |
#endif |
#endif |
|
|
int hex_output; |
#define PRINTHAT (fortran_output?PUTS("**"):PUTS("^")) |
|
|
#ifdef FPRINT |
#ifdef FPRINT |
FILE *asir_out; |
FILE *asir_out; |
#define OUT asir_out |
#define OUT asir_out |
char DFORMAT[BUFSIZ]; |
char DFORMAT[BUFSIZ]; |
|
int hex_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 38 char DFORMAT[BUFSIZ]; |
|
Line 90 char DFORMAT[BUFSIZ]; |
|
#define PRINTUI printui |
#define PRINTUI printui |
#define PRINTGF2MAT printgf2mat |
#define PRINTGF2MAT printgf2mat |
#define PRINTGFMMAT printgfmmat |
#define PRINTGFMMAT printgfmmat |
|
#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 51 char DFORMAT[BUFSIZ]; |
|
Line 105 char DFORMAT[BUFSIZ]; |
|
static char *buf; |
static char *buf; |
#define OUT buf |
#define OUT buf |
extern char DFORMAT[BUFSIZ]; |
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 TAIL while ( *OUT ) OUT++; |
#define PUTS(s) strcat(OUT,s) |
#define PUTS(s) strcat(OUT,s) |
Line 74 extern char DFORMAT[BUFSIZ]; |
|
Line 133 extern char DFORMAT[BUFSIZ]; |
|
#define PRINTUI sprintui |
#define PRINTUI sprintui |
#define PRINTGF2MAT sprintgf2mat |
#define PRINTGF2MAT sprintgf2mat |
#define PRINTGFMMAT sprintgfmmat |
#define PRINTGFMMAT sprintgfmmat |
|
#define PRINTBYTEARRAY sprintbytearray |
|
#define PRINTQUOTE sprintquote |
#define PRINTERR sprinterr |
#define PRINTERR sprinterr |
#define PRINTLF sprintlf |
#define PRINTLF sprintlf |
#define PRINTLOP sprintlop |
#define PRINTLOP sprintlop |
#define PRINTFOP sprintfop |
#define PRINTFOP sprintfop |
#define PRINTEOP sprinteop |
#define PRINTEOP sprinteop |
#define PRINTQOP sprintqop |
#define PRINTQOP sprintqop |
|
#define PRINTUP sprintup |
#endif |
#endif |
|
|
#if defined(THINK_C) |
|
void PRINTEXPR(VL,Obj); |
|
void PRINTNUM(Num); |
|
void PRINTN(N); |
|
void PRINTV(VL,V); |
|
void PRINTP(VL,P); |
|
void PRINTR(VL,R); |
|
void PRINTLIST(VL,LIST); |
|
void PRINTVECT(VL,VECT); |
|
void PRINTMAT(VL,MAT); |
|
void PRINTSTR(STRING); |
|
void PRINTCOMP(VL,COMP); |
|
void PRINTDP(VL,DP); |
|
void PRINTUI(VL,USINT); |
|
void PRINTGF2MAT(VL,GF2MAT); |
|
void PRINTGFMMAT(VL,GFMMAT); |
|
void PRINTERR(VL,ERR); |
|
void PRINTCPLX(C); |
|
void PRINTLM(LM); |
|
void PRINTLF(VL,F); |
|
#else |
|
void PRINTEXPR(); |
void PRINTEXPR(); |
void PRINTNUM(); |
void PRINTNUM(); |
void PRINTN(); |
void PRINTN(); |
|
|
void PRINTUI(); |
void PRINTUI(); |
void PRINTGF2MAT(); |
void PRINTGF2MAT(); |
void PRINTGFMMAT(); |
void PRINTGFMMAT(); |
|
void PRINTBYTEARRAY(); |
|
void PRINTQUOTE(); |
void PRINTERR(); |
void PRINTERR(); |
void PRINTCPLX(); |
void PRINTCPLX(); |
void PRINTLM(); |
void PRINTLM(); |
void PRINTLF(); |
void PRINTLF(); |
#endif |
void PRINTUP2(); |
|
|
#ifdef FPRINT |
#ifdef FPRINT |
void output_init() { |
void output_init() { |
|
|
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 |
|
|
#if defined(THINK_C) |
|
void sor(long *,char,long,long); |
|
#endif |
|
|
|
void printbf(a) |
void printbf(a) |
BF a; |
BF a; |
{ |
{ |
sor(a->body,'g',-1,0); |
sor(a->body,double_output ? 'f' : 'g',-1,0); |
} |
} |
#endif |
#endif |
#endif |
#endif |
|
|
} |
} |
|
|
#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; |
{ |
{ |
/* XXX : not implemented yet */ |
char *str; |
|
char *GENtostr(); |
|
char *GENtostr0(); |
|
|
|
if ( double_output ) { |
|
str = GENtostr0(a->body,myoutbrute); |
|
} else { |
|
str = GENtostr(a->body); |
|
} |
|
TAIL PRINTF(OUT,"%s",str); |
|
free(str); |
} |
} |
#endif |
#endif |
#endif |
#endif |
|
|
PRINTLF(vl,(F)p); break; |
PRINTLF(vl,(F)p); break; |
case O_GFMMAT: |
case O_GFMMAT: |
PRINTGFMMAT(vl,(GFMMAT)p); break; |
PRINTGFMMAT(vl,(GFMMAT)p); break; |
|
case O_BYTEARRAY: |
|
PRINTBYTEARRAY(vl,(BYTEARRAY)p); break; |
|
case O_QUOTE: |
|
PRINTQUOTE(vl,(QUOTE)p); break; |
default: |
default: |
break; |
break; |
} |
} |
|
|
} |
} |
} |
} |
|
|
|
#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 0 |
|
int printmode = PRINTF_E; |
|
#else |
|
int printmode = PRINTF_G; |
|
#endif |
|
|
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,"%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"); |
} |
} |
|
|
} |
} |
PRINTV(vl,v); |
PRINTV(vl,v); |
if ( cmpq(DEG(dc),ONE) ) { |
if ( cmpq(DEG(dc),ONE) ) { |
PUTS("^"); |
PRINTHAT; |
if ( INT(DEG(dc)) && SGN(DEG(dc))>0 ) |
if ( INT(DEG(dc)) && SGN(DEG(dc))>0 ) |
PRINTNUM((Num)DEG(dc)); |
PRINTNUM((Num)DEG(dc)); |
else { |
else { |
|
|
else if ( (vid)v->attr == V_PF ) { |
else if ( (vid)v->attr == V_PF ) { |
pf = ((PFINS)v->priv)->pf; ad = ((PFINS)v->priv)->ad; |
pf = ((PFINS)v->priv)->pf; ad = ((PFINS)v->priv)->ad; |
if ( !strcmp(NAME(pf),"pow") ) { |
if ( !strcmp(NAME(pf),"pow") ) { |
PUTS("("); PRINTR(vl,(R)ad[0].arg); 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 ) { |
} else if ( !pf->argc ) { |
TAIL PRINTF(OUT,"%s",NAME(pf)); |
TAIL PRINTF(OUT,"%s",NAME(pf)); |
|
|
} |
} |
} |
} |
|
|
|
void PRINTBYTEARRAY(vl,array) |
|
VL vl; |
|
BYTEARRAY array; |
|
{ |
|
int len,i; |
|
unsigned int t; |
|
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) |
void PRINTERR(vl,e) |
VL vl; |
VL vl; |
ERR e; |
ERR e; |
|
|
{ |
{ |
int d,i; |
int d,i; |
|
|
if ( !p ) |
if ( !p ) { |
fprintf(asir_out,"0"); |
TAIL PRINTF(OUT,"0"); |
else { |
} else { |
d = degup2(p); |
d = degup2(p); |
fprintf(asir_out,"("); |
TAIL PRINTF(OUT,"("); |
if ( !d ) |
if ( !d ) { |
fprintf(asir_out,"1"); |
TAIL PRINTF(OUT,"1"); |
else if ( d == 1 ) |
} else if ( d == 1 ) { |
fprintf(asir_out,"@"); |
TAIL PRINTF(OUT,"@"); |
else |
} else { |
fprintf(asir_out,"@^%d",d); |
PUTS("@"); PRINTHAT; TAIL PRINTF(OUT,"%d",d); |
|
} |
for ( i = d-1; i >= 0; i-- ) { |
for ( i = d-1; i >= 0; i-- ) { |
if ( p->b[i/BSH] & (1<<(i%BSH)) ) |
if ( p->b[i/BSH] & (1<<(i%BSH)) ) |
if ( !i ) |
if ( !i ) { |
fprintf(asir_out,"+1"); |
TAIL PRINTF(OUT,"+1"); |
else if ( i == 1 ) |
} else if ( i == 1 ) { |
fprintf(asir_out,"+@"); |
TAIL PRINTF(OUT,"+@"); |
else |
} else { |
fprintf(asir_out,"+@^%d",i); |
PUTS("+@"); PRINTHAT; TAIL PRINTF(OUT,"%d",i); |
|
} |
} |
} |
fprintf(asir_out,")"); |
TAIL PRINTF(OUT,")"); |
} |
} |
} |
} |
|
|
|
|
PUTS("*@p"); |
PUTS("*@p"); |
} else { |
} else { |
PRINTNUM(n->c[d]); |
PRINTNUM(n->c[d]); |
TAIL PRINTF(OUT,"*@p^%d",d); |
PUTS("*@p"); PRINTHAT; TAIL PRINTF(OUT,"%d",d); |
} |
} |
for ( i = d-1; i >= 0; i-- ) { |
for ( i = d-1; i >= 0; i-- ) { |
if ( n->c[i] ) { |
if ( n->c[i] ) { |
PUTS("+("); PRINTNUM(n->c[i]); PUTS(")"); |
PUTS("+("); PRINTNUM(n->c[i]); PUTS(")"); |
if ( i >= 2 ) { |
if ( i >= 2 ) { |
TAIL PRINTF(OUT,"*@p^%d",i); |
PUTS("*@p"); PRINTHAT; TAIL PRINTF(OUT,"%d",i); |
} else if ( i == 1 ) |
} else if ( i == 1 ) |
PUTS("*@p"); |
PUTS("*@p"); |
} |
} |