[BACK]Return to hol.k CVS log [TXT][DIR] Up to [local] / OpenXM / src / k097
File: [local] / OpenXM / src / k097 / hol.k (download)

Revision 1.2, Fri Jan 21 03:01:25 2000 UTC (24 years, 3 months ago) by takayama
Branch: MAIN
CVS Tags: maekawa-ipv6, R_1_3_1-2, RELEASE_20000124, RELEASE_1_3_1_13b, RELEASE_1_2_3_12, RELEASE_1_2_3, RELEASE_1_2_2_KNOPPIX_b, RELEASE_1_2_2_KNOPPIX, RELEASE_1_2_2, RELEASE_1_2_1, RELEASE_1_1_3, RELEASE_1_1_2, KNOPPIX_2006, HEAD, DEB_REL_1_2_3-9
Changes since 1.1: +1 -0 lines 

Added the tag $OpenXM$

/* $OpenXM: OpenXM/src/k097/hol.k,v 1.2 2000/01/21 03:01:25 takayama Exp $ */
/* hol.k,  1998, 12/14 */
/* packages for holonomic systems. */
sm1("[(chattr) 0 /npower ] extension ");
sm1(" [(parse) (cohom.sm1) pushfile] extension ");
sm1("[(chattr) 1 /npower] extension ");
sm1(" [(parse) (oxasir.sm1) pushfile] extension ");
/* class Hol extends Object { };
*/
   def Gkz(a,b) {
     local a1,b1;
     sm1(a," { { (integer) dc } map } map /a1 set ", b,
           " {(integer) dc} map  /b1 set [a1 b1] gkz /FunctionValue set ");
   }
   def Rrank(a) {
     sm1(a," rrank /FunctionValue set ");
   }
   def DeRham(a,b) {
     local a1,b1;
     a1 = a; b1 = b;
     if (IsArray(b1)) {
        b1 = Map(b1,"ToString");
        b1 = sm1(b1, " from_records ");
     }
     if (IsPolynomial(b1)) {
        b1 = ToString(b1);
     }
     if (IsPolynomial(a1)) {
        a1 = ToString(a1);
     }
     sm1("[ ", a1, b1, " ] deRham {(universalNumber) dc } map /FunctionValue set ");
   }
   def Annfs(a,b) {
     local a1,b1;
     a1 = a; b1 = b;
     if (IsArray(b1)) {
        b1 = Map(b1,"ToString");
        b1 = sm1(b1, " from_records ");
     }
     if (IsPolynomial(b1)) {
        b1 = ToString(b1);
     }
     if (IsPolynomial(a1)) {
        a1 = ToString(a1);
     }
     sm1("[ ", a1, b1, " ] annfs /FunctionValue set ");
   }
   def Fctr(a) {
     sm1(a," fctr /FunctionValue set ");
   }



HelpAdd(["Hol.",
 ["Annfs:  RingD(\"x,y\"); x=Poly(\"x\"); y=Poly(\"y\"); Annfs(x^3-y^2,[x,y]):   ",
  "Annfs:  Annfs(\"x^3-y^2\",\"x,y\"):   ",
  "Rrank:  Rrank(Gkz([[1,1,1,1],[0,1,3,4]],[1,2])):   ",
  "DeRham:  RingD(\"x,y\"); x=Poly(\"x\"); y=Poly(\"y\"); DeRham(x^3-y^2,[x,y]):   ",
  "Fctr:  Fctr(x^10-1):   ",
  "Gkz       "
]]);

OutputPrompt ;