| version 1.3, 2000/08/21 08:31:38 |
version 1.16, 2003/12/25 02:40:24 |
|
|
| * shall be made on your publication or presentation in any form of the |
* shall be made on your publication or presentation in any form of the |
| * results obtained by use of the SOFTWARE. |
* results obtained by use of the SOFTWARE. |
| * (4) In the event that you modify the SOFTWARE, you shall notify FLL by |
* (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 |
* for such modification or the source code of the modified part of the |
| * SOFTWARE. |
* SOFTWARE. |
| * |
* |
|
|
| * 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/cpexpr.c,v 1.2 2000/02/07 03:14:39 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/io/cpexpr.c,v 1.15 2003/02/14 22:29:15 ohara Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| #include "parse.h" |
#include "parse.h" |
| #include "al.h" |
#include "al.h" |
| #include "base.h" |
#include "base.h" |
| |
#if PARI |
| |
#include "genpari.h" |
| |
#endif |
| |
|
| extern int hex_output,fortran_output; |
extern int hex_output,fortran_output,double_output,real_digit; |
| |
|
| #define PRINTHAT (fortran_output?PUTS("**"):PUTS("^")) |
#define PRINTHAT (fortran_output?PUTS("**"):PUTS("^")) |
| |
|
| Line 76 extern int hex_output,fortran_output; |
|
| Line 79 extern int hex_output,fortran_output; |
|
| #define PRINTUI length_ui |
#define PRINTUI length_ui |
| #define PRINTGF2MAT length_gf2mat |
#define PRINTGF2MAT length_gf2mat |
| #define PRINTGFMMAT length_gfmmat |
#define PRINTGFMMAT length_gfmmat |
| |
#define PRINTBYTEARRAY length_bytearray |
| |
#define PRINTQUOTE length_QUOTE |
| #define PRINTERR length_err |
#define PRINTERR length_err |
| #define PRINTLF length_lf |
#define PRINTLF length_lf |
| #define PRINTLOP length_lop |
#define PRINTLOP length_lop |
| Line 83 extern int hex_output,fortran_output; |
|
| Line 88 extern int hex_output,fortran_output; |
|
| #define PRINTEOP length_eop |
#define PRINTEOP length_eop |
| #define PRINTQOP length_qop |
#define PRINTQOP length_qop |
| #define PRINTUP length_up |
#define PRINTUP length_up |
| |
#define PRINTUM length_um |
| |
#define PRINTSF length_sf |
| |
|
| void PRINTEXPR(); |
void PRINTEXPR(); |
| void PRINTNUM(); |
void PRINTNUM(); |
|
|
| 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(); |
| |
void PRINTUP(); |
| |
void PRINTUM(); |
| void PRINTUP2(); |
void PRINTUP2(); |
| |
void PRINTFOP(); |
| |
void PRINTEOP(); |
| |
void PRINTLOP(); |
| |
void PRINTQOP(); |
| |
void PRINTSF(); |
| |
|
| static int total_length; |
static int total_length; |
| |
|
| int estimate_length(vl,p) |
int estimate_length(VL vl,pointer p) |
| VL vl; |
|
| pointer p; |
|
| { |
{ |
| total_length = 0; |
total_length = 0; |
| PRINTEXPR(vl,p); |
PRINTEXPR(vl,p); |
| return total_length; |
return total_length; |
| } |
} |
| |
|
| #if PARI |
#if defined(PARI) |
| void PRINTBF(a) |
void PRINTBF(BF a) |
| BF a; |
|
| { |
{ |
| char *str; |
char *str; |
| char *GENtostr(); |
char *GENtostr0(); |
| |
void myoutbrute(); |
| |
|
| str = GENtostr(a->body); |
if ( double_output ) { |
| |
str = GENtostr0(a->body,myoutbrute); |
| |
} else { |
| |
str = GENtostr(a->body); |
| |
} |
| total_length += strlen(str); |
total_length += strlen(str); |
| free(str); |
free(str); |
| } |
} |
| #endif |
#endif |
| |
|
| void PRINTEXPR(vl,p) |
void PRINTEXPR(VL vl,pointer p) |
| VL vl; |
|
| pointer p; |
|
| { |
{ |
| if ( !p ) { |
if ( !p ) { |
| total_length++; |
total_length++; |
|
|
| 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; |
| } |
} |
| } |
} |
| |
|
| void PRINTN(n) |
void PRINTN(N n) |
| N n; |
|
| { |
{ |
| double ceil(); |
double ceil(); |
| |
|
|
|
| total_length += (int)(ceil(0.31*((double)(BSH*PL(n))))+1); |
total_length += (int)(ceil(0.31*((double)(BSH*PL(n))))+1); |
| } |
} |
| |
|
| void PRINTNUM(q) |
void PRINTNUM(Num q) |
| Num q; |
|
| { |
{ |
| if ( !q ) { |
if ( !q ) { |
| PUTS("0"); |
PUTS("0"); |
|
|
| } |
} |
| break; |
break; |
| case N_R: |
case N_R: |
| total_length += 20; /* XXX */ |
if ( double_output ) |
| |
total_length += 400+real_digit; /* XXX */ |
| |
else |
| |
total_length += 20+real_digit; /* XXX */ |
| break; |
break; |
| case N_A: |
case N_A: |
| PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")"); |
PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")"); |
| break; |
break; |
| #if PARI |
#if defined(PARI) |
| case N_B: |
case N_B: |
| PRINTBF((BF)q); break; |
PRINTBF((BF)q); break; |
| #endif |
#endif |
|
|
| case N_GFPN: |
case N_GFPN: |
| PRINTUP((UP)(((GFPN)q)->body)); |
PRINTUP((UP)(((GFPN)q)->body)); |
| break; |
break; |
| |
case N_GFS: |
| |
total_length += 13; /* XXX */ |
| |
break; |
| |
case N_GFSN: |
| |
PRINTUM(BDY((GFSN)q)); |
| |
break; |
| default: |
default: |
| break; |
break; |
| } |
} |
| } |
} |
| |
|
| void PRINTCPLX(a) |
void PRINTCPLX(C a) |
| C a; |
|
| { |
{ |
| PUTS("("); |
PUTS("("); |
| if ( a->r ) |
if ( a->r ) |
|
|
| PUTS(")"); |
PUTS(")"); |
| } |
} |
| |
|
| void PRINTP(vl,p) |
void PRINTP(VL vl,P p) |
| VL vl; |
|
| P p; |
|
| { |
{ |
| V v; |
V v; |
| DCP dc; |
DCP dc; |
| int t; |
|
| |
|
| if ( !p ) |
if ( !p ) |
| PUTS("0"); |
PUTS("0"); |
|
|
| |
|
| extern int hideargs; |
extern int hideargs; |
| |
|
| void PRINTV(vl,v) |
void PRINTV(VL vl,V v) |
| VL vl; |
|
| V v; |
|
| { |
{ |
| PF pf; |
PF pf; |
| PFAD ad; |
PFAD ad; |
|
|
| 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(")"); 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 ) |
} else if ( !pf->argc ) |
| PUTS(NAME(pf)); |
PUTS(NAME(pf)); |
| else { |
else { |
|
|
| } |
} |
| } |
} |
| |
|
| void PRINTR(vl,a) |
void PRINTR(VL vl,R a) |
| VL vl; |
|
| R a; |
|
| { |
{ |
| if ( !a ) |
if ( !a ) |
| PUTS("0"); |
PUTS("0"); |
|
|
| } |
} |
| } |
} |
| |
|
| void PRINTVECT(vl,vect) |
void PRINTVECT(VL vl,VECT vect) |
| VL vl; |
|
| VECT vect; |
|
| { |
{ |
| int i; |
int i; |
| pointer *ptr; |
pointer *ptr; |
|
|
| PUTS("]"); |
PUTS("]"); |
| } |
} |
| |
|
| void PRINTMAT(vl,mat) |
void PRINTMAT(VL vl,MAT mat) |
| VL vl; |
|
| MAT mat; |
|
| { |
{ |
| int i,j,r,c; |
int i,j,r,c; |
| pointer *ptr; |
pointer *ptr; |
|
|
| } |
} |
| } |
} |
| |
|
| void PRINTLIST(vl,list) |
void PRINTLIST(VL vl,LIST list) |
| VL vl; |
|
| LIST list; |
|
| { |
{ |
| NODE tnode; |
NODE tnode; |
| |
|
|
|
| PUTS("]"); |
PUTS("]"); |
| } |
} |
| |
|
| void PRINTSTR(str) |
void PRINTSTR(STRING str) |
| STRING str; |
|
| { |
{ |
| char *p; |
char *p; |
| |
|
|
|
| } |
} |
| } |
} |
| |
|
| void PRINTCOMP(vl,c) |
void PRINTCOMP(VL vl,COMP c) |
| VL vl; |
|
| COMP c; |
|
| { |
{ |
| int n,i; |
int n,i; |
| |
|
|
|
| PUTS("}"); |
PUTS("}"); |
| } |
} |
| |
|
| void PRINTDP(vl,d) |
void PRINTDP(VL vl,DP d) |
| VL vl; |
|
| DP d; |
|
| { |
{ |
| int n,i; |
int n,i; |
| MP m; |
MP m; |
|
|
| } |
} |
| } |
} |
| |
|
| void PRINTUI(vl,u) |
void PRINTUI(VL vl,USINT u) |
| VL vl; |
|
| USINT u; |
|
| { |
{ |
| total_length += 10; |
total_length += 10; |
| } |
} |
| |
|
| void PRINTGF2MAT(vl,mat) |
void PRINTGF2MAT(VL vl,GF2MAT mat) |
| VL vl; |
|
| GF2MAT mat; |
|
| { |
{ |
| int row,col,w,i,j,k,m; |
int row,col,w,i,j,k,m; |
| unsigned int t; |
unsigned int t; |
|
|
| } |
} |
| } |
} |
| |
|
| void PRINTGFMMAT(vl,mat) |
void PRINTGFMMAT(VL vl,GFMMAT mat) |
| VL vl; |
|
| GFMMAT mat; |
|
| { |
{ |
| int row,col,i,j; |
int row,col,i,j; |
| unsigned int t; |
|
| unsigned int **b; |
unsigned int **b; |
| |
|
| row = mat->row; |
row = mat->row; |
|
|
| } |
} |
| } |
} |
| |
|
| void PRINTERR(vl,e) |
void PRINTBYTEARRAY(VL vl,BYTEARRAY array) |
| VL vl; |
|
| ERR e; |
|
| { |
{ |
| |
/* |xx xx ... xx| */ |
| |
total_length += 1+3*array->len; |
| |
} |
| |
|
| |
void PRINTQUOTE(VL vl,QUOTE quote) |
| |
{ |
| |
/* <...quoted...> */ |
| |
total_length += 20; |
| |
} |
| |
|
| |
void PRINTERR(VL vl,ERR e) |
| |
{ |
| PUTS("error("); PRINTEXPR(vl,e->body); PUTS(")"); |
PUTS("error("); PRINTEXPR(vl,e->body); PUTS(")"); |
| } |
} |
| |
|
| void PRINTUP2(p) |
void PRINTUP2(UP2 p) |
| UP2 p; |
|
| { |
{ |
| int d,i; |
int d,i; |
| |
|
|
|
| } |
} |
| } |
} |
| |
|
| void PRINTLF(vl,f) |
void PRINTLF(VL vl,F f) |
| VL vl; |
|
| F f; |
|
| { |
{ |
| switch ( FOP(f) ) { |
switch ( FOP(f) ) { |
| case AL_TRUE: |
case AL_TRUE: |
|
|
| } |
} |
| } |
} |
| |
|
| PRINTFOP(vl,f) |
void PRINTFOP(VL vl,F f) |
| VL vl; |
|
| F f; |
|
| { |
{ |
| char *op; |
char *op; |
| NODE n; |
NODE n; |
|
|
| } |
} |
| } |
} |
| |
|
| PRINTEOP(vl,f) |
void PRINTEOP(VL vl,F f) |
| VL vl; |
|
| F f; |
|
| { |
{ |
| oFOP op; |
oFOP op; |
| char *sop; |
char *sop; |
|
|
| PUTS(")"); |
PUTS(")"); |
| } |
} |
| |
|
| PRINTLOP(vl,f) |
void PRINTLOP(VL vl,F f) |
| VL vl; |
|
| F f; |
|
| { |
{ |
| char *op; |
char *op; |
| |
|
|
|
| PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0"); |
PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0"); |
| } |
} |
| |
|
| PRINTQOP(vl,f) |
void PRINTQOP(VL vl,F f) |
| VL vl; |
|
| F f; |
|
| { |
{ |
| char *op; |
char *op; |
| |
|
|
|
| PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")"); |
PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")"); |
| } |
} |
| |
|
| PRINTUP(n) |
void PRINTUP(UP n) |
| UP n; |
|
| { |
{ |
| int i,d; |
int i,d; |
| |
|
|
|
| } |
} |
| PUTS(")"); |
PUTS(")"); |
| } |
} |
| |
} |
| |
|
| |
void PRINTUM(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; total_length += 13; |
| |
} |
| |
for ( i = d-1; i >= 0; i-- ) { |
| |
if ( n->c[i] ) { |
| |
PUTS("+("); PRINTSF(n->c[i]); PUTS(")"); |
| |
if ( i >= 2 ) { |
| |
PUTS("*@s"); PRINTHAT; total_length += 13; |
| |
} else if ( i == 1 ) |
| |
PUTS("*@s"); |
| |
} |
| |
} |
| |
PUTS(")"); |
| |
} |
| |
} |
| |
|
| |
void PRINTSF(unsigned int i) |
| |
{ |
| |
if ( !i ) { |
| |
PUTS("0"); |
| |
} else |
| |
total_length += 15; |
| } |
} |