Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fix 2 d perp #346

Merged
merged 3 commits into from
Jan 24, 2023
Merged

Fix 2 d perp #346

merged 3 commits into from
Jan 24, 2023

Conversation

jshipton
Copy link
Contributor

This fixes the issue with ufl.replace not finding the perp of a function defined on a subcomponent of a mixed function space. This affects shallow water simulations in the plane. The fix is to label (of course!) the Coriolis term with a label called perp, who's value is the function that defines the perp operator. This means that we can delay the application of the perp operator until after the subject of the form has been replaced. Ideally we would find a better fix for this...

…ect is replaced as otherwise replace can no longer find the subject
Copy link
Contributor

@tommbendall tommbendall left a comment

Choose a reason for hiding this comment

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

Even though we'd rather not have to fix this, I think it's very neat and good news that we've understood what's going on!

gusto/equations.py Outdated Show resolved Hide resolved
Copy link
Contributor

@tommbendall tommbendall left a comment

Choose a reason for hiding this comment

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

Thanks, all looks great now. And good to know we're also completing issue #272 !

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants