Return to cv0.sm1 CVS log

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

Annotation of OpenXM/src/kan96xx/Doc/cv0.sm1, Revision 1.1.1.1

1.1       maekawa     1:
                      2: /example1 {
                      3:   [(x,y,a,b,c,d)   %% change here
                      4:     ring_of_differential_operators
                      5:    (Dx,Dy)         %% change here
                      6:   elimination_order 0] swap01 define_ring
                      7:
                      8:   [[(x*Dx-a). (-b).] [(-c). ((x-1)*Dx-d).] [(Dy). (0).] [(0). (Dy).]] %% give equations
                      9:   /ff set
                     10:   ff { { [[$h$. $1$.]] replace } map } map /ff set
                     11:   /ff ff homogenize  def
                     12:   [ff] groebner /ans set
                     13:   ans 0 get {[[$h$. $1$.]] replace} map /gg set
                     14:   (Now, you get the characteristic variety) message
                     15:   (When your result is given as [I_1  I_2 ...], ) message
                     16:   (the characteristic variety is  V(I_1) \cup V(I_2) \cup ... ) message
                     17:   gg characteristic-v print (  ) message ( ) message
                     18: } def
                     19:
                     20:
                     21: /characteristic-v {
                     22:   /arg1 set
                     23:   [/gb  /lps /i /n /ans /maxp /ansp /k] pushVariables
                     24:   [  /gb arg1 def
                     25:      /ans [ ] def
                     26:      /maxp 0 def
                     27:      /lps gb {lpoint} map def
                     28:      0 1 << lps length 1 sub >>
                     29:      {
                     30:        /i set
                     31:        lps i get maxp gt
                     32:        { /maxp lps i get def }
                     33:        {  }
                     34:        ifelse
                     35:      } for
                     36:
                     37:      %%lps print
                     38:      /ans [
                     39:       0 1 maxp { pop [ ]   } for
                     40:      ] def
                     41:
                     42:      gb toVectors /gb set
                     43:
                     44:      0 1 << lps length 1 sub >>
                     45:      {
                     46:        /i set  /k lps i get def
                     47:        /ansp ans k get def
                     48:        << gb i get >> k  get principal /f set
                     49:        /ansp ansp [f] join def
                     50:        ans k ansp put
                     51:      } for
                     52:
                     53:      /arg1 ans def
                     54:   ] pop
                     55:   popVariables
                     56:   arg1
                     57: } def
                     58:
                     59: %%%%%%%%%%%%%%%%%%%%%%%%%%
                     60:
                     61:
                     62: ( ) message
                     63: (cv0.sm1 11/10, 1994. This program computes characteristic varieties.) message
                     64: (Type in example1 ; for a demo.) message
                     65:
                     66:
                     67: