Return to hol.k CVS log

Up to [local] / OpenXM / src / k097

Annotation of OpenXM/src/k097/hol.k, Revision 1.1.1.1

1.1       maekawa     1: /* hol.k,  1998, 12/14 */
                      2: /* packages for holonomic systems. */
                      3: sm1("[(chattr) 0 /npower ] extension ");
                      4: sm1(" [(parse) (cohom.sm1) pushfile] extension ");
                      5: sm1("[(chattr) 1 /npower] extension ");
                      6: sm1(" [(parse) (oxasir.sm1) pushfile] extension ");
                      7: /* class Hol extends Object { };
                      8: */
                      9:    def Gkz(a,b) {
                     10:      local a1,b1;
                     11:      sm1(a," { { (integer) dc } map } map /a1 set ", b,
                     12:            " {(integer) dc} map  /b1 set [a1 b1] gkz /FunctionValue set ");
                     13:    }
                     14:    def Rrank(a) {
                     15:      sm1(a," rrank /FunctionValue set ");
                     16:    }
                     17:    def DeRham(a,b) {
                     18:      local a1,b1;
                     19:      a1 = a; b1 = b;
                     20:      if (IsArray(b1)) {
                     21:         b1 = Map(b1,"ToString");
                     22:         b1 = sm1(b1, " from_records ");
                     23:      }
                     24:      if (IsPolynomial(b1)) {
                     25:         b1 = ToString(b1);
                     26:      }
                     27:      if (IsPolynomial(a1)) {
                     28:         a1 = ToString(a1);
                     29:      }
                     30:      sm1("[ ", a1, b1, " ] deRham {(universalNumber) dc } map /FunctionValue set ");
                     31:    }
                     32:    def Annfs(a,b) {
                     33:      local a1,b1;
                     34:      a1 = a; b1 = b;
                     35:      if (IsArray(b1)) {
                     36:         b1 = Map(b1,"ToString");
                     37:         b1 = sm1(b1, " from_records ");
                     38:      }
                     39:      if (IsPolynomial(b1)) {
                     40:         b1 = ToString(b1);
                     41:      }
                     42:      if (IsPolynomial(a1)) {
                     43:         a1 = ToString(a1);
                     44:      }
                     45:      sm1("[ ", a1, b1, " ] annfs /FunctionValue set ");
                     46:    }
                     47:    def Fctr(a) {
                     48:      sm1(a," fctr /FunctionValue set ");
                     49:    }
                     50:
                     51:
                     52:
                     53: HelpAdd(["Hol.",
                     54:  ["Annfs:  RingD(\"x,y\"); x=Poly(\"x\"); y=Poly(\"y\"); Annfs(x^3-y^2,[x,y]):   ",
                     55:   "Annfs:  Annfs(\"x^3-y^2\",\"x,y\"):   ",
                     56:   "Rrank:  Rrank(Gkz([[1,1,1,1],[0,1,3,4]],[1,2])):   ",
                     57:   "DeRham:  RingD(\"x,y\"); x=Poly(\"x\"); y=Poly(\"y\"); DeRham(x^3-y^2,[x,y]):   ",
                     58:   "Fctr:  Fctr(x^10-1):   ",
                     59:   "Gkz       "
                     60: ]]);
                     61:
                     62: OutputPrompt ;
                     63: