version 1.16, 2001/09/03 07:01:09 |
version 1.21, 2003/02/14 22:29:15 |
|
|
* 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: OpenXM_contrib2/asir2000/io/spexpr.c,v 1.20 2002/09/11 07:23:27 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "al.h" |
#include "al.h" |
|
|
#define PRINTV printv |
#define PRINTV printv |
#define PRINTEXPR printexpr |
#define PRINTEXPR printexpr |
#define PRINTNUM printnum |
#define PRINTNUM printnum |
#define PRINTP printp |
#define PRINTP asir_printp |
#define PRINTR printr |
#define PRINTR printr |
#define PRINTLIST printlist |
#define PRINTLIST printlist |
#define PRINTVECT printvect |
#define PRINTVECT printvect |
Line 124 extern int print_quote; |
|
Line 124 extern int print_quote; |
|
#define PRINTV sprintv |
#define PRINTV sprintv |
#define PRINTEXPR sprintexpr |
#define PRINTEXPR sprintexpr |
#define PRINTNUM sprintnum |
#define PRINTNUM sprintnum |
#define PRINTP sprintp |
#define PRINTP asir_sprintp |
#define PRINTR sprintr |
#define PRINTR sprintr |
#define PRINTLIST sprintlist |
#define PRINTLIST sprintlist |
#define PRINTVECT sprintvect |
#define PRINTVECT sprintvect |
Line 170 void PRINTCPLX(); |
|
Line 170 void PRINTCPLX(); |
|
void PRINTLM(); |
void PRINTLM(); |
void PRINTLF(); |
void PRINTLF(); |
void PRINTUP2(); |
void PRINTUP2(); |
|
void PRINTUP(); |
|
void PRINTUM(); |
|
void PRINTFOP(); |
|
void PRINTEOP(); |
|
void PRINTLOP(); |
|
void PRINTQOP(); |
|
void PRINTSF(); |
|
|
#ifdef FPRINT |
#ifdef FPRINT |
void output_init() { |
void output_init() { |
|
|
{ |
{ |
if ( NUM(p) ) |
if ( NUM(p) ) |
#if defined(INTERVAL) |
#if defined(INTERVAL) |
if ( NID(p) != N_IP && NID(p) != N_ID && NID(p) != N_IT && NID(p) != N_IF |
if ( NID(p) != N_IP && NID(p) != N_IntervalDouble && NID(p) != N_IntervalQuad && NID(p) != N_IntervalBigFloat |
&& compnum(CO,(Num)p,0) < 0 ) |
&& compnum(CO,(Num)p,0) < 0 ) |
#else |
#else |
if ( compnum(CO,(Num)p,0) < 0 ) |
if ( compnum(CO,(Num)p,0) < 0 ) |
|
|
return (mmono(COEF(DC(p)))); |
return (mmono(COEF(DC(p)))); |
} |
} |
|
|
#if PARI |
#if defined(PARI) |
void printbf(a) |
void printbf(a) |
BF a; |
BF a; |
{ |
{ |
|
void sor(); |
|
|
sor(a->body,double_output ? 'f' : 'g',-1,0); |
sor(a->body,double_output ? 'f' : 'g',-1,0); |
} |
} |
#endif |
#endif |
|
|
s[0] = 0; buf = s; |
s[0] = 0; buf = s; |
} |
} |
|
|
#if PARI |
#if defined(PARI) |
#include "genpari.h" |
#include "genpari.h" |
|
|
void myoutbrute(g) |
void myoutbrute(g) |
|
|
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 |
#if defined(INTERVAL) |
#if defined(INTERVAL) |
case N_IP: |
case N_IP: |
case N_IF: |
case N_IntervalBigFloat: |
PUTS("["); |
PUTS("["); |
PRINTNUM(INF((Itv)q)); |
PRINTNUM(INF((Itv)q)); |
PUTS(","); |
PUTS(","); |
PRINTNUM(SUP((Itv)q)); |
PRINTNUM(SUP((Itv)q)); |
PUTS("]"); |
PUTS("]"); |
break; |
break; |
case N_ID: |
case N_IntervalDouble: |
switch (printmode) { |
switch (printmode) { |
case PRINTF_E: |
case PRINTF_E: |
TAIL PRINTF(OUT, "[%.16e,%.16e]",INF((ItvD)q),SUP((ItvD)q)); |
TAIL PRINTF(OUT, "[%.16e,%.16e]",INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
#if defined(ITVDEBUG) |
#if defined(ITVDEBUG) |
printbin(INF((ItvD)q)); |
printbin(INF((IntervalDouble)q)); |
printbin(SUP((ItvD)q)); |
printbin(SUP((IntervalDouble)q)); |
#endif |
#endif |
break; |
break; |
case MID_PRINTF_G: |
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); |
TAIL PRINTF(OUT, "<%g,%g>", (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
break; |
break; |
case MID_PRINTF_E: |
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); |
TAIL PRINTF(OUT, "<%.16e,%.16e>", (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
break; |
break; |
case PRINTF_G: |
case PRINTF_G: |
default: |
default: |
TAIL PRINTF(OUT, "[%g,%g]",INF((ItvD)q),SUP((ItvD)q)); |
TAIL PRINTF(OUT, "[%g,%g]",INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
break; |
break; |
} |
} |
break; |
break; |
|
|
if ( a->i ) { |
if ( a->i ) { |
#if defined(INTERVAL) |
#if defined(INTERVAL) |
if ( a->r && ((compnum(0,a->i,0) > 0) |
if ( a->r && ((compnum(0,a->i,0) > 0) |
|| NID(a->i) == N_IP || NID(a->i) == N_ID |
|| NID(a->i) == N_IP || NID(a->i) == N_IntervalDouble |
|| NID(a->i) == N_IT || NID(a->i) == N_IF) ) |
|| NID(a->i) == N_IntervalQuad || NID(a->i) == N_IntervalBigFloat) ) |
#else |
#else |
if ( a->r && (compnum(0,a->i,0) > 0) ) |
if ( a->r && (compnum(0,a->i,0) > 0) ) |
#endif |
#endif |
|
|
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 ) { |
TAIL PRINTF(OUT,"%s",NAME(pf)); |
TAIL PRINTF(OUT,"%s",NAME(pf)); |
} else { |
} else { |
|
|
GFMMAT mat; |
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; |
|
|
BYTEARRAY array; |
BYTEARRAY array; |
{ |
{ |
int len,i; |
int len,i; |
unsigned int t; |
|
unsigned char *b; |
unsigned char *b; |
|
|
len = array->len; |
len = array->len; |
|
|
} |
} |
} |
} |
|
|
PRINTFOP(vl,f) |
void PRINTFOP(vl,f) |
VL vl; |
VL vl; |
F f; |
F f; |
{ |
{ |
|
|
} |
} |
} |
} |
|
|
PRINTEOP(vl,f) |
void PRINTEOP(vl,f) |
VL vl; |
VL vl; |
F f; |
F f; |
{ |
{ |
|
|
PUTS(")"); |
PUTS(")"); |
} |
} |
|
|
PRINTLOP(vl,f) |
void PRINTLOP(vl,f) |
VL vl; |
VL vl; |
F f; |
F f; |
{ |
{ |
|
|
PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0"); |
PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0"); |
} |
} |
|
|
PRINTQOP(vl,f) |
void PRINTQOP(vl,f) |
VL vl; |
VL vl; |
F f; |
F f; |
{ |
{ |
|
|
PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")"); |
PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")"); |
} |
} |
|
|
PRINTUP(n) |
void PRINTUP(n) |
UP n; |
UP n; |
{ |
{ |
int i,d; |
int i,d; |
|
|
} |
} |
} |
} |
|
|
PRINTUM(n) |
void PRINTUM(n) |
UM n; |
UM n; |
{ |
{ |
int i,d; |
int i,d; |
|
|
} |
} |
} |
} |
|
|
PRINTSF(i) |
void PRINTSF(i) |
unsigned int i; |
unsigned int i; |
{ |
{ |
if ( !i ) { |
if ( !i ) { |