[BACK]Return to _distm.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / engine

Diff for /OpenXM_contrib2/asir2000/engine/_distm.c between version 1.4 and 1.8

version 1.4, 2000/11/07 06:06:39 version 1.8, 2001/03/19 04:02:03
Line 45 
Line 45 
  * 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/_distm.c,v 1.3 2000/08/22 05:04:05 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/engine/_distm.c,v 1.7 2001/01/11 02:26:49 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "inline.h"  #include "inline.h"
   
 #define NV(p) ((p)->nv)  
 #define C(p) ((p)->c)  
 #if 0  
 #define ITOS(p) (((unsigned int)(p))>>1)  
 #define STOI(i) ((P)((((unsigned int)(i))<<1)|1))  
 #else  
 #define ITOS(p) (((unsigned int)(p))&0x7fffffff)  
 #define STOI(i) ((P)((unsigned int)(i)|0x80000000))  
 #endif  
   
 struct cdlm {  
         int c;  
         DL d;  
 };  
   
 extern int (*cmpdl)();  extern int (*cmpdl)();
 extern int do_weyl;  extern int do_weyl;
   
   void dpto_dp();
 void _dptodp();  void _dptodp();
 void adddl_dup();  void _adddl_dup();
 void adddl_destructive();  void adddl_destructive();
 void _mulmdm_dup();  void _mulmdm_dup();
 void _free_dp();  void _free_dp();
Line 85  DP _dp_free_list;
Line 71  DP _dp_free_list;
 DL _dl_free_list;  DL _dl_free_list;
 int current_dl_length;  int current_dl_length;
   
 #define _NEWDL_NOINIT(d,n) if ((n)!= current_dl_length){_dl_free_list=0; current_dl_length=(n);} if(!_dl_free_list)_DL_alloc(); (d)=_dl_free_list; _dl_free_list = *((DL *)_dl_free_list)  void _free_private_storage()
 #define _NEWDL(d,n) if ((n)!= current_dl_length){_dl_free_list=0; current_dl_length=(n);} if(!_dl_free_list)_DL_alloc(); (d)=_dl_free_list; _dl_free_list = *((DL *)_dl_free_list); bzero((d),(((n)+1)*sizeof(int)))  {
 #define _NEWMP(m) if(!_mp_free_list)_MP_alloc(); (m)=_mp_free_list; _mp_free_list = NEXT(_mp_free_list)          _mp_free_list = 0;
 #define _MKDP(n,m,d) if(!_dp_free_list)_DP_alloc(); (d)=_dp_free_list; _dp_free_list = (DP)BDY(_dp_free_list); (d)->nv=(n); BDY(d)=(m)          _dp_free_list = 0;
           _dl_free_list = 0;
           GC_gcollect();
   }
   
 #define _NEXTMP(r,c) \  
 if(!(r)){_NEWMP(r);(c)=(r);}else{_NEWMP(NEXT(c));(c)=NEXT(c);}  
   
 #define _NEXTMP2(r,c,s) \  
 if(!(r)){(c)=(r)=(s);}else{NEXT(c)=(s);(c)=(s);}  
   
 #define FREEDL(m) *((DL *)m)=_dl_free_list; _dl_free_list=(m)  
 #define FREEMP(m) NEXT(m)=_mp_free_list; _mp_free_list=(m)  
 #define FREEDP(m) BDY(m)=(MP)_dp_free_list; _dp_free_list=(m)  
   
   
 void _DL_alloc()  void _DL_alloc()
 {  {
         int *p;          int *p;
Line 109  void _DL_alloc()
Line 87  void _DL_alloc()
   
 /*      fprintf(stderr,"DL_alloc : %d \n",++DL_alloc_count); */  /*      fprintf(stderr,"DL_alloc : %d \n",++DL_alloc_count); */
         dl_len = (current_dl_length+1);          dl_len = (current_dl_length+1);
   #if defined(LONG_IS_64BIT)
           if ( dl_len & 1 )
                   dl_len += 1;
   #endif
         p = (int *)GC_malloc(128*dl_len*sizeof(int));          p = (int *)GC_malloc(128*dl_len*sizeof(int));
         for ( i = 0; i < 128; i++, p += dl_len ) {          for ( i = 0; i < 128; i++, p += dl_len ) {
                 *(DL *)p = _dl_free_list;                  *(DL *)p = _dl_free_list;
Line 167  DP p1,p2,*pr;
Line 149  DP p1,p2,*pr;
                                         if ( t ) {                                          if ( t ) {
                                                 _NEXTMP2(mr0,mr,s); C(mr) = STOI(t);                                                  _NEXTMP2(mr0,mr,s); C(mr) = STOI(t);
                                         } else {                                          } else {
                                                 FREEDL(s->dl); FREEMP(s);                                                  _FREEDL(s->dl); _FREEMP(s);
                                         }                                          }
                                         s = m2; m2 = NEXT(m2); FREEDL(s->dl); FREEMP(s);                                          s = m2; m2 = NEXT(m2); _FREEDL(s->dl); _FREEMP(s);
                                         break;                                          break;
                                 case 1:                                  case 1:
                                         s = m1; m1 = NEXT(m1); _NEXTMP2(mr0,mr,s);                                          s = m1; m1 = NEXT(m1); _NEXTMP2(mr0,mr,s);
Line 196  DP p1,p2,*pr;
Line 178  DP p1,p2,*pr;
                 _MKDP(NV(p1),mr0,*pr);                  _MKDP(NV(p1),mr0,*pr);
                 if ( *pr )                  if ( *pr )
                         (*pr)->sugar = MAX(p1->sugar,p2->sugar);                          (*pr)->sugar = MAX(p1->sugar,p2->sugar);
                 FREEDP(p1); FREEDP(p2);                  _FREEDP(p1); _FREEDP(p2);
         }          }
 }  }
   
Line 464  int rtablen;
Line 446  int rtablen;
 #if 0  #if 0
                 _comm_mulmd_tab(mod,n,rtab,curlen,tab,k+1,tmptab);                  _comm_mulmd_tab(mod,n,rtab,curlen,tab,k+1,tmptab);
                 for ( j = 0; j < curlen; j++ )                  for ( j = 0; j < curlen; j++ )
                         if ( rtab[j].d ) { FREEDL(rtab[j].d); }                          if ( rtab[j].d ) { _FREEDL(rtab[j].d); }
                 for ( j = 0; j <= min; j++ )                  for ( j = 0; j <= min; j++ )
                         if ( tab[j].d ) { FREEDL(tab[j].d); }                          if ( tab[j].d ) { _FREEDL(tab[j].d); }
                 curlen *= k+1;                  curlen *= k+1;
                 bcopy(tmptab,rtab,curlen*sizeof(struct cdlm));                  bcopy(tmptab,rtab,curlen*sizeof(struct cdlm));
 #else  #else
                 _comm_mulmd_tab_destructive(mod,n,rtab,curlen,tab,k+1);                  _comm_mulmd_tab_destructive(mod,n,rtab,curlen,tab,k+1);
                 for ( j = 0; j <= min; j++ )                  for ( j = 0; j <= min; j++ )
                         if ( tab[j].d ) { FREEDL(tab[j].d); }                          if ( tab[j].d ) { _FREEDL(tab[j].d); }
                 curlen *= k+1;                  curlen *= k+1;
 #endif  #endif
         }          }
Line 510  struct cdlm *rt;
Line 492  struct cdlm *rt;
                 for ( i = 0; i < n; i++, p++ ) {                  for ( i = 0; i < n; i++, p++ ) {
                         if ( t[i].c ) {                          if ( t[i].c ) {
                                 p->c = dmar(t[i].c,c,0,mod);                                  p->c = dmar(t[i].c,c,0,mod);
                                 adddl_dup(nv,t[i].d,d,&p->d);                                  _adddl_dup(nv,t[i].d,d,&p->d);
                         }                          }
                 }                  }
         }          }
Line 537  int n1;
Line 519  int n1;
                 for ( i = 0; i < n; i++, p++ ) {                  for ( i = 0; i < n; i++, p++ ) {
                         if ( t[i].c ) {                          if ( t[i].c ) {
                                 p->c = dmar(t[i].c,c,0,mod);                                  p->c = dmar(t[i].c,c,0,mod);
                                 adddl_dup(nv,t[i].d,d,&p->d);                                  _adddl_dup(nv,t[i].d,d,&p->d);
                         }                          }
                 }                  }
         }          }
Line 551  int n1;
Line 533  int n1;
                 }                  }
 }  }
   
 void _dp_red_mod_destructive(p1,p2,mod,rp)  void dlto_dl(d,dr)
 DP p1,p2;  DL d;
 int mod;  DL *dr;
 DP *rp;  
 {  {
         int i,n;          int i,n;
         DL d1,d2,d;          DL t;
         MP m;  
         DP t,s;  
         int c,c1;  
   
         n = p1->nv; d1 = BDY(p1)->dl; d2 = BDY(p2)->dl;          n = current_dl_length;
         _NEWDL(d,n); d->td = d1->td - d2->td;          _NEWDL(t,n); *dr = t;
           t->td = d->td;
         for ( i = 0; i < n; i++ )          for ( i = 0; i < n; i++ )
                 d->d[i] = d1->d[i]-d2->d[i];                  t->d[i] = d->d[i];
         c = invm(ITOS(BDY(p2)->c),mod); c1 = dmar(c,ITOS(BDY(p1)->c),0,mod);  
         _NEWMP(m); m->dl = d; m->c = STOI(mod-c1); NEXT(m) = 0;  
         _MKDP(n,m,s); s->sugar = d->td;  
         _mulmd_dup(mod,s,p2,&t); _free_dp(s);  
         _addmd_destructive(mod,p1,t,rp);  
 }  }
   
 void _dp_sp_mod_dup(p1,p2,mod,rp)  
 DP p1,p2;  
 int mod;  
 DP *rp;  
 {  
         int i,n,td;  
         int *w;  
         DL d1,d2,d;  
         MP m;  
         DP t,s,u;  
   
         n = p1->nv; d1 = BDY(p1)->dl; d2 = BDY(p2)->dl;  
         w = (int *)ALLOCA(n*sizeof(int));  
         for ( i = 0, td = 0; i < n; i++ ) {  
                 w[i] = MAX(d1->d[i],d2->d[i]); td += w[i];  
         }  
         _NEWDL(d,n); d->td = td - d1->td;  
         for ( i = 0; i < n; i++ )  
                 d->d[i] = w[i] - d1->d[i];  
         _NEWMP(m); m->dl = d; m->c = BDY(p2)->c; NEXT(m) = 0;  
         _MKDP(n,m,s); s->sugar = d->td; _mulmd_dup(mod,s,p1,&t); _free_dp(s);  
         _NEWDL(d,n); d->td = td - d2->td;  
         for ( i = 0; i < n; i++ )  
                 d->d[i] = w[i] - d2->d[i];  
         _NEWMP(m); m->dl = d; m->c = STOI(mod - ITOS(BDY(p1)->c)); NEXT(m) = 0;  
         _MKDP(n,m,s); s->sugar = d->td; _mulmd_dup(mod,s,p2,&u); _free_dp(s);  
         _addmd_destructive(mod,t,u,rp);  
 }  
   
 void _dltodl(d,dr)  void _dltodl(d,dr)
 DL d;  DL d;
 DL *dr;  DL *dr;
Line 616  DL *dr;
Line 561  DL *dr;
                 t->d[i] = d->d[i];                  t->d[i] = d->d[i];
 }  }
   
 void adddl_dup(n,d1,d2,dr)  void _adddl_dup(n,d1,d2,dr)
 int n;  int n;
 DL d1,d2;  DL d1,d2;
 DL *dr;  DL *dr;
Line 631  DL *dr;
Line 576  DL *dr;
                 dt->d[i] = d1->d[i]+d2->d[i];                  dt->d[i] = d1->d[i]+d2->d[i];
 }  }
   
 /* d1 += d2 */  
   
 void adddl_destructive(n,d1,d2)  
 int n;  
 DL d1,d2;  
 {  
         DL dt;  
         int i;  
   
         d1->td += d2->td;  
         for ( i = 0; i < n; i++ )  
                 d1->d[i] += d2->d[i];  
 }  
   
 void _free_dlarray(a,n)  void _free_dlarray(a,n)
 DL *a;  DL *a;
 int n;  int n;
 {  {
         int i;          int i;
   
         for ( i = 0; i < n; i++ ) { FREEDL(a[i]); }          for ( i = 0; i < n; i++ ) { _FREEDL(a[i]); }
 }  }
   
 void _free_dp(f)  void _free_dp(f)
Line 663  DP f;
Line 594  DP f;
                 return;                  return;
         m = f->body;          m = f->body;
         while ( m ) {          while ( m ) {
                 s = m; m = NEXT(m); FREEDL(s->dl); FREEMP(s);                  s = m; m = NEXT(m); _FREEDL(s->dl); _FREEMP(s);
         }          }
         FREEDP(f);          _FREEDP(f);
 }  }
   
   void dpto_dp(p,r)
   DP p;
   DP *r;
   {
           MP m,mr0,mr;
           DL t;
   
           if ( !p )
                   *r = 0;
           else {
                   /* XXX : dummy call to set current_dl_length */
                   _NEWDL_NOINIT(t,NV(p));
   
                   for ( m = BDY(p), mr0 = 0; m; m = NEXT(m) ) {
                           _NEXTMP(mr0,mr);
                           dlto_dl(m->dl,&mr->dl);
                           mr->c = m->c;
                   }
                   NEXT(mr) = 0;
                   _MKDP(p->nv,mr0,*r);
                   (*r)->sugar = p->sugar;
           }
   }
   
 void _dptodp(p,r)  void _dptodp(p,r)
 DP p;  DP p;
 DP *r;  DP *r;
Line 688  DP *r;
Line 643  DP *r;
         }          }
 }  }
   
 void _dp_nf_mod_destructive(b,g,ps,mod,full,rp)  /*
 NODE b;   * destructive merge of two list
 DP g;   *
 DP *ps;   * p1, p2 : list of DL
 int mod,full;   * return : a merged list
 DP *rp;   */
   
   NODE _symb_merge(m1,m2,n)
   NODE m1,m2;
   int n;
 {  {
         DP u,p,d,s,t;          NODE top,prev,cur,m,t;
         NODE l;  
         MP m,mr,mrd;  
         int sugar,psugar,n,h_reducible,i;  
   
         if ( !g ) {          if ( !m1 )
                 *rp = 0; return;                  return m2;
         }          else if ( !m2 )
         sugar = g->sugar;                  return m1;
         n = g->nv;          else {
         for ( d = 0; g; ) {                  switch ( (*cmpdl)(n,(DL)BDY(m1),(DL)BDY(m2)) ) {
                 for ( h_reducible = 0, l = b; l; l = NEXT(l) ) {                          case 0:
                         if ( dp_redble(g,p = ps[(int)BDY(l)]) ) {                                  top = m1; _FREEDL((DL)BDY(m2)); m = NEXT(m2);
                                 h_reducible = 1;  
                                 psugar = (BDY(g)->dl->td - BDY(p)->dl->td) + p->sugar;  
                                 _dp_red_mod_destructive(g,p,mod,&u); g = u;  
                                 sugar = MAX(sugar,psugar);  
                                 if ( !g ) {  
                                         if ( d )  
                                                 d->sugar = sugar;  
                                         _dptodp(d,rp); _free_dp(d); return;  
                                 }  
                                 break;                                  break;
                         }                          case 1:
                                   top = m1; m = m2;
                                   break;
                           case -1:
                                   top = m2; m = m1;
                                   break;
                 }                  }
                 if ( !h_reducible ) {                  prev = top; cur = NEXT(top);
                         /* head term is not reducible */                  /* BDY(prev) > BDY(m) always holds */
                         if ( !full ) {                  while ( cur && m ) {
                                 if ( g )                          switch ( (*cmpdl)(n,(DL)BDY(cur),(DL)BDY(m)) ) {
                                         g->sugar = sugar;                                  case 0:
                                 _dptodp(g,rp); _free_dp(g); return;                                          _FREEDL(BDY(m)); m = NEXT(m);
                         } else {                                          prev = cur; cur = NEXT(cur);
                                 m = BDY(g);                                          break;
                                 if ( NEXT(m) ) {                                  case 1:
                                         BDY(g) = NEXT(m); NEXT(m) = 0;                                          t = NEXT(cur); NEXT(cur) = m; m = t;
                                 } else {                                          prev = cur; cur = NEXT(cur);
                                         FREEDP(g); g = 0;                                          break;
                                 }                                  case -1:
                                 if ( d ) {                                          NEXT(prev) = m; m = cur;
                                         for ( mrd = BDY(d); NEXT(mrd); mrd = NEXT(mrd) );                                          prev = NEXT(prev); cur = NEXT(prev);
                                         NEXT(mrd) = m;                                          break;
                                 } else {  
                                         _MKDP(n,m,d);  
                                 }  
                         }                          }
                 }                  }
                   if ( !cur )
                           NEXT(prev) = m;
                   return top;
         }          }
         if ( d )  
                 d->sugar = sugar;  
         _dptodp(d,rp); _free_dp(d);  
 }  }

Legend:
Removed from v.1.4  
changed lines
  Added in v.1.8

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>