Changeset 2ed3a3 in git for Singular/LIB/matrix.lib

Timestamp:

Nov 20, 2008, 6:17:33 PM (15 years ago)

Author:

Motsak Oleksandr <motsak@…>

Branches:

(u'spielwiese', 'ec94ef7a30b928574c0c3daf41f6804dff5f6b69')

Children:

c16f9922fe5d93bd72e4b2440cd6c7cd43889ca6

Parents:

2a03adf5a9e042150758e454f9153509b05e3633

Message:

*motsak: reviewed exterior/symmetric powers.


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

File:

• Singular/LIB/matrix.lib (modified) (9 diffs)

Singular/LIB/matrix.lib

@@ -1 +1 @@
 ///////////////////////////////////////////////////////////////////////////////
-version="$Id: matrix.lib,v 1.43 2008-10-01 15:29:23 Singular Exp $";
+version="$Id: matrix.lib,v 1.44 2008-11-20 17:17:33 motsak Exp $";
 category="Linear Algebra";
 info="
@@ -34 +34 @@
  rm_unitcol(A);         remove unit columns and associated rows of A
  headStand(A);          A[n-i+1,m-j+1]:=A[i,j]
- symmetricBasis(n,k);   ring, with a basis of a k-th symmetric power of an n-dim v.space
- symmetricPower(A,k);   module, k-th symmetric power of a module/matrix A
- exteriorBasis(n,k);    ring, with a basis of a k-th exterior power of an n-dim v.space
- exteriorPower(A,k);    module, k-th exterior power of a module/matrix A 
+ symmetricBasis(n,k[,s]); basis of k-th symmetric power of n-dim v.space
+ exteriorBasis(n,k[,s]); basis of k-th exterior power of n-dim v.space
+ symmetricPower(A,k);   k-th symmetric power of a module/matrix A
+ exteriorPower(A,k);    k-th exterior power of a module/matrix A 
           (parameters in square brackets [] are optional)
 ";
@@ -1074 +1074 @@
 //////////////////////////////////////////////////////////////////////////////
 
-// Symmetric/Exterior powers, thanks to Oleksandr Iena for his persistence ;-)
-
-proc symmetricBasis(int n, int k)
-"USAGE:    symmetricBasis(n, k); n int, k int
-RETURN:   ring: the basis of the k-th symmetric power of a n-dim vector space
-NOTE:     The output consists of a polynomial ring in n variables, together with the ideal called I.
-          The ideal I is the basis of the k-th symmetric power of a n-dim vector space (ordered lexicographically).
-          The char. of the returned ring is the same as of current basis ring or zero.
+// Symmetric/Exterior powers thanks to Oleksandr Iena for his persistence ;-)
+
+proc symmetricBasis(int n, int k, list #)
+"USAGE:    symmetricBasis(n, k[,s]); n int, k int, s string
+RETURN:   poynomial ring containing the ideal \"symBasis\", 
+          being a basis of the k-th symmetric power of an n-dim vector space.
+NOTE:     The output polynomial ring has characteristics 0 and n variables 
+          named \"S(i)\", where the base variable name S is either given by the
+          optional string argument(which must not contain brackets) or equal to
+          "e" by default.
 SEE ALSO: exteriorBasis
 KEYWORDS: symmetric basis
 EXAMPLE:  example symmetricBasis; shows an example"
 {
-  int p = 0;
-  if( nameof(basering) != "basering" )
+//------------------------ handle optional base variable name---------------
+  string S = "e";
+  if( size(#) > 0 )
   {
-    p = char(basering);
+    if( typeof(#[1]) != "string" )
+    {
+      ERROR("Wrong optional argument: must be a string");
+    }
+    S = #[1];
+    if( (find(S, "(") + find(S, ")")) > 0 )
+    {
+      ERROR("Wrong optional argument: must be a string without brackets"); 
+    }
   }
 
-  ring T = (p), (e(1..n)), dp;
-  ideal I = simplify( NF(maxideal(k), std(0)), 1 + 2 + 8 );
-  int N = ncols(I);
-  I = sort(I)[1]; // lex
-  I = I[N..1];
-
-  export I;
-  return(T);
-}
-example
-{ "EXAMPLE:"; echo = 2;
-
-def r = symmetricBasis(2, 3); setring r; r;
-I; // basis of 3rd sym. power of a 2-dim v.s.
-kill r;
-
-ring r = (0),(a, b, c, d), dp; r;
-def g = symmetricBasis(3, 2); setring g; g; I;
-
-kill g, r;
-
-ring r = (32003),(a, b, c, d), dp; r;
-def g = symmetricBasis(4, 2); setring g; g; I;
+//------------------------- create ring container for symmetric power basis-
+  execute("ring @@@SYM_POWER_RING_NAME=(0),("+S+"(1.."+string(n)+")),dp;");
+
+//------------------------- choose symmetric basis -------------------------
+  ideal symBasis = maxideal(k);
+
+//------------------------- export and return      -------------------------
+  export symBasis;
+  return(basering);
+}
+example
+{ "EXAMPLE:"; echo = 2;
+
+// basis of the 3-rd symmetricPower of a 4-dim vector space:
+def R = symmetricBasis(4, 3, "@e"); setring R; 
+R;  // container ring:
+symBasis; // symmetric basis:
 }
 
 //////////////////////////////////////////////////////////////////////////////
 
-proc exteriorBasis(int n, int k)
-"USAGE:    exteriorBasis(n, k); n int, k int
-RETURN:   ring: the basis of the k-th exterior power of a n-dim vector space
-NOTE:     The output consists of a polynomial ring in n variables, together with the ideal called I.
-          The ideal I is the basis of the k-th exterior power of a n-dim vector space (ordered lexicographically).
-          The char. of the returned ring is the same as of current basis ring or zero.
-
+proc exteriorBasis(int n, int k, list #)
+"USAGE:    exteriorBasis(n, k[,s]); n int, k int, s string
+RETURN:   poynomial ring containing the ideal \"extBasis\", 
+          being a basis of the k-th exterior power of an n-dim vector space.
+NOTE:     The output polynomial ring has characteristics 0 and n variables 
+          named \"S(i)\", where the base variable name S is either given by the
+          optional string argument(which must not contain brackets) or equal to
+          "e" by default.
 SEE ALSO: symmetricBasis
 KEYWORDS: exterior basis
 EXAMPLE:  example exteriorBasis; shows an example"
 {
-  int p = 0;
-  if( nameof(basering) != "basering" )
+//------------------------ handle optional base variable name---------------
+  string S = "e";
+  if( size(#) > 0 )
   {
-    p = char(basering);
+    if( typeof(#[1]) != "string" )
+    {
+      ERROR("Wrong optional argument: must be a string");
+    }
+    S = #[1];
+    if( (find(S, "(") + find(S, ")")) > 0 )
+    {
+      ERROR("Wrong optional argument: must be a string without brackets"); 
+    }
   }
-    
-  ring S = (p), (e(1..n)), dp;
+
+//------------------------- create ring container for symmetric power basis-
+  execute("ring @@@EXT_POWER_RING_NAME=(0),("+S+"(1.."+string(n)+")),dp;");
+
+//------------------------- choose exterior basis -------------------------
   def T = SuperCommutative(); setring T;
-  ideal I = simplify( NF(maxideal(k), std(0)), 1 + 2 + 8 );
-  int N = ncols(I);
-  I = sort(I)[1];
-  I = I[N..1 ];
-
-  export I;
-  return(T);
-}
-example
-{ "EXAMPLE:"; echo = 2;
-
-def r = exteriorBasis(2, 3); setring r; r;
-"Basis: ", I; // basis of 3rd sym. power of a 2-dim v.s.
-kill r;
-
-ring r = (0),(a, b, c, d), dp; r;
-def g = exteriorBasis(3, 2); setring g; g; "Basis: ", I;
-
-kill g, r;
-
-ring r = (32003),(a, b, c, d), dp; r;
-def g = exteriorBasis(4, 2); setring g; g; "Basis: ", I;
-}
-
-//////////////////////////////////////////////////////////////////////////////
-
-proc symmetricPower(module A, int k)
-"USAGE:    symmetricPower(A, k); A module, k int
-RETURN:   module: the k-th symmetric power of A
-NOTE:     the chosen bases (ordered lexicographically) and
-          temporary data  may be shown if printlevel is big enough
-SEE ALSO: exteriorPower
-KEYWORDS: symmetric power
-EXAMPLE:  example symmetricPower; shows an example"
-{
-  int p = printlevel - voice + 2;
+  ideal extBasis = simplify( NF(maxideal(k), std(0)), 1 + 2 + 8 );
+
+//------------------------- export and return      -------------------------
+  export extBasis;
+  return(basering);
+}
+example
+{ "EXAMPLE:"; echo = 2;
+// basis of the 3-rd symmetricPower of a 4-dim vector space:
+def r = exteriorBasis(4, 3, "@e"); setring r; 
+r; // container ring:
+extBasis; // exterior basis:
+}
+
+
+static proc chooseSafeVarName(string prefix, string suffix)
+"USAGE: give appropreate prefix for variable names
+RETURN: safe variable name (repeated prefix + suffix)
+"
+{
+  string V = varstr(basering);
+  string S = suffix;
+  while( find(V, S) > 0 )
+  {
+    S = prefix + S;
+  }
+  return(S);
+}
+
+static proc mapPower(int p, module A, int k, def Tn, def Tm)
+"USAGE: by both symmetric- and exterior-Power"
+NOTE: everything over the basering!
+      module A (matrix of the map), int k (power)
+      rings Tn is source- and Tm is image-ring with bases
+          resp. Ink and Imk.
+      M = max dim of Image, N - dim. of source      
+SEE ALSO: symmetricPower, exteriorPower"
+{
   def save = basering;
 
+  int n = nvars(save);
   int M = nrows(A);
   int N = ncols(A);
@@ -1181 +1202 @@
   int i, j;
 
-  /////////////////////////////////////////////////////////////////
-  def Tn = symmetricBasis(N, k); setring Tn;
-  ideal Ink = I;
-  dbprint(p-3, "Temporary Source Ring", basering);
-  dbprint(p, "S^k(Source Basis)", Ink);
-
-  /////////////////////////////////////////////////////////////////
-  def Tm = symmetricBasis(M, k); setring Tm;
-  ideal Imk = I;
-  dbprint(p-3, "Temporary Image Ring", basering);
-  dbprint(p, "S^k(Image Basis)", Imk);
-
-  /////////////////////////////////////////////////////////////////
-  def Rm = save + Tm;
-  setring Rm;
+//------------------------- compute matrix of single images ------------------
+  def Rm = save + Tm;  setring Rm; 
   dbprint(p-2, "Temporary Working Ring", Rm);
 
@@ -1207 +1215 @@
     for( j = M; j > 0; j-- )
     {
-      t = t + A[i][j] * var( nvars(save) + j); // tensor product!!!
+      t = t + A[i][j] * var(n + j);
     }
 
@@ -1213 +1221 @@
   }
 
-  dbprint(p-1, "Matrix of simgle images", B);
-
-  map TMap = Tn, B; ideal C = TMap(Ink); // apply S^k A to Source basis vectors... (Ink)
-  C = NF(C, std(0));
-
+  dbprint(p-1, "Matrix of single images", B);
+
+//------------------------- compute image ---------------------
+  // apply S^k(A): Tn -> Rm  to Source basis vectors Ink:
+  map TMap = Tn, B; 
+
+  ideal C = NF(TMap(Ink), std(0)); 
   dbprint(p-1, "Image Matrix: ", C);
 
-  // and write it in Image basis (Imk)
-
+
+//------------------------- write it in Image basis ---------------------
   ideal Imk = imap(Tm, Imk);
 
@@ -1247 +1257 @@
 
     D[i] = tt;
-//    tt;
   }
 
-  // D; print(D);
-
-
+//------------------------- map it back and return  ---------------------
   setring save;
-
-  // basering;
-
-  module AA = imap(Rm, D);
-
-//  nrows(AA) = Nk; // ????
-//  ncols(AA) = Mk;
-
-//  Nk, Mk;
-
-//  AA[Mk] = AA[Mk] * 1;
-
-  return(AA);
-}
-example
-{ "EXAMPLE:"; echo = 2;
-int save = printlevel; printlevel = 1;
+  return( imap(Rm, D) );
+}
+
+
+
+
+//////////////////////////////////////////////////////////////////////////////
+
+proc symmetricPower(module A, int k)
+"USAGE:    symmetricPower(A, k); A module, k int
+RETURN:   module: the k-th symmetric power of A
+NOTE:     the chosen bases and most of intermediate data will be shown if 
+          printlevel is big enough
+SEE ALSO: exteriorPower
+KEYWORDS: symmetric power
+EXAMPLE:  example symmetricPower; shows an example"
+{
+  int p = printlevel - voice + 2;
+
+  def save = basering;
+
+  int M = nrows(A);
+  int N = ncols(A);
+
+  string S = chooseSafeVarName("@", "@_e");
+
+//------------------------- choose source basis -------------------------
+  def Tn = symmetricBasis(N, k, S); setring Tn;
+  ideal Ink = symBasis;
+  export Ink;
+  dbprint(p-3, "Temporary Source Ring", basering);
+  dbprint(p, "S^k(Source Basis)", Ink);
+
+//------------------------- choose image basis -------------------------
+  def Tm = symmetricBasis(M, k, S); setring Tm;
+  ideal Imk = symBasis;
+  export Imk;
+  dbprint(p-3, "Temporary Image Ring", basering);
+  dbprint(p, "S^k(Image Basis)", Imk);
+
+//------------------------- compute and return S^k(A) in chosen bases --
+  setring save;  
+
+  return(mapPower(p, A, k, Tn, Tm));
+}
+example
+{ "EXAMPLE:"; echo = 2;
 
 ring r = (0),(a, b, c, d), dp; r;
-
-module A = a*gen(1), b * gen(1), c*gen(1), d * gen(1);
-
-print(symmetricPower(A, 2));
-print(symmetricPower(A, 3));
-
-module B = a*gen(1) + c* gen(2), b * gen(1) + d * gen(2);
-
+module B = a*gen(1) + c* gen(2), b * gen(1) + d * gen(2); print(B);
+
+// symmetric power over a commutative K-algebra:
 print(symmetricPower(B, 2));
 print(symmetricPower(B, 3));
 
-kill r;
-ring r = (0),(a, b, c, d), dp;def g = SuperCommutative();setring g;
-
-module A = a*gen(1), b * gen(1), c*gen(1), d * gen(1);
-
-print(symmetricPower(A, 2));
-print(symmetricPower(A, 3));
-
-module B = a*gen(1) + c* gen(2), b * gen(1) + d * gen(2);
-
-print(symmetricPower(B, 2));
-print(symmetricPower(B, 3));
-
-kill r;
-
-printlevel = save;
+// symmetric power over an exterior algebra:
+def g = SuperCommutative(); setring g; g;
+
+module B = a*gen(1) + c* gen(2), b * gen(1) + d * gen(2); print(B);
+
+print(symmetricPower(B, 2)); // much smaller!
+print(symmetricPower(B, 3)); // zero! (over an exterior algebra!)
+
 }
 
@@ -1307 +1331 @@
 "USAGE:    exteriorPower(A, k); A module, k int
 RETURN:   module: the k-th exterior power of A
-NOTE:     the chosen bases and temporary data
-          may be shown if printlevel is big enough.
-          Last rows may be invisible if zero.
+NOTE:     the chosen bases and most of intermediate data will be shown if 
+          printlevel is big enough. Last rows will be invisible if zero.
 SEE ALSO: symmetricPower
 KEYWORDS: exterior power
@@ -1320 +1343 @@
   int N = ncols(A);
 
-  int i, j;
-
-  /////////////////////////////////////////////////////////////////
-  def Tn = exteriorBasis(N, k); setring Tn;
-  ideal Ink = I;
+  string S = chooseSafeVarName("@", "@_e");
+
+//------------------------- choose source basis -------------------------
+  def Tn = exteriorBasis(N, k, S); setring Tn;
+  ideal Ink = extBasis;
+  export Ink;
   dbprint(p-3, "Temporary Source Ring", basering);
   dbprint(p, "E^k(Source Basis)", Ink);
 
-  /////////////////////////////////////////////////////////////////
-  def Tm = exteriorBasis(M, k); setring Tm;
-  ideal Imk = I;
+//------------------------- choose image basis -------------------------
+  def Tm = exteriorBasis(M, k, S); setring Tm;
+  ideal Imk = extBasis;
+  export Imk;
   dbprint(p-3, "Temporary Image Ring", basering);
   dbprint(p, "E^k(Image Basis)", Imk);
 
-
-  def Rm = save + Tm;
-  setring Rm;
-  dbprint(p-2, "Temporary Working Ring", Rm);
-
-  module A = imap(save, A);
-
-  ideal B; poly t;
-
-  for( i = N; i > 0; i-- )
-  {
-    t = 0;
-    for( j = M; j > 0; j-- )
-    {
-      t = t + A[i][j] * var( nvars(save) + j); // tensor product!!!
-    }
-
-    B[i] = t;
-  }
-
-  dbprint(p-1, "Matrix of simgle images", B);
-
-  map TMap = Tn, B; ideal C = TMap(Ink); // apply S^k A to Source basis vectors... (Ink)
-  C = NF(C, std(0));
-
-  dbprint(p-1, "Image Matrix: ", C);
-
-  // and write it in Image basis (Imk)
-
-  ideal Imk = imap(Tm, Imk);
-
-  module D; poly lm; vector tt;
-
-  for( i = ncols(C); i > 0; i-- )
-  {
-    t = C[i];
-    tt = 0;
-
-    while( t != 0 )
-    {
-      lm = leadmonom(t);
-      //    lm;
-      for( j = ncols(Imk); j > 0; j-- )
-      {
-        if( lm / Imk[j] != 0 )
-        {
-          tt = tt + (lead(t) / Imk[j]) * gen(j);
-          break;
-        }
-      }
-      t = t - lead(t);
-    }
-
-    D[i] = tt;
-//    tt;
-  }
-
-  // D; print(D);
-
-
+//------------------------- compute and return E^k(A) in chosen bases --
   setring save;
-
-  // basering;
-
-  module AA = imap(Rm, D);
-
-//  nrows(AA) = Nk; // ????
-//  ncols(AA) = Mk;
-
-//  Nk, Mk;
-
-//  AA[Mk] = AA[Mk] * 1;
-
-  return(AA);
-}
-example
-{ "EXAMPLE:"; echo = 2;
-  int save = printlevel; printlevel = 1;
-
-  ring r = (0),(a, b, c, d), dp; r;
-
-  module A = a*gen(1), b * gen(1), c*gen(1), d * gen(1);
+  return(mapPower(p, A, k, Tn, Tm));
+}
+example
+{ "EXAMPLE:"; echo = 2;
+  ring r = (0),(a, b, c, d, e, f), dp; 
+  r; "base ring:";
+
+  module B = a*gen(1) + c*gen(2) + e*gen(3), 
+             b*gen(1) + d*gen(2) + f*gen(3), 
+                        e*gen(1) + f*gen(3);
+
+  print(B);
+  print(exteriorPower(B, 2)); 
+  print(exteriorPower(B, 3)); 
+
+  def g = SuperCommutative(); setring g; g;
+
+  module A = a*gen(1), b * gen(1), c*gen(2), d * gen(2);
+  print(A);
 
   print(exteriorPower(A, 2));
-  print(exteriorPower(A, 3));
-
-  module B = a*gen(1) + c* gen(2), b * gen(1) + d * gen(2);
+
+  module B = a*gen(1) + c*gen(2) + e*gen(3), 
+             b*gen(1) + d*gen(2) + f*gen(3), 
+                        e*gen(1) + f*gen(3);
+  print(B);
 
   print(exteriorPower(B, 2));
   print(exteriorPower(B, 3));
 
-  kill r;
-  ring r = (0),(a, b, c, d), dp;def g = SuperCommutative();setring g;
-
-  module A = a*gen(1), b * gen(1), c*gen(1), d * gen(1);
-
-  print(exteriorPower(A, 2));
-  print(exteriorPower(A, 3));
-
-  module B = a*gen(1) + c* gen(2), b * gen(1) + d * gen(2);
-
-  print(exteriorPower(B, 2));
-  print(exteriorPower(B, 3));
-
-  kill r;
-
-  printlevel = save;
 }