version 1.20, 2005/04/07 08:33:12 |
version 1.23, 2005/11/16 23:42:54 |
|
|
* 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/parse/arith.c,v 1.19 2005/02/08 18:06:05 saito Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/parse/arith.c,v 1.22 2005/10/03 00:29:11 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
Line 93 struct oAFUNC afunc[] = { |
|
Line 93 struct oAFUNC afunc[] = { |
|
/* O_DPV=22 */ {adddv,subdv,muldv,notdef,notdef,chsgndv,compdv}, |
/* O_DPV=22 */ {adddv,subdv,muldv,notdef,notdef,chsgndv,compdv}, |
/* O_QUOTEARG=23 */ {notdef,notdef,notdef,notdef,notdef,notdef,compqa}, |
/* O_QUOTEARG=23 */ {notdef,notdef,notdef,notdef,notdef,notdef,compqa}, |
/* O_MAT=24 */ {AddMatI,SubMatI,MulMatG,notdef,notdef,ChsgnI,(int(*)())notdef}, |
/* O_MAT=24 */ {AddMatI,SubMatI,MulMatG,notdef,notdef,ChsgnI,(int(*)())notdef}, |
|
/* O_NBP=25 */ {addnbp,subnbp,mulnbp,notdef,pwrnbp,chsgnnbp,(int(*)())notdef}, |
}; |
}; |
|
|
void arf_init() { |
void arf_init() { |
|
|
*r = a; |
*r = a; |
else if ( OID(a) == OID(b) ) |
else if ( OID(a) == OID(b) ) |
(*afunc[OID(a)].add)(vl,a,b,r); |
(*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); |
(*afunc[mid].add)(vl,a,b,r); |
else if ( (mid = MAX(OID(a),OID(b))) == O_DP && MIN(OID(a),OID(b)) <= O_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); |
(*afunc[mid].add)(vl,a,b,r); |
|
|
*r = a; |
*r = a; |
else if ( OID(a) == OID(b) ) |
else if ( OID(a) == OID(b) ) |
(*afunc[OID(a)].sub)(vl,a,b,r); |
(*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); |
(*afunc[mid].sub)(vl,a,b,r); |
else if ( (mid = MAX(OID(a),OID(b))) == O_DP && MIN(OID(a),OID(b)) <= O_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); |
(*afunc[mid].sub)(vl,a,b,r); |
|
|
*r = 0; |
*r = 0; |
} else if ( (aid = OID(a)) == (bid = OID(b)) ) |
} else if ( (aid = OID(a)) == (bid = OID(b)) ) |
(*(afunc[aid].mul))(vl,a,b,r); |
(*(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); |
(*afunc[mid].mul)(vl,a,b,r); |
else { |
else { |
switch ( aid ) { |
switch ( aid ) { |
|
|
else if ( (OID(a) == OID(b)) ) |
else if ( (OID(a) == OID(b)) ) |
(*(afunc[OID(a)].div))(vl,a,b,r); |
(*(afunc[OID(a)].div))(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_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); |
(*afunc[mid].div)(vl,a,b,r); |
else |
else |
notdef(vl,a,b,r); |
notdef(vl,a,b,r); |
|
|
*r = 0; |
*r = 0; |
else |
else |
mkpow(vl,a,e,r); |
mkpow(vl,a,e,r); |
} else if ( OID(a) == O_QUOTE ) |
} else if ( !e ) { |
(*(afunc[O_QUOTE].pwr))(vl,a,e,r); |
|
else if ( !e ) { |
|
if ( OID(a) == O_MAT ) |
if ( OID(a) == O_MAT ) |
(*(afunc[O_MAT].pwr))(vl,a,e,r); |
(*(afunc[O_MAT].pwr))(vl,a,e,r); |
else |
else |
*r = (pointer)ONE; |
*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 ( (OID(a) == O_P) && (SGN((Q)e) < 0) ) { |
if ( allow_laurent ) |
if ( allow_laurent ) |
(*(afunc[O_P].pwr))(vl,a,e,r); |
(*(afunc[O_P].pwr))(vl,a,e,r); |