コード例 #1
0
ファイル: decomp.c プロジェクト: nhatcher/lie
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;
}
コード例 #2
0
ファイル: lr.c プロジェクト: d4g33z/lie
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;
}
コード例 #3
0
ファイル: decomp.c プロジェクト: nhatcher/lie
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;
}