version 1.41, 2006/08/08 02:39:17 |
version 1.45, 2012/12/17 07:20:44 |
|
|
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* |
* |
* $OpenXM: OpenXM_contrib2/asir2000/engine/dist.c,v 1.40 2006/04/16 00:51:13 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/dist.c,v 1.44 2011/03/30 02:43:18 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
|
|
|
|
#define ORD_HOMO_WW_DRL_ZIGZAG 13 |
#define ORD_HOMO_WW_DRL_ZIGZAG 13 |
|
|
int cmpdl_drl_zigzag(), cmpdl_homo_ww_drl_zigzag(); |
int cmpdl_drl_zigzag(), cmpdl_homo_ww_drl_zigzag(); |
|
int cmpdl_top_weight(); |
|
|
int (*cmpdl)()=cmpdl_revgradlex; |
int (*cmpdl)()=cmpdl_revgradlex; |
|
int (*cmpdl_tie_breaker)(); |
int (*primitive_cmpdl[3])() = {cmpdl_revgradlex,cmpdl_gradlex,cmpdl_lex}; |
int (*primitive_cmpdl[3])() = {cmpdl_revgradlex,cmpdl_gradlex,cmpdl_lex}; |
|
|
int do_weyl; |
int do_weyl; |
Line 118 int has_sfcoef_p(P f) |
|
Line 120 int has_sfcoef_p(P f) |
|
} |
} |
} |
} |
|
|
|
extern N *current_top_weight_vector; |
|
static int current_top_weight_len; |
|
|
void initd(struct order_spec *spec) |
void initd(struct order_spec *spec) |
{ |
{ |
|
int len,i; |
|
|
switch ( spec->id ) { |
switch ( spec->id ) { |
case 3: |
case 3: |
cmpdl = cmpdl_composite; |
cmpdl = cmpdl_composite; |
Line 165 void initd(struct order_spec *spec) |
|
Line 172 void initd(struct order_spec *spec) |
|
} |
} |
break; |
break; |
} |
} |
|
if ( current_top_weight_vector ) { |
|
cmpdl_tie_breaker = cmpdl; |
|
cmpdl = cmpdl_top_weight; |
|
for ( len = 0, i = 0; i < spec->nv; i++ ) |
|
if ( current_top_weight_vector[i] ) |
|
len = MAX(PL(current_top_weight_vector[i]),len); |
|
current_top_weight_len = len; |
|
} |
dp_current_spec = spec; |
dp_current_spec = spec; |
} |
} |
|
|
Line 642 void comm_muld(VL vl,DP p1,DP p2,DP *pr) |
|
Line 657 void comm_muld(VL vl,DP p1,DP p2,DP *pr) |
|
l = l1; |
l = l1; |
} |
} |
if ( l > wlen ) { |
if ( l > wlen ) { |
if ( w ) GC_free(w); |
if ( w ) GCFREE(w); |
w = (MP *)MALLOC(l*sizeof(MP)); |
w = (MP *)MALLOC(l*sizeof(MP)); |
wlen = l; |
wlen = l; |
} |
} |
Line 680 void comm_muld_trunc(VL vl,DP p1,DP p2,DL dl,DP *pr) |
|
Line 695 void comm_muld_trunc(VL vl,DP p1,DP p2,DL dl,DP *pr) |
|
l = l1; |
l = l1; |
} |
} |
if ( l > wlen ) { |
if ( l > wlen ) { |
if ( w ) GC_free(w); |
if ( w ) GCFREE(w); |
w = (MP *)MALLOC(l*sizeof(MP)); |
w = (MP *)MALLOC(l*sizeof(MP)); |
wlen = l; |
wlen = l; |
} |
} |
Line 814 void weyl_muld(VL vl,DP p1,DP p2,DP *pr) |
|
Line 829 void weyl_muld(VL vl,DP p1,DP p2,DP *pr) |
|
else { |
else { |
for ( m = BDY(p1), l = 0; m; m = NEXT(m), l++ ); |
for ( m = BDY(p1), l = 0; m; m = NEXT(m), l++ ); |
if ( l > wlen ) { |
if ( l > wlen ) { |
if ( w ) GC_free(w); |
if ( w ) GCFREE(w); |
w = (MP *)MALLOC(l*sizeof(MP)); |
w = (MP *)MALLOC(l*sizeof(MP)); |
wlen = l; |
wlen = l; |
} |
} |
Line 846 void weyl_muldm(VL vl,MP m0,DP p,DP *pr) |
|
Line 861 void weyl_muldm(VL vl,MP m0,DP p,DP *pr) |
|
else { |
else { |
for ( m = BDY(p), l = 0; m; m = NEXT(m), l++ ); |
for ( m = BDY(p), l = 0; m; m = NEXT(m), l++ ); |
if ( l > wlen ) { |
if ( l > wlen ) { |
if ( w ) GC_free(w); |
if ( w ) GCFREE(w); |
w = (MP *)MALLOC(l*sizeof(MP)); |
w = (MP *)MALLOC(l*sizeof(MP)); |
wlen = l; |
wlen = l; |
} |
} |
Line 858 void weyl_muldm(VL vl,MP m0,DP p,DP *pr) |
|
Line 873 void weyl_muldm(VL vl,MP m0,DP p,DP *pr) |
|
for ( i = 0, tlen = 1; i < n2; i++ ) |
for ( i = 0, tlen = 1; i < n2; i++ ) |
tlen *= d0->d[n2+i]+1; |
tlen *= d0->d[n2+i]+1; |
if ( tlen > rtlen ) { |
if ( tlen > rtlen ) { |
if ( tab ) GC_free(tab); |
if ( tab ) GCFREE(tab); |
if ( psum ) GC_free(psum); |
if ( psum ) GCFREE(psum); |
rtlen = tlen; |
rtlen = tlen; |
tab = (struct cdl *)MALLOC(rtlen*sizeof(struct cdl)); |
tab = (struct cdl *)MALLOC(rtlen*sizeof(struct cdl)); |
psum = (MP *)MALLOC(rtlen*sizeof(MP)); |
psum = (MP *)MALLOC(rtlen*sizeof(MP)); |
Line 921 void weyl_mulmm(VL vl,MP m0,MP m1,int n,struct cdl *rt |
|
Line 936 void weyl_mulmm(VL vl,MP m0,MP m1,int n,struct cdl *rt |
|
rtab[0].d = d; |
rtab[0].d = d; |
|
|
if ( rtablen > tmptablen ) { |
if ( rtablen > tmptablen ) { |
if ( tmptab ) GC_free(tmptab); |
if ( tmptab ) GCFREE(tmptab); |
tmptab = (struct cdl *)MALLOC(rtablen*sizeof(struct cdl)); |
tmptab = (struct cdl *)MALLOC(rtablen*sizeof(struct cdl)); |
tmptablen = rtablen; |
tmptablen = rtablen; |
} |
} |
Line 947 void weyl_mulmm(VL vl,MP m0,MP m1,int n,struct cdl *rt |
|
Line 962 void weyl_mulmm(VL vl,MP m0,MP m1,int n,struct cdl *rt |
|
continue; |
continue; |
} |
} |
if ( k+1 > tablen ) { |
if ( k+1 > tablen ) { |
if ( tab ) GC_free(tab); |
if ( tab ) GCFREE(tab); |
if ( ctab ) GC_free(ctab); |
if ( ctab ) GCFREE(ctab); |
tablen = k+1; |
tablen = k+1; |
tab = (struct cdl *)MALLOC(tablen*sizeof(struct cdl)); |
tab = (struct cdl *)MALLOC(tablen*sizeof(struct cdl)); |
ctab = (Q *)MALLOC(tablen*sizeof(Q)); |
ctab = (Q *)MALLOC(tablen*sizeof(Q)); |
Line 1089 void adddl(int n,DL d1,DL d2,DL *dr) |
|
Line 1104 void adddl(int n,DL d1,DL d2,DL *dr) |
|
DL dt; |
DL dt; |
int i; |
int i; |
|
|
if ( !d1->td ) |
*dr = dt = (DL)MALLOC_ATOMIC((n+1)*sizeof(int)); |
*dr = d2; |
dt->td = d1->td + d2->td; |
else if ( !d2->td ) |
for ( i = 0; i < n; i++ ) |
*dr = d1; |
dt->d[i] = d1->d[i]+d2->d[i]; |
else { |
|
*dr = dt = (DL)MALLOC_ATOMIC((n+1)*sizeof(int)); |
|
dt->td = d1->td + d2->td; |
|
for ( i = 0; i < n; i++ ) |
|
dt->d[i] = d1->d[i]+d2->d[i]; |
|
} |
|
} |
} |
|
|
/* d1 += d2 */ |
/* d1 += d2 */ |
Line 1640 int cmpdl_matrix(int n,DL d1,DL d2) |
|
Line 1649 int cmpdl_matrix(int n,DL d1,DL d2) |
|
return 0; |
return 0; |
} |
} |
|
|
|
int cmpdl_top_weight(int n,DL d1,DL d2) |
|
{ |
|
int *w; |
|
N sum,wm,wma,t; |
|
N *a; |
|
struct oN tn; |
|
int len,i,sgn,tsgn; |
|
int *t1,*t2; |
|
|
|
w = (int *)ALLOCA(n*sizeof(int)); |
|
len = current_top_weight_len+3; |
|
t1 = d1->d; t2 = d2->d; |
|
for ( i = 0; i < n; i++ ) w[i] = t1[i]-t2[i]; |
|
sum = (N)W_ALLOC(len); sgn = 0; |
|
wm = (N)W_ALLOC(len); |
|
wma = (N)W_ALLOC(len); |
|
a = current_top_weight_vector; |
|
for ( i = 0; i < n; i++ ) { |
|
if ( !a[i] || !w[i] ) continue; |
|
tn.p = 1; |
|
if ( w[i] > 0 ) { |
|
tn.b[0] = w[i]; tsgn = 1; |
|
} else { |
|
tn.b[0] = -w[i]; tsgn = -1; |
|
} |
|
_muln(a[i],&tn,wm); |
|
if ( !sgn ) { |
|
sgn = tsgn; |
|
t = wm; wm = sum; sum = t; |
|
} else if ( sgn == tsgn ) { |
|
_addn(sum,wm,wma); |
|
if ( !PL(wma) ) |
|
sgn = 0; |
|
t = wma; wma = sum; sum = t; |
|
} else { |
|
sgn *= _subn(sum,wm,wma); |
|
t = wma; wma = sum; sum = t; |
|
} |
|
} |
|
if ( sgn > 0 ) return 1; |
|
else if ( sgn < 0 ) return -1; |
|
else return (*cmpdl_tie_breaker)(n,d1,d2); |
|
} |
|
|
GeoBucket create_bucket() |
GeoBucket create_bucket() |
{ |
{ |
GeoBucket g; |
GeoBucket g; |
Line 1885 P separate_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
Line 1938 P separate_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
NEWNBM(t); NEWNBMBDY(t,d1); t->d = d1; t->c = (P)ONE; |
NEWNBM(t); NEWNBMBDY(t,d1); t->d = d1; t->c = (P)ONE; |
for ( i = 0; i < d1; i++ ) { |
for ( i = 0; i < d1; i++ ) { |
if ( NBM_GET(a->b,i+1) ) NBM_SET(t->b,i); |
if ( NBM_GET(a->b,i+1) ) NBM_SET(t->b,i); |
|
else NBM_CLR(t->b,i); |
|
} |
|
*ar = nbmtonbp(t); |
|
} |
|
|
|
return a->c; |
|
} |
|
|
|
/* a=c*rest*x -> a0= rest*x, ar=rest, at=x */ |
|
|
|
P separate_tail_nbm(NBM a,NBP *a0,NBP *ar,NBP *at) |
|
{ |
|
int i,d,d1; |
|
NBM t; |
|
|
|
if ( !(d=a->d) ) error("separate_tail_nbm : invalid argument"); |
|
|
|
if ( a0 ) { |
|
NEWNBM(t); t->d = a->d; t->b = a->b; t->c = (P)ONE; |
|
*a0 = nbmtonbp(t); |
|
} |
|
|
|
d1 = a->d-1; |
|
if ( at ) { |
|
NEWNBM(t); NEWNBMBDY(t,1); t->d = 1; t->c = (P)ONE; |
|
if ( NBM_GET(a->b,d1) ) NBM_SET(t->b,0); |
|
else NBM_CLR(t->b,0); |
|
*at = nbmtonbp(t); |
|
} |
|
|
|
if ( ar ) { |
|
NEWNBM(t); NEWNBMBDY(t,d1); t->d = d1; t->c = (P)ONE; |
|
for ( i = 0; i < d1; i++ ) { |
|
if ( NBM_GET(a->b,i) ) NBM_SET(t->b,i); |
else NBM_CLR(t->b,i); |
else NBM_CLR(t->b,i); |
} |
} |
*ar = nbmtonbp(t); |
*ar = nbmtonbp(t); |