Return to hol.k CVS log

Up to [local] / OpenXM / src / k097

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

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