Changeset a3a116 in git for Singular/LIB/multigrading.lib
- Timestamp:
- Mar 25, 2011, 8:00:56 PM (12 years ago)
- Branches:
- (u'spielwiese', 'f6c3dc58b0df4bd712574325fe76d0626174ad97')
- Children:
- 611871aebf624c3f38dad1d47584a712f5bb90cf
- Parents:
- 119e27bf4a3415daccdfebbe95c3e15772a8462c
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/multigrading.lib
r119e27b ra3a116 4674 4674 0; 4675 4675 4676 intmat L = preimageLattice(P,B); 4677 4678 print(L); // should be a (3x2)-matrix with columns e.g. [1,1,-1] and [0,3,-2] (the generated lattice should be identical) 4679 4680 kill B, P, L; 4681 4676 // should be a (3x2)-matrix with columns e.g. [1,1,-1] and [0,3,-2] (the generated lattice should be identical) 4677 print(preimageLattice(P,B)); 4682 4678 } 4683 4679 … … 4784 4780 if ( r == n ) 4785 4781 { 4786 intmat U[ n][1]; // U now is the n-dim zero-vector4782 intmat U[1][n]; // U now is the n-dim zero-vector 4787 4783 } 4788 4784 else 4789 4785 { 4790 4786 // we want a matrix with column operations so we transpose 4791 list L = hermiteNormalForm(B, "transform"); //hermite(transpose(B), "transform"); 4787 intmat BB = transpose(B); 4788 list L = hermiteNormalForm(BB, "transform"); 4792 4789 intmat U = transpose(L[2]); 4793 4790 4791 4794 4792 // delete rows 1 to r 4795 4793 intmat Udel[nrows(U) - r][ncols(U)]; … … 4822 4820 4,8; 4823 4821 4824 // should result in a (2x4)-matrix with columns 4825 // [-1, 2], [2, -3], [-1, 0] and [0, 1] 4826 intmat U = projectLattice(B); 4822 // should result in a (2x4)-matrix such that the corresponding lattice is created by 4823 // [-1, 2], [-2, 3], [-1, 0] and [0, 1] 4824 print(projectLattice(B)); 4825 4826 // another example 4827 4828 intmat BB[4][2] = 4829 1,0, 4830 0,1, 4831 0,0, 4832 0,0; 4833 4834 // should result in a (2x4)-matrix such that the corresponding lattice is created by 4835 // [0,0],[0,0],[1,0],[0,1] 4836 print(projectLattice(BB)); 4837 4838 // one more example 4839 4840 intmat BBB[3][4] = 4841 1,0,1,2, 4842 1,1,0,0, 4843 3,0,0,3; 4844 4845 // should result in the (1x3)-matrix that consists of just zeros 4846 print(projectLattice(BBB)); 4827 4847 4828 4848 } … … 4914 4934 8,11; 4915 4935 4916 intmat U = intersectLattices(A,B); 4917 4918 // should result in a (2x4)-matrix with columns 4936 // should result in a (3x2)-matrix with columns 4919 4937 // e.g. [3, 3, 3], [-3, 0, 3] (the lattice should be the same) 4920 print( U);4938 print(intersectLattices(A,B)); 4921 4939 } 4922 4940
Note: See TracChangeset
for help on using the changeset viewer.