[BACK]Return to contiguity CVS log [TXT][DIR] Up to [local] / OpenXM / src / asir-contrib / packages / sample

Diff for /OpenXM/src/asir-contrib/packages/sample/contiguity between version 1.1 and 1.3

version 1.1, 2000/01/21 10:39:44 version 1.3, 2000/01/24 10:44:34
Line 1 
Line 1 
 /* $OpenXM$ */  /* $OpenXM: OpenXM/src/asir-contrib/packages/sample/contiguity,v 1.2 2000/01/22 11:20:00 takayama Exp $ */
   
 /* contiguity */  /* contiguity */
   
Line 23  def contiguity_eliminated(G) {
Line 23  def contiguity_eliminated(G) {
 def contiguity_f4_a_bf() {  def contiguity_f4_a_bf() {
    V = [x,y,a,b,c,cp];     V = [x,y,a,b,c,cp];
    W = [[dx,1,dy,1,x,1,y,1]];     W = [[dx,1,dy,1,x,1,y,1]];
    G = sm1_gb([[dx*(x*dx+c)-(x*dx+y*dy+a)*(x*dx+y*dy+b),     F = [x*dx+y*dy+a,
                 dy*(y*dy+cp)-(x*dx+y*dy+a)*(x*dx+y*dy+b),          sm1_mul(dx,x*dx+c-1,V)-sm1_mul(x*dx+y*dy+a,x*dx+y*dy+b,V),
                 x*dx+y*dy+a],V,W]);          sm1_mul(dy,y*dy+cp-1,V)-sm1_mul(x*dx+y*dy+a,x*dx+y*dy+b,V)];
      G = sm1_gb([F,V,W]);
    G = G[0];     G = G[0];
    B = contiguity_eliminated(G);     B = contiguity_eliminated(G);
    print("b-function for HG is ",0);     print("b-function for HG is ",0);
Line 37  def contiguity_f4_a_bf() {
Line 37  def contiguity_f4_a_bf() {
 def contiguity_f4_a() {  def contiguity_f4_a() {
    V = [x,y,a,b,c,cp];     V = [x,y,a,b,c,cp];
    W = [[dx,1,dy,1,x,1,y,1]];     W = [[dx,1,dy,1,x,1,y,1]];
    F = [x*dx+y*dy+a,     F = [x*dx+y*dy+a,
         dx*(x*dx+c)-(x*dx+y*dy+a)*(x*dx+y*dy+b),          sm1_mul(dx,x*dx+c-1,V)-sm1_mul(x*dx+y*dy+a,x*dx+y*dy+b,V),
         dy*(y*dy+cp)-(x*dx+y*dy+a)*(x*dx+y*dy+b)];          sm1_mul(dy,y*dy+cp-1,V)-sm1_mul(x*dx+y*dy+a,x*dx+y*dy+b,V)];
    G = sm1_syz([F,V,W]);     G = sm1_syz([F,V,W]);
   
    Gb = map(first,G[1][0]);     Gb = map(first,G[1][0]);
Line 48  def contiguity_f4_a() {
Line 48  def contiguity_f4_a() {
    for (I=0; I<length(Gb); I++) {     for (I=0; I<length(Gb); I++) {
      if (contiguity_param_only(Gb[I]) != 0) {       if (contiguity_param_only(Gb[I]) != 0) {
         Ans = append(Ans,[ [Gb[I],C[I][0]] ]);          Ans = append(Ans,[ [Gb[I],C[I][0]] ]);
           /* print(sm1_mul(C[I],F,[x,y,a,b,c,cp])); */ /* check ! */
      }       }
    }     }
    return(Ans);     return(Ans);

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3

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