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

Diff for /OpenXM_contrib2/asir2000/engine/Fgfs.c between version 1.8 and 1.9

version 1.8, 2002/11/01 05:43:35 version 1.9, 2002/11/01 06:47:41
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/engine/Fgfs.c,v 1.7 2002/10/30 08:07:11 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/Fgfs.c,v 1.8 2002/11/01 05:43:35 noro Exp $ */
   
 #include "ca.h"  #include "ca.h"
   
Line 843  void mfctrsf_hensel(VL vl,VL rvl,P f,P pp0,P u0,P v0,P
Line 843  void mfctrsf_hensel(VL vl,VL rvl,P f,P pp0,P u0,P v0,P
         cv = (P *)ALLOCA((n+1)*sizeof(P));          cv = (P *)ALLOCA((n+1)*sizeof(P));
   
         /* ydy = y^dy */          /* ydy = y^dy */
         ydy = C_UMALLOC(dy); COEF(ydy)[dy] = 1;          ydy = C_UMALLOC(dy); DEG(ydy) = dy; COEF(ydy)[dy] = 1;
         setmod_gfsn(ydy);          setmod_gfsn(ydy);
   
         /* (R[y]/(y^dy))[x,X] */          /* (R[y]/(y^dy))[x,X] */
Line 980  void gfsn_univariate_to_sfbm(P f,int dy,BM *r)
Line 980  void gfsn_univariate_to_sfbm(P f,int dy,BM *r)
                 for ( i = DEG(cy); i >= 0; i-- )                  for ( i = DEG(cy); i >= 0; i-- )
                         COEF(COEF(b)[i])[d] = COEF(cy)[i];                          COEF(COEF(b)[i])[d] = COEF(cy)[i];
         }          }
           for ( i = 0; i <= dy; i++ )
                   degum(COEF(b)[i],dx);
         *r = b;          *r = b;
 }  }
   
Line 1021  void poly_to_gfsn_poly_main(P f,V v,P *r)
Line 1023  void poly_to_gfsn_poly_main(P f,V v,P *r)
         GFSN g;          GFSN g;
         DCP dc,dct,dc0;          DCP dc,dct,dc0;
   
         if ( !f || NUM(f) )          if ( !f )
                 *r = f;                  *r = f;
         else if ( VR(f) == v ) {          else if ( NUM(f) || VR(f) == v ) {
                 d = getdeg(v,f);                  d = getdeg(v,f);
                 u = UMALLOC(d);                  u = UMALLOC(d);
                 ptosfum(f,u);                  ptosfum(f,u);
Line 1080  void gfsn_poly_to_poly_main(P f,V v,P *r)
Line 1082  void gfsn_poly_to_poly_main(P f,V v,P *r)
                 MKP(VR(f),dc0,*r);                  MKP(VR(f),dc0,*r);
         }          }
 }  }
   
   void printsfum(UM f)
   {
           int i;
   
           for ( i = DEG(f); i >= 0; i-- ) {
                   printf("+(");
                   printf("%d",IFTOF(COEF(f)[i]));
                   printf(")*y^%d",i);
           }
   }
   
   void printsfbm(BM f)
   {
           int i;
   
           for ( i = DEG(f); i >= 0; i-- ) {
                   printf("+(");
                   printsfum(COEF(f)[i]);
                   printf(")*y^%d",i);
           }
   }
   
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>