version 1.1, 2018/09/19 05:45:07 |
version 1.2, 2018/09/28 08:20:28 |
Line 272 void pwrp(VL vl,P p,Z q,P *pr) |
|
Line 272 void pwrp(VL vl,P p,Z q,P *pr) |
|
} else if ( !INT(q) ) { |
} else if ( !INT(q) ) { |
error("pwrp: can't calculate fractional power."); *pr = 0; |
error("pwrp: can't calculate fractional power."); *pr = 0; |
} else if ( smallz(q) ) { |
} else if ( smallz(q) ) { |
n = QTOS(q); x = (P *)ALLOCA((n+1)*sizeof(pointer)); |
n = ZTOS(q); x = (P *)ALLOCA((n+1)*sizeof(pointer)); |
NEWDC(dct); DEG(dct) = DEG(dc); COEF(dct) = COEF(dc); |
NEWDC(dct); DEG(dct) = DEG(dc); COEF(dct) = COEF(dc); |
NEXT(dct) = 0; MKP(VR(p),dct,t); |
NEXT(dct) = 0; MKP(VR(p),dct,t); |
for ( i = 0, u = (P)One; i < n; i++ ) { |
for ( i = 0, u = (P)One; i < n; i++ ) { |