version 1.1.1.1, 1999/12/03 07:39:11 |
version 1.7, 2000/12/15 05:30:08 |
|
|
/* $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.6 2000/12/13 10:54:09 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; |
|
|
#define TAIL |
#define TAIL |
#define PUTS(s) fputs(s,OUT) |
#define PUTS(s) fputs(s,OUT) |
Line 38 char DFORMAT[BUFSIZ]; |
|
Line 88 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 PRINTERR printerr |
#define PRINTERR printerr |
#define PRINTLF printlf |
#define PRINTLF printlf |
#define PRINTLOP printlop |
#define PRINTLOP printlop |
Line 51 char DFORMAT[BUFSIZ]; |
|
Line 102 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; |
|
|
#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 127 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 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 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() { |
|
|
} |
} |
|
|
#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 |
|
|
void sprintbf(a) |
void sprintbf(a) |
BF a; |
BF a; |
{ |
{ |
/* XXX : not implemented yet */ |
char *str; |
|
char *GENtostr(); |
|
|
|
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; |
default: |
default: |
break; |
break; |
} |
} |
|
|
} |
} |
break; |
break; |
case N_R: |
case N_R: |
TAIL PRINTF(OUT,"%g",BDY((Real)q)); |
TAIL PRINTF(OUT,double_output?"%f":"%g",BDY((Real)q)); |
break; |
break; |
case N_A: |
case N_A: |
PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")"); |
PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")"); |
|
|
} |
} |
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 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"); |
} |
} |