OpenXM/src/asir-contrib/packages/sample/contiguity - diff

Return to contiguity CVS log

Up to [local] / OpenXM / src / asir-contrib / packages / sample

version 1.1, 2000/01/21 10:39:44

version 1.3, 2000/01/24 10:44:34

Line 1

/* $OpenXM$ */

/* $OpenXM: OpenXM/src/asir-contrib/packages/sample/contiguity,v 1.2 2000/01/22 11:20:00 takayama Exp $ */

/* contiguity */

Line 23 def contiguity_eliminated(G) {

def contiguity_f4_a_bf() {

V = [x,y,a,b,c,cp];

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];

B = contiguity_eliminated(G);

print("b-function for HG is ",0);

Line 37 def contiguity_f4_a_bf() {

def contiguity_f4_a() {

V = [x,y,a,b,c,cp];

W = [[dx,1,dy,1,x,1,y,1]];

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]);

Gb = map(first,G[1][0]);

Line 48 def contiguity_f4_a() {

for (I=0; I<length(Gb); I++) {

if (contiguity_param_only(Gb[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);