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

Diff for /OpenXM_contrib2/asir2000/builtin/gf.c between version 1.9 and 1.14

version 1.9, 2001/06/28 08:57:20 version 1.14, 2001/10/09 01:36:06
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/gf.c,v 1.8 2001/06/25 10:01:27 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/gf.c,v 1.13 2001/09/03 07:56:19 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
Line 67  void nullspace_lm(LM **,int,int *);
Line 67  void nullspace_lm(LM **,int,int *);
 void nullspace_gf2n(GF2N **,int,int *);  void nullspace_gf2n(GF2N **,int,int *);
 void nullspace_gfpn(GFPN **,int,int *);  void nullspace_gfpn(GFPN **,int,int *);
 void nullspace_gfs(GFS **,int,int *);  void nullspace_gfs(GFS **,int,int *);
   void nullspace_gfsn(GFSN **,int,int *);
 void null_to_sol(int **,int *,int,int,UM *);  void null_to_sol(int **,int *,int,int,UM *);
   
 void showgfmat(UM **,int);  void showgfmat(UM **,int);
Line 209  NODE *rp;
Line 210  NODE *rp;
         BM fl;          BM fl;
         BM *r;          BM *r;
         VL vl,nvl;          VL vl,nvl;
         int i,fn,dx,dy;          int i,fn,dx,dy,d;
         NODE t,top;          NODE t,top;
         UM fm,hm,q;          UM fm,hm,q;
         UM *gm;          UM *gm;
Line 229  NODE *rp;
Line 230  NODE *rp;
         x = VR((P)BDY(mfl));          x = VR((P)BDY(mfl));
         y = vl->v == x ? vl->next->v : vl->v;          y = vl->v == x ? vl->next->v : vl->v;
   
         for ( i = 0, t = mfl; i < fn; i++, t = NEXT(t) ) {          for ( i = 0, t = mfl, d = 0; i < fn; i++, t = NEXT(t) ) {
                 gm[i] = (pointer)UMALLOC(getdeg(x,(P)BDY(t)));                  gm[i] = (pointer)UMALLOC(getdeg(x,(P)BDY(t)));
                 ptosfum((P)BDY(t),gm[i]);                  ptosfum((P)BDY(t),gm[i]);
                   d += DEG(gm[i]);
         }          }
   
         /* reorder f if necessary */          /* reorder f if necessary */
Line 240  NODE *rp;
Line 242  NODE *rp;
                 vl = nvl; f = g;                  vl = nvl; f = g;
         }          }
         dx = getdeg(x,f);          dx = getdeg(x,f);
           if ( dx != d )
                   error("sfuhensel : product of factors has incompatible degree");
   
         dy = getdeg(y,f);          dy = getdeg(y,f);
         if ( bound < dy+1 ) bound = dy+1;          dy = MAX(dy,bound);
         fl = BMALLOC(dx,bound);          fl = BMALLOC(dx,dy);
         ptosfbm(bound,f,fl);          ptosfbm(dy,f,fl);
         shiftsfbm(bound,fl,FTOIF(CONT(ev)));          if ( ev ) shiftsfbm(fl,FTOIF(CONT(ev)));
   
         /* fm = fl mod y */          /* fm = fl mod y */
         fm = W_UMALLOC(dx);          fm = W_UMALLOC(dx);
Line 265  NODE *rp;
Line 270  NODE *rp;
         r[i] = fl;          r[i] = fl;
   
         for ( i = fn-1, top = 0; i >= 0; i-- ) {          for ( i = fn-1, top = 0; i >= 0; i-- ) {
                 sfbmtop(bound,r[i],x,y,&s);                  sfbmtop(r[i],x,y,&s);
                 reorderp(CO,vl,s,&u);                  reorderp(CO,vl,s,&u);
                 MKNODE(t,u,top); top = t;                  MKNODE(t,u,top); top = t;
         }          }
Line 879  NODE arg;
Line 884  NODE arg;
 LIST *rp;  LIST *rp;
 {  {
         int i,j,n;          int i,j,n;
         Q mod;  
         MAT mat,r;          MAT mat,r;
         VECT u;          VECT u;
         Q q;          Q q;
Line 905  LIST *rp;
Line 909  LIST *rp;
                         nullspace_gfpn((GFPN **)w,n,ind); break;                          nullspace_gfpn((GFPN **)w,n,ind); break;
                 case FF_GFS:                  case FF_GFS:
                         nullspace_gfs((GFS **)w,n,ind); break;                          nullspace_gfs((GFS **)w,n,ind); break;
                   case FF_GFSN:
                           nullspace_gfsn((GFSN **)w,n,ind); break;
                 default:                  default:
                         error("nullspace_ff : current_ff is not set");                          error("nullspace_ff : current_ff is not set");
         }          }
Line 1081  int *ind;
Line 1087  int *ind;
                         chsgngfs(u[j],&h);                          chsgngfs(u[j],&h);
                         for ( s = j; s < n; s++ ) {                          for ( s = j; s < n; s++ ) {
                                 mulgfs(h,t[s],&w); addgfs(w,u[s],&w1); u[s] = w1;                                  mulgfs(h,t[s],&w); addgfs(w,u[s],&w1); u[s] = w1;
                           }
                   }
           }
   }
   
   void nullspace_gfsn(mat,n,ind)
   GFSN **mat;
   int n;
   int *ind;
   {
           int i,j,l,s;
           GFSN w,w1,h,inv;
           GFSN *t,*u;
   
           bzero(ind,n*sizeof(int));
           ind[0] = 0;
   
           for ( i = j = 0; j < n; i++, j++ ) {
                   for ( ; j < n; j++ ) {
                           for ( l = i; l < n; l++ )
                                   if ( mat[l][j] )
                                           break;
                           if ( l < n ) {
                                   t = mat[i]; mat[i] = mat[l]; mat[l] = t; break;
                           } else
                                   ind[j] = 1;
                   }
                   if ( j == n )
                           break;
                   invgfsn(mat[i][j],&inv);
                   for ( s = j, t = mat[i]; s < n; s++ ) {
                           mulgfsn(t[s],inv,&w); t[s] = w;
                   }
                   for ( l = 0; l < n; l++ ) {
                           if ( l == i )
                                   continue;
                           u = mat[l];
                           chsgngfsn(u[j],&h);
                           for ( s = j; s < n; s++ ) {
                                   mulgfsn(h,t[s],&w); addgfsn(w,u[s],&w1); u[s] = w1;
                         }                          }
                 }                  }
         }          }
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.14
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>