Changeset a12b3f in git


Ignore:
Timestamp:
Aug 15, 1997, 10:00:37 AM (27 years ago)
Author:
Kai Krüger <krueger@…>
Branches:
(u'spielwiese', '73c62e0961bcfc8f1a00420b41eec2ea3c0ef6e9')
Children:
a21f1f2001491f06c9e455a96001e48e05143a88
Parents:
9f50b21dd037356d2991e1e8bb5f9f8324c06a76
Message:
* makedbm.lib pre-alpha


git-svn-id: file:///usr/local/Singular/svn/trunk@626 2c84dea3-7e68-4137-9b89-c4e89433aadc
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/makedbm.lib

    r9f50b2 ra12b3f  
    1 proc read_dbm (link l)
     1// $Id: makedbm.lib,v 1.2 1997-08-15 08:00:37 krueger Exp $
     2//=========================================================================
     3//
     4// Please send bugs and comments to krueger@mathematik.uni-kl.de
     5//
     6//=============================================================================
     7
     8LIBRARY:  makedbm.lib     some usefull tools needed by the Arnold-Classifier.
     9
     10   dbm_read(l);          read all entries from a DBM-databaes pointed by l
     11   dbm_getnext(l);       read next entry from a DBM-databaes pointed by l
     12
     13//=============================================================================
     14
     15proc dbm_read (link l)
    216{
    317  string s="";
     
    1024}
    1125
     26//=============================================================================
     27proc dbm_getnext (link l)
     28{
     29  string s="";
     30  s=read(l);
     31  if( s != "" ) { s,"=",read(l,s); }
     32}
     33
     34//=============================================================================
    1235proc create_sing_dbm
    1336{
    1437  link l="DBM:rw NFlist";
    1538  write(l, "A[k]", "x^(k+1)");
    16 write(l, "D[k]", "x2y+y^(k-1)");
    17 write(l, "E[6k]", "x3+y^(3*k+1)+a*x*(y^(2*k+1))");
    18 write(l, "E[6k+1]", "x3+x*(y^(2*k+1))+a*(y^(3*k+2))");
    19 write(l, "E[6k+2]", "x3+y^(3*k+2)+a*x*(y^(2*k+2))");
    20 write(l, "J[k,0]", "x3+b*x2*y^k+y^(3*k)c*x*y^(2*k+1)");
    21 write(l, "J[k,r]", "x3+x2*y^k+a*y^(3*k+r)");
    22 write(l, "X[1,0]", "x4+a*x2y2+y4");
    23 write(l, "X[1,r]", "x4+x2y2+a*y^(4+r)");
    24 write(l, "X[k,0]", "x4+b*x3y^k+a*x2y^(2k) + xy^(3k)");
    25 write(l, "X[k,r]", "x4+a*x3y^k+x2y^(2*k) + b*(y^(4*k+r))");
    26 write(l, "W[12k]", "x4+y^(4*k+1)+a*x*(y^(3*k+1))+c*x2*(y^(2*k+1))");
    27 write(l, "W[12k+1]", "x4+x*(y^(3*k+1))+a*x2*(y^(2*k+1))+c*y^(4*k+2)");
    28 write(l, "W[12k+5]", "x4+x*(y^(3*k+2))+a*x2*(y^(2*k+2))+b*y^(4*k+3)");
    29 write(l, "W[12k+6]", "x4+y^(4*k+3)+a*x*(y^(3*k+3))+b*x2*(y^(2*k+2))");
    30 write(l, "W[k,0]", "x4+b*x2*(y^(2*k+1))+a*x*(y^(3*k+2))+y^(4*k+2)");
    31 write(l, "W[k,r]", "x4+a*x3*(y^(k+1))+x2*(y^(2*k+1))+b*y^(4*k+2+r)");
    32 write(l, "W#[k,2r-1]", "(x2+y^(2*k+1))^2+b*x*(y^(3*k+1+r))+a*y^(4*k+2+r)");
    33 write(l, "W#[k,2r]", "(x2+y^(2*k+1))^2+b*x2*(y^(2*k+1+r))+a*x*(y^(3*k+2+r))");
    34 write(l, "Y[1,r,s]", "x^(4+r)+a*x2*(y2)+y^(4+s)");
    35 write(l, "Y[k,r,s]", "((x+a*y^k)^2 + b*y^(2*k+s))*( x2 + y^(2*k+r))");
    36 write(l, "Z[1,0]", "x3y + x2y3 + xy6 +y7");
    37 write(l, "Z[1,r]", "x3y + x2*(y^3) + a*y^(7+r)");
    38 write(l, "Z[k,r,s]", "(x2+a*x*(y^k)+b*y^(2*k+r))*(x2+y^(2*k+2*r+s))");
    39 write(l, "Z[k,r,0]", "(x+a*(y^k))*(x3+d*x2*(y^(k+r))+c*x*(y^(2*k+2*r+1))+y^(3*k+3*r))");
    40 write(l, "Z[k,12k+6r-1]", "(x+a*(y^k))*(x3+b*x*(y^(2*k+2*r+1))+y^(3*k+3*r+1))");
    41 write(l, "Z[k,12k+6r]", "(x+a*(y^k))*(x3+x*(y^(2*k+2*r+1))+b*y^(3*k+3*r+2))");
    42 write(l, "Z[k,12k+6r+1]", "(x+a*(y^k))*(x3+b*x*(y^(2*k+2*r+2))+y^(3*k+3*r+2))");
    43 write(l, "Z[k,0]", "y*(x3+d*x2*(y^k)+c*x*(y^(2*k+1))+y^(3*k))");
    44 write(l, "Z[k,r]", "y*(x3+x2y^(k+1)+b*(y^(3*k+r+3)))");
    45 write(l, "Z[6k+5]", "y*(x3+b*x*(y^(2*k+1))+y^(3*k+1))");
    46 write(l, "Z[6k+6]", "y*(x3+x*(y^(2*k+1))+b*y^(3*k+2))");
    47 write(l, "Z[6k+7]", "y*(x3+b*x*(y^(2*k+2))+y^(3*k+2))");
    48 write(l, "Q[k,0]", "x3+z2y+b*x2*(y^k)+x*(y^(2*k))");
    49 write(l, "Q[k,r]", "x3+z2y+x2*(y^k)+b*y^(3*k+r)");
    50 write(l, "Q[6k+4]", "x3+z2y+y^(3*k+1)+b*x*(y^(2*k+1))");
    51 write(l, "Q[6k+5]", "x3+z2y+x*(y^(2*k+1))+b*y^(3*k+2)");
    52 write(l, "Q[6k+6]", "x3+z2y+y^(3*k+2)+b*x*(y^(2*k+2))");
    53 write(l, "S[12k-1]", "x2z+z2y+y^(4*k)+a*x*(y^(3*k))+c*z*(y^(2*k+1))");
    54 write(l, "S[12k]", "x2z+z2y+x*(y^(3*k))+c*y^(4*k+1)+a*z*(y^(2*k+1))");
    55 write(l, "S[k,0]", "x2z+z2y+y^(4*k+1)+a*x*(y^(3*k+1))+b*z*(y^(2*k+1))");
    56 write(l, "S[k,r]", "x2z+z2y+x2*(y^(2*k))+a*x3*(y^k)+b*y^(4*k+r+1)");
    57 write(l, "S#[k,2r-1]", "x2z+z2y+z*(y^(2*k+1))+b*x*(y^(3*k+r))+a*(y^(4*k+r+1))");
    58 write(l, "S#[k,2r]", "x2z+z2y+z*(y^(2*k+1))+b*x2*(y^(2*k+r))");
    59 write(l, "S[12k+4]", "x2z+z2y+x*(y^(3*k+1))+a*z*(y^(2*k+2))+b*y^(4*k+2)");
    60 write(l, "S[12k+5]", "x2z+z2y+y^(4*k+2)+a*x*(y^(3*k+2))+b*z*(y^(2*k+2))");
    61 write(l, "U[12k]", "x3+z2x+y^(3*k+1)+a*x*(y^(2*k+1))+b*z*(y^(2*k+1))+d*x2*(y^(k+1))");
    62 write(l, "U[k,2r-1]", "x3+z2x+x*(y^(2*k+1))+a*x2*(y^(k+1))+b*(y^(3*k+r+2))+c*z*(y^(2*k+r+1))");
    63 write(l, "U[k,2r]", "x3+z2x+x*(y^(2*k+1))+a*x2*(y^(k+1))+b*z*(y^(2*k+r+1))+c*z2*(y^(k+r))");
    64 write(l, "U[12k+4]", "x3+z2x+y^(3*k+2)+a*x*(y^(2*k+2))+b*z*(y^(2*k+2))+c*x2*(y^(k+1))");
    65 write(l, "V[1,0]", "x2y+z4+a*z3y+b*z2y2+y3z");
    66 write(l, "V[1,r]", "x2y+z4+b*z3y+z2y2+a*(y^(r+4))");
    67 write(l, "V#[1,2r-1]", "x2y+z3y+a*z2y2+y4+b*x*(z^(r+2))");
    68 write(l, "V#[1,2r]", "x2y+z3y+a*z2y2+y4+b*(z^(r+4))");
    69 write(l, "T[k,r,s]", "x^k+y^r+z^s+xyz");
     39  write(l, "D[k]", "x2y+y^(k-1)");
     40  write(l, "E[6k]", "x3+y^(3*k+1)+a*x*(y^(2*k+1))");
     41  write(l, "E[6k+1]", "x3+x*(y^(2*k+1))+a*(y^(3*k+2))");
     42  write(l, "E[6k+2]", "x3+y^(3*k+2)+a*x*(y^(2*k+2))");
     43  write(l, "J[k,0]", "x3+b*x2*y^k+y^(3*k)c*x*y^(2*k+1)");
     44  write(l, "J[k,r]", "x3+x2*y^k+a*y^(3*k+r)");
     45  write(l, "X[1,0]", "x4+a*x2y2+y4");
     46  write(l, "X[1,r]", "x4+x2y2+a*y^(4+r)");
     47  write(l, "X[k,0]", "x4+b*x3y^k+a*x2y^(2k) + xy^(3k)");
     48  write(l, "X[k,r]", "x4+a*x3y^k+x2y^(2*k) + b*(y^(4*k+r))");
     49  write(l, "W[12k]", "x4+y^(4*k+1)+a*x*(y^(3*k+1))+c*x2*(y^(2*k+1))");
     50  write(l, "W[12k+1]", "x4+x*(y^(3*k+1))+a*x2*(y^(2*k+1))+c*y^(4*k+2)");
     51  write(l, "W[12k+5]", "x4+x*(y^(3*k+2))+a*x2*(y^(2*k+2))+b*y^(4*k+3)");
     52  write(l, "W[12k+6]", "x4+y^(4*k+3)+a*x*(y^(3*k+3))+b*x2*(y^(2*k+2))");
     53  write(l, "W[k,0]", "x4+b*x2*(y^(2*k+1))+a*x*(y^(3*k+2))+y^(4*k+2)");
     54  write(l, "W[k,r]", "x4+a*x3*(y^(k+1))+x2*(y^(2*k+1))+b*y^(4*k+2+r)");
     55  write(l, "W#[k,2r-1]", "(x2+y^(2*k+1))^2+b*x*(y^(3*k+1+r))+a*y^(4*k+2+r)");
     56  write(l, "W#[k,2r]", "(x2+y^(2*k+1))^2+b*x2*(y^(2*k+1+r))+a*x*(y^(3*k+2+r))");
     57  write(l, "Y[1,r,s]", "x^(4+r)+a*x2*(y2)+y^(4+s)");
     58  write(l, "Y[k,r,s]", "((x+a*y^k)^2 + b*y^(2*k+s))*( x2 + y^(2*k+r))");
     59  write(l, "Z[1,0]", "x3y + x2y3 + xy6 +y7");
     60  write(l, "Z[1,r]", "x3y + x2*(y^3) + a*y^(7+r)");
     61  write(l, "Z[k,r,s]", "(x2+a*x*(y^k)+b*y^(2*k+r))*(x2+y^(2*k+2*r+s))");
     62  write(l, "Z[k,r,0]", "(x+a*(y^k))*(x3+d*x2*(y^(k+r))+c*x*(y^(2*k+2*r+1))+y^(3*k+3*r))");
     63  write(l, "Z[k,12k+6r-1]", "(x+a*(y^k))*(x3+b*x*(y^(2*k+2*r+1))+y^(3*k+3*r+1))");
     64  write(l, "Z[k,12k+6r]", "(x+a*(y^k))*(x3+x*(y^(2*k+2*r+1))+b*y^(3*k+3*r+2))");
     65  write(l, "Z[k,12k+6r+1]", "(x+a*(y^k))*(x3+b*x*(y^(2*k+2*r+2))+y^(3*k+3*r+2))");
     66  write(l, "Z[k,0]", "y*(x3+d*x2*(y^k)+c*x*(y^(2*k+1))+y^(3*k))");
     67  write(l, "Z[k,r]", "y*(x3+x2y^(k+1)+b*(y^(3*k+r+3)))");
     68  write(l, "Z[6k+5]", "y*(x3+b*x*(y^(2*k+1))+y^(3*k+1))");
     69  write(l, "Z[6k+6]", "y*(x3+x*(y^(2*k+1))+b*y^(3*k+2))");
     70  write(l, "Z[6k+7]", "y*(x3+b*x*(y^(2*k+2))+y^(3*k+2))");
     71  write(l, "Q[k,0]", "x3+z2y+b*x2*(y^k)+x*(y^(2*k))");
     72  write(l, "Q[k,r]", "x3+z2y+x2*(y^k)+b*y^(3*k+r)");
     73  write(l, "Q[6k+4]", "x3+z2y+y^(3*k+1)+b*x*(y^(2*k+1))");
     74  write(l, "Q[6k+5]", "x3+z2y+x*(y^(2*k+1))+b*y^(3*k+2)");
     75  write(l, "Q[6k+6]", "x3+z2y+y^(3*k+2)+b*x*(y^(2*k+2))");
     76  write(l, "S[12k-1]", "x2z+z2y+y^(4*k)+a*x*(y^(3*k))+c*z*(y^(2*k+1))");
     77  write(l, "S[12k]", "x2z+z2y+x*(y^(3*k))+c*y^(4*k+1)+a*z*(y^(2*k+1))");
     78  write(l, "S[k,0]", "x2z+z2y+y^(4*k+1)+a*x*(y^(3*k+1))+b*z*(y^(2*k+1))");
     79  write(l, "S[k,r]", "x2z+z2y+x2*(y^(2*k))+a*x3*(y^k)+b*y^(4*k+r+1)");
     80  write(l, "S#[k,2r-1]", "x2z+z2y+z*(y^(2*k+1))+b*x*(y^(3*k+r))+a*(y^(4*k+r+1))");
     81  write(l, "S#[k,2r]", "x2z+z2y+z*(y^(2*k+1))+b*x2*(y^(2*k+r))");
     82  write(l, "S[12k+4]", "x2z+z2y+x*(y^(3*k+1))+a*z*(y^(2*k+2))+b*y^(4*k+2)");
     83  write(l, "S[12k+5]", "x2z+z2y+y^(4*k+2)+a*x*(y^(3*k+2))+b*z*(y^(2*k+2))");
     84  write(l, "U[12k]", "x3+z2x+y^(3*k+1)+a*x*(y^(2*k+1))+b*z*(y^(2*k+1))+d*x2*(y^(k+1))");
     85  write(l, "U[k,2r-1]", "x3+z2x+x*(y^(2*k+1))+a*x2*(y^(k+1))+b*(y^(3*k+r+2))+c*z*(y^(2*k+r+1))");
     86  write(l, "U[k,2r]", "x3+z2x+x*(y^(2*k+1))+a*x2*(y^(k+1))+b*z*(y^(2*k+r+1))+c*z2*(y^(k+r))");
     87  write(l, "U[12k+4]", "x3+z2x+y^(3*k+2)+a*x*(y^(2*k+2))+b*z*(y^(2*k+2))+c*x2*(y^(k+1))");
     88  write(l, "V[1,0]", "x2y+z4+a*z3y+b*z2y2+y3z");
     89  write(l, "V[1,r]", "x2y+z4+b*z3y+z2y2+a*(y^(r+4))");
     90  write(l, "V#[1,2r-1]", "x2y+z3y+a*z2y2+y4+b*x*(z^(r+2))");
     91  write(l, "V#[1,2r]", "x2y+z3y+a*z2y2+y4+b*(z^(r+4))");
     92  write(l, "T[k,r,s]", "x^k+y^r+z^s+xyz");
    7093
    7194//////////////////////////////////////////////////////////////////
    7295// DatenFormat: crk=#; Mu=#; MlrCd=#;
    7396string s;
    74 s ="crk=1; Mu=k; MlnCd=k;";
    75 write(l, "I_A[k]", s);
    76 s = "crk=2; Mu=k; MlnCd=1,1,k-3";
    77 write(l, "I_D[k]", s);
    78 s = "crk=2; Mu=6*k; MlnCd=1,2*k,2*k-1";
    79 write(l, "I_E[6k]", s);
    80 s = "crk=2; Mu=6*k+1; MlnCd=1,2*k,2*k";
    81 write(l, "I_E[6k+1]", s);
     97  s ="crk=1; Mu=k; MlnCd=k;";
     98  write(l, "I_A[k]", s);
     99  s = "crk=2; Mu=k; MlnCd=1,1,k-3";
     100  write(l, "I_D[k]", s);
     101  s = "crk=2; Mu=6*k; MlnCd=1,2*k,2*k-1";
     102  write(l, "I_E[6k]", s);
     103  s = "crk=2; Mu=6*k+1; MlnCd=1,2*k,2*k";
     104  write(l, "I_E[6k+1]", s);
    82105//"I_E[6k+1]=", read(l, "I_E[6k+1]");
    83 s = "crk=2; Mu=6*k+2; MlnCd=1,2*k+1,2*k-1";
    84 write(l, "I_E[6k+2]", s);
    85 s = "crk=2; Mu=6*k-2; MlnCd=1,2*k-1,2*k-1";
    86 write(l, "I_J[k,0]", s);
     106  s = "crk=2; Mu=6*k+2; MlnCd=1,2*k+1,2*k-1";
     107  write(l, "I_E[6k+2]", s);
     108  s = "crk=2; Mu=6*k-2; MlnCd=1,2*k-1,2*k-1";
     109  write(l, "I_J[k,0]", s);
    87110//"I_J[k,0]=", read(l, "I_J[k,0]");
    88 s = "crk=2; Mu=6*k-2+r; MlnCd=1,2*k-1,2*k-1+r";
    89 write(l, "I_J[k,r]", s);
     111  s = "crk=2; Mu=6*k-2+r; MlnCd=1,2*k-1,2*k-1+r";
     112  write(l, "I_J[k,r]", s);
    90113//"I_J[k,r]=", read(l, "I_J[k,r]");
    91 s = "crk=2; Mu=9; MlnCd=1,1,1,1,1";
    92 write(l, "I_X[1,0]", s);
    93 s = "crk=2; Mu=9+r; MlnCd=1,1,1,1,1+r";
    94 write(l, "I_X[1,r]", s);
    95 s = "crk=2; Mu=12*k-3; MlnCd=1,1,2*k-1,2*k-1,2*k-1";
    96 write(l, "I_X[k,0]", s);
    97 s = "crk=2; Mu=12*k-3+r; MlnCd=1,1,2*k-1,2*k-1,2*k-1+r";
    98 write(l, "I_X[k,r]", s);
    99 s = "crk=2; Mu=12*k;";
    100 write(l, "I_W[12k]", s);
    101 s = "crk=2; Mu=12*k+1;";
    102 write(l, "I_W[12k+1]", s);
    103 s = "crk=2; Mu=12*k+5;";
    104 write(l, "I_W[12k+5]", s);
    105 s = "crk=2; Mu=12*k+6;";
    106 write(l, "I_W[12k+6]", s);
    107 s = "crk=2; Mu=12*k+3;";
    108 write(l, "I_W[k,0]", s);
    109 s = "crk=2; Mu=12*k+3+r;";
    110 write(l, "I_W[k,r]", s);
    111 s = "crk=2; Mu=12*k+2+2*r;";
    112 write(l, "I_W#[k,2r-1]", s);
    113 s = "crk=2; Mu=12*k+3+2*r;";
    114 write(l, "I_W#[k,2r]", s);
    115 s = "crk=2; Mu=9+r+s;";
    116 write(l, "I_Y[1,r,s]", s);
    117 s = "crk=2; Mu=12*k-3+r+s;";
    118 write(l, "I_Y[k,r,s]", s);
    119 s = "crk=2; Mu=15;";
    120 write(l, "I_Z[1,0]", s);
    121 s = "crk=2; Mu=15+r;";
    122 write(l, "I_Z[1,r]", s);
    123 s = "crk=2; Mu=9+6*k+r;";
    124 write(l, "I_Z[k,r]", s);
    125 s = "crk=2; Mu=12*k+6*r-3;";
    126 write(l, "I_Z[k,r,0]", s);
    127 s = "crk=2; Mu=12*k+6*r+s-3;";
    128 write(l, "I_Z[k,r,s]", s);
    129 s = "crk=2; Mu=12*k+6*r-1;";
    130 write(l, "I_Z[k,12k+6r-1]", s);
    131 s = "crk=2; Mu=12*k+6*r;";
    132 write(l, "I_Z[k,12k+6r]", s);
    133 s = "crk=2; Mu=12*k+6*r+1;";
    134 write(l, "I_Z[k,12k+6r+1]", s);
    135 s = "crk=2; Mu=9+6*k;";
    136 write(l, "I_Z[k,0]", s);
    137 s = "crk=2; Mu=6*(r+1)-1;";
    138 write(l, "I_Z[6k+5]", s);
    139 s = "crk=2; Mu=6*(r+1);";
    140 write(l, "I_Z[6k+6]", s);
    141 s = "crk=2; Mu=6*(r+1)+1;";
    142 write(l, "I_Z[6k+7]", s);
    143 s = "crk=3; Mu=6*k+2;";
    144 write(l, "I_Q[k,0]", s);
    145 s = "crk=3; Mu=6*k+2+r;";
    146 write(l, "I_Q[k,r]", s);
    147 s = "crk=3; Mu=6*k+4;";
    148 write(l, "I_Q[6k+4]", s);
    149 s = "crk=3; Mu=6*k+5;";
    150 write(l, "I_Q[6k+5]", s);
    151 s = "crk=3; Mu=6*k+6;";
    152 write(l, "I_Q[6k+6]", s);
    153 s = "crk=3; Mu=12*k-1;";
    154 write(l, "I_S[12k-1]", s);
    155 s = "crk=3; Mu=12*k;";
    156 write(l, "I_S[12k]", s);
    157 s = "crk=3; Mu=12*k+2;";
    158 write(l, "I_S[k,0]", s);
    159 s = "crk=3; Mu=12*k+2+r;";
    160 write(l, "I_S[k,r]", s);
    161 s = "crk=3; Mu=12*k+2*r+1;";
    162 write(l, "I_S#[k,2r-1]", s);
    163 s = "crk=3; Mu=12*k+2*r+2;";
    164 write(l, "I_S#[k,2r]", s);
    165 s = "crk=3; Mu=12*k+4;";
    166 write(l, "I_S[12k+4]", s);
    167 s = "crk=3; Mu=12*k+5;";
    168 write(l, "I_S[12k+5]", s);
    169 s = "crk=3; Mu=12*k;";
    170 write(l, "I_U[12k]", s);
    171 s = "crk=3; Mu=12*k+4;";
    172 write(l, "I_U[12k+4]", s);
    173 s = "crk=3; Mu=12*k+1+2*r;";
    174 write(l, "I_U[k,2r-1]", s);
    175 s = "crk=3; Mu=12*k+2+2*r;";
    176 write(l, "I_U[k,2r]", s);
    177 s = "crk=3; Mu=15;";
    178 write(l, "I_V[1,0]", s);
    179 s = "crk=3; Mu=15+r;";
    180 write(l, "I_V[1,r]", s);
    181 s = "crk=3; Mu=14+2*r;";
    182 write(l, "I_V#[1,2r-1]", s);
    183 s = "crk=3; Mu=15+2*r;";
    184 write(l, "I_V#[1,2r]", s);
    185 s = "crk=3; Mu=0;";
    186 write(l, "I_T[k,r,s]", s);
    187 read_sing_dbm();
    188 close(l);
     114  s = "crk=2; Mu=9; MlnCd=1,1,1,1,1";
     115  write(l, "I_X[1,0]", s);
     116  s = "crk=2; Mu=9+r; MlnCd=1,1,1,1,1+r";
     117  write(l, "I_X[1,r]", s);
     118  s = "crk=2; Mu=12*k-3; MlnCd=1,1,2*k-1,2*k-1,2*k-1";
     119  write(l, "I_X[k,0]", s);
     120  s = "crk=2; Mu=12*k-3+r; MlnCd=1,1,2*k-1,2*k-1,2*k-1+r";
     121  write(l, "I_X[k,r]", s);
     122  s = "crk=2; Mu=12*k;";
     123  write(l, "I_W[12k]", s);
     124  s = "crk=2; Mu=12*k+1;";
     125  write(l, "I_W[12k+1]", s);
     126  s = "crk=2; Mu=12*k+5;";
     127  write(l, "I_W[12k+5]", s);
     128  s = "crk=2; Mu=12*k+6;";
     129  write(l, "I_W[12k+6]", s);
     130  s = "crk=2; Mu=12*k+3;";
     131  write(l, "I_W[k,0]", s);
     132  s = "crk=2; Mu=12*k+3+r;";
     133  write(l, "I_W[k,r]", s);
     134  s = "crk=2; Mu=12*k+2+2*r;";
     135  write(l, "I_W#[k,2r-1]", s);
     136  s = "crk=2; Mu=12*k+3+2*r;";
     137  write(l, "I_W#[k,2r]", s);
     138  s = "crk=2; Mu=9+r+s;";
     139  write(l, "I_Y[1,r,s]", s);
     140  s = "crk=2; Mu=12*k-3+r+s;";
     141  write(l, "I_Y[k,r,s]", s);
     142  s = "crk=2; Mu=15;";
     143  write(l, "I_Z[1,0]", s);
     144  s = "crk=2; Mu=15+r;";
     145  write(l, "I_Z[1,r]", s);
     146  s = "crk=2; Mu=9+6*k+r;";
     147  write(l, "I_Z[k,r]", s);
     148  s = "crk=2; Mu=12*k+6*r-3;";
     149  write(l, "I_Z[k,r,0]", s);
     150  s = "crk=2; Mu=12*k+6*r+s-3;";
     151  write(l, "I_Z[k,r,s]", s);
     152  s = "crk=2; Mu=12*k+6*r-1;";
     153  write(l, "I_Z[k,12k+6r-1]", s);
     154  s = "crk=2; Mu=12*k+6*r;";
     155  write(l, "I_Z[k,12k+6r]", s);
     156  s = "crk=2; Mu=12*k+6*r+1;";
     157  write(l, "I_Z[k,12k+6r+1]", s);
     158  s = "crk=2; Mu=9+6*k;";
     159  write(l, "I_Z[k,0]", s);
     160  s = "crk=2; Mu=6*(r+1)-1;";
     161  write(l, "I_Z[6k+5]", s);
     162  s = "crk=2; Mu=6*(r+1);";
     163  write(l, "I_Z[6k+6]", s);
     164  s = "crk=2; Mu=6*(r+1)+1;";
     165  write(l, "I_Z[6k+7]", s);
     166  s = "crk=3; Mu=6*k+2;";
     167  write(l, "I_Q[k,0]", s);
     168  s = "crk=3; Mu=6*k+2+r;";
     169  write(l, "I_Q[k,r]", s);
     170  s = "crk=3; Mu=6*k+4;";
     171  write(l, "I_Q[6k+4]", s);
     172  s = "crk=3; Mu=6*k+5;";
     173  write(l, "I_Q[6k+5]", s);
     174  s = "crk=3; Mu=6*k+6;";
     175  write(l, "I_Q[6k+6]", s);
     176  s = "crk=3; Mu=12*k-1;";
     177  write(l, "I_S[12k-1]", s);
     178  s = "crk=3; Mu=12*k;";
     179  write(l, "I_S[12k]", s);
     180  s = "crk=3; Mu=12*k+2;";
     181  write(l, "I_S[k,0]", s);
     182  s = "crk=3; Mu=12*k+2+r;";
     183  write(l, "I_S[k,r]", s);
     184  s = "crk=3; Mu=12*k+2*r+1;";
     185  write(l, "I_S#[k,2r-1]", s);
     186  s = "crk=3; Mu=12*k+2*r+2;";
     187  write(l, "I_S#[k,2r]", s);
     188  s = "crk=3; Mu=12*k+4;";
     189  write(l, "I_S[12k+4]", s);
     190  s = "crk=3; Mu=12*k+5;";
     191  write(l, "I_S[12k+5]", s);
     192  s = "crk=3; Mu=12*k;";
     193  write(l, "I_U[12k]", s);
     194  s = "crk=3; Mu=12*k+4;";
     195  write(l, "I_U[12k+4]", s);
     196  s = "crk=3; Mu=12*k+1+2*r;";
     197  write(l, "I_U[k,2r-1]", s);
     198  s = "crk=3; Mu=12*k+2+2*r;";
     199  write(l, "I_U[k,2r]", s);
     200  s = "crk=3; Mu=15;";
     201  write(l, "I_V[1,0]", s);
     202  s = "crk=3; Mu=15+r;";
     203  write(l, "I_V[1,r]", s);
     204  s = "crk=3; Mu=14+2*r;";
     205  write(l, "I_V#[1,2r-1]", s);
     206  s = "crk=3; Mu=15+2*r;";
     207  write(l, "I_V#[1,2r]", s);
     208  s = "crk=3; Mu=0;";
     209  write(l, "I_T[k,r,s]", s);
     210  close(l);
    189211}
    190212
    191213proc read_sing_dbm
    192214{
    193   link l="DBM:rw NFlist";
     215  link l="DBM: NFlist";
    194216  read(l, "A[k]");
    195217  read(l, "D[k]");
     
    247269  read(l, "V#[1,2r]");
    248270  read(l, "T[k,r,s]");
    249 }
     271  close(l);
     272}
Note: See TracChangeset for help on using the changeset viewer.