Post a reply
Note:If not registered, provide any username. For more comfort, register here.
Message body:
Enter your message here, it may contain no more than 60000 characters. 

:D :) :( :o :shock: :? 8) :lol: :x :P :oops: :cry: :evil: :twisted: :roll: :wink: :!: :?: :idea: :arrow: :| :mrgreen:
Font size:
Font colour
BBCode is ON
[img] is ON
[flash] is OFF
[url] is ON
Smilies are ON
Disable BBCode
Disable smilies
Do not automatically parse URLs
Confirmation of post
To prevent automated posts the board requires you to enter a confirmation code. The code is displayed in the image you should see below. If you are visually impaired or cannot otherwise read this code please contact the %sBoard Administrator%s.
Confirmation code:
Enter the code exactly as it appears. All letters are case insensitive, there is no zero.

Topic review - How to Compute Homology of a Graded Chain Complex
Author Message
  Post subject:  How to Compute Homology of a Graded Chain Complex  Reply with quote
I'm new to singular. I'd like to compute the Betti numbers of a graded k[x,y]-module M arising as a homology module of a chain complex of free graded k[x,y]-modules, and I'm running into difficulties: I'm using the "homology" command to obtain a presentation of M, but this seems to discard information about the graded structure on M. I wonder if my trouble is related to bug #67 in the Singular issue tracker.

In any case, I'd appreciate your advice on how to circumvent this problem. Below is some Singular code which demonstrates the issue in detail:

LIB "homolog.lib";
ring R=2,(x,y),dp;

//d2 will represent a homomorphism of free graded k[x,y]-modules
module d2;
intvec grades_C1=1,1,1;

//d1 will also represent a homomorphism of free graded k[x,y]-modules
module d1;
intvec grades_C0=1,1,1;

module triv=0;
//Compute a presentation for M:=ker(d1)/im(d2)
def M=homology(d2,d1,triv,triv);

//Next command outputs the grade information for the generators of M.
//Note the null output; homology apparently is not returning a graded object, whereas
//the generator for M should live at grade 1.

//Compute a free resolution of M
resolution resM=res(M,0);

//Since the grading was lost for M, the Betti numbers are off by 1.
print(betti(resM), "betti");
Post Posted: Mon Jun 05, 2017 4:27 am

It is currently Thu Mar 22, 2018 12:08 pm
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group