version 1.21, 2020/02/05 04:56:10 |
version 1.24, 2020/10/06 06:31:19 |
|
|
* 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/asir2018/engine/dist.c,v 1.20 2020/02/03 05:51:52 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2018/engine/dist.c,v 1.23 2020/06/19 10:18:13 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
|
|
Line 591 void _addtodl(int n,DL d1,DL d2) |
|
Line 591 void _addtodl(int n,DL d1,DL d2) |
|
d2->d[i] += d1->d[i]; |
d2->d[i] += d1->d[i]; |
} |
} |
|
|
|
void _subfromdl(int n,DL d1,DL d2) |
|
{ |
|
int i; |
|
|
|
d2->td -= d1->td; |
|
for ( i = 0; i < n; i++ ) |
|
d2->d[i] -= d1->d[i]; |
|
} |
|
|
void _copydl(int n,DL d1,DL d2) |
void _copydl(int n,DL d1,DL d2) |
{ |
{ |
int i; |
int i; |
Line 2755 DPM dpm_sp_nf_asir(VECT psv,int i,int j,DPM *nf); |
|
Line 2764 DPM dpm_sp_nf_asir(VECT psv,int i,int j,DPM *nf); |
|
void dpm_sort(DPM p,DPM *r); |
void dpm_sort(DPM p,DPM *r); |
|
|
extern int DP_Multiple; |
extern int DP_Multiple; |
|
extern int DP_Print; |
|
|
void dpm_nf_z(NODE b,DPM g,VECT psv,int full,int multiple,DPM *rp); |
void dpm_nf_z(NODE b,DPM g,VECT psv,int full,int multiple,DPM *rp); |
NODE dpm_sort_list(NODE l); |
NODE dpm_sort_list(NODE l); |
Line 3039 int compdp_lex(DP *a,DP *b) |
|
Line 3049 int compdp_lex(DP *a,DP *b) |
|
return -cmpdl_lex(compdp_nv,BDY(*a)->dl,BDY(*b)->dl); |
return -cmpdl_lex(compdp_nv,BDY(*a)->dl,BDY(*b)->dl); |
} |
} |
|
|
DMMstack_array dpm_schreyer_frame(NODE g) |
DMMstack_array dpm_schreyer_frame(NODE g,int lex) |
{ |
{ |
LIST l; |
LIST l; |
NODE nd,in,b0,b,n1,n2,n3,t; |
NODE nd,in,b0,b,n1,n2,n3,t; |
Line 3087 DMMstack_array dpm_schreyer_frame(NODE g) |
|
Line 3097 DMMstack_array dpm_schreyer_frame(NODE g) |
|
} |
} |
s = s1; |
s = s1; |
level++; |
level++; |
|
if ( DP_Print ) printf("level=%d,len=%d\n",level,n); |
|
|
|
|
/* create new list */ |
/* create new list */ |
MKVECT(psv,n+1); |
MKVECT(psv,n+1); |
ps = (DPM *)BDY(psv); |
ps = (DPM *)BDY(psv); |
Line 3122 DMMstack_array dpm_schreyer_frame(NODE g) |
|
Line 3132 DMMstack_array dpm_schreyer_frame(NODE g) |
|
} |
} |
if ( h ) m[p2] = h; |
if ( h ) m[p2] = h; |
} |
} |
#if 0 |
if ( lex ) { |
// compress m to m2 |
// compress m to m2 |
for ( j = 0, n2 = NEXT(n1); n2; n2 = NEXT(n2) ) { |
for ( j = 0, n2 = NEXT(n1); n2; n2 = NEXT(n2) ) { |
p2 = (long)BDY(n2); |
p2 = (long)BDY(n2); |
if ( m[p2] ) m2[j++] = m[p2]; |
if ( m[p2] ) m2[j++] = m[p2]; |
} |
} |
qsort(m2,j,sizeof(DP),(int (*)(const void *,const void *))compdp_lex); |
qsort(m2,j,sizeof(DP),(int (*)(const void *,const void *))compdp_lex); |
for ( k = 0; k < j; k++ ) { |
for ( k = 0; k < j; k++ ) { |
NEWDMM(dmm); dmm->dl = BDY(m2[k])->dl; dmm->pos = p1; dmm->c = (Obj)ONE; |
NEWDMM(dmm); dmm->dl = BDY(m2[k])->dl; dmm->pos = p1; dmm->c = (Obj)ONE; |
MKDPM(nv,dmm,dpm); |
|
NEXTNODE(b0,b); BDY(b) = (pointer)dpm; |
|
} |
|
#else |
|
for ( n2 = NEXT(n1); n2; n2 = NEXT(n2) ) { |
|
p2 = (long)BDY(n2); |
|
if ( m[p2] ) { |
|
NEWDMM(dmm); dmm->dl = BDY(m[p2])->dl; dmm->pos = p1; dmm->c = (Obj)ONE; |
|
MKDPM(nv,dmm,dpm); |
MKDPM(nv,dmm,dpm); |
NEXTNODE(b0,b); BDY(b) = (pointer)dpm; |
NEXTNODE(b0,b); BDY(b) = (pointer)dpm; |
} |
} |
|
} else { |
|
for ( n2 = NEXT(n1); n2; n2 = NEXT(n2) ) { |
|
p2 = (long)BDY(n2); |
|
if ( m[p2] ) { |
|
NEWDMM(dmm); dmm->dl = BDY(m[p2])->dl; dmm->pos = p1; dmm->c = (Obj)ONE; |
|
MKDPM(nv,dmm,dpm); |
|
NEXTNODE(b0,b); BDY(b) = (pointer)dpm; |
|
} |
|
} |
} |
} |
#endif |
|
} |
} |
} |
} |
if ( !b0 ) { |
if ( !b0 ) { |
Line 3210 int compdmm_schreyer(int n,DMM m1,DMM m2) |
|
Line 3220 int compdmm_schreyer(int n,DMM m1,DMM m2) |
|
else return 0; |
else return 0; |
} |
} |
} |
} |
|
/* XXX */ |
|
return 0; |
} |
} |
|
|
int compdmm_schreyer_old(int n,DMM m1,DMM m2) |
int compdmm_schreyer_old(int n,DMM m1,DMM m2) |
Line 3292 int compdmm(int n,DMM m1,DMM m2) |
|
Line 3304 int compdmm(int n,DMM m1,DMM m2) |
|
else return (*cmpdl)(n,m1->dl,m2->dl); |
else return (*cmpdl)(n,m1->dl,m2->dl); |
default: |
default: |
error("compdmm : invalid dpm_ordtype"); |
error("compdmm : invalid dpm_ordtype"); |
|
return 0; |
} |
} |
} |
} |
|
|