[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.1 and 1.6

version 1.1, 1999/12/03 07:39:07 version 1.6, 2001/06/25 01:35:20
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/asir99/builtin/gf.c,v 1.1.1.1 1999/11/10 08:12:25 noro Exp $ */  /*
    * Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED
    * All rights reserved.
    *
    * FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited,
    * non-exclusive and royalty-free license to use, copy, modify and
    * redistribute, solely for non-commercial and non-profit purposes, the
    * computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and
    * conditions of this Agreement. For the avoidance of doubt, you acquire
    * only a limited right to use the SOFTWARE hereunder, and FLL or any
    * third party developer retains all rights, including but not limited to
    * copyrights, in and to the SOFTWARE.
    *
    * (1) FLL does not grant you a license in any way for commercial
    * purposes. You may use the SOFTWARE only for non-commercial and
    * non-profit purposes only, such as academic, research and internal
    * business use.
    * (2) The SOFTWARE is protected by the Copyright Law of Japan and
    * international copyright treaties. If you make copies of the SOFTWARE,
    * with or without modification, as permitted hereunder, you shall affix
    * to all such copies of the SOFTWARE the above copyright notice.
    * (3) An explicit reference to this SOFTWARE and its copyright owner
    * shall be made on your publication or presentation in any form of the
    * results obtained by use of the SOFTWARE.
    * (4) In the event that you modify the SOFTWARE, you shall notify FLL by
    * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification
    * for such modification or the source code of the modified part of the
    * SOFTWARE.
    *
    * THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL
    * MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND
    * EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS
    * FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES'
    * RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY
    * MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY.
    * UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT,
    * OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY
    * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL
    * DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES
    * ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES
    * FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY
    * DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF
    * SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART
    * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY
    * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
    * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
    *
    * $OpenXM: OpenXM_contrib2/asir2000/builtin/gf.c,v 1.5 2001/05/09 01:41:41 noro Exp $
   */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
   
Line 18  void nullspace(UM **,UM,int,int,int *);
Line 66  void nullspace(UM **,UM,int,int,int *);
 void nullspace_lm(LM **,int,int *);  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 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 25  void pwr_mod(P,P,V,P,int,N,P *);
Line 74  void pwr_mod(P,P,V,P,int,N,P *);
 void rem_mod(P,P,V,P,int,P *);  void rem_mod(P,P,V,P,int,P *);
   
 void Pnullspace(),Pgcda_mod(),Pftest(),Presfmain(),Ppwr_mod(),Puhensel();  void Pnullspace(),Pgcda_mod(),Pftest(),Presfmain(),Ppwr_mod(),Puhensel();
   void Psfuhensel();
   
 void Pnullspace_ff();  void Pnullspace_ff();
   
Line 34  void Plinear_form_to_vect(),Pvect_to_linear_form();
Line 84  void Plinear_form_to_vect(),Pvect_to_linear_form();
 void solve_linear_equation_gf2n(GF2N **,int,int,int *);  void solve_linear_equation_gf2n(GF2N **,int,int,int *);
 void linear_form_to_array(P,VL,int,Num *);  void linear_form_to_array(P,VL,int,Num *);
 void array_to_linear_form(Num *,VL,int,P *);  void array_to_linear_form(Num *,VL,int,P *);
   void sfuhensel(P,V,V,NODE *);
   
 extern int current_ff;  extern int current_ff;
   
Line 46  struct ftab gf_tab[] = {
Line 97  struct ftab gf_tab[] = {
         {"resfmain",Presfmain,4},          {"resfmain",Presfmain,4},
         {"pwr_mod",Ppwr_mod,6},          {"pwr_mod",Ppwr_mod,6},
         {"uhensel",Puhensel,4},          {"uhensel",Puhensel,4},
           {"sfuhensel",Psfuhensel,3},
         {0,0,0},          {0,0,0},
 };  };
   
Line 130  NODE *rp;
Line 182  NODE *rp;
         *rp = top;          *rp = top;
 }  }
   
   void Psfuhensel(arg,rp)
   NODE arg;
   LIST *rp;
   {
           P f;
           V x,y;
           NODE r;
   
           f = (P)ARG0(arg);
           x = VR((P)ARG1(arg));
           y = VR((P)ARG2(arg));
           sfuhensel(f,x,y,&r);
           MKLIST(*rp,r);
   }
   
   void sfuhensel(f,x,y,rp)
   P f;
   V x,y;
   NODE *rp;
   {
           ML lift;
           int i;
           P s,u;
           NODE t,top;
           VL vl,vl1;
   
           NEWVL(vl1); vl1->v = y; NEXT(vl1) = 0;
           NEWVL(vl); vl->v = x; NEXT(vl) = vl1;
           sfhensel(5,f,x,&lift);
           for ( i = lift->n-1, top = 0; i >= 0; i-- ) {
                   sflumtop(lift->bound,lift->c[i],x,y,&s);
                   reorderp(CO,vl,s,&u);
                   MKNODE(t,u,top); top = t;
           }
           *rp = top;
   }
   
 void Presfmain(arg,rp)  void Presfmain(arg,rp)
 NODE arg;  NODE arg;
 LIST *rp;  LIST *rp;
Line 180  P *mfl;
Line 269  P *mfl;
 ML *listp;  ML *listp;
 {  {
         register int i,j;          register int i,j;
         int q,n,bound;          int q,n,bound,inv,lc;
         int *p;          int *p;
         int **pp;          int **pp;
         ML blist,clist,bqlist,cqlist,rlist;          ML blist,clist,bqlist,cqlist,rlist;
         UM *b;          UM *b;
         LUM fl,tl;          LUM fl,tl;
         LUM *l;          LUM *l;
           UM w;
   
           w = W_UMALLOC(UDEG(f));
         blist = MLALLOC(nf); blist->n = nf; blist->mod = mod;          blist = MLALLOC(nf); blist->n = nf; blist->mod = mod;
         for ( i = 0; i < nf; i++ ) {  
           /* c[0] must have lc(f) */
           blist->c[0] = (pointer)UMALLOC(UDEG(mfl[0]));
           ptoum(mod,mfl[0],w);
           inv = invm(w->c[UDEG(mfl[0])],mod);
           lc = rem(NM((Q)LC(f)),mod);
           if ( SGN((Q)LC(f)) < 0 )
                   lc = (mod-lc)%mod;
           lc = dmar(inv,lc,0,mod);
           if ( lc == 1 )
                   copyum(w,blist->c[0]);
           else
                   mulsum(mod,w,lc,blist->c[0]);
   
           /* c[i] (i=1,...,nf-1) must be monic */
           for ( i = 1; i < nf; i++ ) {
                 blist->c[i] = (pointer)UMALLOC(UDEG(mfl[i]));                  blist->c[i] = (pointer)UMALLOC(UDEG(mfl[i]));
                 ptoum(mod,mfl[i],blist->c[i]);                  ptoum(mod,mfl[i],w);
                   inv = invm(w->c[UDEG(mfl[i])],mod);
                   if ( inv == 1 )
                           copyum(w,blist->c[i]);
                   else
                           mulsum(mod,w,inv,blist->c[i]);
         }          }
   
         gcdgen(f,blist,&clist); henprep(f,blist,clist,&bqlist,&cqlist);          gcdgen(f,blist,&clist); henprep(f,blist,clist,&bqlist,&cqlist);
         n = bqlist->n; q = bqlist->mod;          n = bqlist->n; q = bqlist->mod;
         bqlist->bound = cqlist->bound = bound = mignotte(q,f);          bqlist->bound = cqlist->bound = bound = mignotte(q,f);
Line 738  LIST *rp;
Line 850  LIST *rp;
                         nullspace_gf2n((GF2N **)w,n,ind); break;                          nullspace_gf2n((GF2N **)w,n,ind); break;
                 case FF_GFPN:                  case FF_GFPN:
                         nullspace_gfpn((GFPN **)w,n,ind); break;                          nullspace_gfpn((GFPN **)w,n,ind); break;
                   case FF_GFS:
                           nullspace_gfs((GFS **)w,n,ind); break;
                 default:                  default:
                         error("nullspace_ff : current_ff is not set");                          error("nullspace_ff : current_ff is not set");
         }          }
Line 876  int *ind;
Line 990  int *ind;
                 }                  }
         }          }
 }  }
   
   void nullspace_gfs(mat,n,ind)
   GFS **mat;
   int n;
   int *ind;
   {
           int i,j,l,s;
           GFS w,w1,h,inv;
           GFS *t,*u;
           GFS one;
   
           bzero(ind,n*sizeof(int));
           ind[0] = 0;
           mqtogfs(ONEM,&one);
   
           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;
                   divgfs(one,mat[i][j],&inv);
                   for ( s = j, t = mat[i]; s < n; s++ ) {
                           mulgfs(t[s],inv,&w); t[s] = w;
                   }
                   for ( l = 0; l < n; l++ ) {
                           if ( l == i )
                                   continue;
                           u = mat[l];
                           chsgngfs(u[j],&h);
                           for ( s = j; s < n; s++ ) {
                                   mulgfs(h,t[s],&w); addgfs(w,u[s],&w1); u[s] = w1;
                           }
                   }
           }
   }
   
 /* p = a(0)vl[0]+a(1)vl[1]+...+a(m-1)vl[m-1]+a(m) -> array = [a(0) a(1) ... a(m)] */  /* p = a(0)vl[0]+a(1)vl[1]+...+a(m-1)vl[m-1]+a(m) -> array = [a(0) a(1) ... a(m)] */
   
 void linear_form_to_array(p,vl,m,array)  void linear_form_to_array(p,vl,m,array)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>