[BACK]Return to module_syz.rr CVS log [TXT][DIR] Up to [local] / OpenXM / src / asir-contrib / testing / noro

Diff for /OpenXM/src/asir-contrib/testing/noro/module_syz.rr between version 1.3 and 1.4

version 1.3, 2020/01/22 08:30:40 version 1.4, 2020/01/22 22:44:21
Line 1 
Line 1 
 module newsyz;  module newsyz;
   
 localf module_syz;  localf module_syz, module_syz_old;
 localf simplify_syz, icont, mod, remove_cont,ordcheck;  localf simplify_syz, icont, mod, remove_cont,ordcheck;
 localf complsb, complsb_sd, sortlsb, find_pos, find_pos, reduce, lres_setup, dpm_sort1, comp_pos;  localf complsb, complsb_sd, sortlsb, find_pos, find_pos, reduce, lres_setup, dpm_sort1, comp_pos;
 localf fres,minres,sres,minsres,lres, create_base_ord, simplify_k, simplify_by_k, remove_k, remove_k1, extract_nonzero;  localf fres,minres,sres,minsres,lres, create_base_ord, simplify_k, simplify_by_k, remove_k, remove_k1, extract_nonzero;
Line 64  def module_syz(F,V,H,Ord)
Line 64  def module_syz(F,V,H,Ord)
     B = cons(F[I]+dpm_dptodpm(dp_ptod(1,V),N+I+1),B);      B = cons(F[I]+dpm_dptodpm(dp_ptod(1,V),N+I+1),B);
   }    }
   B = reverse(B);    B = reverse(B);
   if ( Weyl )    if ( H >= 2 ) {
     G = nd_weyl_gr(B,V,0,[1,Ord]|dp=1,homo=H);      // finite field
   else if ( F4 ) {      if ( Weyl )
 //    G = nd_f4(B,V,0,[1,Ord]|dp=1,homo=H,nora=1);        G = nd_weyl_gr(B,V,H,[1,Ord]|dp=1);
     Ind = 0;      else if ( F4 )
     while ( 1 ) {        G = nd_f4(B,V,H,[1,Ord]|dp=1);
       G = nd_f4_trace(B,V,H,-lprime(Ind),[1,Ord]|dp=1);      else
       if ( G ) break;        G = nd_gr(B,V,H,[1,Ord]|dp=1);
       else Ind++;    } else {
     }      if ( Weyl )
   } else        G = nd_weyl_gr(B,V,0,[1,Ord]|dp=1,homo=H);
     G = nd_gr(B,V,0,[1,Ord]|dp=1,homo=H);      else if ( F4 ) {
         Ind = 0;
         while ( 1 ) {
           G = nd_f4_trace(B,V,H,-lprime(Ind),[1,Ord]|dp=1);
           if ( G ) break;
           else Ind++;
         }
       } else
         G = nd_gr(B,V,0,[1,Ord]|dp=1,homo=H);
     }
   G0 = []; S0 = []; Gen0 = [];    G0 = []; S0 = []; Gen0 = [];
   for ( T = G; T != []; T = cdr(T) ) {    for ( T = G; T != []; T = cdr(T) ) {
     H = car(T);      H = car(T);
Line 94  def module_syz(F,V,H,Ord)
Line 103  def module_syz(F,V,H,Ord)
     S0 = map(dpm_dtol,S0,V); G0 = map(dpm_dtol,G0,V); Gen0 = map(dpm_dtol,Gen0,V);      S0 = map(dpm_dtol,S0,V); G0 = map(dpm_dtol,G0,V); Gen0 = map(dpm_dtol,Gen0,V);
   }    }
   return [S0,G0,Gen0];    return [S0,G0,Gen0];
   }
   
   def module_syz_old(F,V,H,O)
   {
           Weyl = type(getopt(weyl)) != -1 ? 1 : 0;
           K = length(F);
           if ( type(F[0]) <= 2 ) {
                   for ( T = [], S = F; S != []; S = cdr(S) )
                           T = cons([car(S)],T);
                   F = reverse(T);
           }
           N = length(F[0]);
           B = [];
           for ( I = 0; I < K; I++ ) {
                   E = vector(N+K);
                   for ( J = 0; J < N; J++ ) E[J] = F[I][J];
                   E[N+I] = 1;
                   B = cons(vtol(E),B);
           }
           B = reverse(B);
           if ( H >= 2 ) {
                   if ( Weyl )
                           G = nd_weyl_gr(B,V,H,[1,O]);
                   else
                           G = nd_gr(B,V,H,[1,O]);
           } else {
                   if ( Weyl )
                           G = nd_weyl_gr_trace(B,V,H,-1,[1,O]);
                   else
                           G = nd_gr_trace(B,V,H,-1,[1,O]);
           }
           G0 = []; S0 = []; Gen0 = [];
           for ( T = G; T != []; T = cdr(T) ) {
                   H = car(T);
                   for ( J = 0; J < N; J++ ) if ( H[J] ) break;
                   if ( J == N ) {
                           H1 = vector(K);
                           for ( J = 0; J < K; J++ ) H1[J] = H[N+J];
                           S0 = cons(vtol(H1),S0);
                   } else {
                           H1 = vector(N);
                           for ( J = 0; J < N; J++ ) H1[J] = H[J];
                           G0 = cons(vtol(H1),G0);
                           H1 = vector(K);
                           for ( J = 0; J < K; J++ ) H1[J] = H[N+J];
                           Gen0 = cons(vtol(H1),Gen0);
                   }
           }
           return [S0,G0,Gen0];
 }  }
   
 def fres(F,V,H,O)  def fres(F,V,H,O)
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>