Changeset 4e3305 in git

Timestamp:

Jan 21, 2015, 12:25:22 AM (9 years ago)

Author:

Oleksandr Motsak <motsak@…>

Branches:

(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')

Children:

ca38640ca6928ad815ae7b2bd1f3516f383315affae8e849c66b071ec3b0afcb4779e35629c0d034

Parents:

374d8e1b4eb33ca5848ce22e0d05e58ecc908065

git-author:

Oleksandr Motsak <motsak@mathematik.uni-kl.de>2015-01-21 00:25:22+01:00

git-committer:

Oleksandr Motsak <motsak@mathematik.uni-kl.de>2015-01-21 00:27:24+01:00

Message:

Added constructor for graded objects: grobj

add: example with zero column

File:

• Singular/LIB/gradedModules.lib (modified) (15 diffs)

Singular/LIB/gradedModules.lib

@@ -15 +15 @@
 
     grzero()        presentation of basering(0)^1
+    grobj(M,w)      construct a graded module object 
     grshift(M,d)    shift graded module M by d
     grtwist(r,d)    presentation of basering(d)^r
@@ -304 +305 @@
   // grades & ordering permutation for N.  d[p] should be sorted!
   
-  def N = module(M[p]);  // reorder the starting ideal/module
-
-  attrib( N, "isHomog", w ); // set the grading
-  attrib( N, "rank", size(w) ); // set the grading
+  def N = grobj(module(M[p]), w);  // reorder the starting ideal/module
 
   d = d[p];
@@ -439 +437 @@
  else              { (N,d) = reorder(M, -1); }
 
- kill M; module M = transpose(N);
- attrib( M, "isHomog", -d ); // set the grading
+ kill M; module M = grobj(transpose(N), -d);
 
 // "b";  grview(M);
@@ -447 +444 @@
  // reverse order:
  (N,d) = reorder(M, 1); kill M;
-
-/*
- module M = transpose(N); 
- attrib( M, "isHomog", -d ); // set the grading
-
-// "c";  grview(M);
- 
- kill N,d; module N; intvec d;
- // reverse order:
- (N,d) = reorder(M, -1); kill M;
- module M = transpose(N); 
- attrib( M, "isHomog", -d ); // set the grading
-
-// "d";  grview(M);
- 
- kill N,d; module N; intvec d;
- // reverse order:
- (N,d) = reorder(M, 1);
- 
- kill M;
-*/
 
 // "e"; grview( N );
@@ -524 +500 @@
   if( size(#) > 0 ) { (N,d) = reorder(M, 1, #[1]); }
   else              { (N,d) = reorder(M, 1); }  
-  kill M; module M = transpose(N);
-  attrib( M, "isHomog", -d ); // set the grading
+  kill M; module M = grobj(transpose(N), -d);
 
 // "b";  grview(M);
@@ -533 +508 @@
   (N,d) = reorder(M, -1); kill M;
  
-  module M = transpose(N); 
-  attrib( M, "isHomog", -d ); // set the grading
+  module M = grobj(transpose(N), -d);
 
 // "c";  grview(M);
@@ -611 +585 @@
 {
  module Z = 0;
- attrib(Z,"isHomog",intvec(0));
- return (Z);
+ return ( grobj(Z,intvec(0)) );
 }
 
@@ -653 +626 @@
   module T; T[r] = 0; T = T, module(transpose(B));
   module AB = module(A), transpose(T);
-
-  attrib(AB, "isHomog", c);
-  attrib(AB, "rank", size(c));
-  
-  return(AB);
+  
+  return(grobj(AB, c));
 }
 example
@@ -666 +636 @@
   grview(grpower( grshift(grzero(), 10), 5 ) );
   
-  module A = [x+y, x, 0, 0], [0, x+y, y, 0];
-  attrib(A,"isHomog", intvec(0,0,0,1));
-  attrib(A,"rank", 4);
+  module A = grobj( module([x+y, x, 0, 0], [0, x+y, y, 0]), intvec(0,0,0,1) );  
   grview(A);
   
-  matrix B[3][1] = 0,x,y;
-  attrib(B,"isHomog", intvec(2,1,1));
+  module B = grobj( module([0,x,y]), intvec(15,1,1) );
   grview(B);
 
-  def C = grsum(A,B);
+  module C = grsum(A,B);
 
   print(C);
@@ -681 +648 @@
   grview(C);
 
-  def D = grsum(
+  module D = grsum(
      grsum(grpower(A,2), grtwist(1,1)), 
      grsum(grtwist(1,2), grpower(B,2))
@@ -690 +657 @@
   grview(D);
 
-  def D10 = grshift(D, 10);
+  module D10 = grshift(D, 10);
 
   print(D10);
@@ -696 +663 @@
   grview(D10);
 
-  def DD = grorder(D10);
+  module DD = grorder(D10);
 
   print(DD);
@@ -702 +669 @@
   grview(DD);
 
-  def TTT = grtranspose(DD);
+  module TTT = grtranspose(DD);
 
   print(TTT);
   homog(TTT);
   grview(TTT);
+
+  module F = grobj( module([x,y,0]), intvec(1,1,5) );
+  grview(B);
+  
+  module T = grsum( F, grsum( grtwist(1, 10), B ) );
+  grview(T);
+
 }
 
@@ -715 +689 @@
 {
   intvec a = attrib(M, "isHomog");
-  attrib(M, "isHomog", intvec( a + intvec(d:size(a))) );
-  return (M);
+  return (grobj(M, intvec( a + intvec(d:size(a))) ));
 }
 
@@ -731 +704 @@
 {
   module Z; Z[a] = 0; 
-  attrib(Z, "isHomog", intvec(d:a));
-  attrib(Z, "rank", a);
+  Z = grobj(Z, intvec(d:a));
 
   ASSUME(2, grisequal(Z, grpower( grshift(grzero(), d), a ) )); // optional check
   return(Z); 
-}  
+}
+
+proc grobj(module M, intvec w)
+""
+{
+  attrib(M, "isHomog", w);
+  attrib(M, "rank", size(w));
+  return (M);
+}