=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2000/builtin/round.c,v retrieving revision 1.4 retrieving revision 1.5 diff -u -p -r1.4 -r1.5 --- OpenXM_contrib2/asir2000/builtin/round.c 2005/02/08 17:20:56 1.4 +++ OpenXM_contrib2/asir2000/builtin/round.c 2007/01/30 00:28:26 1.5 @@ -1,208 +1,6 @@ -/* $OpenXM: OpenXM_contrib2/asir2000/builtin/round.c,v 1.3 2003/12/02 07:00:12 noro Exp $ */ - -#undef INTERVAL -#define INTERVAL #include "ca.h" #include "parse.h" -#include "interval.h" -void Padd_pinf(),Padd_minf(); -void Psub_pinf(),Psub_minf(); -void Pmul_pinf(),Pmul_minf(); -void Pdiv_pinf(),Pdiv_minf(); - struct ftab round_tab[] = { - {"add_pinf",Padd_pinf,2}, - {"add_minf",Padd_minf,2}, - {"sub_pinf",Psub_pinf,2}, - {"sub_minf",Psub_minf,2}, - {"mul_pinf",Pmul_pinf,2}, - {"mul_minf",Pmul_minf,2}, - {"div_pinf",Pdiv_pinf,2}, - {"div_minf",Pdiv_minf,2}, {0,0,0}, }; - -void Padd_pinf(NODE arg,Real *rp) -{ - Real a,b,r; - double c; - - a = (Real)ARG0(arg); - b = (Real)ARG1(arg); - if ( !a ) - *rp = b; - else if ( !b ) - *rp = a; - else if ( NID(a) != N_R || NID(b) != N_R ) - error("add_pinf : invalid argument"); - else { - FPPLUSINF - c = BDY(a)+BDY(b); - FPNEAREST - MKReal(c,r); - *rp = r; - } -} - -void Padd_minf(NODE arg,Real *rp) -{ - Real a,b,r; - double c; - - a = (Real)ARG0(arg); - b = (Real)ARG1(arg); - if ( !a ) - *rp = b; - else if ( !b ) - *rp = a; - else if ( NID(a) != N_R || NID(b) != N_R ) - error("add_pinf : invalid argument"); - else { - FPMINUSINF - c = BDY(a)+BDY(b); - FPNEAREST - MKReal(c,r); - *rp = r; - } -} - -void Psub_pinf(NODE arg,Real *rp) -{ - Real a,b,r; - double c; - - a = (Real)ARG0(arg); - b = (Real)ARG1(arg); - if ( !a ) { - if ( !b ) - *rp = 0; - else { - c = -BDY(b); - MKReal(c,r); - *rp = r; - } - } else if ( !b ) - *rp = a; - else if ( NID(a) != N_R || NID(b) != N_R ) - error("sub_pinf : invalid argument"); - else { - FPPLUSINF - c = BDY(a)-BDY(b); - FPNEAREST - MKReal(c,r); - *rp = r; - } -} - -void Psub_minf(NODE arg,Real *rp) -{ - Real a,b,r; - double c; - - a = (Real)ARG0(arg); - b = (Real)ARG1(arg); - if ( !a ) { - if ( !b ) - *rp = 0; - else { - c = -BDY(b); - MKReal(c,r); - *rp = r; - } - } else if ( !b ) - *rp = a; - else if ( NID(a) != N_R || NID(b) != N_R ) - error("sub_minf : invalid argument"); - else { - FPMINUSINF - c = BDY(a)-BDY(b); - FPNEAREST - MKReal(c,r); - *rp = r; - } -} - -void Pmul_pinf(NODE arg,Real *rp) -{ - Real a,b,r; - double c; - - a = (Real)ARG0(arg); - b = (Real)ARG1(arg); - if ( !a || !b ) - *rp = 0; - else if ( NID(a) != N_R || NID(b) != N_R ) - error("mul_pinf : invalid argument"); - else { - FPPLUSINF - c = BDY(a)*BDY(b); - FPNEAREST - MKReal(c,r); - *rp = r; - } -} - -void Pmul_minf(NODE arg,Real *rp) -{ - Real a,b,r; - double c; - - a = (Real)ARG0(arg); - b = (Real)ARG1(arg); - if ( !a || !b ) - *rp = 0; - else if ( NID(a) != N_R || NID(b) != N_R ) - error("mul_minf : invalid argument"); - else { - FPMINUSINF - c = BDY(a)*BDY(b); - FPNEAREST - MKReal(c,r); - *rp = r; - } -} - -void Pdiv_pinf(NODE arg,Real *rp) -{ - Real a,b,r; - double c; - - a = (Real)ARG0(arg); - b = (Real)ARG1(arg); - if ( !b ) - error("div_pinf : division by 0"); - else if ( !a ) - *rp = 0; - else if ( NID(a) != N_R || NID(b) != N_R ) - error("div_pinf : invalid argument"); - else { - FPPLUSINF - c = BDY(a)/BDY(b); - FPNEAREST - MKReal(c,r); - *rp = r; - } -} - -void Pdiv_minf(NODE arg,Real *rp) -{ - Real a,b,r; - double c; - - a = (Real)ARG0(arg); - b = (Real)ARG1(arg); - if ( !b ) - error("div_minf : division by 0"); - else if ( !a ) - *rp = 0; - else if ( NID(a) != N_R || NID(b) != N_R ) - error("div_minf : invalid argument"); - else { - FPMINUSINF - c = BDY(a)/BDY(b); - FPNEAREST - MKReal(c,r); - *rp = r; - } -}