Macaulay2计算ℙ³上商层Q上同调的代码报错求助

错误根源 #

你用A/B直接做商的方式不符合Macaulay2的规则：

1. 做商的两个模必须共享相同的 ambient module，但OO_X^1 (-1)和OO_X^4的扭转不同，不满足这个条件；
2. 正合序列中的商层Q是同态的余核，不是抽象的模商——你完全没定义从$\mathcal{O}(-1)$到$\mathcal{O}^{\oplus4}$的单射映射，这是正合序列的核心部分。

修正后的完整代码 #

下面是可以正确运行的代码，包含映射定义、余核构造和上同调计算：

X = Proj(QQ[x,y,z,w])
-- 定义序列中的两个层
O_minus1 = OO_X(-1)
O_4 = OO_X^4
-- 定义从O(-1)到O^4的单射映射（这里用坐标映射作为示例，可根据你的需求替换）
f = map(O_4, O_minus1, matrix{{x, y, z, w}})
-- 构造商层Q：映射f的余核
Q = cokernel f
-- 计算Q的各阶上同调（ℙ³的上同调非平凡的是H⁰到H³）
for i from 0 to 3 do (
    print("H^" | toString(i) | "(Q): " | toString(cohomology(i, Q)))
)

关键细节解释 #

• 映射定义：map(目标层, 源层, 矩阵)是Macaulay2中定义层同态的标准方式，这里的矩阵{{x,y,z,w}}对应把$\mathcal{O}(-1)$的截面$s$映射为$\mathcal{O}^{\oplus4}$的截面$(x s, y s, z s, w s)$，是一个单射。
• 余核构造：cokernel f就是正合序列中$f: \mathcal{O}(-1) \to \mathcal{O}^{\oplus4}$的商层Q，符合你给出的正合序列定义。
• 上同调计算：cohomology(i, Q)返回第i个上同调群作为QQ上的向量空间，循环遍历0到3阶即可得到所有非平凡的上同调结果。

内容的提问来源于stack exchange，提问作者Protágoras