version 1.99, 2004/09/14 07:23:34 |
version 1.101, 2004/09/14 10:00:26 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.98 2004/08/18 00:17:02 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.100 2004/09/14 09:25:48 noro Exp $ */ |
|
|
#include "nd.h" |
#include "nd.h" |
|
|
Line 3214 ND nd_quo(ND p,ND d) |
|
Line 3214 ND nd_quo(ND p,ND d) |
|
divq(HCQ(p),HCQ(d),&q); |
divq(HCQ(p),HCQ(d),&q); |
chsgnq(q,&CQ(mq)); |
chsgnq(q,&CQ(mq)); |
t = nd_mul_nm_trunc(mq,d,HDL(d)); |
t = nd_mul_nm_trunc(mq,d,HDL(d)); |
|
CQ(mq) = q; |
p = nd_add(0,p,t); |
p = nd_add(0,p,t); |
} |
} |
NEXT(mq) = 0; |
NEXT(mq) = 0; |
Line 4725 void nd_det(MAT f,P *rp) |
|
Line 4726 void nd_det(MAT f,P *rp) |
|
} |
} |
d = mjj; |
d = mjj; |
} |
} |
|
if ( sgn < 0 ) |
|
nd_mul_c_q(d,mone); |
dv = ndtondv(0,d); |
dv = ndtondv(0,d); |
*rp = ndvtop(0,CO,fv,dv); |
*rp = ndvtop(0,CO,fv,dv); |
} |
} |