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

Add constant mode settings to Amg_data #6222

Merged
merged 2 commits into from
Feb 2, 2025

Conversation

quangx
Copy link
Contributor

@quangx quangx commented Jan 31, 2025

Addresses an issue in which constant_mode is created but not used.

#6221

@tjhei

Copy link
Member

@tjhei tjhei left a comment

Choose a reason for hiding this comment

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

Thank you.

Copy link
Member

@tjhei tjhei left a comment

Choose a reason for hiding this comment

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

Looks like iteration numbers are similar and only sometimes better. Maybe not too surprising as meshes are coarse in tests.

Maybe we should run one of the examples with melt transport before/after?

Copy link
Member

@gassmoeller gassmoeller left a comment

Choose a reason for hiding this comment

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

Yes, this looks like a clear bug fix. Can you also add a changelog entry that makes clear we now use the constant modes in the melt solver, just to make users aware results may change?

Do I see this correctly that this change generally reduces the number of iterations in tests, but in some cases increases them?

@@ -525,6 +525,7 @@ namespace aspect

Amg_data.smoother_sweeps = 2;
Amg_data.coarse_type = "symmetric Gauss-Seidel";
Amg_data.constant_modes=constant_modes;
Copy link
Member

Choose a reason for hiding this comment

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

Sorry to nitpick, but just to keep our code style here consistent, please make this into:

Suggested change
Amg_data.constant_modes=constant_modes;
Amg_data.constant_modes = constant_modes;

Copy link
Member

@gassmoeller gassmoeller left a comment

Choose a reason for hiding this comment

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

Thanks for the fix! Ready to merge.

@gassmoeller gassmoeller merged commit cacbbb2 into geodynamics:main Feb 2, 2025
8 checks passed
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