version 1.4, 2000/12/10 03:12:19 |
version 1.8, 2001/01/05 11:14:25 |
|
|
/* $OpenXM: OpenXM/src/k097/slib.k,v 1.3 2000/12/10 02:21:46 takayama Exp $ */ |
/* $OpenXM: OpenXM/src/k097/slib.k,v 1.7 2001/01/04 12:29:31 takayama Exp $ */ |
/* slib.sm1, standard library. */ |
/* slib.sm1, standard library. */ |
/* April 26-- , 1996 */ |
/* April 26-- , 1996 */ |
/* Don't use function names that is already used as a postscipt macro names*/ |
/* Don't use function names that is already used as a postscipt macro names*/ |
Line 37 def Append(f,g) { return(Join(f,[g])); } |
|
Line 37 def Append(f,g) { return(Join(f,[g])); } |
|
|
|
def Length(f) { sm1(f," length (universalNumber) dc /FunctionValue set"); } |
def Length(f) { sm1(f," length (universalNumber) dc /FunctionValue set"); } |
|
|
def Indexed(name,i) { |
|
sm1(name,i," s.Indexed /FunctionValue set "); |
|
} |
|
/* Indexed2("a",2,3) ---> "a[2,3]" */ |
|
def Indexed2(name,i,j) { |
|
sm1(name,i,j," s.Indexed2 /FunctionValue set "); |
|
} |
|
def Transpose(mat) { |
def Transpose(mat) { |
sm1(mat," transpose /FunctionValue set "); |
sm1(mat," transpose /FunctionValue set "); |
} |
} |
Line 78 def LiftStd(F) { /* Print("Input is "); Println(F); */ |
|
Line 71 def LiftStd(F) { /* Print("Input is "); Println(F); */ |
|
[arg1 [(needBack)]] groebner |
[arg1 [(needBack)]] groebner |
/FunctionValue set "); } |
/FunctionValue set "); } |
|
|
|
/* |
def Reduction(f,G) { |
def Reduction(f,G) { |
sm1(f,G," reduction /FunctionValue set "); |
sm1(f,G," reduction /FunctionValue set "); |
} |
} |
|
*/ |
|
def Reduction(f,myset) { |
|
local n, indexTable, set2, i, j, tmp, t_syz,r,rng, |
|
vsize,tt; |
|
vsize = null; |
|
r = GetRing(Poly("1")); /* Save the current ring */ |
|
rng = GetRing(f); |
|
if (Tag(rng) == 0) { |
|
rng = GetRing(myset); |
|
} |
|
if (Tag(rng) != 0) {SetRing(rng);} |
|
|
|
if (IsArray(f)) { |
|
vsize = Length(f); |
|
sm1(" [f] fromVectors 0 get /f set "); |
|
} |
|
|
|
n = Length(myset); |
|
if (n > 0) { |
|
if (IsArray(myset[0])) { |
|
if (vsize != Length(myset[0])) { |
|
Error("Reduction: size mismatch."); |
|
} |
|
sm1(" myset fromVectors /myset set "); |
|
} |
|
} |
|
|
|
indexTable = NewArray(n); |
|
set2 = [ ]; |
|
j = 0; |
|
for (i=0; i<n; i++) { |
|
if (Tag(myset[i]) == 0) { |
|
indexTable[i] = -1; |
|
}else if (myset[i] == Poly("0")) { |
|
indexTable[i] = -1; |
|
}else{ |
|
set2 = Append(set2,myset[i]); |
|
indexTable[i] = j; |
|
j++; |
|
} |
|
} |
|
sm1(" f set2 (gradedPolySet) dc reduction /tmp set "); |
|
t_syz = NewArray(n); |
|
for (i=0; i<n; i++) { |
|
if (indexTable[i] != -1) { |
|
t_syz[i] = tmp[2, indexTable[i]]; |
|
}else{ |
|
t_syz[i] = Poly("0"); |
|
} |
|
} |
|
if (Tag(vsize) != 0) { |
|
tt = tmp[0]; |
|
sm1(" [vsize (integer) dc tt] toVectors /tt set "); |
|
tmp[0] = tt; |
|
} |
|
SetRing(r); |
|
return([tmp[0],tmp[1],t_syz]); |
|
} |
|
|
|
|
|
|
def IntegerToSm1Integer(f) { |
def IntegerToSm1Integer(f) { |
sm1(f, " (integer) dc /FunctionValue set "); |
sm1(f, " (integer) dc /FunctionValue set "); |
} |
} |
Line 94 def RingD(vList,weightMatrix,pp) { |
|
Line 146 def RingD(vList,weightMatrix,pp) { |
|
sm1("[", vList, |
sm1("[", vList, |
"ring_of_differential_operators ( ) elimination_order 0 ] define_ring |
"ring_of_differential_operators ( ) elimination_order 0 ] define_ring |
/tmp set "); |
/tmp set "); |
|
SetRingVariables(); |
return(tmp); |
return(tmp); |
} else ; |
} else ; |
if (argsize == 2) { |
if (argsize == 2) { |
Line 117 def RingD(vList,weightMatrix,pp) { |
|
Line 170 def RingD(vList,weightMatrix,pp) { |
|
ringpp = |
ringpp = |
sm1("[", vList, |
sm1("[", vList, |
"ring_of_differential_operators ", new0, " weight_vector",pp, " ] define_ring"); |
"ring_of_differential_operators ", new0, " weight_vector",pp, " ] define_ring"); |
/* setRingVariables(); It doesn't work. It's a mystery. */ |
SetRingVariables(); |
sm1(" @@@.indexMode.flag.save @@@.indexMode "); |
sm1(" @@@.indexMode.flag.save @@@.indexMode "); |
return( ringpp ); |
return( ringpp ); |
} |
} |
|
|
sm1( "[(N)] system_variable (universalNumber) dc /FunctionValue set "); |
sm1( "[(N)] system_variable (universalNumber) dc /FunctionValue set "); |
} |
} |
|
|
SetRingVariables_Verbose = true; |
SetRingVariables_Verbose = false; |
def SetRingVariables() { |
def SetRingVariables() { |
/* Don't use local variables in this function, |
/* Don't use local variables in this function, |
because we set global variables in this function. |
because we set global variables in this function. |
Line 152 def SetRingVariables() { |
|
Line 205 def SetRingVariables() { |
|
if (SetRingVariables_Verbose ) { |
if (SetRingVariables_Verbose ) { |
Print("SetRingVariables() Setting the global variables : "); |
Print("SetRingVariables() Setting the global variables : "); |
} |
} |
k00setRingVariables(0,sm1( "[(CC)] system_variable (universalNumber) dc ")); |
if (k00setRingVariables(0,sm1( "[(N)] system_variable (universalNumber) dc "))) { |
k00setRingVariables(sm1( "[(C)] system_variable (universalNumber) dc "), |
sm1(" define_ring_variables "); |
sm1( "[(LL)] system_variable (universalNumber) dc ")); |
} |
k00setRingVariables(sm1( "[(L)] system_variable (universalNumber) dc "), |
|
sm1( "[(MM)] system_variable (universalNumber) dc ")); |
|
k00setRingVariables(sm1( "[(M)] system_variable (universalNumber) dc "), |
|
sm1( "[(NN)] system_variable (universalNumber) dc ")); |
|
if (SetRingVariables_Verbose) {Ln();} |
if (SetRingVariables_Verbose) {Ln();} |
} |
} |
|
|
Line 169 def k00AreThereLeftBrace(s) { |
|
Line 218 def k00AreThereLeftBrace(s) { |
|
if (jj != -1) return(true); else return(false); |
if (jj != -1) return(true); else return(false); |
} |
} |
|
|
def void k00setRingVariables(tmp002_p,tmp002_q) { |
def k00setRingVariables(p,q) { |
/* tmp002_ must not be used as variables names. */ |
local v,i; |
local tmp002_i,tmp002_v,tmp002_str; |
for (i = p; i< q; i++) { |
PSfor (tmp002_i=tmp002_p;tmp002_i<tmp002_q;tmp002_i++) { |
v = getxvar(i); |
tmp002_v = getxvar(tmp002_i); |
if (k00AreThereLeftBrace(v)) { |
if (k00AreThereLeftBrace(tmp002_v)) { ; } |
return(false); |
else { |
|
if (SetRingVariables_Verbose) {Print(tmp002_v); Print(" ");} |
|
str = AddString(["/",tmp002_v," $",tmp002_v,"$ (poly) data_conversion def "]); |
|
sm1("[(parse) ",str," ] extension "); |
|
} |
} |
tmp002_v = getdvar(tmp002_i); |
v = getdvar(i); |
if (k00AreThereLeftBrace(tmp002_v)) { ; } |
if (k00AreThereLeftBrace(v)) { |
else { |
return(false); |
if (SetRingVariables_Verbose) {Print(tmp002_v); Print(" ");} |
|
str = AddString(["/",tmp002_v," $",tmp002_v,"$ (poly) data_conversion def "]); |
|
sm1("[(parse) ",str," ] extension "); |
|
} |
} |
} |
} |
|
return(true); |
} |
} |
/* ---------------------------------- */ |
/* ---------------------------------- */ |
|
|
Line 359 HelpAdd(["GetPathName", |
|
Line 402 HelpAdd(["GetPathName", |
|
["GetPathName(s) checks if the file s exists in the current directory or", |
["GetPathName(s) checks if the file s exists in the current directory or", |
"in LOAD_K_PATH. If there exists, it returns the path name (string s)."]]); |
"in LOAD_K_PATH. If there exists, it returns the path name (string s)."]]); |
|
|
|
def Load_sm1(fnames,flag) { |
|
local ppp,n,i,cmd; |
|
if (Boundp(flag)) { |
|
}else{ |
|
n = Length(fnames); |
|
for (i=0; i<n; i++) { |
|
ppp = GetPathName(fnames[i]); |
|
if (Tag(ppp) != 0) { |
|
sm1(" [(parse) ppp pushfile ] extension "); |
|
cmd = AddString(["/",flag," 1 def "]); |
|
sm1(" [(parse) cmd ] extension "); |
|
i=n; /* break; */ |
|
} |
|
} |
|
} |
|
} |
|
|
|
HelpAdd(["Load_sm1", |
|
["Load_sm1(s,flag) loads a sm1 program from s[0], s[1], ....", |
|
"If loading is succeeded, the already-loaded flag is set to true.", |
|
"(list s, string flag)."]]); |
|
|
|
def GetRing(f) { |
|
sm1(" f getRing /FunctionValue set "); |
|
} |
|
|
|
def SetRing(r) { |
|
sm1(" r ring_def "); |
|
} |
|
|
|
def ReParse(a) { |
|
local c; |
|
if (IsArray(a)) { |
|
c = Map(a,"ReParse"); |
|
}else{ |
|
sm1(a," toString . /c set"); |
|
} |
|
return(c); |
|
} |
|
HelpAdd(["ReParse", |
|
["Reparse(obj): ", |
|
"It parses the given object in the current ring." |
|
]]); |
|
|
|
def void Pmat(a) { |
|
sm1(" a pmat "); |
|
} |
|
HelpAdd(["Pmat", |
|
["Pmat(m): ", |
|
"Print the array m in a pretty way." |
|
]]); |