version 1.1, 1999/12/03 07:39:11 |
version 1.3, 2000/08/21 08:31:38 |
|
|
/* $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@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.2 2000/02/07 03:14:39 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; |
|
|
#define TAIL |
#define TAIL |
#define PUTS(s) fputs(s,OUT) |
#define PUTS(s) fputs(s,OUT) |
Line 51 char DFORMAT[BUFSIZ]; |
|
Line 100 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 80 extern char DFORMAT[BUFSIZ]; |
|
Line 131 extern char DFORMAT[BUFSIZ]; |
|
#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(); |
Line 122 void PRINTERR(); |
|
Line 153 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; |
{ |
{ |
|
|
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 |
|
|
} |
} |
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)); |
|
|
{ |
{ |
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"); |
} |
} |