===================================================================
RCS file: /home/cvs/OpenXM/src/kan96xx/Kan/dr.sm1,v
retrieving revision 1.5
retrieving revision 1.16
diff -u -p -r1.5 -r1.16
--- OpenXM/src/kan96xx/Kan/dr.sm1	2000/12/10 03:33:46	1.5
+++ OpenXM/src/kan96xx/Kan/dr.sm1	2003/07/17 12:11:09	1.16
@@ -1,4 +1,4 @@
-% $OpenXM: OpenXM/src/kan96xx/Kan/dr.sm1,v 1.4 2000/06/15 07:33:29 takayama Exp $
+% $OpenXM: OpenXM/src/kan96xx/Kan/dr.sm1,v 1.15 2003/07/17 09:10:54 takayama Exp $
 %% dr.sm1 (Define Ring) 1994/9/25, 26
 %% This file is error clean.
 
@@ -257,6 +257,8 @@
   /arg1 set
   [/vars /n /i /xList /dList /param] pushVariables
   [
+     (This is an obsolete macro. Use ring_of_differential_difference_operators)
+      error 
      (mmLarger) (matrix) switch_function
      (mpMult)   (difference) switch_function
      (red@)     (module1) switch_function
@@ -278,7 +280,44 @@
 } def
 
 
+/ring_of_differential_difference_operators {
+  /arg1 set
+  [/vars /n /i /xList /dList /param /dvar /evar /vars2 ] pushVariables
+  [
+     /vars arg1 def
+     vars tag 6 eq not {
+       ( List is expected as the argument for ring_of_differential_difference_operators ) error
+     } { } ifelse
+     vars 0 get /dvar set
+     vars 1 get /evar set
+     (mmLarger) (matrix) switch_function
+     (mpMult)   (difference) switch_function
+     (red@)     (module1) switch_function
+     (groebner) (standard) switch_function
+     (isSameComponent) (x) switch_function
 
+     [dvar to_records pop] /vars set %[x y z]
+     vars reverse /xList set         %[z y x]
+
+     [evar to_records pop] /vars2 set %[s1 s2]
+
+     vars2 reverse  {@@@.Esymbol 2 1 roll 2 cat_n} map
+     xList 
+     join /xList set   %[Es2 Es1 z y x]
+
+     vars2 reverse
+     vars {@@@.Dsymbol 2 1 roll 2 cat_n} map
+     reverse join /dList set              %[s2 s1 Dz Dy Dx]
+     [(H)] xList join [@@@.esymbol] join /xList set
+     [(h)] dList join [@@@.Esymbol] join /dList set
+     [0 1 1 << vars2 length 1 add >>  << xList length >>
+        1 1 << vars2 length 1 add >> << xList length 1 sub >> ] /param set
+     [ xList dList param ] /arg1 set
+  ] pop
+  popVariables
+  arg1
+} def
+
 /reverse {
   /arg1 set
   arg1 length 1 lt
@@ -1474,6 +1513,13 @@
   (type?) data_conversion  RationalFunctionP eq
 } def
 
+[(isRing)
+ [(obj isRing bool) ]
+] putUsages
+/isRing {
+  (type?) data_conversion  RingP eq
+} def
+
 /toString.tmp {
   /arg1 set
   [/obj /fname] pushVariables
@@ -2276,7 +2322,7 @@ newline
 
 /tolower.aux {
   /arg1 set
-  arg1 64 gt  arg1 96 lt and
+  arg1 64 gt  arg1 91 lt and
   { arg1 32 add }
   { arg1 } ifelse
 } def
@@ -2460,7 +2506,7 @@ newline
 [(ord_w)
  [(ff [v1 w1 v2 w2 ... vm wm] ord_w d)
   (poly ff; string v1; integer w1; ...)
-  (order of ff by the weight vector [w1 w2 ...])
+  (order of the initial of ff by the weight vector [w1 w2 ...])
   (Example: [(x,y) ring_of_polynomials 0] define_ring )
   (          (x^2 y^3-x). [(x) 2 (y) 1] ord_w ::)
  ]
@@ -2488,6 +2534,46 @@ newline
   arg1
 } def
 
+[(ord_w_all)
+ [(ff [v1 w1 v2 w2 ... vm wm] ord_w d)
+  (poly ff; string v1; integer w1; ...)
+  (order of ff by the weight vector [w1 w2 ...])
+  (Example: [(x,y,t) ring_of_polynomials 0] define_ring )
+  (          (x^2 y^3-x-t). [(t) 1 ] ord_w_all ::)
+ ]
+] putUsages 
+/ord_w_all {
+  /arg2 set /arg1 set
+  [/ord_w_all-in /fff /fff-in /www /rrr /iii /ddd /zzz /ddd-tmp] pushVariables
+  [
+    /fff arg1 def
+    /www arg2 def
+    fff (0). eq { /ddd -intInfinity def /ord_w_all.LLL goto} { } ifelse
+    /ddd -intInfinity def
+    fff (ring) dc /rrr set
+    /zzz (0) rrr ,, def
+    fff init /fff-in set
+    fff fff-in sub /fff set
+    {
+     /ddd-tmp 0 def
+     0 2 www length 1 sub {
+       /iii set
+       fff-in << www iii get rrr ,, >> degree
+       << www iii 1 add get >> mul
+       ddd-tmp add /ddd-tmp set
+     } for
+     ddd-tmp ddd gt { /ddd ddd-tmp def }  {  } ifelse
+     fff zzz eq { exit } {  } ifelse
+     fff init /fff-in set
+     fff fff-in sub /fff set
+    } loop
+    /ord_w_all.LLL
+    /arg1 ddd def
+  ] pop
+  popVariables
+  arg1
+} def
+
 [(laplace0)
  [
  (f [v1 ... vn] laplace0 g)
@@ -3562,14 +3648,15 @@ $ [ff ff] fromVectors :: $ 
 ]] putUsages
 
 /define_ring_variables {
-  [/in-define_ring_variables /v /p /v0] pushVariables
+  [/in-define_ring_variables /drv._v /drv._p /drv._v0] pushVariables
+%% You cannot use these names for names for polynomials.
   [
-     /v getVariableNames def
-     /v0 v def
-     v { dup /p set (/) 2 1 roll ( $) p ($. def ) } map cat
-     /v set
-%     v message
-     [(parse) v] extension 
+     /drv._v getVariableNames def
+     /drv._v0 drv._v def
+     drv._v { dup /drv._p set (/) 2 1 roll ( $) drv._p ($. def ) } map cat
+     /drv._v set
+%     drv._v message
+     [(parse) drv._v] extension 
   ] pop
   popVariables
 } def
@@ -3598,6 +3685,89 @@ $ [ff ff] fromVectors :: $ 
   (it returns 1 else it returns 0.)
   $Example: (hoge) boundp ::$
 ]] putUsages
+[(isSubstr) 
+ [
+  (s1 s2 isSubstr pos)
+  (If s1 is a substring of s2, isSubstr returns the position in s2 from which)
+  (s1 is contained in s2.)
+  (If s1 is not a substring of s2, then isSubstr returns -1.)
+ ]
+] putUsages
+/isSubstr {
+  /arg2 set /arg1 set
+  [/in-isSubstr /s1 /s2 /i1 /i2 /n1 /n2
+   /ans /flg
+  ]  pushVariables
+  [
+    /s1 arg1 def
+    /s2 arg2 def
+    s1 (array) dc /s1 set
+    s2 (array) dc /s2 set
+    /n1 s1 length def
+    /n2 s2 length def
+    /ans -1 def
+    0 1 n2 n1 sub {
+      /i2 set
+      /flg 1 def
+      0 1 n1 1 sub {
+        /i1 set
+        s1 i1 get s2 i2 i1 add get eq {
+        } { 
+          /flg 0 def exit
+        } ifelse
+      } for
+      flg {
+        /ans i2 def  
+        /isSubstr.L2 goto
+      } { /ans -1 def } ifelse
+    } for
+    /isSubstr.L2
+    /arg1 ans def
+  ] pop
+  popVariables
+  arg1
+} def
+
+[(execve)
+ [  
+   (command execve)
+   ([arg0 arg1 arg2 ...] execve )
+   (It executes the command by the system call execve.)
+   (cf. system, forkExec)
+ ] 
+] putUsages
+
+/execve {
+  /execve.arg set
+  [(forkExec) execve.arg [ ] 1] extension
+} def
+
+[(beginEcart)
+ [  
+   (beginEcart)
+   (Set the environments for the ecart division algorithm.)
+ ] 
+] putUsages
+
+/beginEcart {
+  (red@) (ecart) switch_function
+  [(Ecart) 1] system_variable
+  [(CheckHomogenization) 0] system_variable
+  [(ReduceLowerTerms) 0] system_variable
+  [(AutoReduce) 0] system_variable
+} def
+
+[(endEcart)
+ [  
+   (endEcart)
+   (End of using the ecart division algorithm.)
+ ] 
+] putUsages
+
+/endEcart {
+  (red@) (standard) switch_function
+  [(Ecart) 0] system_variable
+} def
 
 ;