version 1.17, 2001/10/09 01:36:22 |
version 1.18, 2002/01/08 04:14:39 |
|
|
* 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.17 2001/10/09 01:36:22 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "al.h" |
#include "al.h" |
|
|
{ |
{ |
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 ) |
|
|
#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 |