version 1.18, 2002/01/08 04:14:39 |
version 1.22, 2003/05/29 16:45:01 |
|
|
* 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/spexpr.c,v 1.17 2001/10/09 01:36:22 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/io/spexpr.c,v 1.21 2003/02/14 22:29:15 ohara Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "al.h" |
#include "al.h" |
|
|
#endif |
#endif |
|
|
#define PRINTHAT (fortran_output?PUTS("**"):PUTS("^")) |
#define PRINTHAT (fortran_output?PUTS("**"):PUTS("^")) |
|
extern int prefixmode; |
|
|
#ifdef FPRINT |
#ifdef FPRINT |
FILE *asir_out; |
FILE *asir_out; |
|
|
#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 113 extern int double_output; |
|
Line 114 extern int double_output; |
|
extern int real_digit; |
extern int real_digit; |
extern int print_quote; |
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) |
#define PRINTF sprintf |
#define PRINTF sprintf |
Line 124 extern int print_quote; |
|
Line 126 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 |
|
|
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; |
{ |
{ |
|
|
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_IntervalBigFloat: |
case N_IntervalBigFloat: |
PUTS("["); |
if ( prefixmode == 0 ) { |
PRINTNUM(INF((Itv)q)); |
PUTS("["); |
PUTS(","); |
PRINTNUM(INF((Itv)q)); |
PRINTNUM(SUP((Itv)q)); |
PUTS(","); |
PUTS("]"); |
PRINTNUM(SUP((Itv)q)); |
|
PUTS("]"); |
|
} else { |
|
PUTS("intval("); |
|
PRINTNUM(INF((Itv)q)); |
|
PUTS(","); |
|
PRINTNUM(SUP((Itv)q)); |
|
PUTS(")"); |
|
} |
break; |
break; |
case N_IntervalDouble: |
case N_IntervalDouble: |
switch (printmode) { |
switch (printmode) { |
case PRINTF_E: |
case PRINTF_E: |
TAIL PRINTF(OUT, "[%.16e,%.16e]",INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
if (prefixmode == 0 ) { |
|
TAIL PRINTF(OUT, "[%.16e,%.16e]",INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
|
} else { |
|
TAIL PRINTF(OUT, "intval(%.16e,%.16e)",INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
|
} |
#if defined(ITVDEBUG) |
#if defined(ITVDEBUG) |
printbin(INF((IntervalDouble)q)); |
printbin(INF((IntervalDouble)q)); |
printbin(SUP((IntervalDouble)q)); |
printbin(SUP((IntervalDouble)q)); |
#endif |
#endif |
break; |
break; |
case MID_PRINTF_G: |
case MID_PRINTF_G: |
TAIL PRINTF(OUT, "<%g,%g>", (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
if (prefixmode == 0) { |
|
TAIL PRINTF(OUT, "<%g,%g>", (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
|
} else { |
|
TAIL PRINTF(OUT, "intvalm(%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((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
if (prefixmode == 0) { |
|
TAIL PRINTF(OUT, "<%.16e,%.16e>", (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
|
} else { |
|
TAIL PRINTF(OUT, "intvalm(%.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((IntervalDouble)q),SUP((IntervalDouble)q)); |
if (prefixmode == 0) { |
|
TAIL PRINTF(OUT, "[%g,%g]",INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
|
} else { |
|
TAIL PRINTF(OUT, "intval(%g,%g)",INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
|
} |
break; |
break; |
} |
} |
break; |
break; |
|
|
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 { |
|
|
int i; |
int i; |
pointer *ptr; |
pointer *ptr; |
|
|
PUTS("[ "); |
if ( prefixmode == 0 ) { |
for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) { |
PUTS("[ "); |
PRINTEXPR(vl,ptr[i]); PUTS(" "); |
for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) { |
|
PRINTEXPR(vl,ptr[i]); PUTS(" "); |
|
} |
|
PUTS("]"); |
|
} else { |
|
PUTS("vect("); |
|
for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) { |
|
if ( i != 0 ) PUTS(", "); |
|
PRINTEXPR(vl,ptr[i]); |
|
} |
|
PUTS(")"); |
} |
} |
PUTS("]"); |
|
} |
} |
|
|
void PRINTMAT(vl,mat) |
void PRINTMAT(vl,mat) |
|
|
int i,j,r,c; |
int i,j,r,c; |
pointer *ptr; |
pointer *ptr; |
|
|
for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) { |
if ( prefixmode == 0 ) { |
PUTS("[ "); |
for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) { |
for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) { |
PUTS("[ "); |
PRINTEXPR(vl,ptr[j]); PUTS(" "); |
for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) { |
|
PRINTEXPR(vl,ptr[j]); PUTS(" "); |
|
} |
|
PUTS("]"); |
|
if ( i < r - 1 ) |
|
PUTS("\n"); |
} |
} |
PUTS("]"); |
} else { |
if ( i < r - 1 ) |
PUTS("mat(\n"); |
PUTS("\n"); |
for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) { |
|
if ( i != 0 ) PUTS(",\n"); |
|
PUTS("[ "); |
|
for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) { |
|
PRINTEXPR(vl,ptr[j]); PUTS(" "); |
|
} |
|
PUTS("]"); |
|
} |
|
PUTS(")"); |
} |
} |
} |
} |
|
|