Changeset cb2f4dc in git


Ignore:
Timestamp:
Jan 19, 2015, 9:16:48 PM (9 years ago)
Author:
Oleksandr Motsak <motsak@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
374d8e1b4eb33ca5848ce22e0d05e58ecc908065
Parents:
4f2aed6d338208a9af82e4da005570718c40093f
git-author:
Oleksandr Motsak <motsak@mathematik.uni-kl.de>2015-01-19 21:16:48+01:00
git-committer:
Oleksandr Motsak <motsak@mathematik.uni-kl.de>2015-01-19 21:18:44+01:00
Message:
Starting to fix grdeg/grorder/grtranspose for zero columns + test for such a case
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/gradedModules.lib

    r4f2aed rcb2f4dc  
    237237
    238238  module L = lead(M[1..m]); // leading module-terms for input column vectors
    239   intvec d = deg(L[1..m]); // their degrees
     239  intvec d = deg(L[1..m]); // their degrees 
    240240  intvec c = leadexp(L[1..m])[n]; // their module-components
    241   intvec gr = w[c]; // weights
     241 
     242  w = intvec(-6665, w); // 0?????
     243  intvec gr = w[1+c]; // weights
    242244 
    243245  gr = gr + d; // finally we compute their graded degrees
     
    664666  grview(grpower( grshift(grzero(), 10), 5 ) );
    665667 
    666   module A = [x+y, x, 0 ], [0, x+y, y];
    667   attrib(A,"isHomog", intvec(1,1,1));
     668  module A = [x+y, x, 0, 0], [0, x+y, y, 0];
     669  attrib(A,"isHomog", intvec(0,0,0,1));
     670  attrib(A,"rank", 4);
    668671  grview(A);
    669672 
    670   matrix B[2][2] = z,0, 0,z;
    671   attrib(B,"isHomog", intvec(-2,-3));
     673  matrix B[3][1] = 0,x,y;
     674  attrib(B,"isHomog", intvec(2,1,1));
    672675  grview(B);
    673676
     
    678681  grview(C);
    679682
    680   def D = grsum(grpower(A, 2), grpower(B, 2));
     683  def D = grsum(
     684     grsum(grpower(A,2), grtwist(1,1)),
     685     grsum(grtwist(1,2), grpower(B,2))
     686     );
    681687
    682688  print(D);
    683689  homog(D);
    684   grview(D); 
     690  grview(D);
    685691
    686692  def D10 = grshift(D, 10);
     
    690696  grview(D10);
    691697
    692   "";  def T = grtwist(5, 10);
    693  
    694   print(T);
    695   homog(T);
    696   grview(T);
    697  
    698   "";  def TT = grsum(grtwist(3, 0), grtwist(5,-1));
    699  
    700   print(TT);
    701   homog(TT);
    702   grview(TT);
     698  def DD = grorder(D10);
     699
     700  print(DD);
     701  homog(DD);
     702  grview(DD);
     703
     704  def TTT = grtranspose(DD);
     705
     706  print(TTT);
     707  homog(TTT);
     708  grview(TTT);
    703709}
    704710
Note: See TracChangeset for help on using the changeset viewer.