Lines Matching refs:deg
49 int deg;
90 for(j = x->deg; j >= 0; j--) {
93 x->deg = (j>0)? j : 0;
99 for(j=0; j <= x->deg; j++) {
200 y->deg = x->deg;
201 atoa(x->coe, y->coe, y->deg+1);
208 for(j=0; j <= y->deg; j++) {
224 if(x->deg == 0 && iszer(x->coe[0])) return 1;
243 int d = x->deg, j;
245 if(y->deg > d) d = y->deg;
247 if((j <= x->deg) && (j <= y->deg)) {
252 if((j <= x->deg) && (j > y->deg)) {
258 y->deg = d;
266 int d = x->deg, j;
268 if(y->deg > d) d = y->deg;
270 if((j <= x->deg) && (j <= y->deg)) {
275 if((j <= x->deg) && (j > y->deg)) {
282 y->deg = d;
291 int dega = a->deg, degb = b->deg, deg = dega + degb;
294 for(d=deg; d>=0; d--) {
307 atoa(mcand, b->coe, deg+1);
308 b->deg = deg;
313 atop(giant *a, poly x, int deg)
316 int adeg = abs(deg);
317 x->deg = adeg;
319 if(deg < 0) {
334 for(j=0; j < y->deg; j++) {
343 int deg = x->deg, j, coedigs;
345 g->sign = words*(1 + deg);
346 for(j=0; j <= deg; j++) {
361 xwords = (2*bits + log_2(x->deg+1) + 32 + 15)/16;
362 ywords = (2*bits + log_2(y->deg+1) + 32 + 15)/16;
372 gtog(y->coe[y->deg], globx); /* Save high coeff. */
373 y->deg += x->deg;
374 gtog(globx, y->coe[y->deg]); /* Move high coeff. */
376 mulg(x->coe[x->deg], y->coe[y->deg]); /* resolve high coeff. */
382 words = (2*bits + log_2(y->deg+1) + 32 + 15)/16;
385 gtog(y->coe[y->deg], globx); /* Save high coeff. */
386 y->deg += y->deg;
387 gtog(globx, y->coe[y->deg]); /* Move high coeff. */
389 mulg(y->coe[y->deg], y->coe[y->deg]); /* resolve high coeff. */
397 for(j=0; j <= x->deg; j++) {
402 for(j=0; j <= y->deg; j++) {
413 if(x->deg != y->deg) return 1;
414 for(j=0; j <= x->deg; j++) {
426 int n, degx = x->deg, degy = y->deg;
449 { int j, deg = x->deg;
451 for(j=0; j <= deg/2; j++) {
452 gtog(x->coe[deg-j], tmp);
453 gtog(x->coe[j], x->coe[deg-j]);
460 recipp(poly f, int deg)
463 int lim = deg + 1, prec = 1;
465 sscratch->deg = 0; itog(1, sscratch->coe[0]);
470 f->deg = prec-1;
473 tscratch->deg = prec-1;
478 tscratch->deg = prec-1;
481 sscratch->deg = prec-1;
494 for(j = start; j <= x->deg; j++) {
500 x->deg -= shift;
501 for(j=0; j<= x->deg; j++) {
515 int degx = x->deg, degy = y->deg, d, shift;
518 y->deg = 0;
534 rscratch->deg = d; justifyp(rscratch);
537 rscratch->deg = d; justifyp(rscratch);
541 if(d < degx) { x->deg = d; justifyp(x);}
543 rscratch->deg = d;
546 x->deg = degx; justifyp(x);
551 for(d = y->deg+1; d <= degx-shift; d++) itog(0, y->coe[d]);
552 y->deg = degx - shift;
567 for(j=0; j <= x->deg; j++) {
577 s[0]->deg = 0;
580 s[1]->deg = 0;
582 s[2]->deg = 0;
585 s[3]->deg = 4;
595 s[4]->deg = 6;
612 cubic->deg = 3;
655 gtog(smonic->coe[smonic->deg], aux); /* High coeff. */
658 s[0]->deg = smonic->deg + 1;
659 for(j=0; j <= s[0]->deg; j++) itog(1, s[0]->coe[j]);
764 for(j=0; j <= x->deg; j++) {printf("%d: ",j); gout(x->coe[j]);}
831 txd->deg = 0; itog(1, txd->coe[0]);
832 tyd->deg = 0; itog(1, tyd->coe[0]);
833 txn->deg = 1; itog(0, txn->coe[0]); itog(1, txn->coe[1]);