version 1.47, 2013/12/20 02:02:24 |
version 1.51, 2015/09/24 04:43:13 |
|
|
* 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.46 2013/11/05 02:55:03 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/dist.c,v 1.50 2014/10/10 09:02:25 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
|
|
Line 71 int (*cmpdl)()=cmpdl_revgradlex; |
|
Line 71 int (*cmpdl)()=cmpdl_revgradlex; |
|
int (*cmpdl_tie_breaker)(); |
int (*cmpdl_tie_breaker)(); |
int (*primitive_cmpdl[3])() = {cmpdl_revgradlex,cmpdl_gradlex,cmpdl_lex}; |
int (*primitive_cmpdl[3])() = {cmpdl_revgradlex,cmpdl_gradlex,cmpdl_lex}; |
|
|
|
Obj current_top_weight; |
|
int current_top_weight_len; |
|
|
int do_weyl; |
int do_weyl; |
|
|
int dp_nelim,dp_fcoeffs; |
int dp_nelim,dp_fcoeffs; |
Line 124 int has_sfcoef_p(P f) |
|
Line 127 int has_sfcoef_p(P f) |
|
} |
} |
} |
} |
|
|
extern N *current_top_weight_vector; |
extern Obj nd_top_weight; |
static int current_top_weight_len; |
|
|
|
|
void reset_top_weight() |
|
{ |
|
cmpdl = cmpdl_tie_breaker; |
|
cmpdl_tie_breaker = 0; |
|
nd_top_weight = 0; |
|
current_top_weight = 0; |
|
current_top_weight_len = 0; |
|
} |
|
|
void initd(struct order_spec *spec) |
void initd(struct order_spec *spec) |
{ |
{ |
int len,i; |
int len,i,k,row; |
|
Q **mat; |
|
|
switch ( spec->id ) { |
switch ( spec->id ) { |
case 3: |
case 3: |
Line 176 void initd(struct order_spec *spec) |
|
Line 188 void initd(struct order_spec *spec) |
|
} |
} |
break; |
break; |
} |
} |
if ( current_top_weight_vector ) { |
if ( current_top_weight ) { |
cmpdl_tie_breaker = cmpdl; |
cmpdl_tie_breaker = cmpdl; |
cmpdl = cmpdl_top_weight; |
cmpdl = cmpdl_top_weight; |
for ( len = 0, i = 0; i < spec->nv; i++ ) |
if ( OID(current_top_weight) == O_VECT ) { |
if ( current_top_weight_vector[i] ) |
mat = (Q **)&BDY((VECT)current_top_weight); |
len = MAX(PL(current_top_weight_vector[i]),len); |
row = 1; |
|
} else { |
|
mat = (Q **)BDY((MAT)current_top_weight); |
|
row = ((MAT)current_top_weight)->row; |
|
} |
|
for ( k = 0, len = 0; k < row; k++ ) |
|
for ( i = 0; i < spec->nv; i++ ) |
|
if ( mat[k][i] ) |
|
len = MAX(PL(NM(mat[k][i])),len); |
current_top_weight_len = len; |
current_top_weight_len = len; |
} |
} |
dp_current_spec = spec; |
dp_current_spec = spec; |
Line 845 void weyl_muld(VL vl,DP p1,DP p2,DP *pr) |
|
Line 865 void weyl_muld(VL vl,DP p1,DP p2,DP *pr) |
|
} |
} |
} |
} |
|
|
|
void actm(VL vl,int nv,MP m1,MP m2,DP *pr) |
|
{ |
|
DL d1,d2,d; |
|
int n2,i,j,k; |
|
Q jq,c,c1; |
|
MP m; |
|
P t; |
|
|
|
d1 = m1->dl; |
|
d2 = m2->dl; |
|
for ( i = 0; i < nv; i++ ) |
|
if ( d1->d[i] > d2->d[i] ) { |
|
*pr = 0; return; |
|
} |
|
NEWDL(d,nv); |
|
c = ONE; |
|
for ( i = 0; i < nv; i++ ) { |
|
for ( j = d2->d[i], k = d1->d[i]; k > 0; k--, j-- ) { |
|
STOQ(j,jq); mulq(c,jq,&c1); c = c1; |
|
} |
|
d->d[i] = d2->d[i]-d1->d[i]; |
|
} |
|
mulp(vl,C(m1),C(m2),&t); |
|
NEWMP(m); |
|
mulp(vl,(P)c,t,&C(m)); |
|
m->dl = d; |
|
MKDP(nv,m,*pr); |
|
} |
|
|
|
void weyl_actd(VL vl,DP p1,DP p2,DP *pr) |
|
{ |
|
int n; |
|
MP m1,m2; |
|
DP d,r,s; |
|
|
|
if ( !p1 || !p2 ) *pr = 0; |
|
else { |
|
n = NV(p1); |
|
r = 0; |
|
for ( m1 = BDY(p1); m1; m1 = NEXT(m1) ) |
|
for ( m2 = BDY(p2); m2; m2 = NEXT(m2) ) { |
|
actm(vl,n,m1,m2,&d); |
|
addd(vl,r,d,&s); r = s; |
|
} |
|
*pr = r; |
|
} |
|
} |
|
|
/* monomial * polynomial */ |
/* monomial * polynomial */ |
|
|
void weyl_muldm(VL vl,MP m0,DP p,DP *pr) |
void weyl_muldm(VL vl,MP m0,DP p,DP *pr) |
Line 1656 int cmpdl_top_weight(int n,DL d1,DL d2) |
|
Line 1724 int cmpdl_top_weight(int n,DL d1,DL d2) |
|
{ |
{ |
int *w; |
int *w; |
N sum,wm,wma,t; |
N sum,wm,wma,t; |
N *a; |
Q **mat; |
|
Q *a; |
struct oN tn; |
struct oN tn; |
int len,i,sgn,tsgn; |
int len,i,sgn,tsgn,row,k; |
int *t1,*t2; |
int *t1,*t2; |
|
|
w = (int *)ALLOCA(n*sizeof(int)); |
w = (int *)ALLOCA(n*sizeof(int)); |
Line 1668 int cmpdl_top_weight(int n,DL d1,DL d2) |
|
Line 1737 int cmpdl_top_weight(int n,DL d1,DL d2) |
|
sum = (N)W_ALLOC(len); sgn = 0; |
sum = (N)W_ALLOC(len); sgn = 0; |
wm = (N)W_ALLOC(len); |
wm = (N)W_ALLOC(len); |
wma = (N)W_ALLOC(len); |
wma = (N)W_ALLOC(len); |
a = current_top_weight_vector; |
if ( OID(current_top_weight) == O_VECT ) { |
for ( i = 0; i < n; i++ ) { |
mat = (Q **)&BDY((VECT)current_top_weight); |
if ( !a[i] || !w[i] ) continue; |
row = 1; |
tn.p = 1; |
} else { |
if ( w[i] > 0 ) { |
mat = (Q **)BDY((MAT)current_top_weight); |
tn.b[0] = w[i]; tsgn = 1; |
row = ((MAT)current_top_weight)->row; |
} 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; |
for ( k = 0; k < row; k++ ) { |
else if ( sgn < 0 ) return -1; |
a = mat[k]; |
else return (*cmpdl_tie_breaker)(n,d1,d2); |
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(NM(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; |
|
} |
|
return (*cmpdl_tie_breaker)(n,d1,d2); |
} |
} |
|
|
GeoBucket create_bucket() |
GeoBucket create_bucket() |