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

Diff for /OpenXM_contrib2/asir2000/engine/Hgfs.c between version 1.1 and 1.2

version 1.1, 2001/06/20 09:32:13 version 1.2, 2001/06/21 07:47:02
Line 1 
Line 1 
 /* $OpenXM$ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.1 2001/06/20 09:32:13 noro Exp $ */
   
 #include "ca.h"  #include "ca.h"
   
Line 128  void randsfum(d,p)
Line 128  void randsfum(d,p)
 int d;  int d;
 UM p;  UM p;
 {  {
         unsigned int n;  
         int i;          int i;
   
         n = ((unsigned int)random()) % d; DEG(p) = n; COEF(p)[n] = _onesf();          for ( i = 0; i < d; i++ )
         for ( i = 0; i < (int)n; i++ )  
                 COEF(p)[i] = _randomsf();                  COEF(p)[i] = _randomsf();
           for ( i = d-1; i >= 0 && !COEF(p)[i]; i-- );
           p->d = i;
 }  }
   
 void pwrmodsfum(p,e,f,pr)  void pwrmodsfum(p,e,f,pr)
Line 166  UM p,f,pr;
Line 166  UM p,f,pr;
         }          }
 }  }
   
 void spwrum0sf(m,f,e,r)  void spwrsfum(m,f,e,r)
 UM f,m,r;  UM f,m,r;
 N e;  N e;
 {  {
Line 180  N e;
Line 180  N e;
                 cpyum(f,r);                  cpyum(f,r);
         else {          else {
                 a = divin(e,2,&e1);                  a = divin(e,2,&e1);
                 t = W_UMALLOC(2*DEG(m)); spwrum0sf(m,f,e1,t);                  t = W_UMALLOC(2*DEG(m)); spwrsfum(m,f,e1,t);
                 s = W_UMALLOC(2*DEG(m)); q = W_UMALLOC(2*DEG(m));                  s = W_UMALLOC(2*DEG(m)); q = W_UMALLOC(2*DEG(m));
                 mulsfum(t,t,s); DEG(s) = divsfum(s,m,q);                  mulsfum(t,t,s); DEG(s) = divsfum(s,m,q);
                 if ( a ) {                  if ( a ) {
Line 190  N e;
Line 190  N e;
         }          }
 }  }
   
   void tracemodsfum(m,f,e,r)
   UM f,m,r;
   int e;
   {
           UM t,s,q,u;
           int i;
   
           q = W_UMALLOC(2*DEG(m)+DEG(f)); /* XXX */
           t = W_UMALLOC(2*DEG(m));
           s = W_UMALLOC(2*DEG(m));
           u = W_UMALLOC(2*DEG(m));
           DEG(f) = divsfum(f,m,q);
           cpyum(f,s);
           cpyum(f,t);
           for ( i = 1; i < e; i++ ) {
                   mulsfum(t,t,u);
                   DEG(u) = divsfum(u,m,q); cpyum(u,t);
                   addsfum(t,s,u); cpyum(u,s);
           }
           cpyum(s,r);
   }
   
 void make_qmatsf(p,tab,mp)  void make_qmatsf(p,tab,mp)
 UM p;  UM p;
 UM *tab;  UM *tab;
Line 469  int d;
Line 491  int d;
         UM t,s,u,w,g,o;          UM t,s,u,w,g,o;
         N n1,n2,n3,n4,n5;          N n1,n2,n3,n4,n5;
         UM *b;          UM *b;
         int n,m,i,q;          int n,m,i,q,ed;
   
         if ( DEG(f) == d ) {          if ( DEG(f) == d ) {
                 r[0] = UMALLOC(d); cpyum(f,r[0]);                  r[0] = UMALLOC(d); cpyum(f,r[0]);
Line 483  int d;
Line 505  int d;
                 w = W_UMALLOC(DEG(f)); g = W_UMALLOC(DEG(f));                  w = W_UMALLOC(DEG(f)); g = W_UMALLOC(DEG(f));
                 o = W_UMALLOC(0); DEG(o) = 0; COEF(o)[0] = _onesf();                  o = W_UMALLOC(0); DEG(o) = 0; COEF(o)[0] = _onesf();
                 q = field_order_sf();                  q = field_order_sf();
                 STON(q,n1); pwrn(n1,d,&n2); subn(n2,ONEN,&n3);                  if ( q % 2 ) {
                 STON(2,n4); divsn(n3,n4,&n5);                          STON(q,n1); pwrn(n1,d,&n2); subn(n2,ONEN,&n3);
                           STON(2,n4); divsn(n3,n4,&n5);
                   } else
                           ed = d*extdeg_sf();
                 while ( 1 ) {                  while ( 1 ) {
                         randsfum(2*d,t); spwrum0sf(f,t,n5,s);                          randsfum(2*d,t);
                         subsfum(s,o,u); cpyum(f,w); gcdsfum(w,u,g);                          if ( q % 2 ) {
                                   spwrsfum(f,t,n5,s); subsfum(s,o,u);
                           } else
                                   tracemodsfum(f,t,ed,u);
                           cpyum(f,w);
                           gcdsfum(w,u,g);
                         if ( (DEG(g) >= 1) && (DEG(g) < DEG(f)) ) {                          if ( (DEG(g) >= 1) && (DEG(g) < DEG(f)) ) {
                                 canzassf(g,d,r);                                  canzassf(g,d,r);
                                 cpyum(f,w); divsfum(w,g,s);                                  cpyum(f,w); divsfum(w,g,s);
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>