[BACK]Return to gr CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / lib

Diff for /OpenXM_contrib2/asir2000/lib/gr between version 1.13 and 1.15

version 1.13, 2001/11/19 00:57:13 version 1.15, 2002/06/12 08:19:04
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/lib/gr,v 1.12 2001/11/01 10:00:19 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/lib/gr,v 1.14 2001/11/19 01:40:05 noro Exp $
 */  */
 extern INIT_COUNT,ITOR_FAIL$  extern INIT_COUNT,ITOR_FAIL$
 extern REMOTE_MATRIX,REMOTE_NF,REMOTE_VARS$  extern REMOTE_MATRIX,REMOTE_NF,REMOTE_VARS$
Line 352  def minipoly(G0,V,O,P,V0)
Line 352  def minipoly(G0,V,O,P,V0)
         if ( !zero_dim(hmlist(G0,V,O),V,O) )          if ( !zero_dim(hmlist(G0,V,O),V,O) )
                 error("tolex : ideal is not zero-dimensional!");                  error("tolex : ideal is not zero-dimensional!");
   
           Pin = P;
           P = ptozp(P);
           CP = sdiv(P,Pin);
         G1 = cons(V0-P,G0);          G1 = cons(V0-P,G0);
         O1 = [[0,1],[O,length(V)]];          O1 = [[0,1],[O,length(V)]];
         V1 = cons(V0,V);          V1 = cons(V0,V);
Line 372  def minipoly(G0,V,O,P,V0)
Line 375  def minipoly(G0,V,O,P,V0)
                         TL = cons(V0^J,TL);                          TL = cons(V0^J,TL);
                 NF = gennf(G1,TL,V1,O1,V0,1)[0];                  NF = gennf(G1,TL,V1,O1,V0,1)[0];
                 R = tolex_main(V1,O1,NF,[MP],M,MB);                  R = tolex_main(V1,O1,NF,[MP],M,MB);
                 return R[0];                  return ptozp(subst(R[0],V0,CP*V0));
         }          }
 }  }
   
Line 1575  def gbcheck(B,V,O)
Line 1578  def gbcheck(B,V,O)
 {  {
         dp_ord(O);          dp_ord(O);
         D = map(dp_ptod,B,V);          D = map(dp_ptod,B,V);
         L = dp_gr_checklist(D);          L = dp_gr_checklist(D,length(V));
         DP = L[0]; Plist = L[1];          DP = L[0]; Plist = L[1];
         for ( IL = [], I = size(DP)[0]-1; I >= 0; I-- )          for ( IL = [], I = size(DP)[0]-1; I >= 0; I-- )
                 IL = cons(I,IL);                  IL = cons(I,IL);

Legend:
Removed from v.1.13  
changed lines
  Added in v.1.15

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>