version 1.23, 2003/06/07 16:40:25 |
version 1.24, 2003/06/09 16:18:10 |
|
|
* 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/pexpr.c,v 1.22 2003/05/29 16:45:01 saito Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/io/spexpr.c,v 1.23 2003/06/07 16:40:25 saito Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "al.h" |
#include "al.h" |
|
|
case N_IP: |
case N_IP: |
case N_IntervalBigFloat: |
case N_IntervalBigFloat: |
switch ( outputstyle ) { |
switch ( outputstyle ) { |
case 0: |
|
PUTS("["); |
|
PRINTNUM(INF((Itv)q)); |
|
PUTS(","); |
|
PRINTNUM(SUP((Itv)q)); |
|
PUTS("]"); |
|
break; |
|
case 1: |
case 1: |
PUTS("intval("); |
PUTS("intval("); |
PRINTNUM(INF((Itv)q)); |
PRINTNUM(INF((Itv)q)); |
|
|
PRINTNUM(SUP((Itv)q)); |
PRINTNUM(SUP((Itv)q)); |
PUTS(")"); |
PUTS(")"); |
break; |
break; |
|
case 0: |
default: |
default: |
PUTS("["); |
PUTS("["); |
PRINTNUM(INF((Itv)q)); |
PRINTNUM(INF((Itv)q)); |
|
|
switch (printmode) { |
switch (printmode) { |
case PRINTF_E: |
case PRINTF_E: |
switch ( outputstyle ) { |
switch ( outputstyle ) { |
case 0: |
|
TAIL PRINTF(OUT, "[%.16e,%.16e]", |
|
INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
|
break; |
|
case 1: |
case 1: |
TAIL PRINTF(OUT, "intval(%.16e,%.16e)", |
TAIL PRINTF(OUT, "intval(%.16e,%.16e)", |
INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
break; |
break; |
|
case 0: |
default: |
default: |
TAIL PRINTF(OUT, "[%.16e,%.16e]", |
TAIL PRINTF(OUT, "[%.16e,%.16e]", |
INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
|
|
break; |
break; |
case MID_PRINTF_G: |
case MID_PRINTF_G: |
switch ( outputstyle ) { |
switch ( outputstyle ) { |
case 0: |
|
TAIL PRINTF(OUT, "<%g,%g>", |
|
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5, |
|
(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
|
break; |
|
case 1: |
case 1: |
TAIL PRINTF(OUT, "intvalm(%g,%g)", |
TAIL PRINTF(OUT, "intvalm(%g,%g)", |
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5, |
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5, |
(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
break; |
break; |
|
case 0: |
default: |
default: |
TAIL PRINTF(OUT, "<%g,%g>", |
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: |
switch ( outputstyle ) { |
switch ( outputstyle ) { |
case 0: |
|
TAIL PRINTF(OUT, "<%.16e,%.16e>", |
|
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5, |
|
(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
|
break; |
|
case 1: |
case 1: |
TAIL PRINTF(OUT, "intvalm(%.16e,%.16e)", |
TAIL PRINTF(OUT, "intvalm(%.16e,%.16e)", |
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5, |
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5, |
(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5); |
break; |
break; |
|
case 0: |
default: |
default: |
TAIL PRINTF(OUT, "<%.16e,%.16e>", |
TAIL PRINTF(OUT, "<%.16e,%.16e>", |
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5, |
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5, |
|
|
case PRINTF_G: |
case PRINTF_G: |
default: |
default: |
switch ( outputstyle ) { |
switch ( outputstyle ) { |
case 0: |
|
TAIL PRINTF(OUT, "[%g,%g]", |
|
INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
|
break; |
|
case 1: |
case 1: |
TAIL PRINTF(OUT, "intval(%g,%g)", |
TAIL PRINTF(OUT, "intval(%g,%g)", |
INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
break; |
break; |
|
case 0: |
default: |
default: |
TAIL PRINTF(OUT, "[%g,%g]", |
TAIL PRINTF(OUT, "[%g,%g]", |
INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
INF((IntervalDouble)q),SUP((IntervalDouble)q)); |
|
|
pointer *ptr; |
pointer *ptr; |
|
|
switch ( outputstyle ) { |
switch ( outputstyle ) { |
case 0: |
|
PUTS("[ "); |
|
for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) { |
|
PRINTEXPR(vl,ptr[i]); PUTS(" "); |
|
} |
|
PUTS("]"); |
|
break; |
|
case 1: |
case 1: |
PUTS("vect("); |
PUTS("vect("); |
for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) { |
for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) { |
if ( i != 0 ) PUTS(", "); |
if ( i != 0 ) PUTS(","); |
PRINTEXPR(vl,ptr[i]); |
PRINTEXPR(vl,ptr[i]); |
} |
} |
PUTS(")"); |
PUTS(")"); |
break; |
break; |
|
case 0: |
default: |
default: |
PUTS("[ "); |
PUTS("[ "); |
for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) { |
for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) { |
|
|
pointer *ptr; |
pointer *ptr; |
|
|
switch ( outputstyle ) { |
switch ( outputstyle ) { |
case 0: |
|
for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) { |
|
PUTS("[ "); |
|
for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) { |
|
PRINTEXPR(vl,ptr[j]); PUTS(" "); |
|
} |
|
PUTS("]"); |
|
if ( i < r - 1 ) |
|
PUTS("\n"); |
|
} |
|
break; |
|
case 1: |
case 1: |
PUTS("mat(\n"); |
PUTS("mat(\n"); |
for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) { |
for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) { |
if ( i != 0 ) PUTS(",\n"); |
if ( i != 0 ) PUTS(",\n"); |
PUTS("[ "); |
PUTS("[ "); |
for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) { |
for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) { |
if ( j != 0 ) PUTS(", "); |
if ( j != 0 ) PUTS(","); |
PRINTEXPR(vl,ptr[j]); |
PRINTEXPR(vl,ptr[j]); |
} |
} |
PUTS("]"); |
PUTS(" ]"); |
} |
} |
PUTS(")"); |
PUTS(")"); |
break; |
break; |
|
case 0: |
default: |
default: |
for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) { |
for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) { |
PUTS("[ "); |
PUTS("[ "); |