M2
M2 copied to clipboard
lyubeznikResolution bug
In the simplicial Complexes package if the lyubeznikResolution function is passed a monomialIdeal that has 6 or more minimal generators an index error occurs. It seems the issue stems from the following lines in M2/Macaulay2/packages/SimplicialComplexes/Code.m2
841 MinGens := flatten entries mingens I;
842 MO := opts.MonomialOrder;
843 R := QQ(monoid[vars(0..#(mingens I)-1)]);
In the last line # is applied to the matrix mingens I. This always returns 5, the number of keys in the dictionary used to store matrices. This can be fixed by replacing line 843 with R := QQ(monoid[vars(0..#(MinGens)-1)]); which instead uses the local variable MinGens.
Great catch! Would you like to open a pull request with the fix?