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

Diff for /OpenXM_contrib2/asir2000/builtin/algnum.c between version 1.9 and 1.10

version 1.9, 2004/12/06 09:29:34 version 1.10, 2005/01/23 14:03:47
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/algnum.c,v 1.8 2004/12/06 01:15:18 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/algnum.c,v 1.9 2004/12/06 09:29:34 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
Line 54  void Pdefpoly(), Pnewalg(), Pmainalg(), Palgtorat(), P
Line 54  void Pdefpoly(), Pnewalg(), Pmainalg(), Palgtorat(), P
 void Palg(), Palgv(), Pgetalgtree();  void Palg(), Palgv(), Pgetalgtree();
 void Pinvalg_le();  void Pinvalg_le();
 void Pset_field(),Palgtodalg(),Pdalgtoalg();  void Pset_field(),Palgtodalg(),Pdalgtoalg();
   void Pinv_or_split_dalg();
   
 void mkalg(P,Alg *);  void mkalg(P,Alg *);
 int cmpalgp(P,P);  int cmpalgp(P,P);
Line 68  struct ftab alg_tab[] = {
Line 69  struct ftab alg_tab[] = {
         {"set_field",Pset_field,1},          {"set_field",Pset_field,1},
         {"algtodalg",Palgtodalg,1},          {"algtodalg",Palgtodalg,1},
         {"dalgtoalg",Pdalgtoalg,1},          {"dalgtoalg",Pdalgtoalg,1},
           {"inv_or_split_dalg",Pinv_or_split_dalg,1},
         {"invalg_le",Pinvalg_le,1},          {"invalg_le",Pinvalg_le,1},
         {"defpoly",Pdefpoly,1},          {"defpoly",Pdefpoly,1},
         {"newalg",Pnewalg,1},          {"newalg",Pnewalg,1},
Line 97  void Palgtodalg(NODE arg,DAlg *rp)
Line 99  void Palgtodalg(NODE arg,DAlg *rp)
 void Pdalgtoalg(NODE arg,Alg *rp)  void Pdalgtoalg(NODE arg,Alg *rp)
 {  {
         dalgtoalg((DAlg)ARG0(arg),rp);          dalgtoalg((DAlg)ARG0(arg),rp);
   }
   
   NODE inv_or_split_dalg(DAlg,DAlg *);
   NumberField     get_numberfield();
   
   void Pinv_or_split_dalg(NODE arg,Obj *rp)
   {
           NODE gen,t,nd0,nd;
           LIST list;
           int l,i,j,k,n;
           DP *ps,*ps1,*psw;
           NumberField nf;
           DAlg inv;
           extern struct order_spec *dp_current_spec;
           struct order_spec *current_spec;
   
           gen = inv_or_split_dalg((DAlg)ARG0(arg),&inv);
           if ( !gen )
                   *rp = (Obj)inv;
           else {
                   nf = get_numberfield();
                   current_spec = dp_current_spec; initd(nf->spec);
                   l = length(gen);
                   n = nf->n;
                   ps = nf->ps;
                   psw = (DP *)ALLOCA((n+l)*sizeof(DP));
                   for ( i = j = 0; i < n; i++ ) {
                           for ( t = gen; t; t = NEXT(t) )
                                   if ( dp_redble(ps[i],(DP)BDY(t)) ) break;
                           if ( !t )
                                   psw[j++] = ps[i];
                   }
                   nd0  = 0;
                   /* gen[0] < gen[1] < ... */
                   /* psw[0] > psw[1] > ... */
                   for ( i = j-1, t = gen; i >= 0 && t; ) {
                           NEXTNODE(nd0,nd);
                           if ( compd(CO,psw[i],(DP)BDY(t)) > 0 ) {
                                   BDY(nd) = BDY(t); t = NEXT(t);
                           } else
                                   BDY(nd) = (pointer)psw[i--];
                   }
                   for ( ; i >= 0; i-- ) {
                           NEXTNODE(nd0,nd); BDY(nd) = (pointer)psw[i];
                   }
                   for ( ; t; t = NEXT(t), k++ ) {
                           NEXTNODE(nd0,nd); BDY(nd) = BDY(t);
                   }
                   NEXT(nd) = 0;
                   MKLIST(list,nd0);
                   initd(current_spec);
                   *rp = (Obj)list;
           }
 }  }
   
 void Pnewalg(arg,rp)  void Pnewalg(arg,rp)
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>