File: [local] / OpenXM / src / k097 / lib / minimal / minimal-test.k (download)
Revision 1.1, Wed May 24 15:31:28 2000 UTC (24 years, 1 month ago) by takayama
Branch: MAIN
Buenas Tardes!
resol.c gets an option for debugging to build a schreyer skelton.
NOTE:
The following method of ssh-forwarding works to commit to kerberos.
ssh -L 8000:kerberos.math.kobe-u.ac.jp:22 unkei.math.kobe-u.ac.jp
ssh-takayama (shell file) : ssh -v -p 8000 -l takayama $*
.bashrc : alias rcvs='cvs -d localhost:/usr/cvs'
|
/* $OpenXM: OpenXM/src/k097/lib/minimal/minimal-test.k,v 1.1 2000/05/24 15:31:28 takayama Exp $ */
load["minimal.k"];
def test5() {
local a,b,c,cc,v;
a = Sannfs3_laScala2("x^3-y^2*z^2");
b = a[0];
v = [x,y,z];
c = Skernel(b[0],v);
c = c[0];
sm1_pmat([c,b[1],v]);
Println("-----------------------------------");
cc = sm1_res_div(c,b[1],v);
sm1_pmat(sm1_gb(cc,v));
c = Skernel(b[1],v);
c = c[0];
cc = sm1_res_div(c,b[2],v);
sm1_pmat(sm1_gb(cc,v));
return(a);
}
def test6() {
local a,b,c,cc,v;
a = Sannfs3("x^3-y^2*z^2");
b = a[0];
v = [x,y,z];
c = Skernel(b[0],v);
c = c[0];
sm1_pmat([c,b[1],v]);
Println("-------ker = im for minimal ?---------------------");
cc = sm1_res_div(c,b[1],v);
sm1_pmat(sm1_gb(cc,v));
c = Skernel(b[1],v);
c = c[0];
cc = sm1_res_div(c,b[2],v);
sm1_pmat(sm1_gb(cc,v));
Println("------ ker=im for Schreyer ?------------------");
b = a[3];
c = Skernel(b[0],v);
c = c[0];
sm1_pmat([c,b[1],v]);
cc = sm1_res_div(c,b[1],v);
sm1_pmat(sm1_gb(cc,v));
c = Skernel(b[1],v);
c = c[0];
cc = sm1_res_div(c,b[2],v);
sm1_pmat(sm1_gb(cc,v));
return(a);
}
/* May 23, Tue */
def test7() {
local a,b,c,cc,v;
a = Sannfs3_laScala2("x^3-y^2*z^2");
b = a[0];
v = [x,y,z];
c = Skernel(b[0],v);
c = c[0];
sm1_pmat([c,b[1],v]);
Println("-------ker = im for minimal ?---------------------");
cc = sm1_res_div(c,b[1],v);
sm1_pmat(sm1_gb(cc,v));
c = Skernel(b[1],v);
c = c[0];
cc = sm1_res_div(c,b[2],v);
sm1_pmat(sm1_gb(cc,v));
Println("------ ker=im for Schreyer ?------------------");
b = a[3];
c = Skernel(b[0],v);
c = c[0];
sm1_pmat([c,b[1],v]);
cc = sm1_res_div(c,b[1],v);
sm1_pmat(sm1_gb(cc,v));
c = Skernel(b[1],v);
c = c[0];
cc = sm1_res_div(c,b[2],v);
sm1_pmat(sm1_gb(cc,v));
return(a);
}
def sm1_resol1(p) {
sm1(" p resol1 /FunctionValue set ");
}
def test8() {
local p,pp,ans,b,c,cc,ww,ww2;
f = "x^3-y^2*z^2";
p = Sannfs(f,"x,y,z");
ww = [["x",1,"y",1,"z",1,"Dx",1,"Dy",1,"Dz",1,"h",1],
["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]];
ww2 = [["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]];
sm1(" p 0 get { [(x) (y) (z) (Dx) (Dy) (Dz)] laplace0 } map /p set ");
Sweyl("x,y,z",ww);
pp = Map(p,"Spoly");
/* return(pp); */
/* pp =
[y*Dy-z*Dz , -2*x*Dx-3*y*Dy+1 , 2*x*Dy*Dz^2-3*y*Dx^2 ,
2*x*Dy^2*Dz-3*z*Dx^2 , 2*x*z*Dz^3-3*y^2*Dx^2+4*x*Dz^2 ]
*/
ans = sm1_resol1([pp,"x,y,z",ww]);
/* Schreyer is in ans. */
v = [x,y,z];
b = ans;
Println("------ ker=im for Schreyer ?------------------");
c = Skernel(b[0],v);
c = c[0];
sm1_pmat([c,b[1],v]);
cc = sm1_res_div(c,b[1],v);
sm1_pmat(sm1_gb(cc,v));
c = Skernel(b[1],v);
c = c[0];
cc = sm1_res_div(c,b[2],v);
sm1_pmat(sm1_gb(cc,v));
return(ans);
}
/*
a = test8();
SisComplex(a):
*/
def test8a() {
local p,pp,ans,b,c,cc,ww, ans_all;
f = "x^3-y^2*z^2";
p = Sannfs(f,"x,y,z");
sm1(" p 0 get { [(x) (y) (z) (Dx) (Dy) (Dz)] laplace0 } map /p set ");
ww = [["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]];
/* Removed "x",1, ... ===> It causes an error. I do not know the reason.*/
Sweyl("x,y,z",ww);
pp = Map(p,"Spoly");
/* return(pp); */
/* pp =
[y*Dy-z*Dz , -2*x*Dx-3*y*Dy+1 , 2*x*Dy*Dz^2-3*y*Dx^2 ,
2*x*Dy^2*Dz-3*z*Dx^2 , 2*x*z*Dz^3-3*y^2*Dx^2+4*x*Dz^2 ]
*/
ans_all = Sschreyer(pp);
ans = ans_all[0];
/* ans = sm1_resol1([pp,"x,y,z",ww]); */
/* Schreyer is in ans. */
v = [x,y,z];
b = ans;
Println("------ ker=im for Schreyer ?------------------");
c = Skernel(b[0],v);
c = c[0];
sm1_pmat([c,b[1],v]);
cc = sm1_res_div(c,b[1],v);
sm1_pmat(sm1_gb(cc,v));
c = Skernel(b[1],v);
c = c[0];
cc = sm1_res_div(c,b[2],v);
sm1_pmat(sm1_gb(cc,v));
return(ans);
}
/* Comparing two constructions */
def test9() {
local p,pp,ans,b,c,cc,ww,ww2,ans_all,ans2;
f = "x^3-y^2*z^2";
p = Sannfs(f,"x,y,z");
ww2 = [["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]];
sm1(" p 0 get { [(x) (y) (z) (Dx) (Dy) (Dz)] laplace0 } map /p set ");
Sweyl("x,y,z",ww2);
pp = Map(p,"Spoly");
ans = sm1_resol1([pp,"x,y,z",ww2]);
f = "x^3-y^2*z^2";
p = Sannfs(f,"x,y,z");
sm1(" p 0 get { [(x) (y) (z) (Dx) (Dy) (Dz)] laplace0 } map /p set ");
ww = [["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]];
Sweyl("x,y,z",ww);
pp = Map(p,"Spoly");
ans_all = Sschreyer(pp);
ans2 = ans_all[0];
return([ans,ans2]);
}