OpenXM/src/kan96xx/Kan/dr.sm1 - diff

Return to dr.sm1 CVS log

Up to [local] / OpenXM / src / kan96xx / Kan

Diff for /OpenXM/src/kan96xx/Kan/dr.sm1 between version 1.3 and 1.10

-version 1.3, 2000/02/02 09:07:59
+version 1.10, 2002/11/08 14:18:01
 Line 1
 Line 1
 Line 1
- % $OpenXM: OpenXM/src/kan96xx/Kan/dr.sm1,v 1.2 1999/11/08 00:16:05 takayama Exp $
+ % $OpenXM: OpenXM/src/kan96xx/Kan/dr.sm1,v 1.9 2002/02/12 08:54:55 takayama Exp $
  %% dr.sm1 (Define Ring) 1994/9/25, 26
  %% This file is error clean.
  @@@.quiet {   }
- { (macro package : dr.sm1,   9/26,1995 --- Version 2/2, 2000. ) message } ifelse
+ { (macro package : dr.sm1,   9/26,1995 --- Version 12/10, 2000. ) message } ifelse
  /ctrlC-hook {
  %%% define your own routing in case of error.
 Line 257
 Line 257
 Line 257
    /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
-Line 278
+Line 280
 Line 278
 Line 280
  } 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 << 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
-Line 1474
+Line 1513
 Line 1474
 Line 1513
    (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
-Line 2276  newline
+Line 2322  newline
 Line 2276  newline
 Line 2322  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
-Line 3533  $ [ff ff] fromVectors :: $
+Line 3579  $ [ff ff] fromVectors :: $
 Line 3533  $ [ff ff] fromVectors :: $
 Line 3579  $ [ff ff] fromVectors :: $
  [(nl is the newline character.)
   $Example: [(You can break line) nl (here.)] cat message$
  ]] putUsages
+ /to_int {
+   /arg1 set
+   [/to-int /ob /ans] pushVariables
+   [
+     /ob arg1 def
+     /ans ob def
+     ob isArray {
+       ob {to_int} map /ans set
+       /LLL.to_int goto
+     } {  } ifelse
+     ob isInteger {
+       ob (universalNumber) dc /ans set
+       /LLL.to_int goto
+     } {  } ifelse
+     /LLL.to_int
+     /arg1 ans def
+   ] pop
+   popVariables
+   arg1
+ } def
+ [(to_int)
+ [(obj to_int obj2)
+  (All integers in obj are changed to universalNumber.)
+  (Example: /ff [1 2 [(hello) (0).]] def  ff { tag } map ::)
+  (         ff to_int { tag } map :: )
+ ]] putUsages
+ /define_ring_variables {
+   [/in-define_ring_variables /drv._v /drv._p /drv._v0] pushVariables
+ %% You cannot use these names for names for polynomials.
+   [
+      /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
+ [(define_ring_variables)
+ [(It binds  a variable <<a>> in the current ring to the sm1 variable <<a>>.)
+  (For example, if x is a variable in the current ring, it defines the sm1)
+  (variable x by /x (x) def)
+ ]] putUsages
+ /boundp {
+   /arg1 set
+   [/a /ans] pushVariables
+   [
+     /a arg1 def
+     [(parse) [(/) a ( load tag 0 eq { /ans 0 def } )
+                     (               { /ans 1 def } ifelse )] cat ] extension
+     /arg1 ans def
+   ] pop
+   popVariables
+   arg1
+ } def
+ [(boundp)
+  [( a boundp b)
+   (string a, b is 0 or 1.)
+   (If the variable named << a >> is bounded to a value,)
+   (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
+1 n2 n1 sub {
+       /i2 set
+       /flg 1 def
+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
  ;

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>