version 1.29, 2000/08/22 05:34:06 |
version 1.35, 2007/07/03 22:05:46 |
|
|
/* $OpenXM: OpenXM/src/k097/lib/minimal/minimal.k,v 1.28 2000/08/21 07:45:22 takayama Exp $ */ |
/* $OpenXM: OpenXM/src/k097/lib/minimal/minimal.k,v 1.34 2001/01/05 11:14:28 takayama Exp $ */ |
#define DEBUG 1 |
#define DEBUG 1 |
Sordinary = false; |
Sordinary = false; |
/* If you run this program on openxm version 1.1.2 (FreeBSD), |
/* If you run this program on openxm version 1.1.2 (FreeBSD), |
|
|
|
|
*/ |
*/ |
|
|
|
/* We cannot use load command in the if statement. */ |
|
load("lib/minimal/cohom.k"); |
|
Load_sm1(["k0-tower.sm1","lib/minimal/k0-tower.sm1"],"k0-tower.sm1.loaded"); |
|
Load_sm1(["new.sm1","lib/minimal/new.sm1"],"new.sm1.loaded"); |
|
sm1(" oxNoX "); |
|
|
load("cohom.k"); |
|
def load_tower() { |
|
if (Boundp("k0-tower.sm1.loaded")) { |
|
}else{ |
|
sm1(" [(parse) (k0-tower.sm1) pushfile ] extension "); |
|
sm1(" [(parse) (new.sm1) pushfile ] extension "); |
|
sm1(" /k0-tower.sm1.loaded 1 def "); |
|
} |
|
sm1(" oxNoX "); |
|
} |
|
load_tower(); |
|
SonAutoReduce = true; |
SonAutoReduce = true; |
def Factor(f) { |
def Factor(f) { |
sm1(f, " fctr /FunctionValue set"); |
sm1(f, " fctr /FunctionValue set"); |
|
|
HelpAdd(["Max", |
HelpAdd(["Max", |
["Max(v) returns the maximal element in v."]]); |
["Max(v) returns the maximal element in v."]]); |
|
|
|
def Kernel(f,v) { |
|
local ans; |
|
/* v : string or ring */ |
|
if (Length(Arglist) < 2) { |
|
sm1(" [f] syz /ans set "); |
|
}else{ |
|
sm1(" [f v] syz /ans set "); |
|
} |
|
return(ans); |
|
} |
|
def Syz(f) { |
|
sm1(" [f] syz /FunctionValue set "); |
|
} |
|
HelpAdd(["Kernel", |
|
["Kernel(f) returns the syzygy of f.", |
|
"Return value [b, c]: b is a set of generators of the syzygies of f", |
|
" : c=[gb, backward transformation, syzygy without", |
|
" dehomogenization", |
|
"Example: Weyl(\"x,y\",[[\"x\",-1,\"Dx\",1]]); ", |
|
" s=Kernel([x*Dx+1,Dx^2+x^5]); s[0]:"]]); |
|
/* cf. sm1_syz in cohom.k */ |
|
def Gb(f) { |
|
sm1(" [f] gb /FunctionValue set "); |
|
} |
|
HelpAdd(["Gb", |
|
["Gb(f) returns the Groebner basis of f.", |
|
"cf. Kernel, Weyl."]]); |
|
|
|
|
/* End of standard functions that should be moved to standard libraries. */ |
/* End of standard functions that should be moved to standard libraries. */ |
def test0() { |
def test0() { |
local f; |
local f; |
|
|
} |
} |
|
|
|
|
|
|
def Sweyl(v,w) { |
def Sweyl(v,w) { |
/* extern WeightOfSweyl ; */ |
/* extern WeightOfSweyl ; */ |
local ww,i,n; |
local ww,i,n; |
Line 224 def StoTower() { |
|
Line 246 def StoTower() { |
|
} |
} |
|
|
def SsetTower(tower) { |
def SsetTower(tower) { |
sm1(" [(AvoidTheSameRing)] pushEnv |
sm1(" [(AvoidTheSameRing)] pushEnv \ |
[ [(AvoidTheSameRing) 0] system_variable |
[ [(AvoidTheSameRing) 0] system_variable \ |
[(gbListTower) tower (list) dc] system_variable |
[(gbListTower) tower (list) dc] system_variable \ |
] pop popEnv "); |
] pop popEnv "); |
/* sm1("(hoge) message show_ring "); */ |
/* sm1("(hoge) message show_ring "); */ |
} |
} |
Line 336 def NewPolynomialVector(size) { |
|
Line 358 def NewPolynomialVector(size) { |
|
} |
} |
|
|
def SturnOffHomogenization() { |
def SturnOffHomogenization() { |
sm1(" |
sm1(" \ |
[(Homogenize)] system_variable 1 eq |
[(Homogenize)] system_variable 1 eq \ |
{ Sverbose { |
{ Sverbose { \ |
(Warning: Homogenization and ReduceLowerTerms options are automatically turned off.) message } { } ifelse |
(Warning: Homogenization and ReduceLowerTerms options are automatically turned off.) message } { } ifelse \ |
[(Homogenize) 0] system_variable |
[(Homogenize) 0] system_variable \ |
[(ReduceLowerTerms) 0] system_variable |
[(ReduceLowerTerms) 0] system_variable \ |
} { } ifelse |
} { } ifelse \ |
"); |
"); |
} |
} |
/* NOTE!!! Be careful these changes of global environmental variables. |
/* NOTE!!! Be careful these changes of global environmental variables. |
Line 350 def SturnOffHomogenization() { |
|
Line 372 def SturnOffHomogenization() { |
|
after computation and interruption. August 15, 2000. |
after computation and interruption. August 15, 2000. |
*/ |
*/ |
def SturnOnHomogenization() { |
def SturnOnHomogenization() { |
sm1(" |
sm1(" \ |
[(Homogenize)] system_variable 0 eq |
[(Homogenize)] system_variable 0 eq \ |
{ Sverbose { |
{ Sverbose { \ |
(Warning: Homogenization and ReduceLowerTerms options are automatically turned ON.) message } { } ifelse |
(Warning: Homogenization and ReduceLowerTerms options are automatically turned ON.) message } { } ifelse \ |
[(Homogenize) 1] system_variable |
[(Homogenize) 1] system_variable \ |
[(ReduceLowerTerms) 1] system_variable |
[(ReduceLowerTerms) 1] system_variable \ |
} { } ifelse |
} { } ifelse \ |
"); |
"); |
} |
} |
|
|
Line 1454 HelpAdd(["IsSameIdeal_h", |
|
Line 1476 HelpAdd(["IsSameIdeal_h", |
|
"cf. ReParse" |
"cf. ReParse" |
]]); |
]]); |
|
|
def ReParse(a) { |
/* |
local c; |
Output of S* functions may cause a trouble because it uses Schreyer orders. |
if (IsArray(a)) { |
In this case, use ReParse(). |
c = Map(a,"ReParse"); |
*/ |
}else{ |
|
sm1(a," toString . /c set"); |
|
} |
|
return(c); |
|
} |
|
HelpAdd(["ReParse", |
|
["Reparse(obj): obj", |
|
"It parses the given object in the current ring.", |
|
"Outputs from SlaScala, Sschreyer may cause a trouble in other functions,", |
|
"because it uses the Schreyer order.", |
|
"In this case, ReParse the outputs from these functions.", |
|
"cf. IsExaxt_h" |
|
]]); |
|
|
|
def ScheckIfSchreyer(s) { |
def ScheckIfSchreyer(s) { |
local ss; |
local ss; |