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

Diff for /OpenXM_contrib2/asir2000/engine/H.c between version 1.4 and 1.5

version 1.4, 2001/07/03 01:41:25 version 1.5, 2001/07/04 07:19:20
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/H.c,v 1.3 2000/08/22 05:04:04 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/engine/H.c,v 1.4 2001/07/03 01:41:25 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "inline.h"  #include "inline.h"
Line 418  LUM *gp;
Line 418  LUM *gp;
                 error("henmain2 : cannot happen(extgcd)");                  error("henmain2 : cannot happen(extgcd)");
 #endif  #endif
   
   #if 1
         fprintf(stderr,"bound=%d\n",bound);          fprintf(stderr,"bound=%d\n",bound);
   #endif
         for ( k = 1; k < bound; k++ ) {          for ( k = 1; k < bound; k++ ) {
   #if 1
                 fprintf(stderr,".");                  fprintf(stderr,".");
   #endif
                 /* at this point, f = gk*hk mod y^k */                  /* at this point, f = gk*hk mod y^k */
   
 #if 0  #if 0
Line 495  LUM *gp;
Line 499  LUM *gp;
                 for ( i = 0; i <= DEG(wh1); i++ )                  for ( i = 0; i <= DEG(wh1); i++ )
                         COEF(hk)[i][k] = COEF(wh1)[i];                          COEF(hk)[i][k] = COEF(wh1)[i];
         }          }
   #if 1
         fprintf(stderr,"\n");          fprintf(stderr,"\n");
   #endif
         *gp = gk;          *gp = gk;
         clearlum(n,bound,f);          clearlum(n,bound,f);
         DEG(f) = dh;          DEG(f) = dh;
Line 849  ML bqlist,cqlist,*listp;
Line 855  ML bqlist,cqlist,*listp;
                 for ( j = DEG(b[i]), pp = COEF(l[i]), px = COEF(b[i]); j >= 0; j-- )                  for ( j = DEG(b[i]), pp = COEF(l[i]), px = COEF(b[i]); j >= 0; j-- )
                         pp[j][0] = px[j];                          pp[j][0] = px[j];
         }          }
   #if 0
         fprintf(stderr,"bound=%d\n",bound);          fprintf(stderr,"bound=%d\n",bound);
   #endif
         for ( i = 1; i < bound; i++ ) {          for ( i = 1; i < bound; i++ ) {
   #if 0
                 fprintf(stderr,".");                  fprintf(stderr,".");
   #endif
                 mullum(mod,i+1,l[0],l[1],wb0);                  mullum(mod,i+1,l[0],l[1],wb0);
                 for ( j = 2; j < np; j++ ) {                  for ( j = 2; j < np; j++ ) {
                         mullum(mod,i+1,l[j],wb0,wb1);                          mullum(mod,i+1,l[j],wb0,wb1);
Line 880  ML bqlist,cqlist,*listp;
Line 890  ML bqlist,cqlist,*listp;
                 for ( j = n, px = COEF(wq0); j >= 0; j-- )                  for ( j = n, px = COEF(wq0); j >= 0; j-- )
                         px[j] = 0;                          px[j] = 0;
         }          }
   #if 0
         fprintf(stderr,"\n");          fprintf(stderr,"\n");
   #endif
 }  }
   
 static double M;  static double M;
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>