Changeset 4ead7f in git

Timestamp:

Feb 25, 2010, 8:50:26 PM (13 years ago)

Author:

Motsak Oleksandr <motsak@…>

Branches:

(u'spielwiese', '0d6b7fcd9813a1ca1ed4220cfa2b104b97a0a003')

Children:

2b22df12b5c5cb19e00dce3e4ffe2323361f8170

Parents:

710469e33da8554c00e953aa1f5227161c8a1286

Message:

cleanup, better examples/tests for sheafCohBGG2

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

File:

• Singular/LIB/sheafcoh.lib (modified) (8 diffs)

Singular/LIB/sheafcoh.lib

@@ -600 +600 @@
    ring r=0,x(1..5),dp;
    module M=0;
-   def A=sheafCohBGG(0,-9,4);
+   def A=sheafCohBGG(M,-9,4);
    // cohomology of cotangential bundle on P^3:
    //-------------------------------------------
@@ -615 +615 @@
 static proc showResult( def R, int l, int h )
 {
-  int PL = printlevel - voice + 2;
+  int PL = 1; // printlevel - voice + 2;
 // int PL = printlevel + 1;
 
@@ -867 +867 @@
   };
 
-  printlevel = printlevel + 1;
+//  printlevel = printlevel + 1;
   Betti = showResult(Betti, l, h ); // Show usual form of cohomology table
-  printlevel = printlevel - 1;
+//  printlevel = printlevel - 1;
 
   if(PL > 0)
@@ -902 +902 @@
    echo = 2;
    int pl = printlevel;
+   int l,h, t;
+   
+   //-------------------------------------------
    // cohomology of structure sheaf on P^4:
    //-------------------------------------------
    ring r=32001,x(1..5),dp;
-   module M= getStructureSheaf();
-
-   int t = timer;
-   def A=sheafCohBGG(0,-9,4);
-   timer - t;
-
-   printlevel = voice; A=sheafCohBGG2(0,-9,4); printlevel = pl;
-
+
+   module M= getStructureSheaf(); // OO_P^4
+
+   l = -12; h = 12; // range of twists: l..h
+
+   printlevel = 0; 
+   //////////////////////////////////////////////
+   t = timer; 
+   
+   def A = sheafCoh(M, l, h); // global Ext method:
+   
+   "Time: ", timer - t;
+   //////////////////////////////////////////////
+   t = timer; 
+
+   A = sheafCohBGG(M, l, h);  // BGG method (without optimization):
+
+   "Time: ", timer - t;
+   //////////////////////////////////////////////
+   t = timer; 
+   
+   A = sheafCohBGG2(M, l, h); // BGG method (with optimization)
+   
+   "Time: ", timer - t;
+   //////////////////////////////////////////////
+   printlevel = pl;
+   
+   kill A, r; 
+
+   //-------------------------------------------
    // cohomology of cotangential bundle on P^3:
    //-------------------------------------------
@@ -919 +944 @@
    module M = getCotangentialBundle();
 
-   int t = timer;
-   def B=sheafCohBGG(M,-8,4);
-   timer - t;
-
-   printlevel = voice; B=sheafCohBGG2(M,-8,4); printlevel = pl;
+   l = -12; h = 11; // range of twists: l..h
+
+
+   //////////////////////////////////////////////
+   printlevel = 0; 
+   t = timer; 
+   
+   def B = sheafCoh(M, l, h); // global Ext method:
+   
+   "Time: ", timer - t;
+   //////////////////////////////////////////////
+   t = timer; 
+
+   B = sheafCohBGG(M, l, h);  // BGG method (without optimization):
+
+   "Time: ", timer - t;
+   //////////////////////////////////////////////
+   t = timer; 
+   
+   B = sheafCohBGG2(M, l, h); // BGG method (with optimization)
+   
+   "Time: ", timer - t;
+   //////////////////////////////////////////////
+   printlevel = pl;
 }
 
@@ -1182 +1226 @@
     intvec v = 0;
     attrib(M,"isHomog",v);
-    homog(M);
+//    homog(M);
 
     attrib(M, "isCoker", 1);
 
-    attrib(M);
+//     attrib(M);
     return(M);
   };
@@ -1206 +1250 @@
           export M;
 
-//          keepring @@@@QQ; // This is a bad idea... :(
+//          keepring @@@@QQ; // This is a bad idea... :(?
           return (@@@@QQ);
         }
@@ -1278 +1322 @@
   resolution T1=mres(maxideal(1),3);
   module M=T1[3];
-  attrib(M,"isHomog");
-  homog(M);
-  attrib(M, "isCoker", 1);
-  attrib(M);
+//  attrib(M,"isHomog");
+//  homog(M);
+  attrib(M, "isCoker", 1);  
+  // attrib(M);
   return (M);
 };