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

Diff for /OpenXM_contrib2/asir2000/engine/nd.c between version 1.187 and 1.190

version 1.187, 2010/04/26 01:06:22 version 1.190, 2010/07/14 04:36:59
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.186 2010/04/23 07:35:44 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.189 2010/05/06 09:22:09 noro Exp $ */
   
 #include "nd.h"  #include "nd.h"
   
Line 2316  ND_pairs nd_newpairs( NODE g, int t )
Line 2316  ND_pairs nd_newpairs( NODE g, int t )
                                 if ( t0 >= min && t0 <= max && t >= min && t <= max )                                  if ( t0 >= min && t0 <= max && t >= min && t <= max )
                                         break;                                          break;
                         }                          }
                         if ( nd_gbblock[i] >= 0 ) {                          if ( nd_gbblock[i] >= 0 )
                                 fprintf(stderr,"(%d,%d)",t0,t);  
                                 continue;                                  continue;
                         }  
                 }                  }
         NEXTND_pairs(r0,r);          NEXTND_pairs(r0,r);
         r->i1 = (long)BDY(h);          r->i1 = (long)BDY(h);
Line 6602  void nd_det(int mod,MAT f,P *rp)
Line 6600  void nd_det(int mod,MAT f,P *rp)
                 bucket = create_pbucket();                  bucket = create_pbucket();
                 if ( mi[k] ) {                  if ( mi[k] ) {
                     nmv = BDY(mjj); len = LEN(mjj);                      nmv = BDY(mjj); len = LEN(mjj);
                                         fprintf(stderr,"len=%d\n",len);  
                     for ( a = 0; a < len; a++, NMV_ADV(nmv) ) {                      for ( a = 0; a < len; a++, NMV_ADV(nmv) ) {
                                                 fprintf(stderr,".");  
                         u = ndv_mul_nmv_trunc(mod,nmv,mi[k],DL(BDY(d)));                          u = ndv_mul_nmv_trunc(mod,nmv,mi[k],DL(BDY(d)));
                         add_pbucket(mod,bucket,u);                          add_pbucket(mod,bucket,u);
                                                 if ( !(a%1000) )  
                                                         fprintf(stderr,"%d\n",a);  
                     }                      }
                 }                  }
                 if ( mj[k] && mij ) {                  if ( mj[k] && mij ) {
                     nmv = BDY(mij); len = LEN(mij);                      nmv = BDY(mij); len = LEN(mij);
                                         fprintf(stderr,"len=%d\n",len);  
                     for ( a = 0; a < len; a++, NMV_ADV(nmv) ) {                      for ( a = 0; a < len; a++, NMV_ADV(nmv) ) {
                                                 fprintf(stderr,".");  
                         u = ndv_mul_nmv_trunc(mod,nmv,mj[k],DL(BDY(d)));                          u = ndv_mul_nmv_trunc(mod,nmv,mj[k],DL(BDY(d)));
                         add_pbucket(mod,bucket,u);                          add_pbucket(mod,bucket,u);
                                                 if ( !(a%1000) )  
                                                         fprintf(stderr,"%d\n",a);  
                     }                      }
                 }                  }
                 u = nd_quo(mod,bucket,d);                  u = nd_quo(mod,bucket,d);
Line 6969  void parse_nd_option(NODE opt)
Line 6959  void parse_nd_option(NODE opt)
                         if ( !value || OID(value) != O_LIST )                          if ( !value || OID(value) != O_LIST )
                                 error("nd_* : invalid value for gbblock option");                                  error("nd_* : invalid value for gbblock option");
                         u = BDY((LIST)value);                          u = BDY((LIST)value);
             nd_gbblock = MALLOC(2*length(t)+1);              nd_gbblock = MALLOC((2*length(u)+1)*sizeof(int));
                         for ( i = 0; u; u = NEXT(u) ) {                          for ( i = 0; u; u = NEXT(u) ) {
                                 p = BDY((LIST)BDY(u));                                  p = BDY((LIST)BDY(u));
                                 s = nd_gbblock[i++] = QTOS((Q)BDY(p));                                  s = nd_gbblock[i++] = QTOS((Q)BDY(p));
Legend:
Removed from v.1.187  
changed lines
  Added in v.1.190
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>