| 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++ ) { |