Skip to content

Improvements in WeylAlgebras#4147

Open
mahrud wants to merge 3 commits intoMacaulay2:developmentfrom
mahrud:Dmodules
Open

Improvements in WeylAlgebras#4147
mahrud wants to merge 3 commits intoMacaulay2:developmentfrom
mahrud:Dmodules

Conversation

@mahrud
Copy link
Copy Markdown
Member

@mahrud mahrud commented Feb 25, 2026

  • fixed documentation of WeylAlgebras
  • added minor improvements in WeylAlgebras
  • cached holonomic basis in holonomicRank

Copy link
Copy Markdown
Contributor

@antonleykin antonleykin left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Most of the cnages are straightforward and good.
There is an issue in holonomicRank that should be addressed but I'm not sure the current solution is the right one.

else holRank = numgens source basis redI;
holRank
)
-- Note: we don't compute dim redI, because coefficientRing W might nonzero dimension
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

What is your definition of the holonomic rank when the coefficient ring is not a field?
Can you give a use example?

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Good question! I haven't thought about the case when the coefficient ring is not a field, but as for use case, this is the backbone of the ConnectionMatrices packages, where a modified version of this method is currently being used:
https://github.com/Macaulay2/M2/blob/stable/M2/Macaulay2/packages/ConnectionMatrices/holonomic.m2

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I guess really the question is what does "basis" do when the coefficient ring is not a field, and I think it just gives a generating set? I'm not an expert, is holonomic rank over integers well defined?

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'm not sure how to resolve this.
@mahrud, can you figure out what ConnectionMatrices authors think?
@michaelPerlman @lorinczandras @christineBerkesch , any opinion on generalization of holonomic rank?

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Sorry, I think I'm confused. I thought you were asking out of curiosity. Is anything else in Dmodules packages designed with non-field coefficient rings in mind? The package tutorial explicitly says:

This package is mostly concerned with computations in the Weyl algebra, the ring of differential operators over affine space (over a field of characteristic zero).

I also can't find any algorithms in the SST book that don't assume working over a field either.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants