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

Diff for /OpenXM_contrib2/asir2000/builtin/gr.c between version 1.27 and 1.28

version 1.27, 2001/09/11 08:56:47 version 1.28, 2001/09/13 03:04:27
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/builtin/gr.c,v 1.26 2001/09/11 01:30:31 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/gr.c,v 1.27 2001/09/11 08:56:47 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
Line 150  double get_rtime();
Line 150  double get_rtime();
 void _dpmod_to_vect(DP,DL *,int *);  void _dpmod_to_vect(DP,DL *,int *);
 void dp_to_vect(DP,DL *,Q *);  void dp_to_vect(DP,DL *,Q *);
 NODE dp_dllist(DP f);  NODE dp_dllist(DP f);
   DLBUCKET dp_dllist_bucket(DP f);
 NODE symb_merge(NODE,NODE,int),_symb_merge(NODE,NODE,int);  NODE symb_merge(NODE,NODE,int),_symb_merge(NODE,NODE,int);
   DLBUCKET symb_merge_bucket(DLBUCKET,DLBUCKET,int);
 extern int dp_nelim;  extern int dp_nelim;
 extern int dp_fcoeffs;  extern int dp_fcoeffs;
 static DP *ps,*psm;  static DP *ps,*psm;
Line 304  DP f;
Line 306  DP f;
         return mp0;          return mp0;
 }  }
   
   void print_dlbucket(d,nv)
   DLBUCKET d;
   int nv;
   {
           int i;
           NODE n;
   
           for ( ; d; d = NEXT(d) ) {
                   fprintf(stderr,"td = %d\n",d->td);
                   for ( n = BDY(d); n; n = NEXT(n) ) {
                           fprintf(stderr,"<");
                           for ( i = 0; i < nv; i++ ) {
                                   fprintf(stderr,"%d",((DL)BDY(n))->d[i]);
                                   if ( i != nv-1 )
                                           fprintf(stderr," ");
                           }
                           fprintf(stderr,">");
                   }
                   fprintf(stderr,"\n");
           }
   }
   
   DLBUCKET dp_dllist_bucket(f)
   DP f;
   {
           MP m;
           NODE n,n0;
           DLBUCKET d,d0;
           int td;
   
           if ( !f )
                   return 0;
           d0 = 0;
           m = BDY(f);
           do {
                   NEXTDLBUCKET(d0,d);
                   n0 = 0;
                   d->td = td = m->dl->td;
                   do {
                           NEXTNODE(n0,n);
                           BDY(n) = (pointer)m->dl;
                           m = NEXT(m);
                   } while ( m && m->dl->td == td );
                   NEXT(n) = 0;
                   BDY(d) = n0;
           } while ( m );
           NEXT(d) = 0;
           return d0;
   }
   
 void pdl(f)  void pdl(f)
 NODE f;  NODE f;
 {  {
Line 694  int m;
Line 746  int m;
 {  {
         int i,j,k,nh,row,col,nv;          int i,j,k,nh,row,col,nv;
         NODE r,g,gall;          NODE r,g,gall;
         NODE s,s0;          NODE sb;
           DLBUCKET s,s0,s1;
         DP_pairs d,dm,dr,t;          DP_pairs d,dm,dr,t;
         DP h,nf,f1,f2,f21,f21r,sp,sp1,sd,sdm,tdp;          DP h,nf,f1,f2,f21,f21r,sp,sp1,sd,sdm,tdp;
         MP mp,mp0;          MP mp,mp0;
Line 729  int m;
Line 782  int m;
                         _dp_sp_mod(ps[t->dp1],ps[t->dp2],m,&sp);                          _dp_sp_mod(ps[t->dp1],ps[t->dp2],m,&sp);
                         if ( sp ) {                          if ( sp ) {
                                 MKNODE(bt,sp,blist); blist = bt;                                  MKNODE(bt,sp,blist); blist = bt;
                                 s0 = symb_merge(s0,dp_dllist(sp),nv);                                  s0 = symb_merge_bucket(s0,dp_dllist_bucket(sp),nv);
   /*                              print_dlbucket(s0,nv); */
                         }                          }
                 }                  }
                 /* s0 : all the terms appeared in symbolic reduction */                  /* s0 : all the terms appeared in symbolic reduction */
 #if 0  #if 0
                 for ( s = s0, nred = 0; s; s = NEXT(s) ) {                  for ( s = s0, nred = 0; s; s = NEXT(s) ) {
                         for ( j = psn-1; j >= 0; j-- )                          sb = BDY(s);
                                 if ( _dl_redble(BDY(ps[j])->dl,BDY(s),nv) )                          for ( ; sb; sb = NEXT(sb) ) {
                                         break;                                  for ( j = psn-1; j >= 0; j-- )
                         if ( j >= 0 ) {                                          if ( _dl_redble(BDY(ps[j])->dl,BDY(sb),nv) )
                                 dltod(BDY(s),nv,&tdp);                                                  break;
                                 dp_subd(tdp,ps[j],&sd);                                  if ( j >= 0 ) {
                                 for ( k = 0, i = 0; k < nv; k++ )                                          dltod(BDY(sb),nv,&tdp);
                                         if ( BDY(sd)->dl->d[k] )                                          dp_subd(tdp,ps[j],&sd);
                                                 i++;                                          for ( k = 0, i = 0; k < nv; k++ )
                                 fprintf(stderr,"%c ",i<=1 ? 'o' : 'x');                                                  if ( BDY(sd)->dl->d[k] )
                                 _dp_mod(sd,m,0,&sdm);                                                          i++;
                                 mulmd_dup(m,sdm,ps[j],&f2);                                          fprintf(stderr,"%c ",i<=1 ? 'o' : 'x');
                                 MKNODE(bt,f2,blist); blist = bt;                                          _dp_mod(sd,m,0,&sdm);
                                 s = symb_merge(s,dp_dllist(f2),nv);                                          mulmd_dup(m,sdm,ps[j],&f2);
                                 nred++;                                          MKNODE(bt,f2,blist); blist = bt;
                                           /* merge the highest degree part into sb */
                                           s1 = dp_dllist_bucket(f2);
                                           symb_merge(sb,BDY(s1),nv);
                                           /* merge the rest into s */
                                           symb_merge_bucket(s,NEXT(s1),nv);
                                           nred++;
                                   }
                         }                          }
                 }                  }
 #else  #else
                 for ( s = s0, nred = 0; s; s = NEXT(s) ) {                  for ( s = s0, nred = 0; s; s = NEXT(s) ) {
                         for ( r = gall; r; r = NEXT(r) )                          sb = BDY(s);
                                 if ( _dl_redble(BDY(ps[(int)BDY(r)])->dl,BDY(s),nv) )                          for ( ; sb; sb = NEXT(sb) ) {
                                         break;                                  for ( r = gall; r; r = NEXT(r) )
                         if ( r ) {                                          if ( _dl_redble(BDY(ps[(int)BDY(r)])->dl,BDY(sb),nv) )
                                 dltod(BDY(s),nv,&tdp);                                                  break;
                                 dp_subd(tdp,ps[(int)BDY(r)],&sd);                                  if ( r ) {
                                 _dp_mod(sd,m,0,&sdm);                                          dltod(BDY(sb),nv,&tdp);
                                 mulmd_dup(m,sdm,ps[(int)BDY(r)],&f2);                                          dp_subd(tdp,ps[(int)BDY(r)],&sd);
                                 MKNODE(bt,f2,blist); blist = bt;                                          _dp_mod(sd,m,0,&sdm);
                                 s = symb_merge(s,dp_dllist(f2),nv);                                          mulmd_dup(m,sdm,ps[(int)BDY(r)],&f2);
                                 nred++;                                          MKNODE(bt,f2,blist); blist = bt;
                                           /* merge the highest degree part into sb */
                                           s1 = dp_dllist_bucket(f2);
                                           symb_merge(sb,BDY(s1),nv);
                                           /* merge the rest into s */
                                           symb_merge_bucket(s,NEXT(s1),nv);
                                           nred++;
                                   }
                         }                          }
                 }                  }
 #endif  #endif
Line 781  int m;
Line 849  int m;
                         ht[i] = BDY((DP)BDY(r))->dl;                          ht[i] = BDY((DP)BDY(r))->dl;
   
                 /* col = number of all terms */                  /* col = number of all terms */
                 for ( s = s0, col = 0; s; s = NEXT(s), col++ );                  for ( s = s0, col = 0; s; s = NEXT(s) )
                           for ( sb = BDY(s); sb; sb = NEXT(sb) )
                                   col++;
   
                 /* head terms of all terms */                  /* head terms of all terms */
                 at = (DL *)MALLOC(col*sizeof(DL));                  at = (DL *)MALLOC(col*sizeof(DL));
                 for ( s = s0, i = 0; i < col; s = NEXT(s), i++ )                  for ( s = s0, i = 0; i < col; s = NEXT(s) )
                         at[i] = (DL)BDY(s);                          for ( sb = BDY(s); sb; sb = NEXT(sb), i++ )
                                   at[i] = (DL)BDY(sb);
   
                 /* store coefficients separately in spmat and redmat */                  /* store coefficients separately in spmat and redmat */
                 nsp = row-nred;                  nsp = row-nred;
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.28
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>