=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2000/parse/arith.c,v retrieving revision 1.18 retrieving revision 1.25 diff -u -p -r1.18 -r1.25 --- OpenXM_contrib2/asir2000/parse/arith.c 2004/12/18 16:50:10 1.18 +++ OpenXM_contrib2/asir2000/parse/arith.c 2005/11/27 00:07:06 1.25 @@ -45,7 +45,7 @@ * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. * - * $OpenXM: OpenXM_contrib2/asir2000/parse/arith.c,v 1.17 2004/08/18 00:17:02 noro Exp $ + * $OpenXM: OpenXM_contrib2/asir2000/parse/arith.c,v 1.24 2005/11/25 07:18:32 noro Exp $ */ #include "ca.h" #include "parse.h" @@ -64,6 +64,8 @@ struct oARF arf[6]; ARF addfs, subfs, mulfs, divfs, remfs, pwrfs; void divsdc(); +int compqa(); +int compquote(); struct oAFUNC afunc[] = { /* ???=0 */ {0,0,0,0,0,0,0}, @@ -83,14 +85,15 @@ struct oAFUNC afunc[] = { /* O_F=14 */ {notdef,notdef,notdef,notdef,notdef,notdef,(int(*)())notdef}, /* O_GFMMAT=15 */ {notdef,notdef,notdef,notdef,notdef,notdef,(int(*)())notdef}, /* O_BYTEARRAY=16 */ {notdef,notdef,notdef,notdef,notdef,notdef,compbytearray}, -/* O_QUOTE=17 */ {addquote,subquote,mulquote,divquote,pwrquote,chsgnquote,(int(*)())notdef}, +/* O_QUOTE=17 */ {addquote,subquote,mulquote,divquote,pwrquote,chsgnquote,compquote}, /* O_OPTLIST=18 */ {notdef,notdef,notdef,notdef,notdef,notdef,(int(*)())notdef}, /* O_SYMBOL=19 */ {notdef,notdef,notdef,notdef,notdef,notdef,(int(*)())notdef}, /* O_RANGE=20 */ {notdef,notdef,notdef,notdef,notdef,notdef,(int(*)())notdef}, /* O_TB=21 */ {notdef,notdef,notdef,notdef,notdef,notdef,(int(*)())notdef}, /* O_DPV=22 */ {adddv,subdv,muldv,notdef,notdef,chsgndv,compdv}, -/* ???=23 */ {0,0,0,0,0,0,0}, -/* O_MAT=24 */ {AddMatI,SubMatI,MulMatG,notdef,notdef,ChsgnI,notdef}, +/* O_QUOTEARG=23 */ {notdef,notdef,notdef,notdef,notdef,notdef,compqa}, +/* O_MAT=24 */ {AddMatI,SubMatI,MulMatG,notdef,notdef,ChsgnI,(int(*)())notdef}, +/* O_NBP=25 */ {addnbp,subnbp,mulnbp,notdef,pwrnbp,chsgnnbp,compnbp}, }; void arf_init() { @@ -114,7 +117,7 @@ Obj a,b,*r; *r = a; else if ( OID(a) == OID(b) ) (*afunc[OID(a)].add)(vl,a,b,r); - else if ( (mid = MAX(OID(a),OID(b))) <= O_R ) + else if ( (mid = MAX(OID(a),OID(b))) <= O_R || mid == O_QUOTE ) (*afunc[mid].add)(vl,a,b,r); else if ( (mid = MAX(OID(a),OID(b))) == O_DP && MIN(OID(a),OID(b)) <= O_R ) (*afunc[mid].add)(vl,a,b,r); @@ -137,7 +140,7 @@ Obj a,b,*r; *r = a; else if ( OID(a) == OID(b) ) (*afunc[OID(a)].sub)(vl,a,b,r); - else if ( (mid = MAX(OID(a),OID(b))) <= O_R ) + else if ( (mid = MAX(OID(a),OID(b))) <= O_R || mid == O_QUOTE ) (*afunc[mid].sub)(vl,a,b,r); else if ( (mid = MAX(OID(a),OID(b))) == O_DP && MIN(OID(a),OID(b)) <= O_R ) (*afunc[mid].sub)(vl,a,b,r); @@ -164,11 +167,11 @@ Obj a,b,*r; *r = 0; } else if ( (aid = OID(a)) == (bid = OID(b)) ) (*(afunc[aid].mul))(vl,a,b,r); - else if ( (mid = MAX(aid,bid)) <= O_R ) + else if ( (mid = MAX(aid,bid)) <= O_R || mid == O_QUOTE ) (*afunc[mid].mul)(vl,a,b,r); else { switch ( aid ) { - case O_N: case O_P: + case O_N: case O_P: case O_NBP: (*afunc[mid].mul)(vl,a,b,r); break; case O_R: @@ -223,7 +226,7 @@ Obj a,b,*r; else if ( (OID(a) == OID(b)) ) (*(afunc[OID(a)].div))(vl,a,b,r); else if ( (mid = MAX(OID(a),OID(b))) <= O_R || - (mid == O_MAT) || (mid == O_VECT) || (mid == O_DP) ) + (mid == O_MAT) || (mid == O_VECT) || (mid == O_DP) || mid == O_QUOTE ) (*afunc[mid].div)(vl,a,b,r); else notdef(vl,a,b,r); @@ -258,14 +261,14 @@ Obj a,e,*r; *r = 0; else mkpow(vl,a,e,r); - } else if ( OID(a) == O_QUOTE ) - (*(afunc[O_QUOTE].pwr))(vl,a,e,r); - else if ( !e ) { + } else if ( !e ) { if ( OID(a) == O_MAT ) (*(afunc[O_MAT].pwr))(vl,a,e,r); else *r = (pointer)ONE; - } else if ( (OID(e) <= O_N) && INT(e) ) { + } else if ( OID(a) == O_QUOTE || OID(e) == O_QUOTE ) + (*(afunc[O_QUOTE].pwr))(vl,a,e,r); + else if ( (OID(e) <= O_N) && INT(e) ) { if ( (OID(a) == O_P) && (SGN((Q)e) < 0) ) { if ( allow_laurent ) (*(afunc[O_P].pwr))(vl,a,e,r);