| version 1.9, 2019/09/19 06:29:47 |
version 1.10, 2019/10/11 03:45:57 |
|
|
| * 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.8 2019/09/13 02:04:42 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2018/engine/dist.c,v 1.9 2019/09/19 06:29:47 noro Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| |
|
| Line 3157 int compdpm(VL vl,DPM p1,DPM p2) |
|
| Line 3157 int compdpm(VL vl,DPM p1,DPM p2) |
|
| DPM dpm_eliminate_term(DPM a,DPM p,Obj c,int pos) |
DPM dpm_eliminate_term(DPM a,DPM p,Obj c,int pos) |
| { |
{ |
| MP d0,d; |
MP d0,d; |
| |
DL dl; |
| DMM m; |
DMM m; |
| DP f; |
DP f; |
| DPM a1,p1,r; |
DPM a1,p1,r; |
| Line 3165 DPM dpm_eliminate_term(DPM a,DPM p,Obj c,int pos) |
|
| Line 3166 DPM dpm_eliminate_term(DPM a,DPM p,Obj c,int pos) |
|
| d0 = 0; |
d0 = 0; |
| for ( m = BDY(a); m; m = NEXT(m) ) |
for ( m = BDY(a); m; m = NEXT(m) ) |
| if ( m->pos == pos ) { |
if ( m->pos == pos ) { |
| NEXTMP(d0,d); d->dl = m->dl; arf_chsgn(m->c,&d->c); |
NEXTMP(d0,d); |
| |
arf_chsgn(m->c,&d->c); |
| |
if ( !dp_current_spec ) d->dl = m->dl; |
| |
else { |
| |
NEWDL(dl,NV(a)); |
| |
_copydl(NV(a),m->dl,dl); |
| |
dl->td -= dp_current_spec->module_top_weight[pos-1]; |
| |
d->dl = dl; |
| |
} |
| } |
} |
| if ( d0 ) { |
if ( d0 ) { |
| NEXT(d) = 0; MKDP(NV(a),d0,f); |
NEXT(d) = 0; MKDP(NV(a),d0,f); |
| Line 3198 DPM dpm_compress(DPM p,int *tab) |
|
| Line 3207 DPM dpm_compress(DPM p,int *tab) |
|
| } |
} |
| |
|
| // input : s, s = syz(m) output simplified s, m |
// input : s, s = syz(m) output simplified s, m |
| void dpm_simplify_syz(LIST s,LIST m,LIST *s1,LIST *m1) |
void dpm_simplify_syz(LIST s,LIST m,LIST *s1,LIST *m1,LIST *w1) |
| { |
{ |
| int lm,ls,i,j,k,pos,nv; |
int lm,ls,i,j,k,pos,nv; |
| DPM *am,*as; |
DPM *am,*as; |
| DPM p; |
DPM p; |
| DMM d; |
DMM d; |
| Obj c; |
Obj c; |
| int *tab,*dd; |
Z q; |
| |
int *tab,*dd,*new_w; |
| NODE t,t1; |
NODE t,t1; |
| |
|
| lm = length(BDY(m)); |
lm = length(BDY(m)); |
| Line 3241 void dpm_simplify_syz(LIST s,LIST m,LIST *s1,LIST *m1) |
|
| Line 3251 void dpm_simplify_syz(LIST s,LIST m,LIST *s1,LIST *m1) |
|
| tab = (int *)MALLOC((lm+1)*sizeof(int)); |
tab = (int *)MALLOC((lm+1)*sizeof(int)); |
| for ( j = 0, i = 1; i <= lm; i++ ) { |
for ( j = 0, i = 1; i <= lm; i++ ) { |
| if ( am[i] ) { j++; tab[i] = j; } |
if ( am[i] ) { j++; tab[i] = j; } |
| else tab[i] = 0; |
else { tab[i] = 0; } |
| } |
} |
| t = 0; |
t = 0; |
| for ( i = ls-1; i >= 0; i-- ) |
for ( i = ls-1; i >= 0; i-- ) |
| Line 3250 void dpm_simplify_syz(LIST s,LIST m,LIST *s1,LIST *m1) |
|
| Line 3260 void dpm_simplify_syz(LIST s,LIST m,LIST *s1,LIST *m1) |
|
| MKNODE(t1,(pointer)p,t); t = t1; |
MKNODE(t1,(pointer)p,t); t = t1; |
| } |
} |
| MKLIST(*s1,t); |
MKLIST(*s1,t); |
| |
|
| |
if ( dp_current_spec->module_rank ) { |
| |
new_w = (int *)MALLOC(j*sizeof(int)); |
| |
for ( j = 0, i = 1; i <= lm; i++ ) |
| |
if ( tab[i] ) { new_w[j++] = dp_current_spec->module_top_weight[i-1]; } |
| |
t = 0; |
| |
for ( i = j-1; i >= 0; i-- ) { |
| |
STOZ(new_w[i],q); |
| |
MKNODE(t1,q,t); t = t1; |
| |
} |
| |
} else |
| |
t = 0; |
| |
MKLIST(*w1,t); |
| |
|
| t = 0; |
t = 0; |
| for ( i = lm; i >= 1; i-- ) |
for ( i = lm; i >= 1; i-- ) |
| if ( am[i] ) { |
if ( am[i] ) { |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to dist.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2018 / engine