void wt_init(lie_Index r) { sorted=copypoly(poly_null(r)); pos_acc=mkpoly(ACCMIN,r); neg_acc=mkpoly(ACCMIN,r); pos_acc->nrows=0; neg_acc->nrows=0; }
poly* LR_tensor(poly* p,poly* q) { _index i,j,n=p->ncols; poly* res=poly_null(n); for (i=0; i<p->nrows; ++i) for (j=0; j<q->nrows; ++j) res=Addmul_pol_pol_bin(res,LR_tensor_irr(p->elm[i],q->elm[j],n) ,mult(p->coef[i],q->coef[j])); return res; }
poly* Vdecomp(poly* p) { lie_Index i,r=Lierank(grp); poly* result=poly_null(r); cur_expon=mkintarray(r); /* large enough */ for (i=0; i<p->nrows; ++i) result=Addmul_pol_pol_bin(result,vdecomp_irr(p->elm[i]),p->coef[i]); freearr(cur_expon); return result; }