Changeset 8f4b23 in git for Singular/LIB/sheafcoh.lib
- Timestamp:
- May 19, 2005, 3:32:35 PM (19 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- 5c1cb7f29d5ce62aa91dc921b7cb8f34a92fe883
- Parents:
- cc5e38df80c917c52275aa45b1835bb4f86a0f7f
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/sheafcoh.lib
rcc5e38 r8f4b23 1 1 S/////////////////////////////////////////////////////////////////////////////// 2 version="$Id: sheafcoh.lib,v 1. 4 2005-05-19 08:14:51Singular Exp $";2 version="$Id: sheafcoh.lib,v 1.5 2005-05-19 13:32:35 Singular Exp $"; 3 3 category="Commutative Algebra"; 4 4 info=" … … 383 383 /////////////////////////////////////////////////////////////////////////////// 384 384 385 proc sheafCoh(module M,int l,int h )385 proc sheafCoh(module M,int l,int h,list #) 386 386 "USAGE: sheafCoh(M,l,h); M module, l,h int 387 387 ASSUME: @code{M} is graded and is either 0 or comes assigned with an … … 405 405 NOTE: The procedure is based on local duality as described in [Eisenbud: 406 406 Computing cohomology. In Vasconcelos: Computational methods in 407 commutative algebra and algebraic geometry. Springer (1998)]. 407 commutative algebra and algebraic geometry. Springer (1998)].@* 408 By default, the procedure uses @code{mres} to compute the Ext 409 modules. If called with the additional parameter @code{\"sres\"}, 410 the @code{sres} command is used instead. 408 411 SEE ALSO: dimH, sheafCohBGG 409 412 EXAMPLE: example sheafCoh; shows an example 410 413 " 411 414 { 415 int use_sres; 412 416 if( typeof(attrib(M,"isHomog"))!="intvec" ) { 413 417 if (size(M)==0) { attrib(M,"isHomog",0); } 414 418 else { ERROR("No admissible degree vector assigned"); } 419 } 420 if (size(#)>0) { 421 if (#[1]=="sres") { use_sres=1; } 415 422 } 416 423 int i,j; … … 420 427 int col=h-l+1; 421 428 intmat newBetti[n+1][col]; 422 list L=Ext_R(v,M,1)[2]; // list of GB for Ext_R 429 if (use_sres) { list L=Ext_R(v,M,1,"sres")[2]; } 430 else { list L=Ext_R(v,M,1)[2]; } 423 431 for (i=l; i<=h; i++) { 424 432 N0=L[n+2]; … … 565 573 print(betti(FI),"betti"); 566 574 module F=FI[2]; 567 def A1=sheafCoh(F,-6,6); 568 def A2=sheafCohBGG(F,-6,6); 575 int t=timer; 576 def A1=sheafCoh(F,-8,8); 577 timer-t; 578 t=timer; 579 def A2=sheafCohBGG(F,-8,8); 580 timer-t; 569 581 570 582 LIB "sheafcoh.lib";
Note: See TracChangeset
for help on using the changeset viewer.